ReaxFF reactive molecular dynamics on silicon pentaerythritol tetranitrate crystal validates the mechanism for the colossal sensitivity.

Recently quantum mechanical (QM) calculations on a single Si-PETN (silicon-pentaerythritol tetranitrate) molecule were used to explain its colossal sensitivity observed experimentally in terms of a unique Liu carbon-silyl nitro-ester rearrangement (R3Si-CH2-O-R2→ R3Si-O-CH2-R2). In this paper we expanded the study of Si-PETN from a single molecule to a bulk system by extending the ReaxFF reactive force field to describe similar Si-C-H-O-N systems with parameters optimized to reproduce QM results. The reaction mechanisms and kinetics of thermal decomposition of solid Si-PETN were investigated using ReaxFF reactive molecular dynamics (ReaxFF-RMD) simulations at various temperatures to explore the origin of the high sensitivity. We find that at lower temperatures, the decomposition of Si-PETN is initiated by the Liu carbon-silyl nitro-ester rearrangement forming Si-O bonds which is not observed in PETN. As the reaction proceeds, the exothermicity of Si-O bond formation promotes the onset of NO2 formation from N-OC bond cleavage which does not occur in PETN. At higher temperatures PETN starts to react by the usual mechanisms of NO2 dissociation and HONO elimination; however, Si-PETN remains far more reactive. These results validate the predictions from QM that the significantly increased sensitivity of Si-PETN arises from a unimolecular process involving the unusual Liu rearrangement but not from multi-molecular collisions. It is the very low energy barrier and the high exothermicity of the Si-O bond formation providing energy early in the decomposition process that is responsible.

Development of a ReaxFF reactive force field for ettringite and study of its mechanical failure modes from reactive dynamics simulations.

Ettringite is a hexacalcium aluminate trisulfate hydrate mineral that forms during Portland cement hydration. Its presence plays an important role in controlling the setting rate of the highly reactive aluminate phases in cement paste and has also been associated with severe cracking in cured hardened cement. To understand how it forms and how its properties influence those of hardened cement and concrete, we have developed a first-principles-based ReaxFF reactive force field for Ca/Al/H/O/S. Here, we report on the development of this ReaxFF force field and on its validation and application using reactive molecular dynamics (RMD) simulations to characterize and understand the elastic, plastic, and failure response of ettringite at the atomic scale. The ReaxFF force field was validated by comparing the lattice parameters, pairwise distribution functions, and elastic constants of an ettringite crystal model obtained from RMD simulations with those from experiments. The predicted results are in close agreement with published experimental data. To characterize the atomistic failure modes of ettringite, we performed stress-strain simulations to find that Ca-O bonds are responsible for failure of the calcium sulfate and tricalcium aluminate (C3A) column in ettringite during uniaxial compression and tension and that hydrogen bond re-formation during compression induces an increase in plastic strain beyond the material's stress-strain proportionality limit. These results provide essential insight into understanding the mechanistic role of this mineral in cement and concrete degradation, and the ReaxFF potential developed in this work serves as a fundamental tool to further study the kinetics of hydration in cement and concrete.

ReaxFF-lg: correction of the ReaxFF reactive force field for London dispersion, with applications to the equations of state for energetic materials.

The practical levels of density functional theory (DFT) for solids (LDA, PBE, PW91, B3LYP) are well-known not to account adequately for the London dispersion (van der Waals attraction) so important in molecular solids, leading to equilibrium volumes for molecular crystals ~10-15% too high. The ReaxFF reactive force field is based on fitting such DFT calculations and suffers from the same problem. In the paper we extend ReaxFF by adding a London dispersion term with a form such that it has low gradients (lg) at valence distances leaving the already optimized valence interactions intact but behaves as 1/R(6) for large distances. We derive here these lg corrections to ReaxFF based on the experimental crystal structure data for graphite, polyethylene (PE), carbon dioxide, and nitrogen and for energetic materials: hexahydro-1,3,5-trinitro-1,3,5-s-triazine (RDX), pentaerythritol tetranitrate (PETN), 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), and nitromethane (NM). After this dispersion correction the average error of predicted equilibrium volumes decreases from 18.5 to 4.2% for the above systems. We find that the calculated crystal structures and equation of state with ReaxFF-lg are in good agreement with experimental results. In particular, we examined the phase transition between α-RDX and γ-RDX, finding that ReaxFF-lg leads to excellent agreement for both the pressure and volume of this transition occurring at ~4.8 GPa and ~2.18 g/cm(3) density from ReaxFF-lg vs 3.9 GPa and ~2.21 g/cm(3) from experiment. We expect ReaxFF-lg to improve the descriptions of the phase diagrams for other energetic materials.

Mechanism and kinetics for the initial steps of pyrolysis and combustion of 1,6-dicyclopropane-2,4-hexyne from ReaxFF reactive dynamics.

We report the kinetic analysis and mechanism for the initial steps of pyrolysis and combustion of a new fuel material, 1,6-dicyclopropane-2,4-hexyne, that has enormous heats of pyrolysis and combustion, making it a potential high-energy fuel or fuel additive. These studies employ the ReaxFF force field for reactive dynamics (RD) simulations of both pyrolysis and combustion processes for both unimolecular and multimolecular systems. We find that both pyrolysis and combustion initiate from unimolecular reactions, with entropy-driven reactions being most important in both processes. Pyrolysis initiates with extrusion of an ethylene molecule from the fuel molecule and is followed quickly by isomerization of the fuel molecule, which induces additional radicals that accelerate the pyrolysis process. In the combustion process, we find three distinct mechanisms for the O(2) attack on the fuel molecule: (1) attack on the cyclopropane, ring expanding to form the cyclic peroxide which then decomposes; (2) attack onto the central single bond of the diyne which then fissions to form two C(5)H(5)O radicals; (3) attack on the alkyne-cyclopropane moiety to form a seven-membered ring peroxide which then decomposes. Each of these unimolecular combustion processes releases energy that induces additional radicals to accelerate the combustion process. Here oxygen has major effects both as the radical acceptor and as the radical producer. We extract both the effective activation energy and the effective pre-exponential factor by kinetic analysis of pyrolysis and combustion from these ReaxFF simulations. The low value of the derived effective activation energy (26.18 kcal/mol for pyrolysis and 16.40 kcal/mol for combustion) reveals the high activity of this fuel molecule.