Researcher profile

D. M. Ceperley

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

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

12 published item(s)

preprint2013arXiv

Nuclear Quantum Effects and Nonlocal Exchange-Correlation Functionals Applied to Liquid Hydrogen at High Pressure

Using first-principles molecular dynamics, we study the influence of nuclear quantum effects (NQEs) and nonlocal exchange--correlation density functionals (DFs) near molecular dissociation in liquid hydrogen. NQEs strongly influence intramolecular properties, such as bond stability, and are thus an essential part of the dissociation process. Moreover, by including DFs that account for either the self-interaction error or dispersion interactions, we find a much better description of molecular dissociation and metallization than previous studies based on classical protons and/or local or semi-local DFs. We obtain excellent agreement with experimentally measured optical properties along pre-compressed Hugoniots, and while we still find a first-order liquid--liquid transition at low temperatures, transition pressures are increased by more than 100 GPa.

preprint2011arXiv

The superfluid-insulator transition in the disordered two-dimensional Bose-Hubbard model

We investigate the superfluid-insulator transition in the disordered two-dimensional Bose-Hubbard model through quantum Monte Carlo simulations. The Bose-Hubbard model is studied in the presence of site disorder and the quantum critical point between the Bose-glass and superfluid is determined in both the grand canonical ($μ/U=0.375$ close to $ρ=1$) and canonical ensemble ($ρ=1$ and 0.5). Particular attention is paid to disorder averaging and it is shown that an extremely large number of disorder realizations are needed in order to obtain reliable results. Typically we average over more than $100,000$ disorder realizations. In the grand canonical ensemble we find $Z t_c/U=0.112(1)$ with $μ/U=0.375$, significantly different from previous studies. When compared to the critical point in the absence of disorder ($Z t_c/U=0.2385$), this result confirms previous findings showing that disorder enlarges the superfluid region. At the critical point, in order to estimate universal features, we compute the dynamic conductivity scaling curves.

preprint2011arXiv

Validity of the scattering length approximation in strongly interacting Fermi systems

We investigate the energy spectrum of systems of two, three and four spin-1/2 fermions with short range attractive interactions both exactly, and within the scattering length approximation. The formation of molecular bound states and the ferromagnetic transition of the excited scattering state are examined systematically as a function of the 2-body scattering length. Identification of the upper branch (scattering states) is discussed and a general approach valid for systems with many particles is given. We show that an adiabatic ferromagnetic transition occurs, but at a critical transition point kF a much higher than predicted from previous calculations, almost all of which use the scattering length approximation. In the 4-particle system the discrepancy is a factor of 2. The exact critical interaction strength calculated in the 4-particle system is consistent with that reported by experiment. To make comparisons with the adiabatic transition, we study the quench dynamics of the pairing instability using the eigenstate wavefunctions.

preprint2010arXiv

Electrical conductivity of high-pressure liquid hydrogen by quantum Monte Carlo methods

We compute the electrical conductivity for liquid hydrogen at high pressure using quantum Monte Carlo. The method uses Coupled Electron-Ion Monte Carlo to generate configurations of liquid hydrogen. For each configuration correlated sampling of electrons is performed in order to calculate a set of lowest many-body eigenstates and current-current correlation functions of the system, which are summed over in the many-body Kubo formula to give AC electrical conductivity directly. The extrapolated DC conductivity at 3000 K for several densities shows a liquid semiconductor to liquid-metal transition at high pressure. Our results are in good agreement with shock-wave data.

preprint2009arXiv

Construction of localized wave functions for a disordered optical lattice and analysis of the resulting Hubbard model parameters

We propose a method to construct localized single particle wave functions using imaginary time projection and thereby determine lattice Hamiltonian parameters. We apply the method to a specific disordered potential generated by an optical lattice experiment and calculate for each instance of disorder, the equivalent lattice model parameters. The probability distributions of the Hubbard parameters are then determined. Tests of localization and eigen-energy convergence are examined.

preprint2009arXiv

Hexatic and mesoscopic phases in the 2D quantum Coulomb system

We study the Wigner crystal melting in a two dimensional quantum system of particles interacting via the 1/r Coulomb potential. We use quantum Monte Carlo methods to calculate its phase diagram, locate the Wigner crystal region, and analyze its instabilities towards the liquid phase. We discuss the role of quantum effects in the critical behavior of the system, and compare our numerical results with the classical theory of melting, and the microemulsion theory of frustrated Coulomb systems. We find a Pomeranchuk effect much larger then in solid helium. In addition, we find that the exponent for the algebraic decay of the hexatic phase differs significantly from the Kosterilitz-Thouless theory of melting. We search for the existence of mesoscopic phases and find evidence of metastable bubbles but no mesoscopic phase that is stable in equilibrium.

preprint2006arXiv

