Bonding and structure in dense multi-component molecular mixtures

Superionic Silica-Water and Silica-Hydrogen Compounds in the Deep Interiors of Uranus and Neptune

Silica, water, and hydrogen are known to be the major components of celestial bodies, and have significant influence on the formation and evolution of giant planets, such as Uranus and Neptune. Thus, it is of fundamental importance to investigate their states and possible reactions under the planetary conditions. Here, using advanced crystal structure searches and first-principles calculations in the Si-O-H system, we find that a silica-water compound (SiO_{2})_{2}(H_{2}O) and a silica-hydrogen compound SiO_{2}H_{2} can exist under high pressures above 450 and 650 GPa, respectively. Further simulations reveal that, at high pressure and high temperature conditions corresponding to the interiors of Uranus and Neptune, these compounds exhibit superionic behavior, in which protons diffuse freely like liquid while the silicon and oxygen framework is fixed as solid. Therefore, these superionic silica-water and silica-hydrogen compounds could be regarded as important components of the deep mantle or core of giants, which also provides an alternative origin for their anomalous magnetic fields. These unexpected physical and chemical properties of the most common natural materials at high pressure offer key clues to understand some abstruse issues including demixing and erosion of the core in giant planets, and shed light on building reliable models for solar giants and exoplanets.

First principles reactive simulation for equation of state prediction

The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman-Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. We also demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.

Coexistence of plastic and partially diffusive phases in a helium-methane compound

Helium and methane are major components of giant icy planets and are abundant in the universe. However, helium is the most inert element in the periodic table and methane is one of the most hydrophobic molecules, thus whether they can react with each other is of fundamental importance. Here, our crystal structure searches and first-principles calculations predict that a He3CH4 compound is stable over a wide range of pressures from 55 to 155 GPa and a HeCH4 compound becomes stable around 105 GPa. As nice examples of pure van der Waals crystals, the insertion of helium atoms changes the original packing of pure methane molecules and also largely hinders the polymerization of methane at higher pressures. After analyzing the diffusive properties during the melting of He3CH4 at high pressure and high temperature, in addition to a plastic methane phase, we have discovered an unusual phase which exhibits coexistence of diffusive helium and plastic methane. In addition, the range of the diffusive behavior within the helium-methane phase diagram is found to be much narrower compared to that of previously predicted helium-water compounds. This may be due to the weaker van der Waals interactions between methane molecules compared to those in helium-water compounds, and that the helium-methane compound melts more easily.

Benchmarking the effective one-component plasma model for warm dense neon and krypton within quantum molecular dynamics simulation

The effective one-component plasma (EOCP) model has provided an efficient approach to obtaining many important thermophysical parameters of hot dense matter [J. Clérouin, et al., Phys. Rev. Lett. 116, 115003 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.115003]. In this paper, we perform extensive quantum molecular dynamics (QMD) simulations to determine the equations of state, ionic structures, and ionic transport properties of neon and krypton within the warm dense matter (WDM) regime where the density (ρ) is up to 12 g/cm^{3} and the temperature (T) is up to 100 kK. The simulated data are then used as a benchmark to explicitly evaluate the EOCP and Yukawa models. It is found that, within present ρ-T regime, the EOCP model can excellently reproduce the diffusion and viscosity coefficients of neon and krypton due to the fact that this model defines a system which nearly reproduces the actual physical states of WDM. Therefore, the EOCP model may be a promising alternative approach to reasonably predicting the transport behaviors of matter in WDM regime at lower QMD computational cost. The evaluation of Yukawa model shows that the consideration of the energy level broadening effect in the average atom model is necessary. Finally, with the help of EOCP model, the Stokes-Einstein relationships about neon and krypton are discussed, and fruitful plasma parameters as well as a practical ρ-T-dependent formula of the effective coupling parameter are obtained. These results not only provide valuable information for future theoretical and experimental studies on dense neon and krypton but also reveal the applicability of the EOCP model and the limitation of the Yukawa model in WDM regime and further support the continuing search for a unified description of ionic transport in dense plasma.

Multicomponent mutual diffusion in the warm, dense matter regime

We present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.

Laser-driven shock compression of "synthetic planetary mixtures" of water, ethanol, and ammonia

Water, methane, and ammonia are commonly considered to be the key components of the interiors of Uranus and Neptune. Modelling the planets' internal structure, evolution, and dynamo heavily relies on the properties of the complex mixtures with uncertain exact composition in their deep interiors. Therefore, characterising icy mixtures with varying composition at planetary conditions of several hundred gigapascal and a few thousand Kelvin is crucial to improve our understanding of the ice giants. In this work, pure water, a water-ethanol mixture, and a water-ethanol-ammonia "synthetic planetary mixture" (SPM) have been compressed through laser-driven decaying shocks along their principal Hugoniot curves up to 270, 280, and 260 GPa, respectively. Measured temperatures spanned from 4000 to 25000 K, just above the coldest predicted adiabatic Uranus and Neptune profiles (3000-4000 K) but more similar to those predicted by more recent models including a thermal boundary layer (7000-14000 K). The experiments were performed at the GEKKO XII and LULI2000 laser facilities using standard optical diagnostics (Doppler velocimetry and optical pyrometry) to measure the thermodynamic state and the shock-front reflectivity at two different wavelengths. The results show that water and the mixtures undergo a similar compression path under single shock loading in agreement with Density Functional Theory Molecular Dynamics (DFT-MD) calculations using the Linear Mixing Approximation (LMA). On the contrary, their shock-front reflectivities behave differently by what concerns both the onset pressures and the saturation values, with possible impact on planetary dynamos.

