Dynamic compression of water to conditions in ice giant interiors.

Recent discoveries of water-rich Neptune-like exoplanets require a more detailed understanding of the phase diagram of H2O at pressure-temperature conditions relevant to their planetary interiors. The unusual non-dipolar magnetic fields of ice giant planets, produced by convecting liquid ionic water, are influenced by exotic high-pressure states of H2O-yet the structure of ice in this state is challenging to determine experimentally. Here we present X-ray diffraction evidence of a body-centered cubic (BCC) structured H2O ice at 200 GPa and ~ 5000 K, deemed ice XIX, using the X-ray Free Electron Laser of the Linac Coherent Light Source to probe the structure of the oxygen sub-lattice during dynamic compression. Although several cubic or orthorhombic structures have been predicted to be the stable structure at these conditions, we show this BCC ice phase is stable to multi-Mbar pressures and temperatures near the melt boundary. This suggests variable and increased electrical conductivity to greater depths in ice giant planets that may promote the generation of multipolar magnetic fields.

Measurement and implications of Saturn's gravity field and ring mass.

The interior structure of Saturn, the depth of its winds, and the mass and age of its rings constrain its formation and evolution. In the final phase of the Cassini mission, the spacecraft dived between the planet and its innermost ring, at altitudes of 2600 to 3900 kilometers above the cloud tops. During six of these crossings, a radio link with Earth was monitored to determine the gravitational field of the planet and the mass of its rings. We find that Saturn's gravity deviates from theoretical expectations and requires differential rotation of the atmosphere extending to a depth of at least 9000 kilometers. The total mass of the rings is (1.54 ± 0.49) × 1019 kilograms (0.41 ± 0.13 times that of the moon Mimas), indicating that the rings may have formed 107 to 108 years ago.

Path integral Monte Carlo simulations of warm dense aluminum.

We perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of 0.1-32.4gcm^{-3} and 10^{4}-10^{8} K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multishock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multishock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and 〈Z〉 to reveal the evolution of the plasma structure and ionization behavior.

Jupiter's atmospheric jet streams extend thousands of kilometres deep.

The depth to which Jupiter's observed east-west jet streams extend has been a long-standing question. Resolving this puzzle has been a primary goal for the Juno spacecraft, which has been in orbit around the gas giant since July 2016. Juno's gravitational measurements have revealed that Jupiter's gravitational field is north-south asymmetric, which is a signature of the planet's atmospheric and interior flows. Here we report that the measured odd gravitational harmonics J3, J5, J7 and J9 indicate that the observed jet streams, as they appear at the cloud level, extend down to depths of thousands of kilometres beneath the cloud level, probably to the region of magnetic dissipation at a depth of about 3,000 kilometres. By inverting the measured gravity values into a wind field, we calculate the most likely vertical profile of the deep atmospheric and interior flow, and the latitudinal dependence of its depth. Furthermore, the even gravity harmonics J8 and J10 resulting from this flow profile also match the measurements, when taking into account the contribution of the interior structure. These results indicate that the mass of the dynamical atmosphere is about one per cent of Jupiter's total mass.

A suppression of differential rotation in Jupiter's deep interior.

Jupiter's atmosphere is rotating differentially, with zones and belts rotating at speeds that differ by up to 100 metres per second. Whether this is also true of the gas giant's interior has been unknown, limiting our ability to probe the structure and composition of the planet. The discovery by the Juno spacecraft that Jupiter's gravity field is north-south asymmetric and the determination of its non-zero odd gravitational harmonics J3, J5, J7 and J9 demonstrates that the observed zonal cloud flow must persist to a depth of about 3,000 kilometres from the cloud tops. Here we report an analysis of Jupiter's even gravitational harmonics J4, J6, J8 and J10 as observed by Juno and compared to the predictions of interior models. We find that the deep interior of the planet rotates nearly as a rigid body, with differential rotation decreasing by at least an order of magnitude compared to the atmosphere. Moreover, we find that the atmospheric zonal flow extends to more than 2,000 kilometres and to less than 3,500 kilometres, making it fully consistent with the constraints obtained independently from the odd gravitational harmonics. This depth corresponds to the point at which the electric conductivity becomes large and magnetic drag should suppress differential rotation. Given that electric conductivity is dependent on planetary mass, we expect the outer, differentially rotating region to be at least three times deeper in Saturn and to be shallower in massive giant planets and brown dwarfs.

Measurement of Jupiter's asymmetric gravity field.

