Computational insights into the structure and decomposition behaviors of 2,4,6-triamino-5-nitropyrimidine-1,3-dioxide under high pressure up to 10 GPa.

CONTEXT

Inspired by the recent successful synthesis of the energetic compound 2,4,6-triamino-5-nitropyrimidine-1,3-dioxide (ICM-102), which displayed a good balance between high energy and sensitivity, the response of the structure and decomposition behaviors of ICM-102 to high pressure was systematically investigated using first principle calculations. ICM-102 exhibited a graphite-like layer structure, with the c-axis and the a-axis mainly contributing to the distance between the molecular planes. As the pressure increased from 1 atm to 10 GPa, this distance decreased from 3.166 to 2.689 Ǻ. The hydrogen bonds had the most contribution to the non-covalent interactions within the same molecular planes, resulting in the b-axis discontinuity. However, van der Waals interactions gradually appeared between molecular planes as the pressure increased to 2.5 GPa. Based on the analysis of crystal orbitals, the distribution of π bonds and the Laplacian bond order (LBO), it was determined that the generation mechanism of HO molecules involved the cleavage of N-O (coordinated oxygen atoms), followed by intermolecular hydrogen transfer reactions, and ultimately the formation of HO molecules through competition with H atoms in the amino groups within the same molecules. More importantly, the pressure dependence of LBO values for N-O revealed that high pressure could inhibit the ICM-102 decomposition process due to reinforcing hydrogen bonds and van der Waals interactions. This work will deepen our understanding of the stability of ICM-102 under high pressure and provide a helpful reference for its potential detonation applications.

METHODS

All simulations, including geometry optimization and vibration analysis under quasi-hydrostatic pressure, were conducted using the CP2K code. The PBE function and the Goldk-Teter-Hutter (GTH) pseudopotential with the double-ζ-with-polarization (DZVP) basis set were employed. Additionally, the semiempirical dispersion correction D3 (BJ) was used to account for the intermolecular dispersion force. The simulations were performed under periodic boundary conditions, with a finest grid level cutoff set to 400 Ry for the Γ point. The Broyden-Flecher-Goldfarb-Shanno (BFGS) optimization method was used, with tighter convergence criteria applied for the subsequent calculations of infrared spectra. Finally, the wave-function analysis, such as non-covalent interaction and LBO, was conducted using the Multiwfn and VMD packages.