Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study
Abstract
We develop an energy density matrix that parallels the onebody reduced density matrix (1RDM) for manybody quantum systems. Just as the density matrix gives access to the number density and occupation numbers, the energy density matrix yields the energy density and orbital occupation energies. The eigenvectors of the matrix provide a natural orbital partitioning of the energy density while the eigenvalues comprise a single particle energy spectrum obeying a total energy sum rule. For meanfield systems the energy density matrix recovers the exact spectrum. When correlation becomes important, the occupation energies resemble quasiparticle energies in some respects. We explore the occupation energy spectrum for the finite 3D homogeneous electron gas in the metallic regime and an isolated oxygen atom with ground state quantum Monte Carlo techniques imple mented in the QMCPACK simulation code. The occupation energy spectrum for the homogeneous electron gas can be described by an effective mass below the Fermi level. Above the Fermi level evanescent behavior in the occupation energies is observed in similar fashion to the occupation numbers of the 1RDM. A direct comparison with total energy differences demonstrates a quantita tive connection between the occupation energies and electron addition and removal energies for the electronmore »
 Authors:
 ORNL
 Publication Date:
 Research Org.:
 Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
 Sponsoring Org.:
 USDOE Office of Science (SC)
 OSTI Identifier:
 1148859
 DOE Contract Number:
 DEAC0500OR22725
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Physical Review B; Journal Volume: 90; Journal Issue: 3
 Country of Publication:
 United States
 Language:
 English
 Subject:
 quantum monte carlo; energy density; density matrix
Citation Formats
Krogel, Jaron T, Kim, Jeongnim, and Reboredo, Fernando A. Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study. United States: N. p., 2014.
Web. doi:10.1103/PhysRevB.90.035125.
Krogel, Jaron T, Kim, Jeongnim, & Reboredo, Fernando A. Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study. United States. doi:10.1103/PhysRevB.90.035125.
Krogel, Jaron T, Kim, Jeongnim, and Reboredo, Fernando A. Wed .
"Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study". United States.
doi:10.1103/PhysRevB.90.035125.
@article{osti_1148859,
title = {Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study},
author = {Krogel, Jaron T and Kim, Jeongnim and Reboredo, Fernando A},
abstractNote = {We develop an energy density matrix that parallels the onebody reduced density matrix (1RDM) for manybody quantum systems. Just as the density matrix gives access to the number density and occupation numbers, the energy density matrix yields the energy density and orbital occupation energies. The eigenvectors of the matrix provide a natural orbital partitioning of the energy density while the eigenvalues comprise a single particle energy spectrum obeying a total energy sum rule. For meanfield systems the energy density matrix recovers the exact spectrum. When correlation becomes important, the occupation energies resemble quasiparticle energies in some respects. We explore the occupation energy spectrum for the finite 3D homogeneous electron gas in the metallic regime and an isolated oxygen atom with ground state quantum Monte Carlo techniques imple mented in the QMCPACK simulation code. The occupation energy spectrum for the homogeneous electron gas can be described by an effective mass below the Fermi level. Above the Fermi level evanescent behavior in the occupation energies is observed in similar fashion to the occupation numbers of the 1RDM. A direct comparison with total energy differences demonstrates a quantita tive connection between the occupation energies and electron addition and removal energies for the electron gas. For the oxygen atom, the association between the ground state occupation energies and particle addition and removal energies becomes only qualitative. The energy density matrix provides a new avenue for describing energetics with quantum Monte Carlo methods which have traditionally been limited to total energies.},
doi = {10.1103/PhysRevB.90.035125},
journal = {Physical Review B},
number = 3,
volume = 90,
place = {United States},
year = {Wed Jan 01 00:00:00 EST 2014},
month = {Wed Jan 01 00:00:00 EST 2014}
}

Dynamical cluster quantum Monte Carlo study of the singleparticle spectra of strongly interacting fermion gases
We study the singleparticle spectral function of resonantly interacting fermions in the unitary regime, as described by the threedimensional attractive Hubbard model in the dilute limit. Our approach, based on the dynamical cluster approximation and the maximum entropy method, shows the emergence of a gap with decreasing temperature, as reported in recent coldatom photoemission experiments, for coupling values that span the BECBCS crossover. By comparing the behavior of the spectral function to that of the imaginary time dynamical pairing susceptibility, we attribute the development of the gap to the formation of local bound atom pairs. 
Interaction picture density matrix quantum Monte Carlo
The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the Nbody thermal density matrix and hence provides access to exact properties of manyparticle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a fourelectron system, comparing ourmore »