The gravity harmonics of a fluid, rotating planet can be decomposed into static components arising from solid-body rotation and dynamic components arising from flows. In the absence of internal dynamics, the gravity field is axially and hemispherically symmetric and is dominated by even zonal gravity harmonics J2n that are approximately proportional to qn, where q is the ratio between centrifugal acceleration and gravity at the planet's equator. Any asymmetry in the gravity field is attributed to differential rotation and deep atmospheric flows. The odd harmonics, J3, J5, J7, J9 and higher, are a measure of the depth of the winds in the different zones of the atmosphere. Here we report measurements of Jupiter's gravity harmonics (both even and odd) through precise Doppler tracking of the Juno spacecraft in its polar orbit around Jupiter. We find a north-south asymmetry, which is a signature of atmospheric and interior flows. Analysis of the harmonics, described in two accompanying papers, provides the vertical profile of the winds and precise constraints for the depth of Jupiter's dynamical atmosphere.

First-principles simulations of warm dense lithium fluoride.

We perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of LiF. Our simulations cover a wide density-temperature range of 2.08-15.70gcm^{-3} and 10^{4}-10^{9} K. Since PIMC and DFT-MD accurately treat effects of atomic shell structure, we find a pronounced compression maximum and a shoulder on the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, LEOS, and models. In addition, we compute pair-correlation functions that reveal an evolving plasma structure and ionization process that is driven by thermal and pressure ionization. Finally, we compute electronic density of states of liquid LiF from DFT-MD simulations and find that the electronic gap can remain open with increasing density and temperature to at least 15.7 gcm^{-3}.

First-principles equation of state and electronic properties of warm dense oxygen.

We perform all-electron path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of oxygen. Our simulations cover a wide density-temperature range of 1-100 g cm(-3) and 10(4)-10(9) K. By combining results from PIMC and DFT-MD, we are able to compute pressures and internal energies from first-principles at all temperatures and provide a coherent equation of state. We compare our first-principles calculations with analytic equations of state, which tend to agree for temperatures above 8 × 10(6) K. Pair-correlation functions and the electronic density of states reveal an evolving plasma structure and ionization process that is driven by temperature and density. As we increase the density at constant temperature, we find that the ionization fraction of the 1s state decreases while the other electronic states move towards the continuum. Finally, the computed shock Hugoniot curves show an increase in compression as the first and second shells are ionized.

All-electron path integral Monte Carlo simulations of warm dense matter: application to water and carbon plasmas.

We develop an all-electron path integral Monte Carlo method with free-particle nodes for warm dense matter and apply it to water and carbon plasmas. We thereby extend path integral Monte Carlo studies beyond hydrogen and helium to elements with core electrons. Path integral Monte Carlo results for pressures, internal energies, and pair-correlation functions compare well with density functional theory molecular dynamics calculations at temperatures of (2.5-7.5)×10(5) K, and both methods together form a coherent equation of state over a density-temperature range of 3-12 g/cm(3) and 10(4)-10(9) K.

Strong coupling and degeneracy effects in inertial confinement fusion implosions.

Accurate knowledge about the equation of state (EOS) of deuterium is critical to inertial confinement fusion (ICF). Low-adiabat ICF implosions routinely access strongly coupled and degenerate plasma conditions. Using the path integral Monte Carlo method, we have derived a first-principles EOS (FPEOS) table of deuterium. It is the first ab initio EOS table which completely covers typical ICF implosion trajectory in the density and temperature ranges of ρ=0.002-1596 g/cm3 and T=1.35 eV-5.5 keV. Discrepancies in internal energy and pressure have been found in strongly coupled and degenerate regimes with respect to SESAME EOS. Hydrodynamics simulations of cryogenic ICF implosions using the FPEOS table have indicated significant differences in peak density, areal density (ρR), and neutron yield relative to SESAME simulations.

Fundamental high-pressure calibration from all-electron quantum Monte Carlo calculations.

We develop an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and use it to construct a primary ultra-high-pressure calibration based on the equation of state of cubic boron nitride. We compute the static contribution to the free energy with the QMC method and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. We compute the anharmonic Raman frequency shift with QMC simulations as a function of pressure and temperature, allowing optical pressure calibration. In contrast to present experimental approaches, small systematic errors in the theoretical EOS do not increase with pressure, and no extrapolation is needed. This all-electron method is applicable to first-row solids, providing a new reference for ab initio calculations of solids and benchmarks for pseudopotential accuracy.

Quantum Monte Carlo computations of phase stability, equations of state, and elasticity of high-pressure silica.

