BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Burkhard Militzer

Burkhard Militzer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciastro-ph.EPphysics.plasm-phcond-mat.stat-mechphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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