Ab initio structural dynamics of pure and nitrogen-containing amorphous carbon.

Amorphous carbon (a-C) has attracted considerable interest due to its desirable properties, which are strongly dependent on its structure, density and impurities. Using ab initio molecular dynamics simulations we show that the sp2/sp3 content and underlying structural order of a-C produced via liquid quenching evolve at high temperatures and pressures on sub-nanosecond timescales. Graphite-like densities ([Formula: see text] 2.7 g/cc) favor the formation of layered arrangements characterized by sp2 disordered bonding resembling recently synthesized monolayer amorphous carbon (MAC), while at diamond-like densities ([Formula: see text] 3.3 g/cc) the resulting structures are dominated by disordered tetrahedral sp3 hybridization typical of diamond-like amorphous carbon (DLC). At intermediate densities the system is a highly compressible mixture of coexisting sp2 and sp3 regions that continue to segregate over 10's of picoseconds. The addition of nitrogen (20.3%) (a-CN) generates major system features similar with those of a-C, but has the unexpected effect of reinforcing the thermodynamically disfavored carbon structural motifs at low and high densities, while inhibiting phase separation in the intermediate region. At the same time, no nitrogen elimination from the carbon framework is observed above [Formula: see text] 2.8 g/cc, suggesting that nitrogen impurities are likely to remain embedded in the carbon structures during fast temperature quenches at high pressures.

Chemical evolution in nitrogen shocked beyond the molecular stability limit.

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. This suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

Chemistry-mediated Ostwald ripening in carbon-rich C/O systems at extreme conditions.

There is significant interest in establishing a capability for tailored synthesis of next-generation carbon-based nanomaterials due to their broad range of applications and high degree of tunability. High pressure (e.g., shockwave-driven) synthesis holds promise as an effective discovery method, but experimental challenges preclude elucidating the processes governing nanocarbon production from carbon-rich precursors that could otherwise guide efforts through the prohibitively expansive design space. Here we report findings from large scale atomistically-resolved simulations of carbon condensation from C/O mixtures subjected to extreme pressures and temperatures, made possible by machine-learned reactive interatomic potentials. We find that liquid nanocarbon formation follows classical growth kinetics driven by Ostwald ripening (i.e., growth of large clusters at the expense of shrinking small ones) and obeys dynamical scaling in a process mediated by carbon chemistry in the surrounding reactive fluid. The results provide direct insight into carbon condensation in a representative system and pave the way for its exploration in higher complexity organic materials. They also suggest that simulations using machine-learned interatomic potentials could eventually be employed as in-silico design tools for new nanomaterials.

Submicrosecond Aggregation during Detonation Synthesis of Nanodiamond.

Detonation nanodiamond (DND) is known to form aggregates that significantly reduce their unique nanoscale properties and require postprocessing to separate. How and when DND aggregates is an important question that has not been answered experimentally and could provide the foundation for approaches to limit aggregation. To answer this question, time-resolved small-angle X-ray scattering was performed during the detonation of high-explosives that are expected to condense particulates in the diamond, graphite, and liquid regions of the carbon phase diagram. DND aggregation into low fractal dimension structures could be observed as early as 0.1 µs, along with a separate scattering population also observed from an explosive that produces primarily graphitic products. A counterexample is the case of a high-explosive that produces nano-onions, where no hierarchical scattering was observed for at least 10 µs behind the detonation front. These results suggest that DND aggregation occurs on time scales comparable to particle formation.

Investigating 3,4-bis(3-nitrofurazan-4-yl)furoxan detonation with a rapidly tuned density functional tight binding model.

We describe a machine learning approach to rapidly tune density functional tight binding models for the description of detonation chemistry in organic molecular materials. Resulting models enable simulations on the several 10s of ps scales characteristic to these processes, with "quantum-accuracy." We use this approach to investigate early shock chemistry in 3,4-bis(3-nitrofurazan-4-yl)furoxan, a hydrogen-free energetic material known to form onion-like nanocarbon particulates following detonation. We find that the ensuing chemistry is significantly characterized by the formation of large CxNyOz species, which are likely precursors to the experimentally observed carbon condensates. Beyond utility as a means of investigating detonation chemistry, the present approach can be used to generate quantum-based reference data for the development of full machine-learned interatomic potentials capable of simulation on even greater time and length scales, i.e., for applications where characteristic time scales exceed the reach of methods including Kohn-Sham density functional theory, which are commonly used for reference data generation.

