The aim of this paper is to determine the effect of polymer density, correlated to the comonomer content, and nanosilica addition on the mechanical and Environmental Stress Cracking Resistance (ESCR) characteristics of high-density polyethylene (HDPE). In this regard, five HDPE samples with similar Melt Flow Index (MFI) and molar mass but various densities were acquired from a petrochemical plant. Two polymerization reactors work in series and differ only in the amount of 1-buene comonomer fed to the second reactor. To ascertain the microstructure of the studied samples, GPC and SSA (successive self-nucleation and annealing) analyses were accomplished. All samples resulted having similar characteristics but slightly various SCB/1000 C=7.26-9.74 (SCB=Short Chain Branching). Consequently, meanwhile studied HDPEs reveal similar notched impact and stress at yield values, the tensile modulus, stress-at-break, and elongation-at-break tend to demonstrate different results with the SCB content. More significantly, ESCR characteristic varied considerably with SCB/1000 C extent, so that higher amount of SCB acknowledged advanced ESCR. Notably, blending HDPE sample containing higher amount of SCB/1000 C, with 3 wt.% of chemically modified nanosilica enhanced ESCR characteristic by 40 %. DFT (Density Functional Theory) calculations unveiled the role of the comonomer, quantitatively by binding energies and qualitatively by Non Covalent Interaction (NCI) plots.