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

David M. Ceperley

David M. Ceperley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elcond-mat.mtrl-sciphysics.comp-phcond-mat.quant-gascond-mat.supr-conmath.NA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
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

9 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
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
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
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