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

David M. Ceperley

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

ResearcherUnclaimed source record
cond-mat.str-elphysics.comp-phcond-mat.mtrl-sciphysics.chem-phcond-mat.otherquant-phcond-mat.quant-gascond-mat.supr-concond-mat.stat-mechmath.NAphysics.atom-ph

Catalog footprint

What is connected

28works
11topics
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

28 published item(s)

preprint2022arXiv

Stability and distortion of fcc-LaH$_{10}$ with path-integral molecular dynamics

The synthesis of the high temperature superconductor LaH$_{10}$ requires pressures in excess of 100 GPa, wherein it adopts a face-centered cubic structure. Upon decompression, this structure undergoes a distortion which still supports superconductivity, but with a much lower critical temperature. Previous calculations have shown that quantum and anharmonic effects are necessary to stabilize the cubic structure, but have not resolved the low pressure distortion. Using large scale path-integral molecular dynamics enabled by a machine learned potential, we show that a rhombohedral distortion appears at sufficiently low pressures, even with quantum and anharmonic effects. We also highlight the importance of quantum zero point motion in stabilizing the cubic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born-Oppenheimer approximation

We develop a formalism to accurately account for the renormalization of electronic structure due to quantum and thermal nuclear motions within the Born-Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from Quantum Monte Carlo calculations in either the canonical or grand canonical ensembles. The formalism applies as well to effective single electron theories such as those based on Density Functional Theory. We show that electronic (Bloch) crystal momentum can be restored by marginalizing the total electron-ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born-Oppenheimer approximation. Based on the Kubo-Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and in general differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200K and 250 GPa and discuss the optical absorption profile of hydrogen crystal at 200K and carbon diamond at 297K.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy gap closure of crystalline molecular hydrogen with pressure

We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical Quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by Coupled Electron Ion Monte Carlo. Nuclear zero point effects cause a large reduction in the gap ($\sim 2eV$). \CP{Depending on the structure,} the fundamental indirect gap closes \CP{between 380GPa and} 530GPa for ideal crystals and 330-380GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to $\sim$450\CP{-500} GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Electronic band gaps from Quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical Quantum Monte Carlo simulations. We discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from DFT calculations. Our procedure is applied to solid molecular hydrogen and compared to experiment for carbon and silicon crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quantum Monte Carlo Compton profiles of solid and liquid lithium

We computed the Compton profile of solid and liquid lithium using quantum Monte Carlo (QMC) and compared with recent experimental measurements obtaining good agreement. Importantly, we find it crucial to account for proper core-valence orthogonalization and to address density differences when comparing with experiment. To account for disorder effects, we sampled finite-temperature configurations using molecular dynamics (MD), then performed diffusion Monte Carlo (DMC) simulations on each configuration. We used Slater-Jastrow wavefunctions and grand-canonical twist-averaged boundary conditions. A QMC pseudopotential correction, derived from an all-electron DMC simulation of the perfect crystal was also used. Our calculations provide the first all-electron QMC benchmark for the Compton profile of lithium crystal and pseudopotential-corrected QMC Compton profiles for both the liquid and solid.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Electron localization properties in high pressure hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo

We analyze in detail the electronic properties of high pressure hydrogen around the liquid-liquid phase transition based on Coupled Electron-Ion Monte Carlo calculations. Computing the off-diagonal single particle density matrix and the momentum distribution we discuss localization properties of the electrons. The abrupt changes of these distributions indicate a metal to insulator transition occurring together with the structural transition from the atomic to molecular fluid. We further discuss the electron-proton and electron-electron pair correlation functions, which also change abruptly at the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discovering correlated fermions using quantum Monte Carlo

It has become increasingly feasible to use quantum Monte Carlo (QMC) methods to study correlated fermion systems for realistic Hamiltonians. We give a summary of these techniques targeted at researchers in the field of correlated electrons, focusing on the fundamentals, capabilities, and current status of this technique. The QMC methods often offer the highest accuracy solutions available for systems in the continuum, and, since they address the many-body problem directly, the simulations can be analyzed to obtain insight into the nature of correlated quantum behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Theory of Finite Size Effects for Electronic Quantum Monte Carlo Calculations of Liquids and Solids

Concentrating on zero temperature Quantum Monte Carlo calculations of electronic systems, we give a general description of the theory of finite size extrapolations of energies to the thermodynamic limit based on one and two-body correlation functions. We introduce new effective procedures, such as using the potential and wavefunction split-up into long and short range functions to simplify the method and we discuss how to treat backflow wavefunctions. Then we explicitly test the accuracy of our method to correct finite size errors on example hydrogen and helium many-body systems and show that the finite size bias can be drastically reduced for even small systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Benchmarking Hydrogen-Helium Mixtures with QMC: Energetics, Pressures, and Forces

