Thermal dependence of the hydrated proton and optimal proton transfer in the protonated water hexamer.

Water is a key ingredient for life and plays a central role as solvent in many biochemical reactions. However, the intrinsically quantum nature of the hydrogen nucleus, revealing itself in a large variety of physical manifestations, including proton transfer, gives rise to unexpected phenomena whose description is still elusive. Here we study, by a combination of state-of-the-art quantum Monte Carlo methods and path-integral molecular dynamics, the structure and hydrogen-bond dynamics of the protonated water hexamer, the fundamental unit for the hydrated proton. We report a remarkably low thermal expansion of the hydrogen bond from zero temperature up to 300 K, owing to the presence of short-Zundel configurations, characterised by proton delocalisation and favoured by the synergy of nuclear quantum effects and thermal activation. The hydrogen bond strength progressively weakens above 300 K, when localised Eigen-like configurations become relevant. Our analysis, supported by the instanton statistics of shuttling protons, reveals that the near-room-temperature range from 250 K to 300 K is optimal for proton transfer in the protonated water hexamer.

The trisulfur radical ion S3•- controls platinum transport by hydrothermal fluids.

Platinum group elements (PGE) are considered to be very poorly soluble in aqueous fluids in most natural hydrothermal-magmatic contexts and industrial processes. Here, we combined in situ X-ray absorption spectroscopy and solubility experiments with atomistic and thermodynamic simulations to demonstrate that the trisulfur radical ion S3•- forms very stable and soluble complexes with both PtII and PtIV in sulfur-bearing aqueous solution at elevated temperatures (∼300 °C). These Pt-bearing species enable (re)mobilization, transfer, and focused precipitation of platinum up to 10,000 times more efficiently than any other common inorganic ligand, such as hydroxide, chloride, sulfate, or sulfide. Our results imply a far more important contribution of sulfur-bearing hydrothermal fluids to PGE transfer and accumulation in the Earth's crust than believed previously. This discovery challenges traditional models of PGE economic concentration from silicate and sulfide melts and provides new possibilities for resource prospecting in hydrothermal shallow crust settings. The exceptionally high capacity of the S3•- ion to bind platinum may also offer new routes for PGE selective extraction from ore and hydrothermal synthesis of noble metal nanomaterials.

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.

Formation of nucleobases in a Miller-Urey reducing atmosphere.

The Miller-Urey experiments pioneered modern research on the molecular origins of life, but their actual relevance in this field was later questioned because the gas mixture used in their research is considered too reducing with respect to the most accepted hypotheses for the conditions on primordial Earth. In particular, the production of only amino acids has been taken as evidence of the limited relevance of the results. Here, we report an experimental work, combined with state-of-the-art computational methods, in which both electric discharge and laser-driven plasma impact simulations were carried out in a reducing atmosphere containing NH3 + CO. We show that RNA nucleobases are synthesized in these experiments, strongly supporting the possibility of the emergence of biologically relevant molecules in a reducing atmosphere. The reconstructed synthetic pathways indicate that small radicals and formamide play a crucial role, in agreement with a number of recent experimental and theoretical results.

Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics.

We introduce a novel approach for a fully quantum description of coupled electron-ion systems from first principles. It combines the variational quantum Monte Carlo solution of the electronic part with the path integral formalism for the quantum nuclear dynamics. On the one hand, the path integral molecular dynamics includes nuclear quantum effects by adding a set of fictitious classical particles (beads) aimed at reproducing nuclear quantum fluctuations via a harmonic kinetic term. On the other hand, variational quantum Monte Carlo can provide Born-Oppenheimer potential energy surfaces with a precision comparable to the most-advanced post-Hartree-Fock approaches, and with a favorable scaling with the system size. In order to cope with the intrinsic noise due to the stochastic nature of quantum Monte Carlo methods, we generalize the path integral molecular dynamics using a Langevin thermostat correlated according to the covariance matrix of quantum Monte Carlo nuclear forces. The variational parameters of the quantum Monte Carlo wave function are evolved during the nuclear dynamics, such that the Born-Oppenheimer potential energy surface is unbiased. Statistical errors on the wave function parameters are reduced by resorting to bead grouping average, which we show to be accurate and well-controlled. Our general algorithm relies on a Trotter breakup between the dynamics driven by ionic forces and the one set by the harmonic interbead couplings. The latter is exactly integrated, even in the presence of the Langevin thermostat, thanks to the mapping onto an Ornstein-Uhlenbeck process. This framework turns out to be also very efficient in the case of noiseless (deterministic) ionic forces. The new implementation is validated on the Zundel ion (H5O2+) by direct comparison with standard path integral Langevin dynamics calculations made with a coupled cluster potential energy surface. Nuclear quantum effects are confirmed to be dominant over thermal effects well beyond room temperature, giving the excess proton an increased mobility by quantum tunneling.

Formamide reaction network in gas phase and solution via a unified theoretical approach: Toward a reconciliation of different prebiotic scenarios.