Active learning for robust, high-complexity reactive atomistic simulations.

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. Here, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn-Sham density functional theory in terms of structure, dynamics, and speciation.

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions.

We describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm-3, affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure-high temperature force-field for condensed phase chemistry in organic materials.

Ultrafast shock synthesis of nanocarbon from a liquid precursor.

Carbon nanoallotropes are important nanomaterials with unusual properties and promising applications. High pressure synthesis has the potential to open new avenues for controlling and designing their physical and chemical characteristics for a broad range of uses but it remains little understood due to persistent conceptual and experimental challenges, in addition to fundamental physics and chemistry questions that are still unresolved after many decades. Here we demonstrate sub-nanosecond nanocarbon synthesis through the application of laser-induced shock-waves to a prototypical organic carbon-rich liquid precursor-liquid carbon monoxide. Overlapping large-scale molecular dynamics simulations capture the atomistic details of the nanoparticles' formation and evolution in a reactive environment and identify classical evaporation-condensation as the mechanism governing their growth on these time scales.

Ultrafast Shock-Induced Reactions in Pentaerythritol Tetranitrate Thin Films.

The chemical and physical processes involved in the shock-to-detonation transition of energetic solids are not fully understood due to difficulties in probing the fast dynamics involved in initiation. Here, we employ shock interferometry experiments with sub-20-ps time resolution to study highly textured (110) pentaerythritol tetranitrate (PETN) thin films during the early stages of shock compression using ultrafast laser-driven shock wave methods. We observe evidence of rapid exothermic chemical reactions in the PETN thin films for interface particle velocities above ∼1.05 km/s as indicated by shock velocities and pressures well above the unreacted Hugoniot. The time scale of our experiment suggests that exothermic reactions begin less than 50 ps behind the shock front for these high-density PETN thin films. Thermochemical calculations for partially reacted Hugoniots also support this interpretation. The experimentally observed time scale of reactivity could be used to narrow possible initiation mechanisms.

Effects of pressure on the structure and lattice dynamics of ammonium perchlorate: A combined experimental and theoretical study.

Ammonium perchlorate NH4ClO4 (AP) was studied using synchrotron angle-dispersive X-ray powder diffraction (XRPD) and Raman spectroscopy. A diamond-anvil cell was used to compress AP up to 50 GPa at room temperature (RT). Density functional theory (DFT) calculations were performed to provide further insight and comparison to the experimental data. A high-pressure barite-type structure (Phase II) forms at ≈4 GPa and appears stable up to 40 GPa. Refined atomic coordinates for Phase II are provided, and details for the Phase I → II transition mechanics are outlined. Pressure-dependent enthalpies computed for DFT-optimized crystal structures confirm the Phase I → II transition sequence, and the interpolated transition pressure is in excellent agreement with the experiment. Evidence for additional (underlying) structural modifications include a marked decrease in the Phase II b'-axis compressibility starting at 15 GPa and an unambiguous stress relaxation in the normalized stress-strain response at 36 GPa. Above 47 GPa, XRD Bragg peaks begin to decrease in amplitude and broaden. The apparent loss of crystalline long-range order likely signals the onset of amorphization. Three isostructural modifications were discovered within Phase II via Raman spectroscopy. A revised RT isothermal phase diagram is discussed based on the findings of this study.

Nanocarbon condensation in detonation.

We analyze the definition of the Gibbs free energy of a nanoparticle in a reactive fluid environment, and propose an approach for predicting the size of carbon nanoparticles produced by the detonation of carbon-rich explosives that regards their condensation as a nucleation process and takes into account absolute entropy effects of the cluster population. The results are consistent with experimental observations and indicate that such entropy considerations are important for determining chemical equilibrium states in energetic materials that contain an excess of carbon. The analysis may be useful for other applications that deal with the nucleation of nanoparticles under reactive conditions.

High-pressure X-ray diffraction, Raman, and computational studies of MgCl2 up to 1 Mbar: Extensive pressure stability of the ß-MgCl2 layered structure.