An accurate understanding of the phase diagram of dense hydrogen and helium mixtures is a crucial component in the construction of accurate models of Jupiter, Saturn, and Jovian extrasolar planets. Though DFT based first principles methods have the potential to provide the accuracy and computational efficiency required for this task, recent benchmarking in hydrogen has shown that achieving this accuracy requires a judicious choice of functional, and a quantification of the errors introduced. In this work, we present a quantum Monte Carlo based benchmarking study of a wide range of density functionals for use in hydrogen-helium mixtures at thermodynamic conditions relevant for Jovian planets. Not only do we continue our program of benchmarking energetics and pressures, but we deploy QMC based force estimators and use them to gain insights into how well the local liquid structure is captured by different density functionals. We find that TPSS, BLYP and vdW-DF are the most accurate functionals by most metrics, and that the enthalpy, energy, and pressure errors are very well behaved as a function of helium concentration. Beyond this, we highlight and analyze the major error trends and relative differences exhibited by the major classes of functionals, and estimate the magnitudes of these effects when possible.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Benchmark of Exchange-Correlation Functionals for High Pressure Hydrogen using Quantum Monte Carlo

The ab-initio phase diagram of dense hydrogen is very sensitive to errors in the treatment of electronic correlation. Recently, it has been shown that the choice of the density functional has a large effect on the predicted location of both the liquid-liquid phase transition and the solid insulator-to-metal transition in dense hydrogen. To identify the most accurate functional for dense hydrogen applications, we systematically benchmark some of the most commonly used functionals using Quantum Monte Carlo. By considering several measures of functional accuracy, we conclude that the van der Waals and hybrid functionals significantly out perform LDA and PBE. We support these conclusions by analyzing the impact of functional choice on structural optimization in the molecular solid, and on the location of the liquid-liquid phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Beyond the Born-Oppenheimer approximation with quantum Monte Carlo

In this work we develop tools that enable the study of non-adiabatic effects with variational and diffusion Monte Carlo methods. We introduce a highly accurate wave function ansatz for electron-ion systems that can involve a combination of both fixed and quantum ions. We explicitly calculate the ground state energies of H$_{2}$, LiH, H$_{2}$O and FHF$^{-}$ using fixed-node quantum Monte Carlo with wave function nodes that explicitly depend on the ion positions. The obtained energies implicitly include the effects arising from quantum nuclei and electron-nucleus coupling. We compare our results to the best theoretical and experimental results available and find excellent agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Molecular-Atomic Transition in the Deuterium Hugoniot with Coupled Electron Ion Monte Carlo

We have performed accurate simulations of the Deuterium Hugoniot using Coupled Electron Ion Monte Carlo (CEIMC). Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the principal Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we found a maximum compression of 4.85, roughly 10% larger than previous theoretical predictions and experimental data but still compatible with the latter because of their large uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exchange-correlation energy for the 3D homogeneous electron gas at arbitrary temperature

We fit finite-temperature path integral Monte Carlo calculations of the exchange-correlation energy of the 3D finite-temperature homogeneous electron gas in the warm-dense regime (r_{s} = (3/4πn)^{1/3} a_{B}^{-1} < 40 and Θ= T/T_{F} > 0.0625). In doing so, we construct a Padé approximant which collapses to Debye-Hückel theory in the high-temperature, low-density limit. Likewise, the zero-temperature limit matches the numerical results of ground-state quantum Monte Carlo, as well as analytical results in the high-density limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Quantum Energy Density: Improved Efficiency for Quantum Monte Carlo

We establish a physically meaningful representation of a quantum energy density for use in Quantum Monte Carlo calculations. The energy density operator, defined in terms of Hamiltonian components and density operators, returns the correct Hamiltonian when integrated over a volume containing a cluster of particles. This property is demonstrated for a helium-neon "gas," showing that atomic energies obtained from the energy density correspond to eigenvalues of isolated systems. The formation energies of defects or interfaces are typically calculated as total energy differences. Using a model of delta-doped silicon (where dopant atoms form a thin plane) we show how interfacial energies can be calculated more efficiently with the energy density, since the region of interest is small. We also demonstrate how the energy density correctly transitions to the bulk limit away from the interface where the correct energy is obtainable from a separate total energy calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Transition to the Metallic State in Low Density Hydrogen

