Source author record

David Ceperley

David 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.mtrl-sci cond-mat.other physics.chem-ph physics.comp-ph quant-ph cond-mat.quant-gas cond-mat.stat-mech physics.plasm-ph

Catalog footprint

What is connected

7works

8topics

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

Phonons of metallic hydrogen with quantum Monte Carlo

We describe a simple scheme to perform phonon calculations with quantum Monte Carlo (QMC) methods, and demonstrate it on metallic hydrogen. Because of the energy and length scales of metallic hydrogen, and the statistical noise inherent to QMC methods, the conventional manner of calculating force constants is prohibitively expensive. We show that our alternate approach is nearly 100 times more efficient in resolving the force constants needed to calculate the phonon spectrum in the harmonic approximation. This requires only the calculation of atomic forces, as in the conventional approach, and otherwise little or no programmatic modifications.

preprint2016arXiv

Interpolated wave functions for nonadiabatic simulations with the fixed-node quantum Monte Carlo method

Simulating nonadiabatic effects with many-body wave function approaches is an open field with many challenges. Recent interest has been driven by new algorithmic developments and improved theoretical understanding of properties unique to electron-ion wave functions. Fixed-node diffusion Monte Caro is one technique that has shown promising results for simulating electron-ion systems. In particular, we focus on the CH molecule for which previous results suggested a relatively significant contribution to the energy from nonadiabatic effects. We propose a new wave function ansatz for diatomic systems which involves interpolating the determinant coefficients calculated from configuration interaction methods. We find this to be an improvement beyond previous wave function forms that have been considered. The calculated nonadiabatic contribution to the energy in the CH molecule is reduced compared to our previous results, but still remains the largest among the molecules under consideration.

preprint2015arXiv

How Large are Nonadiabatic Effects in Atomic and Diatomic Systems?

With recent developments in simulating nonadiabatic systems to high accuracy, it has become possible to determine how much energy is attributed to nuclear quantum effects beyond zero-point energy. In this work we calculate the non-relativistic ground-state energies of atomic and molecular systems without the Born-Oppenheimer approximation. For this purpose we utilize the fixed-node diffusion Monte Carlo method, in which the nodes depend on both the electronic and ionic positions. We report ground-state energies for all systems studied, ionization energies for the first-row atoms and atomization energies for the first-row hydrides. We find the ionization energies of the atoms to be nearly independent of the Born-Oppenheimer approximation, within the accuracy of our results. The atomization energies of molecular systems, however, show small effects of the nonadiabatic coupling between electrons and nuclei.

preprint2014arXiv

Metastable Bose-Einstein Condensation in a Strongly Correlated Optical Lattice

We experimentally and theoretically study the peak fraction of a Bose-Einstein condensate loaded into a cubic optical lattice as the lattice potential depth and entropy per particle are varied. This system is well-described by the superfluid regime of the Bose-Hubbard model, which allows for comparison with mean-field theories and exact quantum Monte Carlo (QMC) simulations. Despite correcting for systematic discrepancies between condensate fraction and peak fraction, we discover that the experiment consistently shows the presence of a condensate at temperatures higher than the critical temperature predicted by QMC simulations. This metastability suggests that turning on the lattice potential is non-adiabatic. To confirm this behavior, we compute the timescales for relaxation in this system, and find that equilibration times are comparable with the known heating rates. The similarity of these timescales implies that turning on the lattice potential adiabatically may be impossible. Our results point to the urgent need for a better theoretical and experimental understanding of the timescales for relaxation and adiabaticity in strongly interacting quantum gases, and the importance of model-independent probes of thermometry in optical lattices.

preprint2012arXiv

Electronic Energy Functionals: Levy-Lieb principle within the Ground State Path Integral Quantum Monte Carlo

We propose a theoretical/computational protocol based on the use of the Ground State (GS) Path Integral (PI) Quantum Monte Carlo (QMC) for the calculation of the kinetic and Coulomb energy density for a system of $N$ interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the $N-1$-conditional probability density within the framework of the Levy-Lieb constrained search principle. The consequences for the development of energy functionals within the context of Density Functional Theory (DFT) are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the $N-1$-conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an {\it on-the-fly} determination of the energy functional for any system.

preprint2009arXiv

Equation of state of metallic hydrogen from Coupled Electron-Ion Monte Carlo simulations

We present a study of hydrogen at pressures higher than molecular dissociation using the Coupled Electron-Ion Monte Carlo method. These calculations use the accurate Reptation Quantum Monte Carlo method to estimate the electronic energy and pressure while doing a Monte Carlo simulation of the protons. In addition to presenting simulation results for the equation of state over a large region of phase space, we report the free energy obtained by thermodynamic integration. We find very good agreement with DFT calculations for pressures beyond 600 GPa and densities above $ρ=1.4 g/cm^3$. Both thermodynamic as well as structural properties are accurately reproduced by DFT calculations. This agreement gives a strong support to the different approximations employed in DFT, specifically the approximate exchange-correlation potential and the use of pseudopotentials for the range of densities considered. We find disagreement with chemical models, which suggests a reinvestigation of planetary models, previously constructed using the Saumon-Chabrier-Van Horn equations of state.

preprint1999arXiv

Characterization of the State of Hydrogen

Fermionic path integral Monte Carlo simulations have been applied to study the equilibrium properties of the hydrogen and deuterium in the density and temperature range of 1.6 < rs < 14.0 and 5000K < T < 167000K. We use this technique to determine the phase diagram by identifying the plasma, the molecular, atomic and metallic regime. We explain how one can identify the phases in the path integral formalism and discuss the state of hydrogen for 5 points in the temperature-density plane. Further we will provide arguments for the nature of the transitions between the regimes.