BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Burkhard Militzer

Burkhard Militzer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPcond-mat.mtrl-sciphysics.plasm-phcond-mat.othercond-mat.stat-mechphysics.comp-phphysics.geo-ph

Catalog footprint

What is connected

31works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

Revelations on Jupiter's Formation, Evolution and Interior: Challenges from Juno Results

The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, with a complex internal structure. These new results imply that current models of Jupiter's formation and evolution require a revision. In this paper, we discuss potential formation and evolution paths that can lead to an internal structure model consistent with Juno data, and the constraints they provide. We note that standard core accretion formation models, including the heavy-element enrichment during planetary growth is consistent with an interior that is inhomogeneous with composition gradients in its deep interior. However, such formation models typically predict that this region, which could be interpreted as a primordial dilute core, is confined to about 10% of Jupiter's total mass. In contrast, structure models that fit Juno data imply that this region contains 30% of the mass or more. One way to explain the origin of this extended region is by invoking a relatively long (about 2 Myrs) formation phase where the growing planet accretes gas and planetesimals delaying the runaway gas accretion. Alternatively, Jupiter's fuzzy core could be a result of a giant impact or convection post-formation. These novel scenarios require somewhat special and specific conditions. Clarity on the plausibility of such conditions could come from future high-resolution observations of planet-forming regions around other stars, from the observed and modeled architectures of extrasolar systems with giant planets, and future Juno data obtained during its extended mission.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Probing the icy shell structure of ocean worlds with gravity-topography admittance

