Shock waves in high-energy materials: the initial chemical events in nitramine RDX

Mechanistic Insights into the Evolution of Graphitic Carbon from Sulfur Containing Polar Aromatic Feedstock

In the realm of carbon fiber research, a variety of structural configurations is noted, comprising crystalline, noncrystalline, and semicrystalline forms. Recent investigations into this domain have revealed an array of intriguing phases of carbon, among which amorphous graphite is the most notable for its unique mechanical, thermal, and electrical properties that arise from its inherent topological disorders. In this study, we utilized the ReaxFF molecular dynamics (MD) simulations to investigate the carbonization and graphitization processes involved in the production of amorphous graphite from benzothiophene, a sulfur-containing polar aromatic precursor. We developed C/H/S ReaxFF force field parameters to describe the high-temperature chemistry of benzothiophene. Our investigation reveals the reaction mechanisms, providing critical insights into the underlying chemical processes toward the formation of amorphous graphite and the structural characteristics of the end products. The formation of volatile gaseous molecules and their continuous elimination led to the development of noncontinuous layered graphite structures analogous to amorphous graphite consisting of pentagons, hexagons, and heptagons. These findings offer unprecedented insights into the carbonization and graphitization processes of sulfur-containing heavy-end aromatic feedstock. This knowledge lays the groundwork for advancing synthesis methods and developing amorphous graphite materials with specific properties.

Implementing reactivity in molecular dynamics simulations with harmonic force fields

The simulation of chemical reactions and mechanical properties including failure from atoms to the micrometer scale remains a longstanding challenge in chemistry and materials science. Bottlenecks include computational feasibility, reliability, and cost. We introduce a method for reactive molecular dynamics simulations using a clean replacement of non-reactive classical harmonic bond potentials with reactive, energy-conserving Morse potentials, called the Reactive INTERFACE Force Field (IFF-R). IFF-R is compatible with force fields for organic and inorganic compounds such as IFF, CHARMM, PCFF, OPLS-AA, and AMBER. Bond dissociation is enabled by three interpretable Morse parameters per bond type and zero energy upon disconnect. Use cases for bond breaking in molecules, failure of polymers, carbon nanostructures, proteins, composite materials, and metals are shown. The simulation of bond forming reactions is included via template-based methods. IFF-R maintains the accuracy of the corresponding non-reactive force fields and is about 30 times faster than prior reactive simulation methods.

Size-dependent shock response mechanisms in nanogranular RDX: a reactive molecular dynamics study

Understanding the shock initiation mechanisms of explosives is pivotal for advancing physicochemical theories and enhancing experimental methodologies. This study delves into the size-dependent shock responses of nanogranular hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) through nonequilibrium reactive molecular dynamics simulations. Utilizing the ReaxFF-lg force field, we examine the influence of the particle size on the decomposition dynamics of RDX under varying shock velocities. Our findings reveal that larger particles promote more significant RDX decomposition at lower velocities due to fluid jet formation and gas compression during void collapse. Conversely, smaller particles exhibit a higher average temperature and a faster decomposition rate under high-velocity shocks, attributed to their increased specific surface area. Detailed chemical reaction pathways are analyzed to elucidate the growth and initiation of reactions during shock waves. The results contribute to resolving the discrepancies observed in experimental studies of shocked granular explosives and provide a deeper understanding of the underlying mechanisms governing their behavior. This research offers valuable insights into the design and control of nano- and submicron-sized explosives with tailored sensitivity to external stimuli.

Corrosion-Resistant Ultrathin Cu Film Deposited on N-Doped Amorphous Carbon Film Substrate and Its Use for Crumpleable Circuit Board

Copper (Cu) is widely used as an industrial electrode due to its high electrical conductivity, mechanical properties, and cost-effectiveness. However, Cu is susceptible to corrosion, which degrades device performance over time. Although various methods (alloying, physical passivation, surface treatment, etc.) are introduced to address the corrosion issue, they can cause decreased conductivity or vertical insulation. Here, using the nitrogen-doped amorphous carbon (a-C:N) thin film is proposed as a substrate on which Cu is directly deposited. This simple method significantly inhibits corrosion of ultrathin Cu (<20 nm) films in humid conditions, enabling the fabrication of ultrathin electronic circuit boards without corrosion under ambient conditions. This study investigates the origin of corrosion resistance through comprehensive microscopic/spectroscopic characterizations and density-functional theory (DFT) calculations: i) diffusion of Cu atoms into the a-C:N driven by stable C-Cu-N bond formation, ii) diffusion of N atoms from the a-C:N to the Cu layer heading the top surface, which is the thermodynamically preferred location for N, and iii) the doped N atoms in Cu layer suppress the inclusion of O into the Cu lattice. By leveraging the ultrathinness and deformability of the circuit board, a transparent electrode and a crumpleable LED lighting device are demonstrated.

Curing and Cross-Linking Processes in the Poly(3,3-bis-azidomethyl oxetane)-tetrahydrofuran/Toluene Diisocyanate/Trimethylolpropane System: A Density Functional Theory and Accelerated ReaxFF Molecular Dynamics Investigation

