Anisotropy in surface-initiated melting of the triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene: A molecular dynamics study

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB)

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. We present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.

A Hotspot's Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material's microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. The origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

A theoretical investigation into the cooperativity effect on the TNT melting point under external electric field

External electric field has been regarded as an effective tool to induce the variation of melting points of molecular crystals. The melting point of 2,4,6-trinitrotoluene (TNT) was calculated by molecular dynamics simulations under external electric field, and the electric field effects on the cooperativity effects of the ternary (TNT)₃ were investigated at the M06-2X/6-311+G(d) and ωB97X-D/6-311++G(2d,p) levels. The results show that the melting points are decreased while the intermolecular interactions are strengthened under the external electric fields, suggesting that the intermolecular interactions cannot be used to explain the decreased melting points. A deduction based on the cooperativity effect is put forward: the enhanced cooperativity effects create the more serious defects in the melting process of the molecular crystal under the external electric fields, and simultaneously the local order parameters are decreased, leading to the decreased melting point. Thus, the cooperativity effect stemmed from the intermolecular C-H∙∙∙O H-bonding interactions controls the change of TNT melting point under the external electric field. Employing the information-theoretic approach (ITA), the origin of the cooperativity effects on the melting points of molecular crystal is revealed. This study opens a new way to challenge the problems involving the melting points for the molecular crystal under the external electric fields. However, note that above deduction needs to be improved; after all, the simple (TNT)₃ model cannot replace the crystal structure.

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic P1̅ phase. First-principles calculations were performed with Perdew-Burke-Ernzerhof (PBE) and PBE0 functionals including thermal effects using the P1̅ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4-10 GPa compared with experimental measurements indexed with P1̅ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

High Explosive Ignition through Chemically Activated Nanoscale Shear Bands

Shock initiation and detonation of high explosives is considered to be controlled through hot spots, which are local regions of elevated temperature that accelerate chemical reactions. Using classical molecular dynamics, we predict the formation of nanoscale shear bands through plastic failure in shocked 1,3,5-triamino-2,4,6-trinitrobenzene high explosive crystal. By scale bridging with quantum-based molecular dynamics, we show that shear bands exhibit lower reaction barriers. While shear bands quickly cool, they remain chemically activated and support increased reaction rates without the local heating typically evoked by the hot spot paradigm. We describe this phenomenon as chemical activation through shear banding.

First-Principles-Based Force Field for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a highly promising energetic material (EM) with high safety. Understanding its microscopic response mechanisms within the external stimulus is meaningful for the design of EMs. In order to comprehend the complicated phenomena, it is necessary to employ molecular simulation methods to investigate the response mechanisms with the force field (FF) at an atomic level. In this work, we developed a tailored FF for LLM-105 based on first-principles calculations. The validity of the FF was evaluated by molecular dynamics simulations. The structural parameters of LLM-105 predicted by FF are in good agreement with the experimental values, such as lattice constant, bond length, bond angle, dihedral angle and center of mass, and so forth. Moreover, the FF possesses good performance to describe the structural response on pressure accurately. In general, our work not only builds a balanced FF in gas and condensed phases, but also provides a useful tool to study the properties about LLM-105 at a large scale, which is helpful to improve the understanding about the balance between energy and safety in EMs.