Electronic transport coefficients from density functional theory across the plasma plane.

We investigate the thermopower and Lorenz number of hydrogen with Kohn-Sham density functional theory (DFT) across the plasma plane toward the near-classical limit, i.e., weakly degenerate and weakly coupled states. Our results are in concordance with certain limiting values for the Lorentz plasma, a model system which only considers electron-ion scattering. Thereby, we clearly show that the widely used method of calculating transport properties via the Kubo-Greenwood (KG) formalism does not capture electron-electron scattering processes. Our discussion also addresses the inadequateness of assuming a Drude-like frequency behavior for the conductivity of nondegenerate plasmas by revisiting the relaxation time approximation within kinetic theory.

Measuring the melting curve of iron at super-Earth core conditions.

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth's interior. We used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth's inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth's mass have the longest dynamos, which provide important shielding against cosmic radiation.

A silicate dynamo in the early Earth.

The Earth's magnetic field has operated for at least 3.4 billion years, yet how the ancient field was produced is still unknown. The core in the early Earth was surrounded by a molten silicate layer, a basal magma ocean that may have survived for more than one billion years. Here we use density functional theory-based molecular dynamics simulations to predict the electrical conductivity of silicate liquid at the conditions of the basal magma ocean: 100-140 GPa, and 4000-6000 K. We find that the electrical conductivity exceeds 10,000 S/m, more than 100 times that measured in silicate liquids at low pressure and temperature. The magnetic Reynolds number computed from our results exceeds the threshold for dynamo activity and the magnetic field strength is similar to that observed in the Archean paleomagnetic record. We therefore conclude that the Archean field was produced by the basal magma ocean.

Comment on "Insulator-metal transition in dense fluid deuterium".

Celliers et al (Reports, 17 August 2018, p. 677), in an attempt to reconcile differences in inferred metallization pressures, provide an alternative temperature analysis of the Knudson et al experiments (Reports, 26 June 2015, p. 1455). We show that this reanalysis implies an anomalously low specific heat for the metallic fluid that is clearly inconsistent with first-principles calculations.

Electrical conductivity of SiO2 at extreme conditions and planetary dynamos.

Ab intio molecular dynamics simulations show that the electrical conductivity of liquid SiO2 is semimetallic at the conditions of the deep molten mantle of early Earth and super-Earths, raising the possibility of silicate dynamos in these bodies. Whereas the electrical conductivity increases uniformly with increasing temperature, it depends nonmonotonically on compression. At very high pressure, the electrical conductivity decreases on compression, opposite to the behavior of many materials. We show that this behavior is caused by a novel compression mechanism: the development of broken charge ordering, and its influence on the electronic band gap.

Density-functional calculations of transport properties in the nondegenerate limit and the role of electron-electron scattering.

We compute electrical and thermal conductivities of hydrogen plasmas in the nondegenerate regime using Kohn-Sham density functional theory (DFT) and an application of the Kubo-Greenwood response formula, and demonstrate that for thermal conductivity, the mean-field treatment of the electron-electron (e-e) interaction therein is insufficient to reproduce the weak-coupling limit obtained by plasma kinetic theories. An explicit e-e scattering correction to the DFT is posited by appealing to Matthiessen's Rule and the results of our computations of conductivities with the quantum Lenard-Balescu (QLB) equation. Further motivation of our correction is provided by an argument arising from the Zubarev quantum kinetic theory approach. Significant emphasis is placed on our efforts to produce properly converged results for plasma transport using Kohn-Sham DFT, so that an accurate assessment of the importance and efficacy of our e-e scattering corrections to the thermal conductivity can be made.

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.

Shock Response and Phase Transitions of MgO at Planetary Impact Conditions.

The moon-forming impact and the subsequent evolution of the proto-Earth is strongly dependent on the properties of materials at the extreme conditions generated by this violent collision. We examine the high pressure behavior of MgO, one of the dominant constituents in Earth's mantle, using high-precision, plate impact shock compression experiments performed on Sandia National Laboratories' Z Machine and extensive quantum calculations using density functional theory (DFT) and quantum Monte Carlo (QMC) methods. The combined data span from ambient conditions to 1.2 TPa and 42 000 K, showing solid-solid and solid-liquid phase boundaries. Furthermore our results indicate that under impact the solid and liquid phases coexist for more than 100 GPa, pushing complete melting to pressures in excess of 600 GPa. The high pressure required for complete shock melting has implications for a broad range of planetary collision events.

Ab initio calculation of the ion feature in x-ray Thomson scattering.

The spectrum of x-ray Thomson scattering is proportional to the dynamic structure factor. An important contribution is the ion feature which describes elastic scattering of x rays off electrons. We apply an ab initio method for the calculation of the form factor of bound electrons, the slope of the screening cloud of free electrons, and the ion-ion structure factor in warm dense beryllium. With the presented method we can calculate the ion feature from first principles. These results will facilitate a better understanding of x-ray scattering in warm dense matter and an accurate measurement of ion temperatures which would allow determining nonequilibrium conditions, e.g., along shock propagation.

