Prebiotic chemistry within a simple impacting icy mixture

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.

A Computational Quantum-Based Perspective on the Molecular Origins of Life's Building Blocks

The search for the chemical origins of life represents a long-standing and continuously debated enigma. Despite its exceptional complexity, in the last decades the field has experienced a revival, also owing to the exponential growth of the computing power allowing for efficiently simulating the behavior of matter-including its quantum nature-under disparate conditions found, e.g., on the primordial Earth and on Earth-like planetary systems (i.e., exoplanets). In this minireview, we focus on some advanced computational methods capable of efficiently solving the Schro¨dinger equation at different levels of approximation (i.e., density functional theory)-such as ab initio molecular dynamics-and which are capable to realistically simulate the behavior of matter under the action of energy sources available in prebiotic contexts. In addition, recently developed metadynamics methods coupled with first-principles simulations are here reviewed and exploited to answer to old enigmas and to propose novel scenarios in the exponentially growing research field embedding the study of the chemical origins of life.

New Signatures of Bio-Molecular Complexity in the Hypervelocity Impact Ejecta of Icy Moon Analogues

Impact delivery of prebiotic compounds to the early Earth from an impacting comet is considered to be one of the possible ways by which prebiotic molecules arrived on the Earth. Given the ubiquity of impact features observed on all planetary bodies, bolide impacts may be a common source of organics on other planetary bodies both in our own and other solar systems. Biomolecules such as amino acids have been detected on comets and are known to be synthesized due to impact-induced shock processing. Here we report the results of a set of hypervelocity impact experiments where we shocked icy mixtures of amino acids mimicking the icy surface of planetary bodies with high-speed projectiles using a two-stage light gas gun and analyzed the ejecta material after impact. Electron microscopic observations of the ejecta have shown the presence of macroscale structures with long polypeptide chains revealed from LCMS analysis. These results suggest a pathway in which impact on cometary ices containing building blocks of life can lead to the synthesis of material architectures that could have played a role in the emergence of life on the Earth and which may be applied to other planetary bodies as well.

Chance and Necessity in the Evolution of Matter to Life: A Comprehensive Hypothesis

Specialists in several branches of life sciences are trying to solve, piece by piece, the immensely complex puzzle of the origin of life. Some parts of the puzzle seem to appear with a rather high degree of clarity, while others remain totally obscure. We cannot be sure that life emerged only on our Earth, but we believe that the presence of large amounts of water in its liquid state is absolutely essential for the emergence and evolution of living matter. We can also assume that the latter exploits everywhere the same light elements, mainly C, H, O, N, S, and P, and somehow manipulates the same simple monomeric and polymeric organic compounds, such as alpha-amino acids, carbohydrates, nucleic bases, and surface-active carboxylic acids. The author contributes to the field by stating that all fundamental particles of our matter are “homochiral” and predominantly produce in an absolute asymmetric synthesis amino acids of L-configuration and carbohydrates of D-series. Another important point is that free atmospheric oxygen mainly stems from the photolysis of water molecules by cosmic irradiation and is not necessarily bound to living organisms on the planet.

Pickup and reactions of molecules on clusters relevant for atmospheric and interstellar processes

In this perspective, we review experiments with molecules picked up on large clusters in molecular beams with the focus on the processes in atmospheric and interstellar chemistry. First, we concentrate on the pickup itself, and we discuss the pickup cross sections. We measure the uptake of different atmospheric molecules on mixed nitric acid-water clusters and determine the accommodation coefficients relevant for aerosol formation in the Earth's atmosphere. Then the coagulation of the adsorbed molecules on the clusters is investigated. In the second part of this perspective, we review examples of different processes triggered by UV-photons or electrons in the clusters with embedded molecules. We start with the photodissociation of hydrogen halides and Freon CF₂Cl₂ on ice nanoparticles in connection with the polar stratospheric ozone depletion. Next, we mention reactions following the excitation and ionization of the molecules adsorbed on clusters. The first ionization-triggered reaction observed between two different molecules picked up on the cluster was the proton transfer between methanol and formic acid deposited on large argon clusters. Finally, negative ion reactions after slow electron attachment are illustrated by two examples: mixed nitric acid-water clusters, and hydrogen peroxide deposited on large Ar_N and (H₂O)_N clusters. The selected examples are discussed from the perspective of the atmospheric and interstellar chemistry, and several future directions are proposed.

Shock Processing of Amino Acids Leading to Complex Structures-Implications to the Origin of Life

