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

B. Militzer

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciastro-ph.EPphysics.plasm-phastro-ph.SR

Catalog footprint

What is connected

12works
4topics
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

12 published item(s)

preprint2022arXiv

Jupiter's inhomogeneous envelope

While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-elements distribution in Jupiter's envelope. In this paper, we assemble the most comprehensive and diverse collection of Jupiter interior models to date and use it to study the distribution of heavy elements in the planet's envelope. We apply a Bayesian statistical approach to our interior model calculations, reproducing the Juno gravitational and atmospheric measurements and constraints from the deep zonal flows. Our results show that the gravity constraints lead to a deep entropy of Jupiter corresponding to a 1 bar temperature 5-15 K higher than traditionally assumed. We also find that uncertainties in the equation of state are crucial when determining the amount of heavy elements in Jupiter's interior. Our models put an upper limit to the inner compact core of Jupiter of 7 Earth masses, independently on the structure model (with or without dilute core) and the equation of state considered. Furthermore, we robustly demonstrate that Jupiter's envelope is inhomogenous, with a heavy-element enrichment in the interior relative to the outer envelope. This implies that heavy element enrichment continued through the gas accretion phase, with important implications for the formation of giant planets in our solar system and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Preliminary Jupiter Model

In anticipation of new observational results for Jupiter's axial moment of inertia and gravitational zonal harmonic coefficients from the forthcoming Juno orbiter, we present a number of preliminary Jupiter interior models. We combine results from ab initio computer simulations of hydrogen-helium mixtures, including immiscibility calculations, with a new nonperturbative calculation of Jupiter's zonal harmonic coefficients, to derive a self-consistent model for the planet's external gravity and moment of inertia. We assume helium rain modified the interior temperature and composition profiles. Our calculation predicts zonal harmonic values to which measurements can be compared. Although some models fit the observed (pre-Juno) second- and fourth-order zonal harmonics to within their error bars, our preferred reference model predicts a fourth-order zonal harmonic whose absolute value lies above the pre-Juno error bars. This model has a dense core of about 12 Earth masses, and a hydrogen-helium-rich envelope with approximately 3 times solar metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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$\times$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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First-Principles Equation of State Calculations of Warm Dense Nitrogen

Using path integral Monte Carlo (PIMC) and density functional molecular dynamics (DFT-MD) simulation methods, we compute a coherent equation of state (EOS) of nitrogen that spans the liquid, warm dense matter (WDM), and plasma regimes. Simulations cover a wide range of density-temperature space, $1.5-13.9$~g$~$cm$^{-3}$ and $10^3-10^9$~K. In the molecular dissociation regime, we extend the pressure-temperature phase diagram beyond previous studies, providing dissociation and Hugoniot curves in good agreement with experiments and previous DFT-MD work. Analysis of pair-correlation functions and the electronic density of states in the WDM regime reveals an evolving plasma structure and ionization process that is driven by temperature and pressure. Our Hugoniot curves display a sharp change in slope in the dissociation regime and feature two compression maxima as the K and L shells are ionized in the WDM regime, which have some significant differences from the predictions of plasma models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First-Principles Prediction of the Softening of the Silicon Shock Hugoniot Curve

Shock compression of silicon (Si) under extremely high pressures (>100 Mbar) was investigated by using two first-principles methods of orbital-free molecular dynamics (OFMD) and path integral Monte Carlo (PIMC). While pressures from the two methods agree very well, PIMC predicts a second compression maximum because of 1s electron ionization that is absent in OFMD calculations since Thomas-Fermi-based theories lack shell structure. The Kohn-Sham density functional theory is used to calculate the equation of state (EOS) of warm dense silicon for low-pressure loadings (P < 100 Mbar). Combining these first-principles EOS results, the principal shock Hugoniot curve of silicon for pressures varying from 1 Mbar to above 10 Gbar was derived. We find that silicon is 20% or more softer than what was predicted by widely-used EOS models. Existing high-pressure experimental data (P = 1 - 2 Mbar) seem to indicate this softening behavior of Si, which calls for future strong-shock experiments (P > 10 Mbar) to benchmark our results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ab Initio Equation of State for Hydrogen-Helium Mixtures with Recalibration of the Giant-Planet Mass-Radius Relation

Using density functional molecular dynamics simulations, we determine the equation of state for hydrogen-helium mixtures spanning density-temperature conditions typical of giant planet interiors, ~0.2-9 g/cc and 1000-80000 K for a typical helium mass fraction of 0.245. In addition to computing internal energy and pressure, we determine the entropy using an ab initio thermodynamic integration technique. A comprehensive equation of state (EOS) table with 391 density-temperature points is constructed and the results are presented in form of two-dimensional free energy fit for interpolation. Deviations between our ab initio EOS and the semi-analytical EOS model by Saumon and Chabrier are analyzed in detail, and we use the results for initial revision of the inferred thermal state of giant planets with known values for mass and radius. Changes are most pronounced for planets in the Jupiter mass range and below. We present a revision to the mass-radius relationship which makes the hottest exoplanets increase in radius by ~0.2 Jupiter radii at fixed entropy and for masses greater than ~0.5 Jupiter mass. This change is large enough to have possible implications for some discrepant "inflated giant exoplanets".

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bonding and Electronic Properties of Ice at High Pressure