The physical and chemical properties of solid propellant are influenced by the composition and structure of the binder, with its network structure being formed through curing and cross-linking reactions. Therefore, understanding the mechanisms of these reactions is crucial. In this study, we investigated the curing and cross-linking mechanisms of poly(3,3-bis-azidomethyl oxetane)-tetrahydrofuran (PBT), toluene diisocyanate (TDI), and trimethylolpropane (TMP) using a combination of density functional theory (DFT) calculations and accelerated ReaxFF molecular dynamics (MD) simulations. DFT calculations revealed that the steric effect of the -CH3 group in TDI exerts a significant influence on the curing reaction between TDI and PBT. Additionally, in the cross-linking process, the energy barrier for TDI reacting with TMP was found to be much lower than that for TDI reacting with the PBT-TDI intermediate. Subsequently, we conducted competing reaction processes of TMP/TDI-PBT-TDI cross-linking and TDI-PBT-TDI self-cross-linking using accelerated MD simulations within the fitted ReaxFF framework. The results showed that the successful frequency of TMP/TDI-PBT-TDI cross-linking was substantially higher than that of TDI-PBT-TDI self-cross-linking, consistent with the energy barrier results from DFT calculations. These findings deepen our understanding of the curing and cross-linking mechanisms of the PBT system, providing valuable insights for the optimization and design of solid propellants.

Molecular Dynamics Simulation of Thermal Nonequilibrium and Chemical Reaction Processes in Hydrogen Combustion

Using reactive force field (ReaxFF) and molecular dynamics simulation, we investigate the combustion process of hydrogen-oxygen systems in initial thermal nonequilibrium states with different translational and rovibrational temperatures for oxygen. The system studied in this work contains 300 oxygen molecules and 700 hydrogen molecules with a density of 7 times the air density. For this system, the characteristic relaxation times of oxygen and hydrogen vibrational energies are 0.173 and 0.249 ns, respectively. 0.6% of hydrogen undergoes a chemical reaction with oxygen during the thermal nonequilibrium relaxation stage. For the distribution of translational energy and vibrational energy of oxygen in the thermal nonequilibrium state, the maximum mean error of the statistical distribution in the simulation and the Boltzmann distribution at temperature calculated from the average kinetic energy of molecules is about 2.25 × 10-5. At the same time, it was observed in the simulation that many-body interactions play a certain role in the combustion process. Furthermore, we compare the ignition time and temperature rise behavior of different combustion mechanisms and molecular dynamics simulations starting from the thermal equilibrium state. These results will provide meaningful references for the construction of thermal nonequilibrium combustion chemical reaction mechanisms.

Enhancing ReaxFF for molecular dynamics simulations of lithium-ion batteries: an interactive reparameterization protocol

Lithium-ion batteries (LIBs) have become an essential technology for the green economy transition, as they are widely used in portable electronics, electric vehicles, and renewable energy systems. The solid-electrolyte interphase (SEI) is a key component for the correct operation, performance, and safety of LIBs. The SEI arises from the initial thermal metastability of the anode-electrolyte interface, and the resulting electrolyte reduction products stabilize the interface by forming an electrochemical buffer window. This article aims to make a first-but important-step towards enhancing the parametrization of a widely-used reactive force field (ReaxFF) for accurate molecular dynamics (MD) simulations of SEI components in LIBs. To this end, we focus on Lithium Fluoride (LiF), an inorganic salt of great interest due to its beneficial properties in the passivation layer. The protocol relies heavily on various Python libraries designed to work with atomistic simulations allowing robust automation of all the reparameterization steps. The proposed set of configurations, and the resulting dataset, allow the new ReaxFF to recover the solid nature of the inorganic salt and improve the mass transport properties prediction from MD simulation. The optimized ReaxFF surpasses the previously available force field by accurately adjusting the diffusivity of lithium in the solid lattice, resulting in a two-order-of-magnitude improvement in its prediction at room temperature. However, our comprehensive investigation of the simulation shows the strong sensitivity of the ReaxFF to the training set, making its ability to interpolate the potential energy surface challenging. Consequently, the current formulation of ReaxFF can be effectively employed to model specific and well-defined phenomena by utilizing the proposed interactive reparameterization protocol to construct the dataset. Overall, this work represents a significant initial step towards refining ReaxFF for precise reactive MD simulations, shedding light on the challenges and limitations of ReaxFF force field parametrization. The demonstrated limitations emphasize the potential for developing more versatile and advanced force fields to upscale ab initio simulation through our interactive reparameterization protocol, enabling more accurate and comprehensive MD simulations in the future.

Significant thermal rectification induced by phonon mismatch of functional groups in a single-molecule junction

Single-molecule junctions (SMJs) may bring exotic physical effects. In this work, a significant thermal rectification effect is observed in a cross-dimensional system, comprising a diamond, a single-molecule junction, and a carbon nanotube (CNT). The molecular dynamics simulations indicate that the interfacial thermal resistance varies with the direction of heat flow, the orientation of the crystal planes of the diamond, and the length of the CNT. We find that the thermal rectification ratio escalates with the length of the CNT, achieving a peak value of 730% with the CNT length of 200 nm. A detailed analysis of phonon vibrations suggests that the primary cause of thermal rectification is the mismatched vibrations between the biphenyl and carbonyl groups. This discovery may offer theoretical insights for both the experimental exploration and practical application of SMJs in efficient thermal management strategy for high power and highly integrated chips.

Effects of Nanoparticle Size on the Thermal Decomposition Mechanisms of 3,5-Diamino-6-hydroxy-2-oxide-4-nitropyrimidone through ReaxFF Large-Scale Molecular Dynamics Simulations

