Spin parameter optimization for spin-polarized extended tight-binding methods

We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential parameters for a given reference dataset, allowing for a nuanced understanding of their impact on diverse molecular properties. In the case of the W4-11 dataset, substantial error reduction was achieved, demonstrating the potential of the optimization. Transferability of the optimized spin-polarization constants over different properties, however, is limited, as we demonstrate by applying the optimized parameters on a set of singlet-triplet gaps in carbenes. Further studies on ionization potentials and electron affinities highlight some inherent limitations of current extended tight-binding methods that can not be resolved by simple parameter optimization. We conclude that the significantly improved accuracy strongly encourages the present re-optimization of the spin-polarization constants, whereas the limited transferability motivates a property-specific optimization strategy.

Size effect on the pyrolysis of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) nanoparticles: a ReaxFF molecular dynamics study

The nanoscale form of the typical insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) exhibits the capability to improve the energy performances while maintaining low sensitivity compared to raw TATB. Investigating the particle size effect on the intrinsic pyrolysis mechanisms facilitates the selection of the TATB particle size in applications to ensure efficient energy release and high safety levels. However, the intrinsic mechanism of this effect remains unclear. This study focuses on pyrolysis as a prerequisite behavior for energy release, employing reactive molecular dynamics simulations to investigate the pyrolysis of TATB nanoparticles with different sizes, aiming to explore the qualitative changes in thermal properties at the atomic level. Results demonstrate that with increasing particle size, the decomposition rate of TATB decreases. Smaller particles exhibit a propensity towards dehydrogenation and C-NO2 bond cleavage reactions. However, larger nano-TATB particles demonstrate a preference for dimerization, which results in the formation of clusters with greater polymerization and increased stability. The highly polymerized clusters are stable under thermal stimulation, inhibiting further decomposition of TATB. These insights reveal the mechanism underlying the qualitative change in the energy performance of TATB nanoparticles at the atomic level.

The Structure Properties of Carbon Materials Formed in 2,4,6-Triamino-1,3,5-Trinitrobenzene Detonation: A Theoretical Insight for Nucleation of Diamond-like Carbon

The structure and properties of nano-carbon materials formed in explosives detonation are always a challenge, not only for the designing and manufacturing of these materials but also for clearly understanding the detonation performance of explosives. Herein, we study the dynamic evolution process of condensed-phase carbon involved in 2,4,6-Triamino-1,3,5-trinitrobenzene (TATB) detonation using the quantum-based molecular dynamics method. Various carbon structures such as, graphene-like, diamond-like, and "diaphite", are obtained under different pressures. The transition from a C sp2- to a sp3-hybrid, driven by the conversion of a hexatomic to a non-hexatomic ring, is detected under high pressure. A tightly bound nucleation mechanism for diamond-like carbon dominated by a graphene-like carbon layer is uncovered. The graphene-like layer is readily constructed at the early stage, which would connect with surrounding carbon atoms or fragments to form the tetrahedral structure, with a high fraction of sp3-hybridized carbon. After that, the deformed carbon layers further coalesce with each other by bonding between carbon atoms within the five-member ring, to form the diamond-like nucleus. The complex "diaphite" configuration is detected during the diamond-like carbon nucleation, which illustrates that the nucleation and growth of detonation nano-diamond would accompany the intergrowth of graphene-like layers.

Raman spectra and vibrational properties of FOX-7 under pressure and temperature: First-principles calculations

FOX-7 (1,1-diamino-2,2-dinitroethene) as one of the widely studied insensitive high explosives exists five polymorphs (α, β, γ, α', ε) whose crystal structures have been determined by XRD (X-rays Diffraction) and which are investigated by a density functional theory (DFT) approach in this work. The calculation results show that the GGA PBE-D2 method can reproduce the experimental crystal structure of FOX-7 polymorphs better. The calculated Raman spectra of FOX-7 polymorphs were compared in detail and fully with the experimental Raman spectra data and it was found that the calculated Raman spectra frequencies have an overall red-shift in middle band (800-1700 cm^-1), and that the maximum deviation does not exceed 4 % (The maximum point is the mode of CC in plane bending). The high-temperature phase transform path (α → β → γ) and the high-pressure phase transform path (α → α'→ε) can be well represented in the computational Raman spectra. In addition, crystal structure of ε-FOX-7 was performed up to 70 GPa to probe Raman spectra and vibrational properties. The results showed that the NH₂ Raman shift is jittering with pressure (not smooth compared to other vibrational modes) and NH₂ anti-symmetry-stretching appears red-shifted. The vibration of hydrogen mixes in all of other vibrational modes. This work shows that the dispersion-corrected GGA PBE method can reproduce the experimental structure, vibrational properties and Raman spectra very well.

Enhanced Effect of an External Electric Field on NH3BH3 Dehydrogenation: an AIMD Study for Thermolysis

