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.

Development of modified embedded-atom model and molecular dynamics simulation of cesium

Molecular Dynamics Simulations for Melting Temperatures of Caf2

Melting Curves And Entropy Of Fusion Of Body-Centered Cubic Tungsten Under Pressure

Structure and Thermodynamic Properties of Beo: Empirical Corrections in the Quasiharmonic Approximation

Ab Initio Calculation of Lattice Dynamics and Thermodynamic Properties of Beryllium

Development of Lattice Inversion Modified Embedded Atom Method and Its Applications

Modified embedded-atom method interatomic potentials for the Mg-Al alloy system

Cohesive zone method based multi particle finite element simulation of compaction densification process of Al and NaCl laminar composite powders

Embedded-atom-method interatomic potentials from lattice inversion

Recent Efforts in Modelling Combined Rate, Size and Thermal Effects on Single Crystal Strength

Molecular dynamics simulation of CaSiO3 perovskite on decompression and re-compression

Modelling Silica using MACE-MP-0 Machine Learnt Interatomic Potentials

MEAM interatomic potential for thermodynamic and mechanical properties of lithium allotropes

An Improved Modified Embedded-Atom Method Potential to Fit the Properties of Silicon at High Temperature

A MD Simulation Study of the Coefficients of Thermal Expansion for β-HMX Crystal

Theoretical study on intermolecular interactions and thermodynamic properties of LLM-105

Amorphous Zirconia-doped Tantala modeling and simulations using explicit multi-element spectral neighbor analysis machine learning potentials (EME-SNAP)

First-Principles Calculations of Elastic and Thermodynamic Properties of Cubic Cdte

Structure phase transformation and equation of state of cerium metal under pressures up to 51 GPa

Modified embedded-atom method interatomic potentials for Ti and Zr

Thermomechanical Properties of Non-Stoichiometric Gadolinium Doped Ceria by Molecular Dynamics Simulations