ReaxFF-lg molecular dynamics method was employed to simulate the decomposition processes of IHEM-1 nanoparticles at high temperatures. The findings indicate that the initial decomposition paths of the nanoparticles with different sizes at varying temperatures are similar, where the bimolecular polymerization reaction occurred first. Particle size has little effect on the initial decomposition pathway, whereas there are differences in the numbers of the species during the decomposition and their evolution trends. The formation of the hydroxyl radicals is the dominant decomposition mechanism with the highest reaction frequency. The degradation rate of the IHEM-1 molecules gradually increases with the increasing temperature. The IHEM-1 nanoparticles with smaller sizes exhibit greater decomposition rate constants. The activation energies for the decomposition are lower than the reported experimental values of bulk explosives, which suggests a higher sensitivity.

Multimodal E-Textile Enabled by One-Step Maskless Patterning of Femtosecond-Laser-Induced Graphene on Nonwoven, Knit, and Woven Textiles

Personal wearable devices are considered important in advanced healthcare, military, and sports applications. Among them, e-textiles are the best candidates because of their intrinsic conformability without any additional device installation. However, e-textile manufacturing to date has a high process complexity and low design flexibility. Here, we report the direct laser writing of e-textiles by converting raw Kevlar textiles to electrically conductive laser-induced graphene (LIG) via femtosecond laser pulses in ambient air. The resulting LIG has high electrical conductivity and chemical reliability with a low sheet resistance of 2.86 Ω/□. Wearable multimodal e-textile sensors and supercapacitors are realized on different types of Kevlar textiles, including nonwoven, knit, and woven structures, by considering their structural textile characteristics. The nonwoven textile exhibits high mechanical stability, making it suitable for applications in temperature sensors and micro-supercapacitors. On the other hand, the knit textile possesses inherent spring-like stretchability, enabling its use in the fabrication of strain sensors for human motion detection. Additionally, the woven textile offers special sensitive pressure-sensing networks between the warp and weft parts, making it suitable for the fabrication of bending sensors used in detecting human voices. This direct laser synthesis of arbitrarily patterned LIGs from various textile structures could result in the facile realization of wearable electronic sensors and energy storage.

Reactive molecular dynamics simulations on the decomposition process of 1,3,5-trinitro-1,3,5-triazine crystal under high temperatures and pressure

CONTEXT

Reactive molecular dynamics simulations were performed to study the decomposition processes of 1,3,5-trinitro-1,3,5-triazine (RDX) crystal under high temperatures (2100, 2400, 2700, and 3000 K) and detonation pressure (34.5 GPa) and 0 GPa. It is found that the initial decomposition paths of RDX under different temperatures coupled with detonation pressure are similar, which is due to the N-NO2 bond breakage to release NO2. The formation rates of N2 and H2O are significantly affected by temperature, while those of CO2 are less influenced. The C atoms finally formed C clusters. As the temperature rises, the decomposition speeds up, indicating that the high temperature accelerates the decomposition. Applying pressure can reduce the reaction energy barrier and accelerate the decomposition.

METHODS

The RDX model was constructed using the Materials Studio 7.0 package. All MD simulations were performed based on the ReaxFF force field in the LAMMPS software package, and the crystals were visualized using the OVITO software package. The time step was 0.1 fs, and the total MD simulation time was 200 ps. DFT calculations were carried out at the B3LYP/6-311G(d,p) level using the Gaussian 09 package.

Deep Potential Molecular Dynamics Study of Chapman-Jouguet Detonation Events of Energetic Materials

Detonation of energetic materials (EMs) is of great importance for military applications, while the understanding of detailed events and mechanisms for the detonation process is scarce. In this study, the first deep neural network potential NNP_Shock for molecular dynamics (MD) simulation of shock-induced detonation of EMs was generated based on a deep potential model, providing DFT accuracy but 106 times the computational efficiency. On this basis, we employ our deep potential to perform MD simulations of shock-induced detonation of high-performance EM material 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20, C6H6N12O12) and obtain the theoretical Chapman-Jouguet (C-J) detonation velocities and pressures directly by multiscale shock technique (MSST) for the first time, which are in good agreement with experiment. In addition, the Hugoniot curves and initial reaction mechanisms were successfully obtained. Therefore, the NNP_Shock potential is competent in research of the detonation performance and shock sensitivity of CL-20, and the method can also be transplanted to studies of other EMs.

Charge equilibration model with shielded long-range Coulomb for reactive molecular dynamics simulations

Atomic description of electrochemical systems requires reactive interaction potential to explicitly describe the chemistry between atoms and molecules and the evolving charge distribution and polarization effects. Calculating Coulomb electrostatic interactions and polarization effects requires a better estimate of the partial charge distribution in molecular systems. However, models such as reactive force fields and charge equilibration (QEq) include Coulomb interactions up to a short-distance cutoff for better computational speeds. Ignoring long-distance electrostatic interaction affects the ability to describe electrochemistry in large systems. We studied the long-range Coulomb effects among charged particles and extended the QEq method to include long-range effects. By this extension, we anticipate a proper account of Coulomb interactions in reactive molecular dynamics simulations. We validate the approach by computing charges on a series of metal-organic frameworks and some simple systems. Results are compared to regular QEq and quantum mechanics calculations. The study shows slightly overestimated charge values in the regular QEq approach. Moreover, our method was combined with Ewald summation to compute forces and evaluate the long-range effects of simple capacitor configurations. There were noticeable differences between the calculated charges with/without long-range Coulomb interactions. The difference, which may have originated from the long-range influence on the capacitor ions, makes the Ewald method a better descriptor of Coulomb electrostatics for charged electrodes. The approach explored in this study enabled the atomic description of electrochemical systems with realistic electrolyte thickness while accounting for the electrostatic effects of charged electrodes throughout the dielectric layer in devices like batteries and emerging solid-state memory.

