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

Unravelling the prebiotic origins of the simplest α-ketoacids in cometary ices: a computational investigation

We have employed the ab initio nanoreactor (AINR) and DFT calculations to explore how the soft impact of comets entering early earth's dense atmosphere could induce chemical reactions in trapped interstellar ice components, leading to the origin of glyoxylic and pyruvic acids the simplest α-ketoacids essential for prebiotic metabolic cycles.

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. Our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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.

Glycine amino acid transformation under impacts by small solar system bodies, simulated via high-pressure torsion method

Impacts by small solar system bodies (meteoroids, asteroids, comets and transitional objects) are characterized by a combination of energy dynamics and chemical modification on both terrestrial and small solar system bodies. In this context, the discovery of glycine amino acid in meteorites and comets has led to a hypothesis that impacts by astronomical bodies could contribute to delivery and polymerization of amino acids in the early Earth to generate proteins as essential molecules for life. Besides the possibility of abiotic polymerization of glycine, its decomposition by impacts could generate reactive groups to form other essential organic biomolecules. In this study, the high-pressure torsion (HPT) method, as a new platform for simulation of impacts by small solar system bodies, was applied to glycine. In comparison with high-pressure shock experiments, the HPT method simultaneously introduces high pressure and deformation strain. It was found that glycine was not polymerized in the experimental condition assayed, but partially decomposed to ethanol under pressures of 1 and 6 GPa and shear strains of < 120 m/m. The detection of ethanol implies the inherent availability of remaining nitrogen-containing groups, which can incorporate to the formation of other organic molecules at the impact site. In addition, this finding highlights a possibility of the origin of ethanol previously detected in comets.

Synthesis of CxNy-rich polycyclic oligomers from primeval monomers in aqueous media

A multichannel, non-thermolytic and efficient pathway is described toward the formation of functionalized carbon nitride-like oligomers, starting from readily available cyanamide and glyoxal (in ratios >2), in aqueous media under mild conditions. Such oligomers can be isolated as stable solids that result from structures involving cyanamide self-additions along with structures formally derived from the condensation of cyanamide, dicyandiamide or melamine with glyoxal, leading occasionally to oxygen-containing units. The oligomeric aggregates have masses up to 500 u, as inferred from mass spectra analyses, and their formation can be rationalized in terms of polyadditions of cyanamide (up to 10-mer) and glyoxal. The latter is not only a willing reaction partner, but also promotes facile condensation by enhancing the reactivity of nitrile fragments and inducing a significant lowering of the energy barriers. This mechanistic surmise is also supported by DFT calculations of the early condensation steps. As a result, melamine/triazine-type structures are obtained in aquatic environments under much milder conditions than those usually required by other synthetic procedures. Moreover, our results also help unveil the abiotic processes affording complex organic matter on celestial bodies and early earth.

Davis Computational Spectroscopy Workflow-From Structure to Spectra

We describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ∼100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

Oxidative β-Cleavage of Fused Cyclobutanols Leading to Hydrofuran-Fused Polycyclic Aromatic Compounds

Treatment of aryl-fused bicyclo[4.2.0]octanols with an oxidant such as phenyliodine diacetate (PIDA) or hypochlorous acid gave dihydrofuran-containing polycyclic aromatic compounds by selective β-cleavage of the cyclobutanol moiety. Mechanistic studies suggest that the oxygen atom of the hydrofuran ring is incorporated from the hydroxy group of the substrate via intramolecular addition. The oxidative transformation should serve as a new method to prepare functionalized polycyclic aromatic compounds.

A molecular beam and computational study on the barrierless gas phase formation of (iso)quinoline in low temperature extraterrestrial environments

Despite remarkable progress toward the understanding of the formation pathways leading to polycyclic aromatic hydrocarbons (PAHs) in combustion systems and in deep space, the complex reaction pathways leading to nitrogen-substituted PAHs (NPAHs) at low temperatures of molecular clouds and hydrocarbon-rich, nitrogen-containing atmospheres of planets and their moons like Titan have remained largely obscure. Here, we demonstrate through laboratory experiments and computations that the simplest prototype of NPAHs - quinoline and isoquinoline (C₉H₇N) - can be synthesized via rapid and de-facto barrier-less reactions involving o-, m- and p-pyridinyl radicals (C₅H₄N˙) with vinylacetylene (C₄H₄) under low-temperature conditions.

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. We have focused our workflow on TiH₂ as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

Minimization of Gibbs Energy in High-Pressure Multiphase, Multicomponent Mixtures through Particle Swarm Optimization

We present a global optimization method to construct phase boundaries in multicomponent mixtures by minimizing the Gibbs energy. The minimization method is, in essence, an extension of the Maxwell construction procedure that is used in single-component systems. For a given temperature, pressure, and overall mixture composition, it reveals the mole fractions of the thermodynamically stable phases and the composition of these phases. Our approach is based on particle swarm optimization (PSO), which is a gradient-free, stochastic method. It is not reliant on good initial guesses for the phase fractions and compositions, which is an important requirement for the high-pressure applications considered in this study because data on phase boundaries at high pressures tend to be extremely limited. One practical use of this method is to create equation-of-state tables needed by continuum-scale, multiphysics codes that are ubiquitous in high-pressure science. Currently, there does not exist a method to generate such tables that rigorously account for changes in phase boundaries due to mixing. We have done extensive testing to demonstrate that PSO can reliably determine the Gibbs energy minimum and can capture nontrivial features like eutectic and peritectic temperatures to produce coherent phase diagrams. As part of our testing, we have developed a PSO-based Helmholtz-energy minimization procedure that we have used to cross-check the results of the Gibbs energy minimization. We conclude with a critique of our approach and provide suggestions for future work, including a PSO-based entropy-maximization method that would enable the aforementioned continuum codes to perform on-the-fly, phase-equilibria calculations of multicomponent mixtures.

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.

