Transferable Reactive Force Fields: Extensions of ReaxFF-lg to Nitromethane

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.

Reactive molecular dynamics study on the thermal decomposition reaction of a triple-base solid propellant

The study of the combustion property of newly designed propellant by means of computational simulation is an efficient pathway for assessment and could avoid exposure to hazardous chemicals. An RDX-modified triple-base solid propellant formula was proposed in this study. Reactive molecular dynamics simulations employing ReaxFF-lg force field were performed to explore the thermal decomposition property of the propellant for a variety of temperatures. The reaction kinetics of the system and major ingredients were analyzed, and the apparent decomposition activation energies were calculated. The population of decomposition intermediates and products is thoroughly investigated. H₂O is consumed at high temperatures indicating a water-gas reaction that could reduce carbon clusters during the combustion of solid propellant. The water-gas reaction, as well as the population of H₂ at high temperature, points out the way of adjusting the formula of the propellant, which is adding fuel and oxidizer to improve combustion temperature and oxygen balance.

Anatomies for the thermal decomposition behavior and product rule of 5,5′-dinitro-2H,2H′-3,3′-bi-1,2,4-triazole

High-performance energetic materials are mainly used in the military, aerospace industry and chemical fields. The ordinary technology of producing energetic materials cannot avoid the domination of its unique needs. At present, revealing the underlying mechanism of the formation of high-energy materials is of great significance for improving their quality characteristics. We pay special attention to the decomposition and reactive molecular dynamics (RMD) simulation of 5,5′-dinitro-2H,2H′-3,3′-bi-1,2,4-triazole (DNBT). Various forms were captured in the simulation, and the form is determined by the temperature of the initial reactant. By observing the heating pattern and morphological changes under the initial thermal equilibrium, interesting temperature jumps were found in 325 K and 350 K. Observation of continuous heating (simulated temperatures are 2600 K, 2900 K, 3200 K and 3500 K) shows that DNBT has the maximum heating rate at 3500 K. In addition, N₂ occupies this dominant position in the product, moreover, N₂ and NO₂ respectively dominate the gas phase products during the initial heating process. According to the transition state analysis results of the intermediates, we found 4 interesting intermediate products, which were determined by high frequency reaction under the 4 simulated temperatures and performed with transition state calculations. It shows that the selection of reactant temperature and its activity is the key to orderly decomposition of DNBT. It is expected that these findings will be widely used in comprehensive decomposition devices and to improve the concept of learning military and industrial technology.

The performance and behavior of DNBT under RMD simulation at high temperature (2600 K, 2900 K, 3200 K and 3500 K).

Structural Rearrangement of Energetic Materials under an External Electric Field: A Case Study of Nitromethane

As a significant stimulus, the external electric field (EEF) can change the decomposition mechanism and energy release of energetic materials (EMs). Hence, understanding the response of EMs to an EEF is greatly meaningful for their safe usage. Herein, the structural arrangement, a crucial factor in the impact sensitivity and detonation performance of EMs, under the EEF ranging from 0.0 to 0.5 V/Å was investigated via molecular dynamics simulation. Nitromethane (NM) was taken as a case study due to the simple structure. The simulation results show that there exists a critical EEF strength between 0.2 and 0.3 V/Å, which can induce the transition of NM molecules from relatively disordered distribution to solidlike ordered and compacted arrangement with a large density. In this ordered structure, NM dipoles are aligned in a head-to-tail pattern parallel to the EEF direction because of the favored dipole-dipole interactions and weak C-H···O hydrogen bonds. As the EEF strength is enhanced, the potential energy and cohesive energy density of the NM system gradually decrease and increase, respectively, indicative of high thermodynamics stability of ordered arrangement. The results reported here also shed light on the potential of the EEF to induce the nucleation and crystallization to explore new polymorphs of EMs.

Quantum mechanical force fields for condensed phase molecular simulations

Molecular simulations are powerful tools for providing atomic-level details into complex chemical and physical processes that occur in the condensed phase. For strongly interacting systems where quantum many-body effects are known to play an important role, density-functional methods are often used to provide the model with the potential energy used to drive dynamics. These methods, however, suffer from two major drawbacks. First, they are often too computationally intensive to practically apply to large systems over long time scales, limiting their scope of application. Second, there remain challenges for these models to obtain the necessary level of accuracy for weak non-bonded interactions to obtain quantitative accuracy for a wide range of condensed phase properties. Quantum mechanical force fields (QMFFs) provide a potential solution to both of these limitations. In this review, we address recent advances in the development of QMFFs for condensed phase simulations. In particular, we examine the development of QMFF models using both approximate and ab initio density-functional models, the treatment of short-ranged non-bonded and long-ranged electrostatic interactions, and stability issues in molecular dynamics calculations. Example calculations are provided for crystalline systems, liquid water, and ionic liquids. We conclude with a perspective for emerging challenges and future research directions.

Photofragmentation of Tetranitromethane: Spin-Unrestricted Time-Dependent Excited-State Molecular Dynamics

In this study, the photofragmentation dynamics of tetranitromethane (TNM) is explored by a spin-unrestricted time-dependent excited-state molecular dynamics (u-TDESMD) algorithm based on Rabi oscillations and principles similar to trajectory surface hopping, with a midintensity field approximation. The leading order process is represented by the molecule undergoing cyclic excitations and de-excitations. During excitation cycles, the nuclear kinetic energy is accumulated to overcome the dissociation barriers in the reactant and a sequence of intermediates. The dissociation pathway includes the ejection of NO₂ groups followed by the formation of NO and CO. The simulated mass spectra at the ab initio level, based on the bond length in possible fragments, are extracted from simulation trajectories. The recently developed methodology has the potential to model and monitor photoreactions with open-shell intermediates and radicals.