The Crystallization of Disordered Materials under Shock Is Governed by Their Network Topology

When the shock load is applied, materials experience incredibly high temperature and pressure conditions on picosecond timescales, usually accompanied by remarkable physical or chemical phenomena. Understanding the underlying physics that governs the kinetics of shocked materials is of great importance for both physics and materials science. Here, combining experiment and large-scale molecular dynamics simulation, the ultrafast nanoscale crystal nucleation process in shocked soda-lime silicate glass is investigated. By adopting topological constraints theory, this study finds that the propensity of nucleation is governed by the connectivity of the atomic network. The densification of local networks, which appears once the crystal starts to grow, results in the underconstrained shell around the crystal and prevents further crystallization. These results shed light on the nanoscale crystallization mechanism of shocked materials from the viewpoint of topological constraint theory.

High-pressure and temperature neural network reactive force field for energetic materials

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive compared with electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network-based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical responses of energetic materials at extreme conditions. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and shock-loading conditions. This new potential shows improved accuracy over the current state-of-the-art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock-loading conditions.

Initial Response of Pentaerythritol Tetranitrate (PETN) under the Coupling Effect of Preheating, Shock and Defect via the Molecular Dynamics Simulations with the Multiscale Shock Technique Method

The initial response of PETN under the coupling of preheating, impact and defects was simulated by Multiscale Shock Technique (MSST) method and molecular dynamics. The temperature change of PETN during impact compression can be divided into three stages: (1) the elastoplastic change of the system caused by initial compression; (2) part of PETN decomposes and releases energy to raise temperature; (3) a secondary chemical reaction occurs, resulting in rapid temperature rise. Under the given conditions, a higher initial preheating temperature will lead to faster decomposition of PETN; The existence of defects will accelerate the decomposition of PETN molecules; Coupling the highest preheating temperature with defects will lead to the fastest decomposition of PETN molecules, while in the defect-free PETN system with a preheating temperature of 300 K, the decomposition of PETN molecules is the slowest. For the case of U_s = 8 km·s^-1, the effect of defects on the initial PETN reaction is greater than the initial preheating temperature; When the impact velocity is greater than 9 km·s^-1, the impact velocity is an important factor affecting the decomposition of PETN molecules. For U_s = 10 km·s^-1, NO₂ is the main initial product in the defective PETN crystal, while in the perfect PETN crystal, it is the combination of NO₂ and HONO. The chemical reaction kinetics analysis shows that the preheating temperature and defects will accelerate the decomposition of PETN. The higher the preheating temperature, the faster the decomposition of PETN. For the case of U_s = 7 km·s^-1, 8 km·s^-1 and 9 km·s^-1, the existence of defects will increase the decomposition rate by more than 50% regardless of the initial preheating temperature. In the case of Us = 10 km·s^-1, the improvement of decomposition rate by defects is not as significant as the initial preheating temperature.

Decomposition mechanism of 1,3,5-trinitro-2,4,6-trinitroaminobenzene under thermal and shock stimuli using ReaxFF molecular dynamics simulations

To obtain atomic-level insights into the decomposition behavior of 1,3,5-trinitro-2,4,6-trinitroaminobenzene (TNTNB) under different stimulations, this study applied reactive molecular dynamics simulations to illustrate the effects of thermal and shock stimuli on the TNTNB crystal. The results show that the initial decomposition of the TNTNB crystal under both thermal and shock stimuli starts with the breakage of the N-NO₂ bond. However, the C₆ ring in TNTNB undergoes structural rearrangement to form a C₃-C₅ bicyclic structure at a constant high temperature. Then, the C₃ and C₅ rings break in turn. The main final products of TNTNB under shock are N₂, CO₂, and H₂O, while NO,  N_2, H₂O and CO are formed instead at 1 atm under a constant high temperature. Pressure is the main reason for this difference. High pressure promotes the complete oxidation of the reactants.

Effects of Different Guests on Pyrolysis Mechanism of α-CL-20/Guest at High Temperatures by Reactive Molecular Dynamics Simulations at High Temperatures