Solid atomic hydrogen is one of the simplest systems to undergo a metal-insulator transition. Near the transition, the electronic degrees of freedom become strongly correlated and their description provides a difficult challenge for theoretical methods. As a result, the order and density of the phase transition are still subject to debate. In this work we use diffusion quantum Monte Carlo to benchmark the transition between paramagnetic and anti-ferromagnetic body centered cubic atomic hydrogen in its ground state. We locate the density of the transition by computing the equation of state for these two phases and identify the phase transition order by computing the band gap near the phase transition. These benchmark results show that the phase transition is continuous and occurs at a Wigner-Seitz radius of $r_s=2.27(3) a_0$. We compare our results to previously reported density functional theory, Hedin's GW approximation, and dynamical mean field theory results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Few-body reference data for multicomponent formalisms: Light nuclei molecules

We present full quantum statistical energetics of some electron-light nuclei systems. This is accomplished with the path integral Monte Carlo method. The effects on energetics arising from the change in the nuclear mass are studied. The obtained results may serve as reference data for the multicomponent density functional theory calculations of light nuclei system. In addition, the results reported here will enable better fitting of todays electron-nuclear energy functionals, for which the description of light nuclei is most challenging, in particular.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas

We perform calculations of the {3D} finite-temperature homogeneous electron gas (HEG) in the warm-dense regime ({r_{s} \equiv (3/4πn)^{1/3}a_{B}^{- 1} = 1.0- 40.0} and {Θ\equiv T/T_{F} = 0.0625- 8.0}) using restricted path integral Monte Carlo (RPIMC). Precise energies, pair correlation functions, and structure factors are obtained. For all densities, we find a significant discrepancy between the ground state parameterized local density approximation (LDA) and our results around {T_{F}} . These results can be used as a benchmark for improved functionals, as well as input for orbital-free DFT formulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Revealing the Condensate and Non-Condensate Distributions in the Inhomogeneous Bose-Hubbard Model

We calculate the condensate fraction and the condensate and non-condensate spatial and momentum distribution of the Bose-Hubbard model in a trap. From our results, it is evident that using approximate distributions can lead to erroneous experimental estimates of the condensate. Strong interactions cause the condensate to develop pedestal-like structures around the central peak that can be mistaken as non-condensate atoms. Near the transition temperature, the peak itself can include a significant non-condensate component. Using distributions generated from QMC simulations, experiments can map their measurements for higher accuracy in identifying phase transitions and temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite-size analysis of the Fermi liquid properties of the homogeneous electron gas

We analyze the extrapolation to the thermodynamic limit of Fermi liquid properties of the homogeneous electron gas in two and three dimensions. Using field theory, we explicitly calculate finite-size effects of the total energy, the renormalization factor, and the effective mass at the Fermi surface within the random phase approximation (RPA) and discuss the validity for general metallic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Free energy methods in Coupled Electron Ion Monte Carlo

Recent progress in simulation methodologies and in computer power allow first principle simulations of condensed systems with Born-Oppenheimer electronic energies obtained by Quantum Monte Carlo methods. Computing free energies and therefore getting a quantitative determination of phase diagrams is one step more demanding in terms of computer resources. In this paper we derive a general relation to compute the free energy of an ab-initio model with Reptation Quantum Monte Carlo (RQMC) energies from the knowledge of the free energy of the same ab-initio model in which the electronic energies are computed by the less demanding but less accurate Variational Monte Carlo (VMC) method. Moreover we devise a procedure to correct transition lines based on the use of the new relation. In order to illustrate the procedure, we consider the liquid-liquid phase transition in hydrogen, a first order transition between a lower pressure, molecular and insulating phase and a higher pressure, partially dissociated and conducting phase. We provide new results along the T = 600K isotherm across the phase transition and find good agreement between the transition pressure and specific volumes at coexistence for the model with RQMC accuracy between the prediction of our procedure and the values that can be directly inferred from the observed plateau in the pressure-volume curve along the isotherm. This work paves the way for future use of VMC in first principle simulations of high pressure hydrogen, an essential simplification when considering larger system sizes or quantum proton effects by Path Integral Monte Carlo methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-Temperature Superconductivity in Atomic Metallic Hydrogen

Superconductivity in the recently proposed ground-state structures of atomic metallic hydrogen is investigated over the pressure range 500 GPa to 3.5 TPa. Near molecular dissociation, the electron--phonon coupling $λ$ and renormalized Coulomb repulsion are similar to the molecular phase. A continuous increase in the critical temperature $T_c$ with pressure is therefore expected, to $\jmmapprox 356$K near 500 GPa. As the atomic phase stabilizes with increasing pressure, $λ$ increases, causing $T_c$ to approach $\jmmapprox 481$K near 700 GPa. At the first atomic--atomic structural phase transformation ($\jmmapprox 1$ -- 1.5 TPa), a discontinuous jump in $λ$ occurs, causing a significant increase in $T_c$ of up to 764K.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved Scaling for Quantum Monte Carlo on Insulators