Magnesium chloride (MgCl2) with the rhombohedral layered CdCl2-type structure (α-MgCl2) has been studied experimentally using synchrotron angle-dispersive powder x-ray diffraction and Raman spectroscopy using a diamond-anvil cell up to 100 GPa at room temperature and theoretically using first-principles density functional calculations. The results reveal a pressure-induced second-order structural phase transition to a hexagonal layered CdI2-type structure (ß-MgCl2) at 0.7 GPa: the stacking sequence of the Cl anions are altered resulting in a reduction of the c-axis length. Theoretical calculations confirm this phase transition sequence and the calculated transition pressure is in excellent agreement with the experiment. Lattice dynamics calculations also reproduce the experimental Raman spectra measured for the ambient and high-pressure phase. According to our experimental results MgCl2 remains in a 2D layered phase up to 100 GPa and further, the 6-fold coordination of Mg cations is retained. Theoretical calculations of relative enthalpy suggest that this extensive pressure stability is due to a low enthalpy of the layered structure ruling out kinetic barrier effects. This observation is unusual, as it contradicts with the general structural behavior of highly compressed AB2 compounds.

Equations of state of anhydrous AlF3 and AlI3: Modeling of extreme condition halide chemistry.

Pressure dependent angle-dispersive x-ray powder diffraction measurements of alpha-phase aluminum trifluoride (α-AlF3) and separately, aluminum triiodide (AlI3) were conducted using a diamond-anvil cell. Results at 295 K extend to 50 GPa. The equations of state of AlF3 and AlI3 were determined through refinements of collected x-ray diffraction patterns. The respective bulk moduli and corresponding pressure derivatives are reported for multiple orders of the Birch-Murnaghan (B-M), finite-strain (F-f), and higher pressure finite-strain (G-g) EOS analysis models. Aluminum trifluoride exhibits an apparent isostructural phase transition at approximately 12 GPa. Aluminum triiodide also undergoes a second-order atomic rearrangement: applied stress transformed a monoclinically distorted face centered cubic (fcc) structure into a standard fcc structural arrangement of iodine atoms. Results from semi-empirical thermochemical computations of energetic materials formulated with fluorine containing reactants were obtained with the aim of predicting the yield of halogenated products.

A simulation assessment of the thermodynamics of dense ion-dipole mixtures with polarization.

Molecular dynamics (MD) simulations are employed to ascertain the relative importance of various electrostatic interaction contributions, including induction interactions, to the thermodynamics of dense, hot ion-dipole mixtures. In the absence of polarization, we find that an MD-constrained free energy term accounting for the ion-dipole interactions, combined with well tested ionic and dipolar contributions, yields a simple, fairly accurate free energy form that may be a better option for describing the thermodynamics of such mixtures than the mean spherical approximation (MSA). Polarization contributions induced by the presence of permanent dipoles and ions are found to be additive to a good approximation, simplifying the thermodynamic modeling. We suggest simple free energy corrections that account for these two effects, based in part on standard perturbative treatments and partly on comparisons with MD simulation. Even though the proposed approximations likely need further study, they provide a first quantitative assessment of polarization contributions at high densities and temperatures and may serve as a guide for future modeling efforts.

Ultrafast shock compression of an oxygen-balanced mixture of nitromethane and hydrogen peroxide.

We apply ultrafast optical interferometry to measure the Hugoniot of an oxygen-balanced mixture of nitromethane and hydrogen peroxide (NM/HP) and compare with Hugoniot data for pure nitromethane (NM) and a 90% hydrogen peroxide/water mixture (HP), as well as theoretical predictions. We observe a 2.1% percent mean pairwise difference between the measured shockwave speed (at the measured piston speed) in unreacted NM/HP and the corresponding "universal" liquid Hugoniot, which is larger than the average standard deviation of our data, 1.4%. Unlike the Hugoniots of both HP and NM, in which measured shock speeds deviate to values greater than the unreacted Hugoniot for piston speeds larger than the respective reaction thresholds, in the NM/HP mixture we observe shock speed deviations to values lower than the unreacted Hugoniot well below the von Neumann pressure (≈28 GPa). Although the trend should reverse for high enough piston speeds, the initial behavior is unexpected. Possible explanations range from mixing effects to a complex index of refraction in the reacted solution. If this is indeed a signature of chemical initiation, it would suggest that the process may not be kinetically limited (on a ~100 ps time scale) between the initiation threshold and the von Neumann pressure.

Nitrogen oxides as a chemistry trap in detonating oxygen-rich materials.

