A 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) /1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)-isopropanol (IPA) interfacial model was constructed to investigate the effect of temperature on cocrystal morphology. A constant volume and temperature molecular dynamics (NVT-MD) simulation was performed on the interfacial model at various temperatures (295-355 K, 20 K intervals). The surface electrostatic potential (ESP) of the CL-20/HMX cocrystal structure and IPA molecule were studied by the B3LYP method at 6-311++G (d, p) level. The surface energies, polarities, adsorption energy, mass density distribution, radial distribution function (RDF), mean square displacement (MSD) and relative changes of attachment energy were analyzed. The results show that polarities of (1 0 0) and (0 1 1) cocrystal surfaces may be more negative and affected by IPA solvent. The adsorption energy per area indicates that growth of the (1 0-2) face in IPA conditions may be more limited, while the (1 0 0) face tends to grow more freely. MSD and diffusion coefficient (D) analyses demonstrated that IPA molecules gather more easily on the cocrystal surface at lower temperatures, and hence have a larger effect on the growth of cocrystal faces. RDF analysis shows that, with the increasing of temperature, the strength of hydrogen bond interactions between cocrystal and solvent becomes stronger, being highest at 335 K for the (1 0 0) and (0 1 1) interfacial models. Results of relative changes of modified attachment energy show that (1 0 0) and (0 1 1) faces tends to be larger than other faces. Moreover, the predicted morphologies at 295 and 355 K are consistent with experimental values, proving that the CL-20/HMX-IPA interfacial model is a reasonable one for this study. Graphical Abstract Construction of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) /1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)-isopropanol (IPA) interfacial model, analysis, and morphology prediction of cocrystal.