The host-guest inclusion strategy has the potential to surpass the limitations of energy density and suboptimal performances of single explosives. The guest molecules can not only enhance the detonation performance of host explosives but also can enhance their stability. Therefore, a deep analysis of the role of guest influence on the pyrolysis decomposition of the host-guest explosive is necessary. The whole decomposition reaction stage of CL-20/H₂O, CL-20/CO₂, CL-20/N₂O, CL-20/NH₂OH was calculated by ReaxFF-MD. The incorporation of CO₂, N₂O and NH₂OH significantly increase the energy levels of CL-20. However, different guests have little influence on the initial decomposition paths of CL-20. The E_a1 and E_a2 values of CL-20/CO₂, CL-20/N₂O, CL-20/NH₂OH systems are higher than the CL-20/H₂O system. Clearly, incorporation of CO₂, N₂O, NH₂OH can inhibit the initial decomposition and intermediate decomposition stage of CL-20/H₂O. Guest molecules become heavily involved in the reaction and influence on the reaction rates. k₁ of CL-20/N₂O and CL-20/NH₂OH systems are significantly larger than that of CL-20/H₂O at high temperatures. k₁ of CL-20/CO₂ system is very complex, which can be affected deeply by temperatures. k₂ of the CL-20/CO₂, CL-20/N₂O systems is significantly smaller than that of CL-20/H₂O at high temperatures. k₂ of CL-20/NH₂OH system shows little difference at high temperatures. For the CL-20/CO₂ system, the k₃ value of CO₂ is slightly higher than that for CL-20/H₂O, CL-20/N₂O, CL-20/NH₂OH systems, while the k₃ values of N₂ and H₂O are slightly smaller than that for the CL-20/H₂O, CL-20/N₂O, CL-20/NH₂OH systems. For the CL-20/N₂O system, the k₃ value of CO₂ is slightly smaller than that for CL-20/H₂O, CL-20/CO₂, CL-20/NH₂OH systems. For the CL-20/NH₂OH system, the k₃ value of H₂O is slightly larger than that for CL-20/H₂O, CL-20/CO₂, CL-20/N₂O systems. These mechanisms revealed that CO₂, N₂O and NH₂OH molecules inhibit the early stages of the initial decomposition of CL-20 and play an important role for the decomposition subsequently.

A coarse-grain reactive model of RDX: Molecular resolution at the μm scale

Predictive models for the thermal, chemical, and mechanical response of high explosives at extreme conditions are important for investigating their performance and safety. We introduce a particle-based, reactive model of 1,3,5-trinitro-1,3,5-triazinane (RDX) with molecular resolution utilizing generalized energy-conserving dissipative particle dynamics with reactions. The model is parameterized with respect to the data from atomistic molecular dynamics simulations as well as from quantum mechanical calculations, thus bridging atomic processes to the mesoscales, including microstructures and defects. It accurately captures the response of RDX under a range of thermal loading conditions compared to atomistic simulations. In addition, the Hugoniot response of the CG model in the overdriven regime reasonably matches atomistic simulations and experiments. Exploiting the model's high computational efficiency, we investigate mesoscale systems involving millions of molecules and characterize size-dependent criticality of hotspots in RDX. The combination of accuracy and computational efficiency of our reactive model provides a tool for investigation of mesoscale phenomena, such as the role of microstructures and defects in the shock-to-deflagration transition, through particle-based simulation.

Evaluating the performance of ReaxFF potentials for sp2 carbon systems (graphene, carbon nanotubes, fullerenes) and a new ReaxFF potential

We study the performance of eleven reactive force fields (ReaxFF), which can be used to study sp2 carbon systems. Among them a new hybrid ReaxFF is proposed combining two others and introducing two different types of C atoms. The advantages of that potential are discussed. We analyze the behavior of ReaxFFs with respect to 1) the structural and mechanical properties of graphene, its response to strain and phonon dispersion relation; 2) the energetics of (n, 0) and (n, n) carbon nanotubes (CNTs), their mechanical properties and response to strain up to fracture; 3) the energetics of the icosahedral C60 fullerene and the 40 C40 fullerene isomers. Seven of them provide not very realistic predictions for graphene, which made us focusing on the remaining, which provide reasonable results for 1) the structure, energy and phonon band structure of graphene, 2) the energetics of CNTs versus their diameter and 3) the energy of C60 and the trend of the energy of the C40 fullerene isomers versus their pentagon adjacencies, in accordance with density functional theory (DFT) calculations and/or experimental data. Moreover, the predicted fracture strain, ultimate tensile strength and strain values of CNTs are inside the range of experimental values, although overestimated with respect to DFT. However, they underestimate the Young's modulus, overestimate the Poisson's ratio of both graphene and CNTs and they display anomalous behavior of the stress - strain and Poisson's ratio - strain curves, whose origin needs further investigation.

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

Regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

Ab initio neural network MD simulation of thermal decomposition of a high energy material CL-20/TNT

CL-20 (2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane, also known as HNIW) is one of the most powerful energetic materials. However, its high sensitivity to environmental stimuli greatly reduces its safety and severely limits its application. In this work, ab initio based neural network potential (NNP) energy surfaces for both β-CL-20 and CL-20/TNT co-crystals were constructed. To accurately simulate the thermal decomposition processes of these two crystal systems, reactive molecular dynamics simulations based on the NNPs were performed. Many important intermediate species and their associated reaction paths during the decomposition had been identified in the simulations and the direct results on detonation temperatures of both systems were provided. The simulations also showed clearly that 2,4,6-trinitrotoluene (TNT) molecules in the co-crystal act as a buffer to slow down the chain reactions triggered by nitrogen dioxide and this effect is more significant at lower temperatures. Specifically, the addition of TNT molecules in the CL-20/TNT co-crystal introduces intermolecular hydrogen bonds between CL-20 and TNT molecules in the system, thereby increasing the thermal stability of the co-crystal. The current reactive molecular dynamics simulation is performed based on the NNP which helps in accelerating the speed of ab initio molecular dynamics (AIMD) simulation by more than 3 orders of magnitude while preserving the accuracy of density functional theory (DFT) calculations. This enabled us to perform longer-time simulations at more realistic temperatures that traditional AIMD methods cannot achieve. With the advantage of the NNP in its powerful fitting ability and transferability, the NNP-based MD simulation can be widely applied to energetic material systems.