The properties of water ice at megabar pressure are characterized with ab initio computer simulations. The focus lies on the metallic Cmcm phase and its insulating distorted analogue with Pnma symmetry. Both phases were recently predicted to occur at 15.5 and 12.5 Mbar respectively [Phys. Rev. Lett. 105 (2010) 195701]. The Fermi surface of the Cmcm phase is analyzed and possibility of Fermi nesting to occur is discussed. The Wannier orbital are computed for the Pnma structure and compared to ice X. While ice X shows typical sp3 hybridization, in the Pnma structure, the orbitals are deformed and no longer all aligned with the hydrogen bonds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

FPEOS: A First-Principles Equation of State Table of Deuterium for Inertial Confinement Fusion Applications

Understanding and designing inertial confinement fusion (ICF) implosions through radiation-hydrodynamics simulations rely on the accurate knowledge of the equation of state (EOS) of the deuterium and tritium fuels. To minimize the drive energy for ignition, the imploding shell of DT fuel must be kept as cold as possible. Such low-adiabat ICF implosions can access to coupled and degenerate plasma conditions, in which the analytical or chemical EOS models become inaccurate. Using the path-integral Monte Carlo (PIMC) simulations we have derived a first-principles EOS (FPEOS) table of deuterium that covers typical ICF fuel conditions at densities ranging from 0.002 to 1596 g/cm3 and temperatures of 1.35 eV to 5.5 keV. We report the internal energy and the pressure, and discuss the structure of the plasma in terms of pair-correlation functions. When compared with the widely used SESAME table and the revised Kerley03 table, discrepancies in the internal energy and in the pressure are identified for moderately coupled and degenerate plasma conditions. In contrast to the SESAME table, the revised Kerley03 table is in better agreement with our FPEOS results over a wide range of densities and temperatures. Although subtle differences still exist for lower temperatures (T < 10 eV) and moderate densities (1 to 10 g/cc), hydrodynamics simulations of cryogenic ICF implosions using the FPEOS table and the Kerley03 table have resulted in similar results for the peak density, areal density (ρR), and neutron yield, which are significantly different from the SESAME simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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 the equation of state of cubic boron nitride(c-BN). We compute the static contribution to the free energy with QMC, and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. Furthermore, we compute the anharmonic Raman frequency shift with QMC as a function of pressure and temperature, allowing optical pressure calibration in table-top experiments. 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 methodology is generally applicable to first-row solids, and can be used to provide a new reference for ab initio calculations of solids and to benchmark pseudopotential accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo for minerals at high pressure: Phase stability, equations of state, and elasticity of silica

Silica is an abundant component of the Earth whose crystalline polymorphs play key roles in its structure and dynamics. As the simplest silicates, understanding pure silica is a prerequisite to understanding the rocky part of the Earth, its majority. First principle density functional theory (DFT) methods have often been used to accurately predict properties of silicates. 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 Schrodinger's equation. Using ground state QMC plus phonons within the quasiharmonic approximation from 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 most well-constrained equations of state and phase boundaries available for silica. QMC indicates a transition to the most dense alpha-PbO2 structure above the core-insulating D" layer, suggesting the absence of significant free silica in the bulk lower mantle, which has been assumed but never proven. We also find an accurate shear elastic constant and its geophysically important softening with pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Frontiers of the physics of dense plasmas and planetary interiors: experiments, theory, applications

Recent developments of dynamic x-ray characterization experiments of dense matter are reviewed, with particular emphasis on conditions relevant to interiors of terrestrial and gas giant planets. These studies include characterization of compressed states of matter in light elements by x-ray scattering and imaging of shocked iron by radiography. Several applications of this work are examined. These include the structure of massive "Super Earth" terrestrial planets around other stars, the 40 known extrasolar gas giants with measured masses and radii, and Jupiter itself, which serves as the benchmark for giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Hot, Dense Helium

Two first-principles simulation techniques, path integral Monte Carlo (PIMC) and density functional molecular dynamics (DFT-MD), are applied to study hot, dense helium in the density-temperature range of 0.387 - 5.35 g/cc and 500 K - 1.28x10^8 K. One coherent equation of state (EOS) is derived by combining DFT-MD data at lower temperatures with PIMC results at higher temperatures. Good agreement between both techniques is found in an intermediate temperature range. For the highest temperatures, the PIMC results converge to the Debye-Hueckel limiting law. In order derive the entropy, a thermodynamically consistent free energy fit is introduced that reproduces the internal energies and pressure derived from the first-principles simulations. The equation of state is presented in form of a table as well as a fit and is compared with chemical models. In addition, the structure of the fluid is analyzed using pair correlation functions. Shock Hugoniot curves are compared with recent laser shock wave experiments.

0 citations0 reviewsNo code linkedOpen access