Growth model interpretation of planet size distribution.

The radii and orbital periods of 4,000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. However, an important question about the composition of planets ranging from 2 to 4 Earth radii (R⊕) still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multicomponent, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 M⊕, if half-ice and half-rock by mass, have radii of 2.5 R⊕, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2- to 4-R⊕ range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are "water worlds."

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.

Shock compression of a fifth period element: liquid xenon to 840 GPa.

Current equation of state (EOS) models for xenon show substantial differences in the Hugoniot above 100 GPa, prompting the need for an improved understanding of xenon's behavior at extreme conditions. We performed shock compression experiments on liquid xenon to determine the Hugoniot up to 840 GPa, using these results to validate density functional theory (DFT) simulations. Despite the nearly fivefold compression, we find that the limiting Thomas-Fermi theory, exact in the high density limit, does not accurately describe the system. Combining the experimental data and DFT calculations, we developed a free-energy-based, multiphase EOS capable of describing xenon over a wide range of pressures and temperatures.

AM05 Density Functional Applied to the Water Molecule, Dimer, and Bulk Liquid.

We compare results for water obtained with the AM05 exchange-correlation density functional ( Armiento, R.; Mattsson, A. E. Phys. Rev. B 2005, 72, 085108 ) with those obtained with five other pure functionals: LDA, PBE, PBEsol, RPBE, and BLYP. For liquid water, AM05 yields an O-O pair correlation function that is more structured than the ones of PBE and BLYP, which, in turn, are more structured than the one of RPBE. However, LDA and PBEsol yields more structured water than AM05. We show that AM05 yields a H2O dimer binding energy of 4.9 kcal/mol. The result is thus within 0.15 kcal/mol of CCSD(T) level theory (5.02 ± 0.05 kcal/mol). We confirm that accuracy in the water dimer binding energy is not a strong indicator for the fidelity of the resulting structure of liquid water.

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.

On "the complete basis set limit" and plane-wave methods in first-principles simulations of water.

Water structure, measured by the height of the first peak in oxygen-oxygen radial distributions, is converged with respect to plane-wave basis energy cutoffs for ab initio molecular dynamics simulations, confirming the reliability of plane-wave methods.

The AM05 density functional applied to solids.

We show that the AM05 functional [Armiento and Mattsson, Phys. Rev. B 72, 085108 (2005)] has the same excellent performance for solids as the hybrid density functionals tested in Paier et al. [J. Chem. Phys. 124, 154709 (2006); 125, 249901 (2006)]. This confirms the original finding that AM05 performs exceptionally well for solids and surfaces. Hartree-Fock hybrid calculations are typically an order of magnitude slower than local or semilocal density functionals such as AM05, which is of a regular semilocal generalized gradient approximation form. The performance of AM05 is on average found to be superior to selecting the best of local density approximation and PBE for each solid. By comparing data from several different electronic-structure codes, we have determined that the numerical errors in this study are equal to or smaller than the corresponding experimental uncertainties.

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.

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.

Self-diffusion rates in Al from combined first-principles and model-potential calculations.

Monovacancy diffusion alone dominates over diffusion due to divacancies and interstitials in Al for all temperatures up to the melting point. Deviations from a single Arrhenius dependence are due to anharmonicity. The conclusion is based on a combination of theoretical methods, from density functional theory to thermodynamic integration, without fitting to experimental data. The calculated diffusion rate agrees with experimental data over 11 orders of magnitude.