Mechanically Robust Irradiation, Atomic Oxygen, and Static-Durable CrOx/CuNi Coatings on Kapton Serving as Space Station Solar Cell Arrays

The polymers that served for solar cell arrays are constantly subject to various hazards, such as atomic oxygen (AO), ion irradiation, or electrostatic discharge (ESD) events. To address these issues, we fabricated and sifted CrO_0.16/CuNi-coated Kapton with a gradient structure with the goal of reaching an equilibrium between AO durability and resistance. The resulting material exhibits an impressively low E_y of 6.61 × 10^-26 cm³ atom^-1, 2.20% of which was detected as pristine Kapton. Self-evolution of the CrO_0.16 coating under 525.4 displacement per atom (dpa) Fe⁺ ion irradiation indicated that it can still maintain a good state of ultrafine nanocrystalline in addition to local amorphization. Its AO-based degradation and irradiation evolution are demonstrated by molecular dynamics (MD) simulations. It is mechanically robust enough to endure the cyclic folding treatments attributed to its gradient structure fabrication. Moreover, the CrO_0.16/CuNi-coated Kapton exhibits alleviated electrostatic accumulation capability and sufficient conductivity. Our strategy has promising potential for creating surface protection on flexible polymers operating in the low Earth orbit (LEO).

Reactive Molecular Dynamics Calculation and Ignition Delay Test of the Mixture of an Additive and 2-Azido-N,N-dimethylethanamine with Dinitrogen Tetroxide

In order to shorten the ignition delay of 2-azido-N,N-dimethylethanamine (DMAZ) and dinitrogen tetroxide (NTO), four amines [tert-butylamine, pyrrole, N,N,N',N'-tetramethyl ethylenediamine (TMEDA), and diethylenetriamine (DABH)] with a mass fraction of 5% were added to DMAZ, and the potential energy change and the product change during the reaction of the mixture of an additive and DMAZ with NTO were analyzed by Reactive molecular dynamics (ReaxFF MD) calculation. Then, the ignition delay of the mixture of the additive and DMAZ as well as pure DMAZ with NTO was measured by a drop experiment with a photoelectric sensor and high-speed camera. The results show that the addition of pyrrole greatly reduced the time to reach the maximum system energy and greatly increased the rate of HNO₂ formation. The dripping of the fuel was approximately a uniform linear motion, and the expression was y = 43.13 + 7.16x. The ignition delay time recorded by the camera was in good agreement with that of the optical signal. The measured ignition delay time for DMAZ with NTO was 261.5 ms. The mixture of pyrrole and DMAZ with NTO had the shortest ignition delay time of 100 ms, and the proportion of shortening the ignition delay time was the largest. The results of the droplet experiment were consistent with those of ReaxFF MD simulation, indicating that HNO₂ plays an important role in the ignition delay, that is, the formation rate of HNO₂ is positively correlated with the ignition delay.

Hot spot formation and initial chemical reaction of PETN containing nanoscale spherical voids under high shock loading

The initial response process of PETN containing nanoscale spherical cavities under impact loading was investigated using the ReaxFF-lg force field combined with the molecular dynamic method. The impact-induced void collapse process, hot spot formation and growth, and chemical reaction processes were determined. The hot spot formation goes through four stages: (1) overall temperature rise due to initial impact compression; (2) temperature rise on the upper surface of the void caused by local plastic deformation; (3) rapid temperature rise caused by molecules entering the interior of the void colliding with the downstream surface of the void; and (4) thermal diffusion between the hot spot and the surrounding region, resulting in a decrease in the temperature of the center of the hot spot and a slow increase in the temperature of the neighboring region. With weak impact, the void shape remains basically symmetric during the void collapse, and the void collapse is mainly caused by local plastic deformation. A strong impact will lead to a more intense material focusing. The void collapse caused by strong impact has a greater effect on the heating of the surrounding material, and the secondary compression formed by the collision between particles makes the hot spot area expand and the central region of the hot spot evolve into an approximate triangular cone. NO2 is produced in large quantities as the initial product during the void collapse to form the hot spot, indicating that the void activates the chemical reactivity of the PETN crystal.

Molecular Forcefield Methods for Describing Energetic Molecular Crystals: A Review

Energetic molecular crystals are widely applied for military and civilian purposes, and molecular forcefields (FF) are indispensable for treating the microscopic issues therein. This article reviews the three types of molecular FFs that are applied widely for describing energetic crystals-classic FFs, consistent FFs, and reactive FFs (ReaxFF). The basic principle of each type of FF is briefed and compared, with the application introduced, predicting polymorph, morphology, thermodynamics, vibration spectra, thermal property, mechanics, and reactivity. Finally, the advantages and disadvantages of these FFs are summarized, and some directions of future development are suggested.

Pressure and Polymorph Dependent Thermal Decomposition Mechanism of Molecular Materials: A Case of 1,3,5,-Trinitro-1,3,5,-triazine

