Source author record

K. P. Driver

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

5 published item(s)

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.

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.

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.

preprint2010arXiv

Phase transformation in Si from semiconducting diamond to metallic beta-Sn phase in QMC and DFT under hydrostatic and anisotropic stress

Silicon undergoes a phase transition from the semiconducting diamond phase to the metallic beta-Sn phase under pressure. We use quantum Monte Carlo calculations to predict the transformation pressure and compare the results to density functional calculations employing the LDA, PBE, PW91, WC, AM05, PBEsol and HSE06 exchange-correlation functionals. Diffusion Monte Carlo predicts a transition pressure of 14.0 +- 1.0 GPa slightly above the experimentally observed transition pressure range of 11.3 to 12.6 GPa. The HSE06 hybrid functional predicts a transition pressure of 12.4 GPa in excellent agreement with experiments. Exchange-correlation functionals using the local-density approximation and generalized-gradient approximations result in transition pressures ranging from 3.5 to 10.0 GPa, well below the experimental values. The transition pressure is sensitive to stress anisotropy. Anisotropy in the stress along any of the cubic axes of the diamond phase of silicon lowers the equilibrium transition pressure and may explain the discrepancy between the various experimental values as well as the small overestimate of the quantum Monte Carlo transition pressure.

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.