Increasing experimental and theoretical evidence points to formamide as a possible hub in the complex network of prebiotic chemical reactions leading from simple precursors like H2, H2O, N2, NH3, CO, and CO2 to key biological molecules like proteins, nucleic acids, and sugars. We present an in-depth computational study of the formation and decomposition reaction channels of formamide by means of ab initio molecular dynamics. To this aim we introduce a new theoretical method combining the metadynamics sampling scheme with a general purpose topological formulation of collective variables able to track a wide range of different reaction mechanisms. Our approach is flexible enough to discover multiple pathways and intermediates starting from minimal insight on the systems, and it allows passing in a seamless way from reactions in gas phase to reactions in liquid phase, with the solvent active role fully taken into account. We obtain crucial new insight into the interplay of the different formamide reaction channels and into environment effects on pathways and barriers. In particular, our results indicate a similar stability of formamide and hydrogen cyanide in solution as well as their relatively facile interconversion, thus reconciling experiments and theory and, possibly, two different and competing prebiotic scenarios. Moreover, although not explicitly sought, formic acid/ammonium formate is produced as an important formamide decomposition byproduct in solution.

Sulfur radical species form gold deposits on Earth.

Current models of the formation and distribution of gold deposits on Earth are based on the long-standing paradigm that hydrogen sulfide and chloride are the ligands responsible for gold mobilization and precipitation by fluids across the lithosphere. Here we challenge this view by demonstrating, using in situ X-ray absorption spectroscopy and solubility measurements, coupled with molecular dynamics and thermodynamic simulations, that sulfur radical species, such as the trisulfur ion S3(-), form very stable and soluble complexes with Au(+) in aqueous solution at elevated temperatures (>250 °C) and pressures (>100 bar). These species enable extraction, transport, and focused precipitation of gold by sulfur-rich fluids 10-100 times more efficiently than sulfide and chloride only. As a result, S3(-) exerts an important control on the source, concentration, and distribution of gold in its major economic deposits from magmatic, hydrothermal, and metamorphic settings. The growth and decay of S3(-) during the fluid generation and evolution is one of the key factors that determine the fate of gold in the lithosphere.

Strong electric fields at a prototypical oxide/water interface probed by ab initio molecular dynamics: MgO(001).

We report a density-functional theory (DFT)-based study of the interface of bulk water with a prototypical oxide surface, MgO(001), and focus our study on the often-overlooked surface electric field. In particular, we observe that the bare MgO(001) surface, although charge-neutral and defectless, has an intense electric field on the Å scale. The MgO(001) surface covered with 1 water monolayer (1 ML) is investigated via a supercell accounting for the experimentally-observed (2 × 3) reconstruction, stable at ambient temperature, and in which two out of six water molecules are dissociated. This 1 ML-hydrated surface is also found to have a high, albeit short-ranged, normal component of the field. Finally, the oxide/water interface is studied via room-temperature ab initio molecular dynamics (AIMD) using 34 H2O molecules between two MgO(001) surfaces. To our best knowledge this is the first AIMD study of the MgO(001)/liquid water interface in which all atoms are treated using DFT and including several layers above the first adsorbed layer. We observe that the surface electric field, averaged over the AIMD trajectories, is still very strong on the fully-wet surface, peaking at about 3 V Å(-1). Even in the presence of bulk-like water, the structure of the first layer in contact with the surface remains similar to the (2 × 3)-reconstructed ice ad-layer on MgO(001). Moreover, we observe proton exchange within the first layer, and between the first and second layers - indeed, the O-O distances close to the surface are found to be distributed towards shorter distances, a property which has been shown to directly promote proton transfer.

Effect of salt on the H-bond symmetrization in ice.

The richness of the phase diagram of water reduces drastically at very high pressures where only two molecular phases, proton-disordered ice VII and proton-ordered ice VIII, are known. Both phases transform to the centered hydrogen bond atomic phase ice X above about 60 GPa, i.e., at pressures experienced in the interior of large ice bodies in the universe, such as Saturn and Neptune, where nonmolecular ice is thought to be the most abundant phase of water. In this work, we investigate, by Raman spectroscopy up to megabar pressures and ab initio simulations, how the transformation of ice VII in ice X is affected by the presence of salt inclusions in the ice lattice. Considerable amounts of salt can be included in ice VII structure under pressure via rock-ice interaction at depth and processes occurring during planetary accretion. Our study reveals that the presence of salt hinders proton order and hydrogen bond symmetrization, and pushes ice VII to ice X transformation to higher and higher pressures as the concentration of salt is increased.

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.

Six-fold-coordinated phosphorus by oxygen in AlPO4 quartz homeotype under high pressure.

AlPO4 belongs to the berlinite quartz homeotype family, which has been the subject of intense high-pressure research triggered by the supposed existence of reversible pressure-induced amorphization. X-ray diffraction experiments, complemented with ab initio calculations, demonstrate the existence of two high-pressure crystalline polymorphs and show that AlPO4 shares the same two-stage densification mechanism as silica. In the first step, a compact hexagonal sublattice of oxygen atoms is formed. In the second step, the cations redistribute in the interstices giving rise to a monoclinic distorted CaCl2 phase. The most outstanding feature of the phase is that phosphorous becomes six-fold coordinated by oxygen, adopting a configuration unknown so far in solid-state science. This finding opens possibilities in the high-pressure chemistry of phosphorus. The close relationship of AlPO4 with silica suggests the existence of completely unexplored families of compounds analogous to those of six-fold-coordinated silicates but based on PO6.