Silica (SiO(2)) is an abundant component of the Earth whose crystalline polymorphs play key roles in its structure and dynamics. First principle density functional theory (DFT) methods have often been used to accurately predict properties of silicates, but fundamental failures occur. Such failures occur even in silica, the simplest silicate, and understanding pure silica is a prerequisite to understanding the rocky part of the Earth. Here, we study silica with quantum Monte Carlo (QMC), which until now was not computationally possible for such complex materials, and find that QMC overcomes the failures of DFT. QMC is a benchmark method that does not rely on density functionals but rather explicitly treats the electrons and their interactions via a stochastic solution of Schrödinger's equation. Using ground-state QMC plus phonons within the quasiharmonic approximation of density functional perturbation theory, we obtain the thermal pressure and equations of state of silica phases up to Earth's core-mantle boundary. Our results provide the best constrained equations of state and phase boundaries available for silica. QMC indicates a transition to the dense alpha-PbO(2) structure above the core-insulating D" layer, but the absence of a seismic signature suggests the transition does not contribute significantly to global seismic discontinuities in the lower mantle. However, the transition could still provide seismic signals from deeply subducted oceanic crust. We also find an accurate shear elastic constant for stishovite and its geophysically important softening with pressure.

First-principles studies of the metallization and the equation of state of solid helium.

The insulator-to-metal transition in hcp solid helium at high pressure is studied with first-principles simulations. Diffusion quantum Monte Carlo (DMC) calculations predict that the band gap closes at a density of 21.3 g/cm;{3} and a pressure of 25.7 terapascals, which is 20% higher in density and 40% higher in pressure than predicted by density functional calculations based on the generalized gradient approximation (GGA). The metallization density derived from GW calculations is found to be in very close agreement with DMC predictions. The zero point motion of the nuclei had no effect on the metallization density within the accuracy of the calculation. Finally, fit functions for the equation of state are presented, and the magnitude of the electronic correlation effects left out of the GGA approximation are discussed.

First principles calculations of shock compressed fluid helium.

The properties of hot dense helium at megabar pressures are studied with two first principles computer simulation techniques: path integral Monte Carlo simulation and density functional molecular dynamics. The simulations predict that the compressibility of helium is substantially increased by electronic excitations that are present in the hot fluid at thermodynamic equilibrium. A maximum compression ratio of 5.24(4)-fold the initial density was predicted for 360 GPa and 150,000 K. This result distinguishes helium from deuterium, for which simulations predicted a maximum compression ratio of 4.3(1). Hugoniot curves for statically precompressed samples are also discussed.

Dense plasma effects on nuclear reaction rates.

Plasma density effects can cause an exponential change in charged particle nuclear reaction rates important in stellar evolution. Reaction rates in dense plasma, with emphasis on quantum aspects, are examined here using path integral Monte Carlo calculations. Quantum mechanics causes a reduction in the many body enhancement of the reaction rate compared to the value for a classical system. This can be attributed to the "quantum smearing" of the short range Coulomb interaction resulting in reduced repulsion between the reacting pair and surrounding particles. Electron screening and ion exchange effects are also examined, with screening reducing and exchange slightly increasing the many body enhancement.

Path integral Monte Carlo simulation of the low-density hydrogen plasma.

Restricted path integral Monte Carlo simulations are used to calculate the equilibrium properties of hydrogen in the density and temperature range of 9.83 x 10(-4)

Calculation of a deuterium double shock Hugoniot from ab initio simulations.

We calculate the equation of state of dense deuterium with two ab initio simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, in the density range of 0.67 < or = rho < or = 1.60 g cm(-3). We derive the double shock Hugoniot and compare with the recent laser-driven double shock wave experiments by Mostovych et al. [Phys. Rev. Lett. 85, 3870 (2000)]. We find excellent agreement between the two types of microscopic simulations, but a significant discrepancy with the laser-driven shock measurements.

Variational density matrix method for warm, condensed matter: application to dense hydrogen

A variational principle for optimizing thermal density matrices is introduced. As a first application, the variational many-body density matrix is written as a determinant of one-body density matrices, which are approximated by Gaussians with the mean, width, and amplitude as variational parameters. The method is illustrated for the particle in an external field problem, the hydrogen molecule and dense hydrogen where the molecular, the dissociated, and the plasma regime are described. Structural and thermodynamic properties (energy, equation of state, and shock Hugoniot) are presented.

Path integral monte carlo calculation of the deuterium hugoniot

Restricted path integral Monte Carlo simulations have been used to calculate the equilibrium properties of deuterium for two densities: 0.674 and 0.838 g cm(-3) ( r(s) = 2.00 and 1.86) in the temperature range of 10(5)