Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1307.5794 · Physical Review B

Quantum Monte Carlo study of the three-dimensional spin-polarized homogeneous electron gas

G. G. Spink, Richard J. Needs, N. D. Drummond

Published Aug 22, 2013• 95 citations
arXivHTMLPDF

Abstract

We have studied the spin-polarized three-dimensional homogeneous electron gas using the diffusion quantum Monte Carlo method, with trial wave functions including backflow and three-body correlations in the Jastrow factor, and we have used twist averaging to reduce finite-size effects. Calculations of the pair-correlation function, including the on-top pair density, as well as the structure factor and the total energy, are reported for systems of 118 electrons in the density range ${r}_{\mathrm{s}}=0.5$--20 a.u., and for spin polarizations of 0, 0.34, 0.66, and 1. We consider the spin resolution of the pair-correlation function and structure factor, and the energy of spin polarization. We show that a control variate method can reduce the variance when twist averaging, and we have achieved higher accuracy and lower noise than earlier quantum Monte Carlo studies.

DOI
10.1103/physrevb.88.085121

Related papers

Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid

N. D. Drummond, Richard J. Needs

We have used the variational and diffusion quantum Monte Carlo methods to calculate the energy, pair-correlation function, static structure factor, and momentum density of the ground state of the two-dimensional (2D) homogeneous electron gas. We have used highly accurate Slater-Jastrow-backflow trial wave functions and twist averaging to reduce finite-size effects where applicable. We compare our results with others in the literature and construct a local-density-approximation exchange-correlation functional for 2D systems.

arXiv:1002.2122•Physical Review B•Feb 13, 2009•40 citations

Correlation energy of the spin-polarized electron liquid studied using quantum Monte Carlo simulations

Sam Azadi, N. D. Drummond, Sam M. Vinko

Variational and diffusion quantum Monte Carlo (VMC and DMC) methods with Slater-Jastrow-backflow trial wave functions are used to study the spin-polarized three-dimensional uniform electron fluid. We report ground state VMC and DMC energies in the density range $0.5\ensuremath{\le}{r}_{\text{s}}\ensuremath{\le}20$. Finite-size errors are corrected using canonical-ensemble twist-averaged boundary conditions and extrapolation of the twist-averaged energy per particle calculated at three system sizes $(N=113, 259, \mathrm{and} 387)$ to the thermodynamic limit of infinite system size. The DMC energies in the thermodynamic limit are used to parametrize a local spin density approximation correlation function for inhomogeneous electron systems.

arXiv:2306.13426•Physical review. B./Physical review. B•Sep 15, 2023•2 citations

Application of quantum Monte Carlo methods to homogeneous electron and electron-hole systems

G. G. Spink

The properties of the macroscopic world around us, and of which we are a part, are largely determined by the low energy, collective behaviour of many interacting particles, including the nuclei and, especially, the electrons present. Although the fundamental laws governing the behaviour of these many-body systems are believed to be known in principle, the practical solution of the equations of quantum mechanics remains a challenging area of research. This thesis is concerned with the application of quantum Monte Carlo methods to two model systems: the spin-polarised homogeneous electron gas, and a hole-doped electron gas. Electronic structure theory is briefly reviewed before discussing in more detail the quantum Monte Carlo methods used in this thesis. A study of the three-dimensional spin-polarised homogeneous electron gas (HEG) is then reported, where the relatively new technique of twist averaging is investigated in detail and accurate energies and pair correlation functions are obtained over densities $r_s = 0.5 – 20$ a.u. and the full range of spin-polarisation, allowing comparison with the Perdew-Zunger interpolation scheme used in local spin density approximation exchange-correlation functionals. Following this, an impurity is added to the electron gas in the form of a positively charged hole, and the interaction is studied. Relaxation energies, pair correlation functions and momentum densities are reported. Trion formation is observed over a range of carrier densities and electron-hole mass ratios in agreement with experiment. Isolated trions are also studied, where the diffusion Monte Carlo method is exact. Methodological innovations developed while carrying out this work are discussed, including a variance reduction technique for twist-averaged calculations and a new trial wave function for impurity-in-HEG calculations.