Many-body wavefunctions for normal liquid $^3$He

We present new trial wave-functions which include 3-body correlations into the backflow coordinates and a 4-body symmetric potential. We show that our wavefunctions lower the energy enough to stabilize the ground state energies of normal liquid $^3$He in the unpolarized state at all pressures in agreement with experiment; however, quantitative discrepancies remain. Further, we include strong spin coupling into the Fermi liquid by adapting pairing wave functions. We demonstrate a new, numerically stable method to evaluate pairing functions which is also useful for Path Integrals calculations at low, but non-zero temperatures.

preprint2006arXiv

Theory of Small Para-Hydrogen Clusters: Magic Numbers and Superfluid Sizes

The interplay between magic number stabilities and superfluidity of small para-hydrogen clusters with sizes $N = 5$ to 40 and temperatures $0.5 K \leq T \leq 4.5 $K is explored with classical and quantum Path Integral Monte Carlo calculations. Clusters with $N < 26$ and T $\leq 1.5 K$ have large superfluid fractions even at the stable magic numbers 13, 19, and 23. In larger clusters, superfluidity is quenched especially at the magic numbers 23, 26, 29, 32, and 37 while below 1 K, superfluidity is recovered for the pairs $(27,28)$, $(30,31)$, and $(35,36)$. For all clusters superfluidity is localized at the surface and correlates with long exchange cycles involving loosely bound surface molecules.

preprint2003arXiv

Multi-spin exchange model the near melting transition of the 2D Wigner crystal

The low temperature properties of fermionic solids are governed by spin exchanges. Near the melting transition, the spin-exchange energy increases as well as the relative contribution of large loops. In this paper, we check the convergence of the multi-spin exchange model and the validity of the Thouless theory near the melting of the Wigner crystal in two dimensions. Exchange energies are computed using Path Integral Monte Carlo for loop sizes up to 8, at $r_s=40$, 50 and 75. The data are then fitted to a geometric model. Then the exchange energies are extrapolated to larger sizes in order to evaluate their contributions to the leading term of the magnetic susceptibility and the specific heat. These results are used to check the convergence of the multi-spin exchange model.

preprint2000arXiv

Electronic structure of self-assembled quantum dots: comparison between density functional theory and diffusion quantum Monte Carlo

We have calculated the exchange, correlation, and total electronic energy of a realistic InAs self-assembled quantum dot embedded in a GaAs matrix as a function of the number of electrons in the dot. The many-body interactions have been treated using the local spin density approximation (LSDA) to density functional theory (DFT) and diffusion quantum Monte Carlo (DMC), so that we may quantify the error introduced by LSDA. The comparison shows that the LSDA errors are about 1-2 meV per electron for the system considered. These errors are small enough to justify the use of LSDA calculations to test models of self-assembled dots against current experimental measurements.

preprint2000arXiv

Quantum Monte Carlo treatment of elastic exciton-exciton scattering

We calculate cross sections for low energy elastic exciton-exciton scattering within the effective mass approximation. Unlike previous theoretical approaches, we give a complete, non-perturbative treatment of the four-particle scattering problem. Diffusion Monte Carlo is used to calculate the essentially exact energies of scattering states, from which phase shifts are determined. For the case of equal-mass electrons and holes, which is equivalent to positronium-positronium scattering, we find a_s = 2.1 a_x for scattering of singlet-excitons and a_s= 1.5 a_x for triplet-excitons, where a_x is the excitonic radius. The spin dependence of the cross sections arises from the spatial exchange symmetry of the scattering wavefunctions. A significant triplet-triplet to singlet-singlet scattering process is found, which is similar to reported effects in recent experiments and theory for excitons in quantum wells. We also show that the scattering length can change sign and diverge for some values of the mass ratio m_h/m_e, an effect not seen in previous perturbative treatments.

preprint1999arXiv

Path Integral Monte Carlo Simulations for Fermion Systems: Pairing in the Electron-Hole Plasma

We review the path integral method wherein quantum systems are mapped with Feynman&#39;s path integrals onto a classical system of &#34;ring-polymers&#34; and then simulated with the Monte Carlo technique. Bose or Fermi statistics correspond to possible &#34;cross-linking&#34; of polymers. As proposed by Feynman, superfluidity and Bose condensation result from macroscopic exchange of bosons. To map fermions onto a positive probability distribution, one must restrict the paths to lie in regions where the fermion density matrix is positive. We discuss a recent application to the two-component electron-hole plasma. At low temperature excitons and bi-excitons form. We have used nodal surfaces incorporating paired fermions and see evidence of a Bose condensation in the energy, specific heat and superfluid density. In the restricted path integral picture, pairing appears as intertwined electron-hole paths. Bose condensation occurs when these intertwined paths wind around the periodic boundaries.