A Hotspot's Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material's microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. The origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. We adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

Extended Lagrangian Born-Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

Extended Lagrangian Born-Oppenheimer molecular dynamics (XL-BOMD) [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] is formulated for orbital-free Hohenberg-Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born-Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. These backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si-O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.

Mechanochemical synthesis of glycine oligomers in a virtual rotational diamond anvil cell

Mechanochemistry of glycine under compression and shear at room temperature is predicted using quantum-based molecular dynamics (QMD) and a simulation design based on rotational diamond anvil cell (RDAC) experiments. Ensembles of high throughput semiempirical density functional tight binding (DFTB) simulations are used to identify chemical trends and bounds for glycine chemistry during rapid shear under compressive loads of up to 15.6 GPa. Significant chemistry is found to occur during compressive shear above 10 GPa. Recovered products consist of small molecules such as water, structural analogs to glycine, heterocyclic molecules, large oligomers, and polypeptides including the simplest polypeptide glycylglycine at up to 4% mass fraction. The population and size of oligomers generally increases with pressure. A number of oligomeric polypeptide precursors and intermediates are also identified that consist of two or three glycine monomers linked together through C-C, C-N, and/or C-O bridges. Even larger oligomers also form that contain peptide C-N bonds and exhibit branched structures. Many of the product molecules exhibit one or more chiral centers. Our simulations demonstrate that athermal mechanical compressive shearing of glycine is a plausible prebiotic route to forming polypeptides.

An Isosymmetric High-Pressure Phase Transition in α-Glycylglycine: A Combined Experimental and Theoretical Study

We investigated the effects of hydrostatic pressure on α-glycylglycine (α-digly) using a combined experimental and theoretical approach. The results of powder X-ray diffraction show a change in compressibility of the axes above 6.7 GPa, but also indicate that the structure remains in the same monoclinic space group, suggesting an isosymmetric phase transition. A noticeable change in the Raman spectra between 6 and 7.5 GPa further supports the observed phase transition. First-principles-based calculations combined with the crystal structure prediction code USPEX predict a number of possible polymorphs at high pressure. An orthorhombic structure with a bent peptide backbone is the lowest enthalpy polymorph above 6.4 GPa; however, it is not consistent with experimental observations. A second monoclinic structure isosymmetric to α-digly, α'-digly, is predicted to become more stable above 11.4 GPa. The partial atomic charges in α'-digly differ from α-digly, and the molecule is bent, possibly indicating different reactivity of α'-digly. The similarity in the lattice parameters predicted from calculations and the axial changes observed experimentally support that the α'-digly phase is likely observed at high pressure. A possible explanation for the isosymmetric phase transition is discussed in terms of relaxing strained hydrogen bonding interactions. Such combined experimental and modeling efforts provide atomic-level insight into how pressure-driven conformational changes alter hydrogen-bonding networks in complicated molecular crystals.

Organic Molecules: Is It Possible to Distinguish Aromatics from Aliphatics Collected by Space Missions in High Speed Impacts?

A prime site of astrobiological interest within the Solar System is the interior ocean of Enceladus. This ocean has already been shown to contain organic molecules and is thought to have the conditions necessary for more complex organic biomolecules to emerge and potentially even life itself. This sub-surface ocean has been accessed by Cassini, an unmanned spacecraft that interacted with the water plumes ejected naturally from Enceladus. The encounter speed with these plumes and their contents was 5 km s−1 and above. Encounters at such speeds allow analysis of vaporised material from submicron-sized particles within the plume, but sampling micron-sized particles remains an open question. The latter particles can impact metal targets exposed on the exterior of future spacecraft, producing impact craters lined with impactor residue, which can then be analysed. Although there is considerable literature on how mineral grains behave in such high-speed impacts, and also on the relationship between the crater residue and the original grain composition, far less is known regarding the behaviour of organic particles. Here we consider a deceptively simple yet fundamental scientific question: for impacts at speeds of around 5–6 kms−1 would the impactor residue alone be sufficient to enable us to recognise the signature conferred by organic particles? Furthermore, would it be possible to identify the organic molecules involved, or at least distinguish between aromatic and aliphatic chemical structures? For polystyrene (aromatic-rich) and poly (methyl methacrylate) (solely aliphatic) latex particles impinging at around 5 km s−1 onto metal targets, we found that sufficient residue is retained at the impact site to permit identification of a carbon-rich projectile, but not of the particular molecules involved, nor is it currently possible to discriminate between aromatic-rich and solely aliphatic particles. This suggests that an alternative analytical method to simple impacts on metal targets is required to enable successful collection of organic samples in a fly-by Enceladus mission.