Periodic DFT study of structural transformations of cocrystal CL-20/MMI under high pressure

J Mol Model. 2024 Apr 5;30(5):124. doi: 10.1007/s00894-024-05918-z.

Abstract

Context: The crystal and molecular structure, electronic properties, optical parameters, and elastic properties of a 1:2 hexanitrohexaazaisowurtzitane (CL-20)/2-mercapto-1-methylimidazole (MMI) cocrystal under 0 ~ 100 GPa hydrostatic pressure were calculated. The results show that the cocrystal CL-20/MMI undergoes three structural transitions at 72 GPa, 95 GPa, and 97 GPa, respectively, and the structural transition occurs in the part of the MMI compound. Structural mutations formed new bonds S1-S2, C2-C7, and N1C5 at 72GPa, 95 GPa, and 97 GPa, respectively. Similarly, the formation of new bonds is confirmed on the basis of an analysis of the changes in lattice constants, cell volumes, and partial densities of states (PDOS) for S1, S2, C2, C7, N1, and C3 at the corresponding pressures. The optical parameters show that the pressure makes the peaks of various optical parameters of CL-20/MMI larger, and the optical activity is enhanced. The optical parameters also confirm the structural mutation of CL-20/MMI under the corresponding pressure.

Method: CL-20/MMI was calculated by using the first-principles norm-conservative pseudopotential based on density functional theory (DFT) in the CASTEP software package. For the optimization results, the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method is selected to optimize the geometry of the cocrystal in the range of 0-100 GPa. GGA/PBE (Perdew-Burke-Ernzerhof) was selected to relax the cocrystal CL-20/MMI fully without constraints at atmospheric pressure. The sampling scheme in the Brillouin zone [10] is the Monkhorst-Pack scheme, and the number of k-point grids was 2 × 2 × 2. By contrast, this study will use the LDA method to calculate.

Keywords: CL-20/MMI; Density functional theory; High pressure; Optical properties; Theoretical study.