The building blocks of life, amino acids, are believed to have been synthesized in the extreme conditions that prevail in space, starting from simple molecules containing hydrogen, carbon, oxygen and nitrogen. However, the fate and role of amino acids when they are subjected to similar processes largely remain unexplored. Here we report, for the first time, that shock processed amino acids tend to form complex agglomerate structures. Such structures are formed on timescales of about 2 ms due to impact induced shock heating and subsequent cooling. This discovery suggests that the building blocks of life could have self-assembled not just on Earth but on other planetary bodies as a result of impact events. Our study also provides further experimental evidence for the ‘threads’ observed in meteorites being due to assemblages of (bio)molecules arising from impact-induced shocks.

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.

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.

Ab initio molecular dynamics study of prebiotic production processes of organic compounds at meteorite impacts on ocean

Recent experiments concerning prebiotic materials syntheses suggest that the iron-bearing meteorite impacts on ocean during Late Heavy Bombardment provided abundant organic compounds associated with biomolecules such as amino acids and nucleobases. However, the molecular mechanism of a series of chemical reactions to produce such compounds is not well understood. In this study, we simulate the shock compression state of a meteorite impact for a model system composed of CO₂ , H₂ O, and metallic iron slab by ab initio molecular dynamics combined with multiscale shock technique, and clarify possible elementary reaction processes up to production of organic compounds. The reactions included not only pathways similar to the Fischer-Tropsch process known as an important hydrocarbon synthesis in many planetary processes but also those resulting in production of a carboxylic acid. It is also found that bicarbonate ions formed from CO₂ and H₂ O participated in some forms in most of these observed elementary reaction processes. These findings would deepen the understanding of the full range of chemical reactions that could occur in the meteorite impact events. © 2018 Wiley Periodicals, Inc.

Generating Converged Accurate Free Energy Surfaces for Chemical Reactions with a Force-Matched Semiempirical Model

We demonstrate the capability of creating robust density functional tight binding (DFTB) models for chemical reactivity in prebiotic mixtures through force matching to short time scale quantum free energy estimates. Molecular dynamics using density functional theory (DFT) is a highly accurate approach to generate free energy surfaces for chemical reactions, but the extreme computational cost often limits the time scales and range of thermodynamic states that can feasibly be studied. In contrast, DFTB is a semiempirical quantum method that affords up to a thousandfold reduction in cost and can recover DFT-level accuracy. Here, we show that a force-matched DFTB model for aqueous glycine condensation reactions yields free energy surfaces that are consistent with experimental observations of reaction energetics. Convergence analysis reveals that multiple nanoseconds of combined trajectory are needed to reach a steady-fluctuating free energy estimate for glycine condensation. Predictive accuracy of force-matched DFTB is demonstrated by direct comparison to DFT, with the two approaches yielding surfaces with large regions that differ by only a few kcal mol^-1.

Meteorite impacts on ancient oceans opened up multiple NH3 production pathways

A recent series of shock experiments by Nakazawa et al. starting in 2005 (e.g. [Nakazawa et al., Earth Planet. Sci. Lett., 2005, 235, 356]) suggested that meteorite impacts on ancient oceans would have yielded a considerable amount of NH₃ to the early Earth from atmospheric N₂ and oceanic H₂O through reduction by meteoritic iron. To clarify the mechanisms, we imitated the impact events by performing multi-scale shock technique-based ab initio molecular dynamics in the framework of density functional theory in combination with multi-scale shock technique (MSST) simulations. Our previous simulations with impact energies close to that of the experiments revealed picosecond-order rapid NH₃ production during shock compression [Shimamura et al., Sci. Rep., 2016, 6, 38952]. It was also shown that the reduction of N₂ took place with an associative mechanism as seen in the catalysis of nitrogenase enzymes. In this study, we performed an MSST-AIMD simulation to investigate the production by meteorite impacts with higher energies, which are closer to the expected values on the early Earth. It was found that the amount of NH₃ produced further increased. We also found that the increased NH₃ production is due to the emergence of multiple reaction mechanisms at increased impact energies. We elucidated that the reduction of N₂ was not only attributed to the associative mechanism but also to a dissociative mechanism as seen in the Haber-Bosch process and to a mechanism through a hydrazinium ion. The emergence of these multiple production mechanisms capable of providing a large amount of NH₃ would support the suggestions from recent experiments much more strongly than was previously believed, i.e., shock-induced NH₃ production played a key role in the origin of life on Earth.

Clustering of Uracil Molecules on Ice Nanoparticles