1,3,5,-Trinitro-1,3,5,-triazine (RDX) serves as an important energetic material and is widely used as various solid propellants and explosives. Understanding the thermal decomposition behaviors of various polymorphs of RDX at high pressure and high temperature is significantly important for safe storage and handling. The present work reveals the early thermal decay mechanisms of two polymorphs (α- and ε-forms) of RDX at high pressure and high temperature by ReaxFF reactive molecular dynamic simulations and climbing image nudged elastic band (CI-NEB) static calculations. It is found that the thermal decomposition rate has positive and negative effects on the pressure for α- and ε-RDX, respectively. This difference originates from the difference of pressure effect on the intermolecular H transfer of the two polymorphs, as we confirm that the bimolecular H transfer rather than the NO₂ partition initiates the decay with a significantly lower energy barrier therein. This finding may be beneficial to understand the pressure and polymorph dependent effect on the decay of RDX and to develop a kinetic model for the combustion of solid RDX.

Molecular dynamic insight into octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) and the nano-HMX decomposition mechanism

The entire decomposition reaction process of a 30 Å HMX nanoparticle at 2000 K by ReaxFF molecular dynamics.

Reaction mechanism of aluminum nanoparticles in explosives under high temperature and high pressure by shock loading

Aluminum nanoparticles (ANPs) can greatly improving the power of explosives. However, the rapid reaction mechanism of ANPs under the simultaneous high temperature and high pressure by shock loading is not...

Theoretical insight into density and stability differences of RDX, HMX and CL-20

In this work, density and stability differences of RDX, HMX and CL-20 are exploited and addressed through static calculations from views of monomolecular parameters, intermolecular interactions (by the proposed BEC method) and crystal packing.

Combined ReaxFF and Ab Initio MD Simulations of Brown Coal Oxidation and Coal–Water Interactions

In this manuscript, we use a combination of Car–Parrinello molecular dynamics (CPMD) and ReaxFF reactive molecular dynamics (ReaxFF-MD) simulations to study the brown coal–water interactions and coal oxidation. Our Car–Parrinello molecular dynamics simulation results reveal that hydrogen bonds dominate the water adsorption process, and oxygen-containing functional groups such as carboxyl play an important role in the interaction between brown coal and water. The discrepancy in hydrogen bonds formation between our simulation results by ab initio molecular dynamics (CPMD) and that by ReaxFF-MD indicates that the ReaxFF force field is not capable of accurately describing the diffusive behaviors of water on lignite at low temperatures. The oxidations of brown coal for both fuel rich and fuel lean conditions at various temperatures were investigated using ReaxFF-MD simulations through which the generation rates of major products were obtained. In addition, it was observed that the density decrease significantly enhances the generation of gaseous products due to the entropy gain by reducing system density. Although the ReaxFF-MD simulation of complete coal combustion process is limited to high temperatures, the combined CPMD and ReaxFF-MD simulations allow us to examine the correlation between water adsorption on brown coal and the initial stage of coal oxidation.

Molecular Dynamics Simulations of the Thermal Decomposition of 3,4-Bis(3-nitrofurazan-4-yl)furoxan

When stimulated, for example, by a high temperature, the physical and chemical properties of energetic materials (EMs) may change, and, in turn, their overall performance is affected. Therefore, thermal stability is crucial for EMs, especially the thermal dynamic behavior. In the past decade, significant efforts have been made to study the thermal dynamic behavior of 3,4-bis(3-nitrofurazan-4-yl)furoxan (DNTF), one of the new high-energy-density materials (HEDMs). However, the thermal decomposition mechanism of DNTF is still not specific or comprehensive. In this study, the self-consistent-charge density-functional tight-binding method was combined with molecular dynamics (MD) simulations to reveal the differences in the thermal decomposition of DNTF under four heating conditions. The O-N (O) bond would fracture first during DNTF initial thermal decomposition at medium and low temperatures, thus triggering the cracking of the whole structure. At 2000 and 2500 K, NO₂ loss on outer ring I is the fastest initial thermal decomposition pathway, and it determines that the decomposition mechanism is different from that of a medium-low temperature. NO₂ is found to be the most active intermediate product; large molecular fragments, such as C₂N₂O, are found for the first time. Hopefully, these results could provide some insights into the decomposition mechanism of new HEDMs.

The Operating Cycle of NO Adsorption and Desorption in Pd-Chabazite for Passive NOx Adsorbers

Pd-doped chabazite (Pd/CHA) offers unique opportunities to adsorb and desorb NO_x in the target temperature range for application as a passive NO_x adsorber (PNA). The ability of Pd/CHA to trap NO_x emissions at low temperatures (<200 °C) is facilitated by the binding of NO_x species at various Pd sites available in the CHA framework. Density functional theory (DFT) simulations are performed to understand Pd speciation in CHA and the interaction of NO with Pd/CHA to explain the mechanisms of NO adsorption, oxidation, and desorption processes. The calculations are used to elucidate the important role of Pd¹⁺ cationic species, anchored at 6MR-3NN, in providing a strong (E_b = -272 kJ/mol) NO adsorption site in Pd/CHA. For NO release, the redox transformation of Pd species comes into play and Pd¹⁺ species are suggested to transform into cationic Pd²⁺, [PdOH]⁺, or [Pd-O-Pd]²⁺ species, all of which show significantly reduced NO binding (-116, -153, and -117 kJ/mol, respectively) as compared to Pd¹⁺. This enables NO desorption at the operating temperature of a downstream catalyst for subsequent catalytic reduction.