How to improve the dehydrogenation properties of ammonia borane (AB, NH₃BH₃) is always a challenge for its practical application in hydrogen storage. In this study, we reveal the enhanced effect of an external electric field (E _ext) on AB dehydrogenation by means of the ab initio molecular dynamics method. The molecular rotation induced by an electrostatic force can facilitate the formation of the H-N···B-H framework, which would aggregate into poly-BN species and further suppress the generation of the volatile byproducts. Meanwhile, the dihydrogen bond (N-H^δ+···^δ-H-B) is favorably formed under E _ext, and the interaction between relevant H atoms is enhanced, leading to a faster H₂ liberation. Correspondingly, the apparent activation energy for AB dissociation is greatly reduced from 18.42 to around 15 kcal·mol^-1 with the application of an electric field, while that for H₂ formation decreases from 20.4 to about 16 kcal·mol^-1. In the whole process, the cleavage of the B-H bond is more favorable than that of the N-H bond, no matter whether the application of E _ext. Our results give a deep insight into a positive effect of an electric field on AB dehydrogenation, which would provide an important inspiration for hydrogen storage in industry applications.

Anisotropic Reaction Properties for Different HMX/HTPB Composites: A Theoretical Study of Shock Decomposition

Plastic-bonded explosives (PBXs) consisting of explosive grains and a polymer binder are commonly synthesized to improve mechanical properties and reduce sensitivity, but their intrinsic chemical behaviors while subjected to stress are not sufficiently understood yet. Here, we construct three composites of β-HMX bonded with the HTPB binder to investigate the reaction characteristics under shock loading using the quantum-based molecular dynamics method. Six typical interactions between HMX and HTPB molecules are detected when the system is subjected to pressure. Although the initial electron structure is modified by the impurity states from HTPB, the metallization process for HMX does not significantly change. The shock decompositions of HMX/HTPB along the (100) and (010) surface are initiated by molecular ring dissociation and hydrogen transfer. The initial oxidations of C and H within HTPB possess advantages. As for the (001) surface, the dissociation is started with alkyl dehydrogenation oxidation, and a stronger hydrogen transfer from HTPB to HMX is detected during the following process. Furthermore, considerable fragment aggregation is observed, which mainly derives from the formation of new C-C and C-N bonds under high pressure. The effect of cluster evolution on the progression of the following reaction is further studied by analyzing the bonded structure and displacement rate.

Photolysis of the Insensitive Explosive 1,3,5-Triamino-2,4,6-trinitrobenzene (TATB)

The explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is of particular interest due to its extreme insensitivity to impact, shock and heat, while providing a good detonation velocity. To determine its fate under environmental conditions, TATB powder was irradiated with simulated sunlight and, in water, under UV light at 254 nm. The hydrolysis of particles submerged in neutral and alkaline solutions was also examined. We found that, by changing experimental conditions (e.g., light source, and mass and physical state of TATB), the intermediates and final products were slightly different. Mono-benzofurazan was the major transformation product in both irradiation systems. Two minor transformation products, the aci-nitro form of TATB and 3,5-diamino-2,4,6-trinitrophenol, were detected under solar light, while 1,3,5-triamino-2-nitroso-4,6-dinitrobenzene, 1,3,5-triamino-2,4-dinitrobenzene and mono-benzofuroxan were produced under UV light. The product identified as 3,5-diamino-2,4,6-trinitrophenol was identical to the one formed in the dark under alkaline conditions (pH 13) and in water incubated at either 50 °C or aged at ambient conditions. Interestingly, when only a few milligrams of TATB were irradiated with simulated sunlight, the aci-isomer and mono-benzofurazan derivative were detected; however, the hydrolysis product 3,5-diamino-2,4,6-trinitrophenol formed only much later in the absence of light. This suggests that the water released from TATB to form mono-benzofurazan was trapped in the interstitial space between the TATB layers and slowly hydrolyzed the relatively stable aci-nitro intermediate to 3,5-diamino-2,4,6-trinitrophenol. This environmentally relevant discovery provides data on the fate of TATB in surface environments exposed to sunlight, which can transform the insoluble substrate into more soluble and corrosive derivatives, such as 3,5-diamino-2,4,6-trinitrophenol, and that some hydrolytic transformation can continue even without light.

A quantum-based molecular dynamics study of the ICM-102/HNO3 host-guest reaction at high temperatures

The contradiction between energy and safety of explosives is better balanced by the host-guest inclusion strategy. Understanding the reaction mechanism of the host-guest explosive is necessary. To deeply analyze the role of the small guest molecules in the host-guest system, a quantum-based molecular dynamics method was used to calculate the initial decomposition reaction of the new host-guest explosive ICM-102/HNO3 against the pure ICM-102 at several high temperatures. The incorporation of HNO3 had no significant influence on the initial decomposition step of ICM-102. Conversely, the earliest intramolecular hydrogen transfer reaction is delayed partly because the H and O atoms of HNO3 connect with the O and H atoms of ICM-102, respectively. As the reaction proceeds, guest molecules get heavily involved in the reaction and increase the reaction rate. The generation rate and quantity of the small oxidizing molecules in the final product were increased significantly in the ICM-102/HNO3 system. These mechanisms revealed that HNO3 molecules inhibit the early stages of the initial decomposition of ICM-102 to some extent, and play an important role in accelerating the decomposition subsequently.