•Dec 22, 2017•0 citations

Low-density phase diagram of the three-dimensional electron gas

Sam Azadi, N. D. Drummond

Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of the three-dimensional homogeneous electron gas at very low density. Fermi fluid and body-centered cubic Wigner crystal ground-state energies are determined using Slater-Jastrow-backflow and Slater-Jastrow many-body wave functions at different densities and spin polarizations in finite simulation cells. Finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Unlike previous studies, our results show that the electron gas undergoes a first-order quantum phase transition directly from a paramagnetic fluid to a body-centered cubic crystal at density parameter ${r}_{\text{s}}=86.6(7)$, with no region of stability for an itinerant ferromagnetic fluid.

arXiv:2201.08743•Physical review. B./Physical review. B•Jun 23, 2022•31 citations

Correlation energy of the paramagnetic electron gas at the thermodynamic limit

Sam Azadi, N. D. Drummond, Sam M. Vinko

The variational and diffusion quantum Monte Carlo methods are used to calculate the correlation energy of the paramagnetic three-dimensional homogeneous electron gas at intermediate to high density. Ground-state energies in finite cells are determined using Slater-Jastrow-backflow trial wave functions, and finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Our correlation energies in the thermodynamic limit are more accurate than previous results. The present diffusion quantum Monte Carlo energies together with our recently reported [Phys. Rev. B 105, 245135 (2022)] results at low density, are used to parametrize the correlation energy of the electron gas using a functional form that satisfies the exact asymptotic behavior at high density.

arXiv:2209.10227•Physical review. B./Physical review. B•Mar 15, 2023•9 citations

Effects of backflow correlation in the three-dimensional electron gas: Quantum Monte Carlo study

Yongkyung Kwon, David M. Ceperley, Richard M. Martin

The correlation energy of the homogeneous three-dimensional interacting electron gas is calculated using the variational and fixed-node diffusion Monte Carlo methods, with trial functions that include backflow and three-body correlations. In the high-density regime ${(r}_{s}<~5)$ the effects of backflow dominate over those due to three-body correlations, but the relative importance of the latter increases as the density decreases. Since the backflow correlations vary the nodes of the trial function, this leads to improved energies in the fixed-node diffusion Monte Carlo calculations. The effects are comparable to those found for the two-dimensional electron gas, leading to much improved variational energies and fixed-node diffusion energies similar to the released-node energies of Ceperley and Alder.