The structure of the icy shells of ocean worlds is important for understanding the stability of their underlying oceans as it controls the rate at which heat can be transported outward and radiated to space. Future spacecraft exploration of the ocean worlds (e.g., by NASA's Europa Clipper mission) will allow for higher-resolution measurements of gravity and shape than currently available. In this paper, we study the sensitivity of gravity-topography admittance to the structure of icy shells in preparation for future data analysis. An analytical viscous relaxation model is used to predict admittance spectra given different shell structures determined by the temperature-dependent viscosity of a tidally heated, conductive shell. We apply these methods to the ocean worlds of Europa and Enceladus. We find that admittance is sensitive to the mechanisms of topography support at different wavelengths and estimate the required gravity performance to resolve transitions between these mechanisms. We find that the Airy isostatic model is unable to accurately describe admittance universally across all wavelengths when the shell thickness is a significant fraction of body's radius. Our models suggest that measurements of admittance at low spherical harmonic degrees are more sensitive to thick shells with high tidal dissipation, and may complement ice-penetrating radar measurements in constraining shell thickness. Finally, we find that admittance may be used to constrain the tidal dissipation within the icy shell, which would be complementary to a more demanding measurement of the tidal phase lag.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Equation of State of Hot, Dense Magnesium Derived with First-PrinciplesComputer Simulations

Using two first-principles computer simulation techniques, path integral Monte-Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium in the regime of warm dense matter, with densities ranging from 0.43 to 86.11~g/cm$^3$~and temperatures from 20,000 K to $5\times10^8$~K. These conditions are relevant for the interiors of giant planets and stars as well as for shock compression measurements and inertial confinement fusion experiments. We study ionization mechanisms and electronic structure of magnesium as a function of density and temperature. We show that the L shell electrons 2s and 2p energy bands merge at high density. This results into a gradual ionization of the L-shell with increasing density and temperature. In this regard, Mg differs from MgO, which is also reflected in the shape of its principal shock Hugoniot curve. For Mg, we predict a single broad pressure-temperature region where the shock compression ratio is approximately 4.9. Mg thus differs from Si and Al plasma that exhibit two well-separated compression maxima on the Hugoniot curve for L and K shell ionizations. Finally we study multiple shocks and effects of preheat and precompression.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Equilibrium Tidal Response of Jupiter: Detectability by Juno

An observation of Jupiter's tidal response is anticipated for the on-going Juno spacecraft mission. We combine self-consistent, numerical models of Jupiter's equilibrium tidal response with observed Doppler shifts from the Juno gravity science experiment to test the sensitivity of the spacecraft to tides raised by the Galilean satellites and the Sun. The concentric Maclaurin spheroid (CMS) method finds the equilibrium shape and gravity field of a rotating, liquid planet with the tide raised by a satellite, expanded in Love numbers ($k_{nm}$). We present improvements to CMS theory that eliminate an unphysical center of mass offset and study in detail the convergence behavior of the CMS approach. We demonstrate that the dependence of $k_{nm}$ with orbital distance is important when considering the combined tidal response for Jupiter. Conversely, the details of the interior structure have a negligible influence on $k_{nm}$, for models that match the zonal harmonics $J_2$, $J_4$ and $J_6$, already measured to high precision by Juno. As the mission continues, improved coverage of Jupiter's gravity field at different phases of Io's orbit is expected to yield an observed value for the degree-2 Love number ($k_{22}$) and potentially select higher--degree $k_{nm}$. We present a test of the sensitivity of the Juno Doppler signal to the calculated $k_{nm}$, which suggests the detectability of $k_{33}$, $k_{42}$ and $k_{31}$, in addition to $k_{22}$. A mismatch of a robust Juno observation with the remarkably small range in calculated Io equilibrium $k_{22}=0.58976\pm0.0001$, would indicate a heretofore uncharacterized dynamic contribution to the tides.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

First-Principles Equation of State Database for Warm Dense Matter Computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al and Si as well as the compounds LiF, B4C, BN, CH4, CH2, C2H3, CH, C2H, MgO, and MgSiO3. For all these materials, we provide the pressure and internal energy over a density-temperature range from ~0.5 to 50 g/cc and from ~10^4 to 10^9 K, which are based on ~5000 different first-principles simulations. We compute isobars, adiabats and shock Hugoniot curves in the regime of L and K shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. We derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H2O, H2O2, Al2O3, CO, and CO2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their endmembers. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the supplementary material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnesium Oxide at Extreme Temperatures and Pressures Studied with First-Principles Simulations

We combine two first-principles computer simulation techniques, path integral Monte-Carlo and density functional theory molecular dynamics, to determine the equation of state of magnesium oxide in the regime of warm dense matter, with densities ranging from 0.35 to 71~g$\,$cm$^{-3}$ and temperatures from 10,000 K to $5\times10^8$~K. These conditions are relevant for the interiors of giant planets and stars as well as for shock wave compression measurements and inertial confinement fusion experiments. We study the electronic structure of MgO and the ionization mechanisms as a function of density and temperature. We show that the L-shell orbitals of magnesium and oxygen hybridize at high density. This results into a gradual ionization of the L-shell with increasing density and temperature. In this regard, MgO behaves differently from pure oxygen, which is reflected in the shape of the MgO principal shock Hugoniot curve. The curve of oxygen shows two compression maxima, while that of MgO shows only one. We predict a maximum compression ratio of 4.66 to occur for a temperature of 6.73 $\times 10^7$ K. Finally we study how multiple shocks and ramp waves can be used to cover a large range of densities and temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Warm, Dense MgSiO$_3$

In order to provide a comprehensive theoretical description of MgSiO$_3$ at extreme conditions, we combine results from path integral Monte Carlo (PIMC) and density functional molecular dynamics simulations (DFT-MD) and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10$^4$ to 10$^8$ K and from 0.321 to 64.2 g$\,$cm$^{-3}$ (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. Finally we analyze the heat capacity and structural properties of the liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Prediction of Chlorine and Fluorine Crystal Structures at High Pressure Using Symmetry Driven Structure Search with Geometric Constraints

The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which has made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Our algorithm extends the symmetry driven structure search method [Phys. Rev. B 98 (2018) 174107] by adding constraints for the bond length and the number of atoms in a molecule, while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range from 10--4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm$\bar{3}$m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to an Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P$4_2$/mmc symmetry at 2500 GPa, and finally into its cubic analogue form with Pm$\bar{3}$n symmetry at 3000 GPa.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-temperature miscibility of iron and rock during terrestrial planet formation

The accretion of a terrestrial body and differentiation of its silicate/oxide mantle from iron core provide abundant energy for heating its interior to temperatures much higher than the present day Earth. The consequences of differentiation on the structure and composition of planets are typically addressed considering only the interaction of molten iron with an immiscible `rocky' phase. We demonstrate that mixing in a representative system of liquid or solid MgO and liquid iron to a single homogeneous liquid occurs at sufficiently low temperature to be present in the aftermath of a giant impact. Applying the thermodynamic integration technique to density functional theory molecular dynamics simulations, we determine the solvus closure temperature for the Fe-MgO system for pressures up to 400 GPa. Solvus closure occurs at $\sim$4000 K at low pressure, and has a weak positive pressure dependence, such that its gradient with respect to depth is less steep than an adiabatic temperature profile. This predicts a new mode of core-mantle differentiation following the most energetic giant impacts, with exsolution of iron from the mixture beginning in the outer layers of the planet. We demonstrate that high-temperature equilibration results in delivery of nominally insoluble Mg-rich material to the early core. Since MgO is the least soluble major mantle component in iron at low temperatures, these results may represent an upper bound on temperature for mixing in terrestrial planets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Concentric Maclaurin Spheroid method with tides and a rotational enhancement of Saturn's tidal response

We extend to three dimensions the Concentric Maclaurin Spheroid method for obtaining the self-consistent shape and gravitational field of a rotating liquid planet, to include a tidal potential from a satellite. We exhibit, for the first time, the important effect of the planetary rotation rate on tidal response of gas giants. Simulations of planets with fast rotation rates like those of Jupiter and Saturn, exhibit significant changes in calculated tidal love numbers $k_{nm}$ when compared with non-rotating bodies. A test model of Saturn fitted to observed zonal gravitational multipole harmonics yields $k_2=0.413$, consistent with a recent observational determination from Cassini astrometry data (Lainey et al., 2016). The calculated love number is robust under reasonable assumptions of interior rotation rate, satellite parameters, and details of Saturn's interior structure. The method is benchmarked against several published test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The properties of heavy elements in giant planet envelopes

The core accretion model for giant planet formation suggests a two layer picture for the initial structure of Jovian planets, with heavy elements in a dense core and a thick H-He envelope. Late planetesimal accretion and core erosion could potentially enrich the H-He envelope in heavy elements, which is supported by the three-fold solar metallicity that was measured in Jupiter's atmosphere by the Galileo entry probe. In order to reproduce the observed gravitational moments of Jupiter and Saturn, models for their interiors include heavy elements, $Z$, in various proportions. However, their effect on the equation of state of the hydrogen-helium mixtures has not been investigated beyond the ideal mixing approximation. In this article, we report results from \textit{ab initio} simulations of fully interacting H-He-$Z$ mixtures in order to characterize their equation of state and to analyze possible consequences for the interior structure and evolution of giant planets. Considering C, N, O, Si, Fe, MgO and SiO$_2$, we show that the behavior of heavy elements in H-He mixtures may still be represented by an ideal mixture if the effective volumes and internal energies are chosen appropriately. In the case of oxygen, we also compute the effect on the entropy. We find the resulting changes in the temperature-pressure profile to be small. A homogeneous distribution of 2\% oxygen by mass changes the temperature in Jupiter's interior by only 80 K.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tidal Response of Preliminary Jupiter Model

In anticipation of improved observational data for Jupiter's gravitational field from the Juno spacecraft, we predict the static tidal response for a variety of Jupiter interior models based on ab initio computer simulations of hydrogen-helium mixtures. We calculate hydrostatic-equilibrium gravity terms using the non-perturbative concentric Maclaurin Spheroid (CMS) method that eliminates lengthy expansions used in the theory of figures. Our method captures terms arising from the coupled tidal and rotational perturbations, which we find to be important for a rapidly-rotating planet like Jupiter. Our predicted static tidal Love number $k_2 = 0.5900$ is $\sim$10\% larger than previous estimates. The value is, as expected, highly correlated with the zonal harmonic coefficient $J_2$, and is thus nearly constant when plausible changes are made to interior structure while holding $J_2$ fixed at the observed value. We note that the predicted static $k_2$ might change due to Jupiter's dynamical response to the Galilean moons, and find reasons to argue that the change may be detectable, although we do not present here a theory of dynamical tides for highly oblate Jovian planets. An accurate model of Jupiter's tidal response will be essential for interpreting Juno observations and identifying tidal signals from effects of other interior dynamics in Jupiter's gravitational field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Understanding Jupiter's Interior

This article provides an overview of how models of giant planet interiors are constructed. We review measurements from past space missions that provide constraints for the interior structure of Jupiter. We discuss typical three-layer interior models that consist of a dense central core and an inner metallic and an outer molecular hydrogen-helium layer. These models rely heavily on experiments, analytical theory, and first-principle computer simulations of hydrogen and helium to understand their behavior up to the extreme pressures ~10 Mbar and temperatures ~10,000 K. We review the various equations of state used in Jupiter models and compare them with shock wave experiments. We discuss the possibility of helium rain, core erosion and double diffusive convection may have important consequences for the structure and evolution of giant planets. In July 2016 the Juno spacecraft entered orbit around Jupiter, promising high-precision measurements of the gravitational field that will allow us to test our understanding of gas giant interiors better than ever before.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Development of Path Integral Monte Carlo Simulations with Localized Nodal Surfaces for Second-Row Elements

We extend the applicability range of fermionic path integral Monte Carlo simulations to heavier elements and lower temperatures by introducing various localized nodal surfaces. Hartree-Fock nodes yield the most accurate prediction for pressure and internal energy that we combine with the results from density functional molecular dynamics simulations to obtain a consistent equation of state for hot, dense silicon under plasma conditions and in the regime of warm dense matter (2.3-18.6 g/cm, 5.0*10^5 - 1.3*10^8 K). The shock Hugoniot curve is derived and the structure of the fluid is characterized with various pair correlation functions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-pressure, temperature elasticity of Fe- and Al-bearing MgSiO3: implications for the Earth's lower mantle

Fe and Al are two of the most important rock-forming elements other than Mg, Si, and O. Their presence in the lower mantle's most abundant minerals, MgSiO_3 bridgmanite, MgSiO_3 post-perovskite and MgO periclase, alters their elastic properties. However, knowledge on the thermoelasticity of Fe- and Al-bearing MgSiO_3 bridgmanite, and post-perovskite is scarce. In this study, we perform ab initio molecular dynamics to calculate the elastic and seismic properties of pure, Fe^{3+}- and Fe^{2+}-, and Al^{3+}-bearing MgSiO_3 perovskite and post-perovskite, over a wide range of pressures, temperatures, and Fe/Al compositions. Our results show that a mineral assemblage resembling pyrolite fits a 1D seismological model well, down to, at least, a few hundred kilometers above the core-mantle boundary, i.e. the top of the D'' region. In D'', a similar composition is still an excellent fit to the average velocities and fairly approximate to the density. We also implement polycrystal plasticity with a geodynamic model to predict resulting seismic anisotropy, and find post-perovskite with predominant (001) slip across all compositions agrees best with seismic observations in the D''.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrogen-Water Mixtures in Giant Planet Interiors Studied with Ab Initio Simulations

We study water-hydrogen mixtures under planetary interior conditions using ab initio molecular dynamics simulations. We determine the thermodynamic properties of various water-hydrogen mixing ratios at temperatures of 2000 and 6000 K for pressures of a few tens of GPa. These conditions are relevant for ice giant planets and for the outer envelope of the gas giants. We find that at 2000 K the mixture is in a molecular regime, while at 6000 K the dissociation of hydrogen and water is important and affects the thermodynamic properties. We study the structure of the liquid and analyze the radial distribution function. We provide estimates for the transport properties, diffusion and viscosity, based on autocorrelation functions. We obtained viscosity estimates of the order of a few tenths of mPa.s for the conditions under consideration. These results are relevant for dynamo simulations of ice giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Miscibility calculations for water and hydrogen in giant planets

We present results from ab initio simulations of liquid water-hydrogen mixtures in the range from 2 to 70 GPa and from 1000 to 6000 K, covering conditions in the interiors of ice giant planets and parts of the outer envelope of gas giant planets. In addition to computing the pressure and the internal energy, we derive the Gibbs free energy by performing a thermodynamic integration. For all conditions under consideration, our simulations predict hydrogen and water to mix in all proportions. The thermodynamic behavior of the mixture can be well described with an ideal mixing approximation. We suggest a substantial fraction of water and hydrogen in giant planets may occur in homogeneously mixed form rather than in separate layers. The extend of mixing depends on the planet's interior dynamics and its conditions of formation, in particular on how much hydrogen was present when icy planetesimals were delivered. Based on our results, we do not predict water-hydrogen mixtures to phase separate during any stage of the evolution of giant planets. We also show that the hydrogen content of an exoplanet is much higher if the mixed interior is assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equations of state and stability of MgSiO$_3$ perovskite and post-perovskite phases from quantum Monte Carlo simulations

We have performed quantum Monte Carlo (QMC) simulations and density functional theory (DFT) calculations to study the equations of state of MgSiO$_3$ perovskite (Pv) and post-perovskite (PPv), up to the pressure and temperature conditions of the base of Earth's lower mantle. The ground state energies were derived using QMC and the temperature dependent Helmholtz free energies were calculated within the quasi-harmonic approximation and density functional perturbation theory. The equations of state for both phases of MgSiO$_3$ agree well with experiments, and better than those from generalized gradient approximation (GGA) calculations. The Pv-PPv phase boundary calculated from our QMC equations of states is also consistent with experiments, and better than previous LDA calculations. We discuss the implications for double crossing of the Pv-PPv boundary in the Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A multiphase equation of state for carbon addressing high pressures and temperatures

We present a 5-phase equation of state (EOS) for elemental carbon. The phases considered are: diamond, BC8, simple-cubic, simple-hexagonal, and the liquid/plasma state. The solid phase free energies are constrained by density functional theory (DFT) calculations. Vibrational contributions to the free energy of each solid phase are treated within the quasiharmonic framework. The liquid free energy model is constrained by fitting to a combination of DFT molecular dynamics performed over the range 10,000 K < T < 100,000 K, and path integral quantum Monte Carlo calculations for T > 100,000 K (both for ρ between 3 and 12 g/cc, with select higher-ρ DFT calculations as well). The liquid free energy model includes an atom-in-jellium approach to account for the effects of ionization due to temperature and pressure in the plasma state, and an ion-thermal model which includes the approach to the ideal gas limit. The precise manner in which the ideal gas limit is reached is greatly constrained by both the highest temperature DFT data and the path integral data, forcing us to discard an ion-thermal model we had used previously in favor of a new one. Predictions are made for the principal Hugoniot and the room-temperature isotherm, and comparisons are made to recent experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solubility of Iron in Metallic Hydrogen and Stability of Dense Cores in Giant Planets

The formation of the giant planets in our solar system, and likely a majority of giant exoplanets, is commonly explained by the accretion of nebular hydrogen and helium onto a large core of terrestrial-like composition. The fate of this core has important consequences for the evolution of the interior structure of the planet. It has recently been shown that H2O, MgO and SiO2 dissolve in liquid metallic hydrogen at high temperature and pressure. In this study, we perform ab initio calculations to study the solubility of an innermost metallic core. We find dissolution of iron to be strongly favored above 2000 K over the entire pressure range (0.4-4 TPa) considered. We compare with and summarize the results for solubilities on other probable core constituents. The calculations imply that giant planet cores are in thermodynamic disequilibrium with surrounding layers, promoting erosion and redistribution of heavy elements. Differences in solubility behavior between iron and rock may influence evolution of interiors, particularly for Saturn-mass planets. Understanding the distribution of iron and other heavy elements in gas giants may be relevant in understanding mass-radius relationships, as well as deviations in transport properties from pure hydrogen-helium mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Superionic to superionic phase change in water: consequences for the interiors of Uranus and Neptune

Using density functional molecular dynamics free energy calculations, we show that the body-centered-cubic phase of superionic ice previously believed to be the only phase is in fact thermodynamically unstable compared to a novel phase with oxygen positions in fcc lattice sites. The novel phase has a lower proton mobility than the bc phase and may exhibit a higher melting temperature. We predict a transition between the two phases at a pressure of 1 +/- 0.5 Mbar, with potential consequences for the interiors of ice giants such as Uranus and Neptune.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab Initio Investigation of a Possible Liquid-Liquid Phase Transition in MgSiO3 at Megabar Pressures

We perform density functional molecular dynamics simulations of liquid and solid MgSiO3 in the pressure range of 120-1600 GPa and for temperatures up to 20000 K in order to provide new insight into the nature of the liquid-liquid phase transition that was recently predicted on the basis of decaying laser shock wave experiments [Phys. Rev. Lett. 108 (2012) 065701]. However, our simulations did not show any signature of a phase transition in the liquid phase. We derive the equation of state for the liquid and solid phases and compute the shock Hugoniot curves. We discuss different thermodynamic functions and by explore alternative interpretations of the experimental findings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab Initio Simulations of Hot, Dense Methane During Shock Experiments

Using density functional theory molecular dynamics simulations, we predict shock Hugoniot curves of precompressed methane up to 75000 K for initial densities ranging from 0.35 to 0.70 g/cc. At 4000 K, we observe the transformation into a metallic, polymeric state consisting of long hydrocarbon chains. These chains persist when the sample is quenched to 300 K, leading to an increase in shock compression. At 6000 K, the sample transforms into a plasma composed of many, short-lived chemical species. We conclude by discussing implications for the interiors of Uranus and Neptune and analyzing the possibility of creating a superionic state of methane in high pressure experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 (PIMC) method with free-particle nodes for warm dense matter and apply it to water and carbon plasmas. We thereby extend PIMC studies beyond hydrogen and helium to elements with core electrons. PIMC pressures, internal energies, and pair-correlation functions compare well with density functional theory molecular dynamics (DFT-MD) at temperatures of (2.5-7.5)$\times10^5$ K and both methods together form a coherent equation of state (EOS) over a density-temperature range of 3--12 g/cm$^3$ and 10$^4$--10$^9$ K.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

H4O and other hydrogen-oxygen compounds at giant-planet core pressures

Water and hydrogen at high pressure make up a substantial fraction of the interiors of giant planets. Using ab initio random structure search methods we investigate the ground-state crystal structures of water, hydrogen, and hydrogen-oxygen compounds. We find that, at pressures beyond 14 Mbar, excess hydrogen is incorporated into the ice phase to form a novel structure with H4O stoichiometry. We also predict two new ground state structures, P2_1/m and I4/mmm, for post-C2/m water ice.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Erosion of icy cores in giant gas planets

Using ab initio simulations we investigate whether water ice is stable in the cores of giant planets, or whether it dissolves into the layer of metallic hydrogen above. By Gibbs free energy calculations we find that for pressures between 10 and 40 Mbar the ice-hydrogen interface is unstable at temperatures above approximately 3000 K, far below the temperature of the core-mantle boundaries in Jupiter and Saturn that are of the order of 10000 K. This implies that the cores of solar and extrasolar giant planets are at least partially eroded.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rocky core solubility in Jupiter and giant exoplanets

Gas giants are believed to form by the accretion of hydrogen-helium gas around an initial protocore of rock and ice. The question of whether the rocky parts of the core dissolve into the fluid H-He layers following formation has significant implications for planetary structure and evolution. Here we use ab initio calculations to study rock solubility in fluid hydrogen, choosing MgO as a representative example of planetary rocky materials, and find MgO to be highly soluble in H for temperatures in excess of approximately 10000 K, implying significant redistribution of rocky core material in Jupiter and larger exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

New Phases of Water Ice Predicted at Megabar Pressures

Based on density functional calculations we predict water ice to attain two new crystal structures with Pbca and Cmcm symmetry at 7.6 and 15.5 Mbar, respectively. The known high pressure ice phases VII, VIII, X, and Pbcm as well as the Pbca phase are all insulating and composed of two interpenetrating hydrogen bonded networks, but the Cmcm structure is metallic and consists of corrugated sheets of H and O atoms. The H atoms are squeezed into octahedral positions between next-nearest O atoms while they occupy tetrahedral positions between nearest O atoms in the ice X, Pbcm, and Pbca phases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sequestration of noble gases in giant planet interiors

The Galileo probe showed that Jupiter's atmosphere is severely depleted in neon compared to protosolar values. We show, via ab initio simulations of the partitioning of neon between hydrogen and helium phases, that the observed depletion can be explained by the sequestration of neon into helium-rich droplets within the postulated hydrogen-helium immiscibility layer of the planet's interior. We also demonstrate that this mechanism will not affect argon, explaining the observed lack of depletion of this gas. This provides strong indirect evidence for hydrogen-helium immiscibility in Jupiter.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Correlations in Hot Dense Helium

Hot dense helium is studied with first-principles computer simulations. By combining path integral Monte Carlo and density functional molecular dynamics, a large temperature and density interval ranging from 1000 to 1000000 K and 0.4 to 5.4 g/cc becomes accessible to first-principles simulations and the changes in the structure of dense hot fluids can be investigated. The focus of this article are pair correlation functions between nuclei, between electrons, and between electrons and nuclei. The density and temperature dependence of these correlation functions is analyzed in order to describe the structure of the dense fluid helium at extreme conditions.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Characterization of the State of Hydrogen

Fermionic path integral Monte Carlo simulations have been applied to study the equilibrium properties of the hydrogen and deuterium in the density and temperature range of 1.6 < rs < 14.0 and 5000K < T < 167000K. We use this technique to determine the phase diagram by identifying the plasma, the molecular, atomic and metallic regime. We explain how one can identify the phases in the path integral formalism and discuss the state of hydrogen for 5 points in the temperature-density plane. Further we will provide arguments for the nature of the transitions between the regimes.

0 citations0 reviewsNo code linkedOpen access