Quantum mechanical corrections to simulated shock Hugoniot temperatures

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 C_xN_yO_z 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.

Synthesis of functionalized nitrogen-containing polycyclic aromatic hydrocarbons and other prebiotic compounds in impacting glycine solutions

Proteinogenic amino acids can be produced on or delivered to a planet via impacting abiotic sources and consequently were likely present before the emergence of life on Earth. However, the role that these materials played in prebiotic scenarios remains an open question, in part because little is known about the survivability and reactivity of astrophysical organic compounds upon impact with a planetary surface. To this end, we use a force-matched semi-empirical quantum simulation method to study impacts of aqueous proteinogenic amino acids at conditions reaching 48 GPa and 3000 K. Here, we probe a relatively unstudied mechanism for prebiotic synthesis where sudden heating and pressurization causes condensation of complex carbon-rich structures from mixtures of glycine, the simplest protein-forming amino acid. These carbon-containing clusters are stable on short timescales and undergo a fundamental structural transition upon expansion and cooling from predominantly sp3-bonded tetrahedral-like moieties to those that are more sp2-bonded and planar. The recovered sp2-bonded structures include large nitrogen containing polycyclic aromatic hydrocarbons (NPAHs) with a number of different functional groups and embedded bonded regions akin to oligo-peptides. A number of small organic molecules with prebiotic relevance are also predicted to form. This work presents an alternate route to gas-phase synthesis for the formation of NPAHs of high complexity and highlights the significance of both the thermodynamic path and local chemical self-assembly in forming prebiotic species during shock synthesis. Our results help determine the role of comets and other celestial bodies in both the delivery and synthesis of potentially significant life building compounds on early Earth.

Dynamic formation of nanodiamond precursors from the decomposition of carbon suboxide (C3O2) under extreme conditions-A ReaxFF study

We use molecular dynamics simulations with the ReaxFF-lg potential to model the high pressure pyrolysis of carbon suboxide (C₃O₂) in mixture with argon as a pressure bath. We show that the reactive simulations catch the experimental behavior of the low-pressure detonation of C₃O₂ (around 10 bars in shock tube experiments) and allow extrapolations to the high-pressure range of solid-state explosive detonation (up to 60 GPa). While at low pressure carbonaceous nanostructures are formed through the aggregation of species such as carbon dimers C₂, it appears that the high pressure deeply modifies the process, with the aggregation of growing C_xO_y heterostructures, in which the oxygen amount is driven by the pressure and the temperature. Pressures in the order of 60 GPa lead to high oxygen ratios, which prevent carbon atoms to get four carbon neighbors (the first condition to get a diamond structure). But a pressure lowering leads to a substantial carbon enrichment through CO₂/CO release and facilitates the formation of pure sp³-carbon phases where diamond precursors can form. These results give new insights on the conditions leading to nanodiamonds during the detonation of carbon-rich high explosives.

Ab initio calculation of thermodynamic potentials and entropies for superionic water

We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII and X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. Differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed.

Prediction of the Chapman–Jouguet chemical equilibrium state in a detonation wave from first principles based reactive molecular dynamics

This Rx2CJ first principle based protocol for predicting the CJ state provides the matching point between atomistic reaction dynamic simulations and the macroscopic properties of detonation, and can be used as a measure of performance for in silico synthesis and characterization of new materials.

A density functional tight binding model with an extended basis set and three-body repulsion for hydrogen under extreme thermodynamic conditions

We present a new DFTB-p3b density functional tight binding model for hydrogen at extremely high pressures and temperatures, which includes a polarizable basis set (p) and a three-body environmentally dependent repulsive potential (3b). We find that use of an extended basis set is necessary under dissociated liquid conditions to account for the substantial p-orbital character of the electronic states around the Fermi energy. The repulsive energy is determined through comparison to cold curve pressures computed from density functional theory (DFT) for the hexagonal close-packed solid, as well as pressures from thermally equilibrated DFT-MD simulations of the liquid phase. In particular, we observe improved agreement in our DFTB-p3b model with previous theoretical and experimental results for the shock Hugoniot of hydrogen up to 100 GPa and 25000 K, compared to a standard DFTB model using pairwise interactions and an s-orbital basis set, only. The DFTB-p3b approach discussed here provides a general method to extend the DFTB method for a wide variety of materials over a significantly larger range of thermodynamic conditions than previously possible.

Equation of state and phase diagram of ammonia at high pressures from ab initio simulations

We present an equation of state as well as a phase diagram of ammonia at high pressures and high temperatures derived from ab initio molecular dynamics simulations. The predicted phases of ammonia are characterized by analyzing diffusion coefficients and structural properties. Both the phase diagram and the subsequently computed Hugoniot curves are compared to experimental results. Furthermore, we discuss two methods that allow us to take into account nuclear quantum effects, which are of considerable importance in molecular fluids. Our data cover pressures up to 330 GPa and a temperature range from 500 K to 10,000 K. This regime is of great interest for interior models of the giant planets Uranus and Neptune, which contain, besides water and methane, significant amounts of ammonia.

Prebiotic chemistry within a simple impacting icy mixture

We present results of prebiotic organic synthesis in shock compressed mixtures of simple ices from quantum molecular dynamics (MD) simulations extended to close to equilibrium time scales. Given the likelihood of an inhospitable prebiotic atmosphere on early Earth, it is possible that impact processes of comets or other icy bodies were a source of prebiotic chemical compounds on the primitive planet. We observe that moderate shock pressures and temperatures within a CO2-rich icy mixture (36 GPa and 2800 K) produce a number of nitrogen containing heterocycles, which dissociate to form functionalized aromatic hydrocarbons upon expansion and cooling to ambient conditions. In contrast, higher shock conditions (48-60 GPa, 3700-4800 K) resulted in the synthesis of long carbon-chain molecules, CH4, and formaldehyde. All shock compression simulations at these conditions have produced significant quantities of simple C-N bonded compounds such as HCN, HNC, and HNCO upon expansion and cooling to ambient conditions. Our results elucidate a mechanism for impact synthesis of prebiotic molecules at realistic impact conditions that is independent of external constraints such as the presence of a catalyst, illuminating UV radiation, or pre-existing conditions on a planet.

New Developments in the Physical Chemistry of Shock Compression

This review discusses new developments in shock compression science with a focus on molecular media. Some basic features of shock and detonation waves, nonlinear excitations that can produce extreme states of high temperature and high pressure, are described. Methods of generating and detecting shock waves are reviewed, especially those using tabletop lasers that can be interfaced with advanced molecular diagnostics. Newer compression methods such as shockless compression and precompression shock that generate states of cold dense molecular matter are discussed. Shock compression creates a metallic form of hydrogen, melts diamond, and makes water a superionic liquid with unique catalytic properties. Our understanding of detonations at the molecular level has improved a great deal as a result of advanced nonequilibrium molecular simulations. Experimental measurements of detailed molecular behavior behind a detonation front might be available soon using femtosecond lasers to produce nanoscale simulated detonation fronts.