Despite decades of research, the chemical processes and states of matter that govern the behavior of energetic materials under detonation conditions are not well understood, including the molecular-level processes that determine decomposition kinetics and energy release. Oxygen content is often employed as a simple and intuitive guide to the development and practical use of explosives, but its effect on detonation chemistry remains little studied, especially for the case of oxygen overabundance. To this end, we have conducted quantum molecular dynamics (QMD) simulations of zero oxygen balance and oxygen-rich mixtures of hydrogen peroxide and nitromethane under detonation-like conditions to near-equilibrium time scales. We find excellent agreement between our extrapolated chemical equilibrium properties and those from thermochemical models for the zero oxygen balance mixture. In contrast, for the oxygen-rich mixture, we observe the formation of nitrogen oxide intermediates, particularly nitrate ions (NO3), that effectively act as an oxygen/nitrogen "trap" by precluding the formation of the equilibrium products N2 and CO2. Our results could have implications for the design and modeling of oxygen-rich energetics in common military and industrial use.

Ultrafast shock initiation of exothermic chemistry in hydrogen peroxide.

We report observations of shock compressed, unreacted hydrogen peroxide at pressures up to the von Neumann pressure for a steady detonation wave, using ultrafast laser-driven shock wave methods. At higher laser drive energy we find evidence of exothermic chemical reactivity occurring in less than 100 ps after the arrival of the shock wave in the sample. The results are consistent with our MD simulations and analysis and suggest that reactivity in hydrogen peroxide is initiated on a sub-100 ps time scale under conditions found just subsequent to the lead shock in a steady detonation wave.

Experimental measurement of speeds of sound in dense supercritical carbon monoxide and development of a high-pressure, high-temperature equation of state.

We report the adiabatic sound speeds for supercritical fluid carbon monoxide along two isotherms, from 0.17 to 2.13 GPa at 297 K and from 0.31 to 3.2 GPa at 600 K. The carbon monoxide was confined in a resistively heated diamond-anvil cell, and the sound speed measurements were conducted in situ using a recently reported variant of the photoacoustic light scattering effect. The measured sound speeds were then used to parametrize a single site dipolar exponential-6 intermolecular potential for carbon monoxide. PρT thermodynamic states, sound speeds, and shock Hugoniots were calculated using the newly parametrized intermolecular potential and compared to previously reported experimental results. Additionally, we generated an analytical equation of state for carbon monoxide by fitting to a grid of calculated PρT states over a range of 0.1-10 GPa and 150-2000 K. A 2% mean variation was found between computed high-pressure solid-phase densities and measured data-a surprising result for a spherical interaction potential. We further computed a rotationally dependent fluid to ß-solid phase boundary; results signal the relative magnitude of short-range rotational disorder under conditions that span existing phase boundary measurements.

Thermodynamics and diffusion in size-symmetric and asymmetric dense electrolytes.

MD simulation results for model size-symmetric and asymmetric electrolytes at high densities and temperatures (well outside the liquid-gas coexistence region) are generated and analyzed focusing on thermodynamic and diffusion properties. An extension of the mean spherical approximation for electrolytes originally derived for charged hard sphere fluids is adapted to these systems by exploiting the separation of short range and Coulomb interaction contributions intrinsic to these theoretical models and is found to perform well for predicting equation of state quantities. The diffusion coefficients of these electrolytes can also be reasonably well predicted using entropy scaling ideas suitably adapted to charged systems and mixtures. Thus, this approach may provide an avenue for studying dense electrolytes or complex molecular systems containing charged species at high pressures and temperatures.

Dissociative melting of ice VII at high pressure.

We have used x-ray diffraction to determine the structure factor of water along its melting line to a static pressure of 57 GPa (570 kbar) and a temperature of more than 1500 K, conditions which correspond to the lower mantle of the Earth, and the interiors of Neptune and Uranus up to a depth of 7000 km. We have also performed corresponding first principles and classical molecular dynamics simulations. Above a pressure of 4 GPa the O-O structure factor is found to be very close to that of a simple soft sphere liquid, thus permitting us to determine the density of liquid water near the melting line. By comparing these results with the density of ice, also determined in this study, we find that the enthalpy of fusion (DeltaH(f)) increases enormously along the melting line, reaching approximately 120 kJ/mole at 40 GPa (compared to 6 kJ/mole at 0 GPa), thus revealing significant molecular dissociation of water upon melting. We speculate that an extended two-phase region could occur in planetary processes involving the adiabatic compression of water.