Combining Kohn-Sham and orbital-free density-functional theory for Hugoniot calculations to extreme pressures.

The shock Hugoniot for lithium 6 deuteride ((6)LiD) was calculated via first principles using Kohn-Sham density-functional theory molecular dynamics (KSMD) for temperatures of 0.5-25 eV. The upper limit of 25 eV represents a practical limit where KSMD is no longer computationally feasible due to the number of electronic bands which are required to be populated. To push the Hugoniot calculations to higher temperatures we make use of orbital-free density-functional theory molecular dynamics (OFMD). Thomas-Fermi-Dirac-based OFMD gives a poor description of the electronic structure at low temperatures so the initial state is not well defined. We propose a method of bootstrapping the Hugoniot from OFMD to the Hugoniot from KSMD between 10 and 20 eV, where the two methods are in agreement. The combination of KSMD and OFMD allows construction of a first-principles Hugoniot from the initial state to 1000 eV. Theoretical shock-compression results are in good agreement with available experimental data and exhibit the appropriate high-temperature limits. We show that a unified KSMD-OFMD Hugoniot can be used to assess the quality of the existing equation-of-state (EOS) models and inform better EOS models based on justifiable physics.

First-principles calculation of entropy for liquid metals.

We demonstrate the accurate calculation of entropies and free energies for a variety of liquid metals using an extension of the two-phase thermodynamic (2PT) model based on a decomposition of the velocity autocorrelation function into gas-like (hard sphere) and solid-like (harmonic) subsystems. The hard sphere model for the gas-like component is shown to give systematically high entropies for liquid metals as a direct result of the unphysical Lorentzian high-frequency tail. Using a memory function framework we derive a generally applicable velocity autocorrelation and frequency spectrum for the diffusive component which recovers the low-frequency (long-time) behavior of the hard sphere model while providing for realistic short-time coherence and high-frequency tails to the spectrum. This approach provides a significant increase in the accuracy of the calculated entropies for liquid metals and is compared to ambient pressure data for liquid sodium, aluminum, gallium, tin, and iron. The use of this method for the determination of melt boundaries is demonstrated with a calculation of the high-pressure bcc melt boundary for sodium. With the significantly improved accuracy available with the memory function treatment for softer interatomic potentials, the 2PT model for entropy calculations should find broader application in high energy density science, warm dense matter, planetary science, geophysics, and material science.

Complex behavior of fluid lithium under extreme conditions.

Lithium is a prototypical simple metal at standard conditions which is well described within the nearly free electron model. However, by changing the density towards expanded or compressed states, the electrical conductivity shows strong and partly unexpected variations. We have performed quantum molecular dynamics simulations for fluid lithium for a wide range of densities and temperatures in order to derive the equation of state, the electrical conductivity, and information about structural and electronic changes along the expansion or compression. Based on these results, we can give a consistent description of the electrical conductivity from the nonmetallic expanded fluid up to the degenerate electron liquid at high densities.

Quantum molecular dynamics simulations for the nonmetal-to-metal transition in fluid helium.

We have performed quantum molecular dynamics simulations for dense helium to study the nonmetal-to-metal transition at high pressures. We present new results for the equation of state and the Hugoniot curve in the warm dense matter region. The optical conductivity is calculated via the Kubo-Greenwood formula from which the dc conductivity is derived. The nonmetal-to-metal transition is identified at about 1 g/cm(3). We compare with experimental results as well as with other theoretical approaches, especially with predictions of chemical models.

Optimized effective potential from a correlated wave function: optimized effective potential-generalized valence bond (OEP-GVB).

The optimized effective potential (OEP) method allows orbital-dependent functionals to be used in density functional theory. Traditionally the orbital-dependent functional of interest has been the Hartree-Fock energy, leading to exact exchange density functional theory. Here we present results that use a generalized valence-bond (GVB) wave function, a multiconfigurational wave function that includes static correlation and dissociates to the proper limits. We demonstrate the effectiveness of the OEP-GVB method by showing the dissociation of H2 and the excitation spectrum of He.

Phase diagram and electrical conductivity of high energy-density water from density functional theory.

The electrical conductivity and structure of water between 2000-70,000 K and 0.1-3.7 g/cm3 is studied by finite temperature density functional theory (DFT). Proton conduction is investigated quantitatively by analyzing diffusion, the pair-correlation function, and Wannier center locations, while the electronic conduction is calculated in the Kubo-Greenwood formalism. The conductivity formulation is valid across three phase transitions (molecular liquid, ionic liquid, superionic, electronic liquid). Above 100 GPa the superionic phase directly borders an electronically conducting fluid, not an insulating ionic fluid, as previously concluded. For simulations of high energy-density systems to be quantitative, we conclude that finite temperature DFT should be employed.