Source author record

Mandy Bethkenhagen

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

ResearcherUnclaimed source record

physics.comp-ph physics.chem-ph physics.plasm-ph astro-ph.SR cond-mat.mtrl-sci cond-mat.stat-mech

Catalog footprint

What is connected

4works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Same and interconvertible high-pressure ice phases

Most experimentally known high-pressure ice phases have a body-centred cubic (bcc) oxygen lattice. Our atomistic simulations show that, amongst these bcc ice phases, ices VII, VII' and X are the same thermodynamic phase under different conditions, whereas superionic ice VII'' has a first-order phase boundary with ice VII'. Moreover, at about 300 GPa, ice X transforms into the Pbcm phase with a sharp structural change but no apparent activation barrier, whilst at higher pressures the barrier gradually increases. Our study thus clarifies the phase behaviour of the high-pressure insulating ices and reveals peculiar solid-solid transition mechanisms not known in other systems.

preprint2020arXiv

Carbon ionization at Gbar pressures: an ab initio perspective on astrophysical high-density plasmas

A realistic description of partially-ionized matter in extreme thermodynamic states is critical to model the interior and evolution of the multiplicity of high-density astrophysical objects. Current predictions of its essential property, the ionization degree, rely widely on analytical approximations that have been challenged recently by a series of experiments. Here, we propose a novel ab initio approach to calculate the ionization degree directly from the dynamic electrical conductivity using the Thomas-Reiche-Kuhn sum rule. This Density Functional Theory framework captures genuinely the condensed matter nature and quantum effects typical for strongly-correlated plasmas. We demonstrate this new capability for carbon and hydrocarbon, which most notably serve as ablator materials in inertial confinement fusion experiments aiming at recreating stellar conditions. We find a significantly higher carbon ionization degree than predicted by commonly used models, yet validating the qualitative behavior of the average atom model Purgatorio. Additionally, we find the carbon ionization state to remain unchanged in the environment of fully-ionized hydrogen. Our results will not only serve as benchmark for traditional models, but more importantly provide an experimentally accessible quantity in the form of the electrical conductivity.

preprint2020arXiv

Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann-Feynman stress tensor for $\mathcal{O}(N)$ Kohn-Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi-Dirac distribution becomes smoother and density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to $\mathcal{O}(N)$ methods. We then describe its evaluation within the $\mathcal{O}(N)$ infinite-cell Clenshaw-Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic as well as insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results, and show convergence with respect to mesh size to established planewave results. We employ the new formulation to compute the viscosity of hydrogen at a million kelvin from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

preprint2019arXiv

High temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from 6 to 3 quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. This parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.