Theoretical study on intermolecular interactions and thermodynamic properties of LLM-105
Fan Jun-Yu,Gao Nan,Wang Peng-Ju,Su Yan,,,,
DOI: https://doi.org/10.7498/aps.73.20231696
IF: 0.906
2024-01-01
Acta Physica Sinica
Abstract:2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a typical high-energy and low-sensitivity energetic material (EM), which has excellent detonation performance and thermal stability. Based on quasi-harmonic approximation, the dispersion corrected density functional theory was used to study the intermolecular interactions and thermodynamic properties of energetic LLM-105 crystal. By introducing the zero-point energy and temperature effect corrections, PBE-D3 dispersion correction scheme could significantly improve the calculation accuracy of structural parameters at experimental temperature (294 K). The temperature dependent lattice parameters and thermal expansion coefficients exhibit strong anisotropy, especially the thermal expansivity in b -axis orientation (intermolecular layers) is much higher than that in the ac plane (intramolecular layers). Through Hirshfeld surface and fingerprint analysis, it was found that the intermolecular interactions of LLM-105 are mainly O···H hydrogen bonding interactions. The change of intermolecular interactions would result in the rotation of nitro group, which could contribute to form new hydrogen-bonding interaction pattern. Mulliken population analysis shows that the bond order of C−NO 2 bond is more sensitive to the change of temperature, which may be the trigger bond for the high-temperature decomposition reaction of LLM-105. The fundamental thermodynamic properties of EMs can not only provide key parameters for mesoscopic or macroscopic thermodynamic simulations, but also gain theoretical insights into the temperature effects of EMs. Specific heat capacity reflects the amount of heat to be supplied to heat the matter and it is important to the risk assessment of EMs during storage or when exposed to external thermal stimuli. Herein, the basic thermodynamic parameters, such as heat capacity, entropy, bulk modulus and elastic constants under different conditions are predicted. Among them, the calculated heat capacity and entropy describe the nonlinear behaviors within the temperature range of 0−500 K, and the calculated isobaric heat capacity C p ( T ) is in good agreement with the available experimental measurements. The elasticity of materials describes the macroscopic response of crystals to external forces, and the bulk modulus B 0 of molecular crystals can be determined through the equation of state, which is an important parameter for evaluating material stiffness. The bulk modulus under adiabatic conditions is in reasonable agreement with experimental values, and the evolution of bulk modulus with temperature reflects the softening behaviors of LLM-105 at temperature. Furthermore, the complete set of second-order elastic constants (SOECs) of LLM-105 is calculated and 13 independent SOECs ( C 11 , C 12 , C 13 , C 15 , C 22 , C 23 , C 25 , C 33 , C 35 , C 44 , C 46 , C 55 , C 66 ) are predicted. With the increasing temperature, all elastic constants gradually decrease due to the weaken of intermolecular interactions of LLM-105. Overall, these results would fundamentally provide a deep understanding of temperature effects and serve as a reference for the experimental measurement of the thermodynamic parameters of EMs.
physics, multidisciplinary