We generate a molecular beam of ice nanoparticles (H₂O)_N, N̅ ≈ 130-220, which picks up several individual gas phase uracil (U) or 5-bromouracil (BrU) molecules. The mass spectra of the doped nanoparticles prove that the uracil and bromouracil molecules coagulate to clusters on the ice nanoparticles. Calculations of U and BrU monomers and dimers on the ice nanoparticles provide theoretical support for the cluster formation. The (U)_mH⁺ and (BrU)_mH⁺ intensity dependencies on m extracted from the mass spectra suggest a smaller tendency of BrU to coagulate compared to U, which is substantiated by a lower mobility of bromouracil on the ice surface. The hydrated U_m·(H₂O)_nH⁺ series are also reported and discussed. On the basis of comparison with the previous experiments, we suggest that the observed propensity for aggregation on ice nanoparticles is a more general trend for biomolecules forming strong hydrogen bonds. This, together with their mobility, leads to their coagulation on ice nanoparticles which is an important aspect for astrochemistry.

Meteorite Impact-Induced Rapid NH3 Production on Early Earth: Ab Initio Molecular Dynamics Simulation

NH₃ is an essential molecule as a nitrogen source for prebiotic amino acid syntheses such as the Strecker reaction. Previous shock experiments demonstrated that meteorite impacts on ancient oceans would have provided a considerable amount of NH₃ from atmospheric N₂ and oceanic H₂O through reduction by meteoritic iron. However, specific production mechanisms remain unclear, and impact velocities employed in the experiments were substantially lower than typical impact velocities of meteorites on the early Earth. Here, to investigate the issues from the atomistic viewpoint, we performed multi-scale shock technique-based ab initio molecular dynamics simulations. The results revealed a rapid production of NH₃ within several picoseconds after the shock, indicating that shocks with greater impact velocities would provide further increase in the yield of NH₃. Meanwhile, the picosecond-order production makes one expect that the important nitrogen source precursors of amino acids were obtained immediately after the impact. It was also observed that the reduction of N₂ proceeded according to an associative mechanism, rather than a dissociative mechanism as in the Haber-Bosch process.

Biomolecule Analogues 2-Hydroxypyridine and 2-Pyridone Base Pairing on Ice Nanoparticles

Ice nanoparticles (H2O)N, N ≈ 450 generated in a molecular beam experiment pick up individual gas phase molecules of 2-hydroxypyridine and 2-pyridone (HP) evaporated in a pickup cell at temperatures between 298 and 343 K. The mass spectra of the doped nanoparticles show evidence for generation of clusters of adsorbed molecules (HP)n up to n = 8. The clusters are ionized either by 70 eV electrons or by two photons at 315 nm (3.94 eV). The two ionization methods yield different spectra, and their comparison provides an insight into the neutral cluster composition, ionization and intracluster ion-molecule reactions, and cluster fragmentation. Quite a few molecules were reported not to coagulate on ice nanoparticles previously. The (HP)n cluster generation on ice nanoparticles represents the first evidence for coagulating of molecules and cluster formation on free ice nanoparticles. For comparison, we investigate the coagulation of HP molecules picked up on large clusters ArN, N ≈ 205, and also (HP)n clusters generated in supersonic expansions with Ar buffer gas. This comparison points to a propensity for the (HP)2 dimer generation on ice nanoparticles. This shows the feasibility of base pairing for model of biological molecules on free ice nanoparticles. This result is important for hypotheses of the biomolecule synthesis on ice grains in the space. We support our findings by theoretical calculations that show, among others, the HP dimer structures on water clusters.

Graphene and nano-diamond synthesis in expansions of molten liquid carbon

Despite their widespread use in high-pressure experiments, little is known about the physical and chemical properties of carbon-containing materials as they expand and cool to ambient conditions. As a result, interpretation of experiments can rely on use of unconstrained models with poor accuracy for the ensuing equation of state properties and final chemical products. To this end, we use quantum simulations to study the free expansion and cooling of carbon from metallic liquid states achieved during shock compression. Expansions from three different sets of shock conditions yielded of a variety of chain and ring structures. We then quantify the relative amounts of graphite-like and diamond-like particles formed during cooling and equilibration. We observe that for all cases, graphene sheets are the majority product formed with more extreme initial conditions producing increasingly larger amounts of diamond particles. Our results can address key needs for future meso-scale models of experiments, where knowledge of material properties and chemical end products can have a pronounced effect on interpreting experimental observables.

Miller experiments in atomistic computer simulations

The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.

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.