Quantum Monte Carlo (QMC) methods are often used to calculate properties of many body quantum systems. The main cost of many QMC methods, for example the variational Monte Carlo (VMC) method, is in constructing a sequence of Slater matrices and computing the ratios of determinants for successive Slater matrices. Recent work has improved the scaling of constructing Slater matrices for insulators so that the cost of constructing Slater matrices in these systems is now linear in the number of particles, whereas computing determinant ratios remains cubic in the number of particles. With the long term aim of simulating much larger systems, we improve the scaling of computing the determinant ratios in the VMC method for simulating insulators by using preconditioned iterative solvers. The main contribution of this paper is the development of a method to efficiently compute for the Slater matrices a sequence of preconditioners that make the iterative solver converge rapidly. This involves cheap preconditioner updates, an effective reordering strategy, and a cheap method to monitor instability of ILUTP preconditioners. Using the resulting preconditioned iterative solvers to compute determinant ratios of consecutive Slater matrices reduces the scaling of QMC algorithms from O(n^3) per sweep to roughly O(n^2), where n is the number of particles, and a sweep is a sequence of n steps, each attempting to move a distinct particle. We demonstrate experimentally that we can achieve the improved scaling without increasing statistical errors. Our results show that preconditioned iterative solvers can dramatically reduce the cost of VMC for large(r) systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The momentum distribution of the homogeneous electron gas

We calculate the off-diagonal density matrix of the homogeneous electron gas at zero temperature using unbiased Reptation Monte Carlo for various densities and extrapolate the momentum distribution, and the kinetic and potential energies to the thermodynamic limit. Our results on the renormalization factor allows us to validate approximate G_0W_0 calculations concerning quasiparticle properties over a broad density region (1 <= r_s <= 10) and show that near the Fermi surface, vertex corrections and self-consistency aspects almost cancel each other out.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Itinerant ferromagnetism in a Fermi gas with contact interaction: Magnetic properties in a dilute Hubbard model

Ground state properties of the repulsive Hubbard model on a cubic lattice are investigated by means of the auxiliary-field quantum Monte Carlo method. We focus on low-density systems with varying on-site interaction $U/t$, as a model relevant to recent experiments on itinerant ferromagnetism in a dilute Fermi gas with contact interaction. Twist-average boundary conditions are used to eliminate open-shell effects and large lattice sizes are studied to reduce finite-size effects. The sign problem is controlled by a generalized constrained path approximation. We find no ferromagnetic phase transition in this model. The ground-state correlations are consistent with those of a paramagnetic Fermi liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Momentum Distribution and Renormalization Factor in Sodium and the Electron Gas

We present experimental and theoretical results on the momentum distribution and the quasiparticle renormalization factor in sodium. From an x-ray Compton-profile measurement of the valence-electron momentum density, we derive its discontinuity at the Fermi wavevector. This yields an accurate measure of the renormalization factor that we compare with quantum Monte Carlo and GW calculations performed both on crystalline sodium and on the homogeneous electron gas. Our calculated results are in good agreement with the experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Zero-Temperature Structures of Atomic Metallic Hydrogen

Ab initio random structure searching with density functional theory was used to determine the zero-temperature structures of atomic metallic hydrogen from 500 GPa to 5 TPa. Including zero point motion in the harmonic approximation, we estimate that molecular hydrogen dissociates into a monatomic body-centered tetragonal structure near 500 GPa (r_s = 1.225), which then remains stable to 2.5 TPa (r_s = 0.969). At higher pressures, hydrogen stabilizes in an ...ABCABC... planar structure that is remarkably similar to the ground state of lithium, which compresses to the face-centered cubic lattice beyond 5 TPa (r_s < 0.86). At this level of theory, our results provide a complete ab initio description of the atomic metallic structures of hydrogen, resolving one of the most fundamental and long outstanding issues concerning the structures of the elements.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Renormalization factor and effective mass of the two-dimensional electron gas

We calculate the momentum distribution of the Fermi liquid phase of the homogeneous, two-dimensional electron gas. We show that, close to the Fermi surface, the momentum distribution of a finite system with $N$ electrons approaches its thermodynamic limit slowly, with leading order corrections scaling as $N^{-1/4}$. These corrections dominate the extrapolation of the renormalization factor, $Z$, and the single particle effective mass, $m^*$, to the infinite system size. We show how convergence can be improved analytically. In the range $1 \le r_s \le 10$, we get a lower renormalization factor $Z$ and a higher effective mass, $m^*>m$, compared to the perturbative RPA values.

0 citations0 reviewsNo code linkedOpen access