Reaction Mechanism and Energetics of Decomposition of Tetrakis(1,3-dimethyltetrazol-5-imidoperchloratomanganese(II)) from Quantum-Mechanics-based Reactive Dynamics

Energetic materials (EMs) are central to construction, space exploration, and defense, but over the past 100 years, their capabilities have improved only minimally as they approach the CHNO energetic ceiling, the maximum energy density possible for EMs based on molecular carbon-hydrogen-nitrogen-oxygen compounds. To breach this ceiling, we experimentally explored redox-frustrated hybrid energetic materials (RFH EMs) in which metal atoms covalently connect a strongly reducing fuel ligand (e.g., tetrazole) to a strong oxidizer (e.g., ClO₄). In this Article, we examine the reaction mechanisms involved in the thermal decomposition of an RFH EM, [Mn(Me₂TzN)(ClO₄]₄ (3, Tz = tetrazole). We use quantum-mechanical molecular reaction dynamics simulations to uncover the atomistic reaction mechanisms underlying this decomposition. We discover a novel initiation mechanism involving oxygen atom transfer from perchlorate to manganese, generating energy that promotes the fission of tetrazole into chemically stable species such as diazomethane, diazenes, triazenes, and methyl azides, which further undergo exothermic decomposition to finally form stable N₂, H₂O, CO, CO₂, Mn-based clusters, and additional incompletely combusted products.

A Review of Recent Developments in Molecular Dynamics Simulations of the Photoelectrochemical Water Splitting Process

In this review, we provide a short overview of the Molecular Dynamics (MD) method and how it can be used to model the water splitting process in photoelectrochemical hydrogen production. We cover classical non-reactive and reactive MD techniques as well as multiscale extensions combining classical MD with quantum chemical and continuum methods. Selected examples of MD investigations of various aqueous semiconductor interfaces with a special focus on TiO2 are discussed. Finally, we identify gaps in the current state-of-the-art where further developments will be needed for better utilization of MD techniques in the field of water splitting.

Mechanism of the improvement of the energy of host-guest explosives by incorporation of small guest molecules: HNO3 and H2O2 promoted C-N bond cleavage of the ring of ICM-102

Host–guest materials exhibit great potential applications as an insensitive high-energy–density explosive and low characteristic signal solid propellant. To investigate the mechanism of the improvement of the energy of host–guest explosives by guest molecules, ReaxFF-lg reactive molecular dynamics simulations were performed to calculate the thermal decomposition reactions of the host–guest explosives systems ICM-102/HNO₃, ICM-102/H₂O₂, and pure ICM-102 under different constant high temperatures and different heating rates. Incorporation of guest molecules significantly increased the energy level of the host–guest system. However, the initial reaction path of the ICM-102 molecule was not changed by the guest molecules. The guest molecules did not initially participate in the host molecule reaction. After a period of time, the H₂O₂ and HNO₃ guest molecules promoted cleavage of the C–N bond of the ICM-102 ring. Stronger oxidation and higher oxygen content resulted in the guest molecules more obviously accelerating destruction of the ICM-102 ring structure. The guest molecules accelerated the initial endothermic reaction of ICM-102, but they played a more important role in the intermediate exothermic reaction stage: incorporation of guest molecules (HNO₃ and H₂O₂) greatly improved the heat release and exothermic reaction rate. Although the energies of the host–guest systems were clearly improved by incorporation of guest molecules, the guest molecules had little effect on the thermal stabilities of the systems.

ParAMS: Parameter Optimization for Atomistic and Molecular Simulations

This work introduces ParAMS-a versatile Python package that aims to make parametrization workflows in computational chemistry and physics more accessible, transparent, and reproducible. We demonstrate how ParAMS facilitates the parameter optimization for potential energy surface (PES) models, which can otherwise be a tedious specialist task. Because of the package's modular structure, various functionality can be easily combined to implement a diversity of parameter optimization protocols. For example, the choice of PES model and the parameter optimization algorithm can be selected independently. An illustration of ParAMS' strengths is provided in two case studies: (i) a density functional-based tight binding (DFTB) repulsive potential for the inorganic ionic crystal ZnO and (ii) a ReaxFF force field for the simulation of organic disulfides.

How Small Molecules Affect the Thermo-Oxidative Aging Mechanism of Polypropylene: A Reactive Molecular Dynamics Study

Understanding the aging mechanism of polypropylene (PP) is fundamental for the fabrication and application of PP-based materials. In this paper, we present our study in which we first used reactive molecular dynamics (RMD) simulations to explore the thermo-oxidative aging of PP in the presence of acetic acid or acetone. We studied the effects of temperature and oxygen on the aging process and discussed the formation pathways of typical small molecule products (H₂, CO, CO₂, CH₄, C₂H₄, and C₂H₆). The effect of two infection agents, acetic acid and acetone, on the aging reaction was analyzed emphatically. The simulation results showed that acetone has a weak impact on accelerating the aging process, while acetic acid has a significant effect, consistent with previous experimental studies. By tracking the simulation trajectories, both acetic acid and acetone produced small active free radicals to further react with other fragment products, thus accelerating the aging process. The first reaction step of acetic acid is often the shedding of the H atom on the hydroxyl group, while the reaction of acetone is often the shedding of the H atom or the methyl. The latter requires higher energy at lower temperatures. This is why the acceleration effect of acetone for the thermo-oxidative aging of PP was not so significant compared to acetic acid in the experimental temperature (383.15 K).