Structural determination of neutral Co n clusters (n = 4-10,13) through IR-UV two-color vibrational spectroscopy and DFT calculations

We recorded IR spectra for neutral cobalt clusters via two-color IR-UV ionization, using the Free Electron Laser for intracavity experiments (FELICE). Well-resolved IR spectra are presented for [Formula: see text] (n  =  4-10, 13) and analyzed with the help of Density Functional Theory calculations using two different correlation exchange functionals: the revisited Tao-Perdew-Staroverov-Scuseria (revTPSS) and the frequently used Perdew-Burke-Ernzerhof (PBE) approaches. Although we have not performed an extensive structure search, we tentatively assign the spectra for all cluster sizes except for n  =  7, and n  =  10. We find that neither of the two functionals chosen clearly outperforms the other in predicting IR spectra, and that relatively low scaling factors of 0.82 (PBE) and 0.8 (revTPSS) are required. In contrast to the magnetic moments, the calculated electric dipole moments fluctuate strongly as a function of cluster size and could therefore be used as an indirect probe to the cluster structure.

Quantum molecular dynamics study on the proton exchange, ionic structures, and transport properties of warm dense hydrogen-deuterium mixtures

Comprehensive knowledge of physical properties such as equation of state (EOS), proton exchange, dynamic structures, diffusion coefficients, and viscosities of hydrogen-deuterium mixtures with densities from 0.1 to 5 g/cm^{3} and temperatures from 1 to 50 kK has been presented via quantum molecular dynamics (QMD) simulations. The existing multi-shock experimental EOS provides an important benchmark to evaluate exchange-correlation functionals. The comparison of simulations with experiments indicates that a nonlocal van der Waals density functional (vdW-DF1) produces excellent results. Fraction analysis of molecules using a weighted integral over pair distribution functions was performed. A dissociation diagram together with a boundary where the proton exchange (H_{2}+D_{2}⇌2HD) occurs was generated, which shows evidence that the HD molecules form as the H_{2} and D_{2} molecules are almost 50% dissociated. The mechanism of proton exchange can be interpreted as a process of dissociation followed by recombination. The ionic structures at extreme conditions were analyzed by the effective coordination number model. High-order cluster, circle, and chain structures can be founded in the strongly coupled warm dense regime. The present QMD diffusion coefficient and viscosity can be used to benchmark two analytical one-component plasma (OCP) models: the Coulomb and Yukawa OCP models.

Stabilization of ammonia-rich hydrate inside icy planets

The interior structure of the giant ice planets Uranus and Neptune, but also of newly discovered exoplanets, is loosely constrained, because limited observational data can be satisfied with various interior models. Although it is known that their mantles comprise large amounts of water, ammonia, and methane ices, it is unclear how these organize themselves within the planets-as homogeneous mixtures, with continuous concentration gradients, or as well-separated layers of specific composition. While individual ices have been studied in great detail under pressure, the properties of their mixtures are much less explored. We show here, using first-principles calculations, that the 2:1 ammonia hydrate, (H₂O)(NH₃)₂, is stabilized at icy planet mantle conditions due to a remarkable structural evolution. Above 65 GPa, we predict it will transform from a hydrogen-bonded molecular solid into a fully ionic phase O^2-([Formula: see text])₂, where all water molecules are completely deprotonated, an unexpected bonding phenomenon not seen before. Ammonia hemihydrate is stable in a sequence of ionic phases up to 500 GPa, pressures found deep within Neptune-like planets, and thus at higher pressures than any other ammonia-water mixture. This suggests it precipitates out of any ammonia-water mixture at sufficiently high pressures and thus forms an important component of icy planets.

Alternative first-principles calculation of entropy for liquids

We present an alternative method for interpreting the velocity autocorrelation function (VACF) of a fluid with application to extracting the entropy in a manner similar to the methods developed by Lin et al. [J. Chem. Phys. 119, 11792 (2003)]JCPSA60021-960610.1063/1.1624057 and improved upon by Desjarlais [Phys. Rev. E 88, 062145 (2013)]PLEEE81539-375510.1103/PhysRevE.88.062145. The liquid VACF is decomposed into two components, one gas and one solid, and each contribution's entropic portion is calculated. However, we fit both the gas and solid portions of the VACF in the time domain. This approach is applied to a single-component liquid (a two-phase model of liquid Al at the melt line) and two different two-component systems: a superionic-to-superionic (bcc to fcc) phase transition in H_{2}O at high temperatures and pressures and a metastable liquid state of MgO. For all three examples, comparisons to existing results in the literature demonstrate the validity of our alternative.

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.

Superionic Phases of the 1:1 Water-Ammonia Mixture

We report four structures for the 1:1 water-ammonia mixture showing superionic behavior at high temperature with the space groups P4/nmm, Ima2, Pma2, and Pm, which have been identified from evolutionary random structure search calculations at 0 K. Analyzing the respective pair distribution functions and diffusive properties the superionic phase is found to be stable in a temperature range between 1000 and 6000 K for pressures up to 800 GPa. We propose a high-pressure phase diagram of the water-ammonia mixture for the first time and compare the self-diffusion coefficients in the mixture to the ones found in water and ammonia. Finally, possible implications on the interior structure of the giant planets Uranus and Neptune are discussed.