Treatment failure in patients with liver hepatocellular carcinoma (LIHC) is primarily caused by tumor progression and therapy resistance. Tumor immunity plays a crucial role in regulating the homeostasis of cells through the process of programmed cell death (PCD). However, the expression profile and clinical significance of PCD-related genes in LIHC require further investigation. In this study, we analyzed twelve commonly observed PCD patterns to construct a prognostic model. We collected RNA-seq data, genomics, and clinical information from TCGA-LIHC and GSE14520 cohorts to validate the prognostic gene signature. We discovered 75 PCD-related differentially expressed genes (DEGs) with prognostic significance in LIHC. Using these genes, we constructed a PCD-related score (PCDscore) with an 11-gene signature through LASSO COX regression analysis. Validation in the GSE14520 cohort demonstrated that LIHC patients with high PCDscore had poorer prognoses. Unsupervised clustering based on the 11 model genes revealed 3 molecular subtypes of LIHC with distinct prognoses. By incorporating PCDscore with clinical features, we constructed a highly predictive nomogram. Additionally, PCDscore was correlated with immune checkpoint genes and immune cell infiltration. LIHC patients with high PCDscore exhibited sensitivity to common chemotherapy drugs (such as cisplatin and docetaxel). To summarize, our study developed a novel PCDscore model that comprehensively analyzed different cell death modes, providing an accurate prediction of clinical prognosis and drug sensitivity for LIHC patients.