arXiv:cond-mat/9803092•Physical review. B, Condensed matter•Sep 15, 1998•236 citations
We have studied the spin-polarized three-dimensional homogeneous electron gas using the diffusion quantum Monte Carlo method, with trial wave functions including backflow and three-body correlations in the Jastrow factor, and we have used twist averaging to reduce finite-size effects. Calculations of the pair-correlation function, including the on-top pair density, as well as the structure factor and the total energy, are reported for systems of 118 electrons in the density range rs=0.5{r}_{\mathrm{s}}=0.5rs​=0.5--20 a.u., and for spin polarizations of 0, 0.34, 0.66, and 1. We consider the spin resolution of the pair-correlation function and structure factor, and the energy of spin polarization. We show that a control variate method can reduce the variance when twist averaging, and we have achieved higher accuracy and lower noise than earlier quantum Monte Carlo studies.
We have used the variational and diffusion quantum Monte Carlo methods to calculate the energy, pair-correlation function, static structure factor, and momentum density of the ground state of the two-dimensional (2D) homogeneous electron gas. We have used highly accurate Slater-Jastrow-backflow trial wave functions and twist averaging to reduce finite-size effects where applicable. We compare our results with others in the literature and construct a local-density-approximation exchange-correlation functional for 2D systems.
Variational and diffusion quantum Monte Carlo (VMC and DMC) methods with Slater-Jastrow-backflow trial wave functions are used to study the spin-polarized three-dimensional uniform electron fluid. We report ground state VMC and DMC energies in the density range $0.5\ensuremath{\le}{r}_{\text{s}}\ensuremath{\le}20$. Finite-size errors are corrected using canonical-ensemble twist-averaged boundary conditions and extrapolation of the twist-averaged energy per particle calculated at three system sizes (N=113,259,and387)(N=113, 259, \mathrm{and} 387)(N=113,259,and387) to the thermodynamic limit of infinite system size. The DMC energies in the thermodynamic limit are used to parametrize a local spin density approximation correlation function for inhomogeneous electron systems.
The properties of the macroscopic world around us, and of which we are a part, are largely determined by the low energy, collective behaviour of many interacting particles, including the nuclei and, especially, the electrons present. Although the fundamental laws governing the behaviour of these many-body systems are believed to be known in principle, the practical solution of the equations of quantum mechanics remains a challenging area of research. This thesis is concerned with the application of quantum Monte Carlo methods to two model systems: the spin-polarised homogeneous electron gas, and a hole-doped electron gas. Electronic structure theory is briefly reviewed before discussing in more detail the quantum Monte Carlo methods used in this thesis. A study of the three-dimensional spin-polarised homogeneous electron gas (HEG) is then reported, where the relatively new technique of twist averaging is investigated in detail and accurate energies and pair correlation functions are obtained over densities rs=0.5–20r_s = 0.5 – 20rs​=0.5–20 a.u. and the full range of spin-polarisation, allowing comparison with the Perdew-Zunger interpolation scheme used in local spin density approximation exchange-correlation functionals. Following this, an impurity is added to the electron gas in the form of a positively charged hole, and the interaction is studied. Relaxation energies, pair correlation functions and momentum densities are reported. Trion formation is observed over a range of carrier densities and electron-hole mass ratios in agreement with experiment. Isolated trions are also studied, where the diffusion Monte Carlo method is exact. Methodological innovations developed while carrying out this work are discussed, including a variance reduction technique for twist-averaged calculations and a new trial wave function for impurity-in-HEG calculations.
Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of the three-dimensional homogeneous electron gas at very low density. Fermi fluid and body-centered cubic Wigner crystal ground-state energies are determined using Slater-Jastrow-backflow and Slater-Jastrow many-body wave functions at different densities and spin polarizations in finite simulation cells. Finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Unlike previous studies, our results show that the electron gas undergoes a first-order quantum phase transition directly from a paramagnetic fluid to a body-centered cubic crystal at density parameter rs=86.6(7){r}_{\text{s}}=86.6(7)rs​=86.6(7), with no region of stability for an itinerant ferromagnetic fluid.
The variational and diffusion quantum Monte Carlo methods are used to calculate the correlation energy of the paramagnetic three-dimensional homogeneous electron gas at intermediate to high density. Ground-state energies in finite cells are determined using Slater-Jastrow-backflow trial wave functions, and finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Our correlation energies in the thermodynamic limit are more accurate than previous results. The present diffusion quantum Monte Carlo energies together with our recently reported [Phys. Rev. B 105, 245135 (2022)] results at low density, are used to parametrize the correlation energy of the electron gas using a functional form that satisfies the exact asymptotic behavior at high density.
The correlation energy of the homogeneous three-dimensional interacting electron gas is calculated using the variational and fixed-node diffusion Monte Carlo methods, with trial functions that include backflow and three-body correlations. In the high-density regime (rs< 5){(r}_{s}<~5)(rs​< 5) the effects of backflow dominate over those due to three-body correlations, but the relative importance of the latter increases as the density decreases. Since the backflow correlations vary the nodes of the trial function, this leads to improved energies in the fixed-node diffusion Monte Carlo calculations. The effects are comparable to those found for the two-dimensional electron gas, leading to much improved variational energies and fixed-node diffusion energies similar to the released-node energies of Ceperley and Alder.