We present a new method for realizing the adiabatic connection approach in density functional theory, which is based on combining accurate variational quantum Monte Carlo calculations with a constrained optimization of the ground state many-body wavefunction for different values of the Coulomb coupling constant. We use the method to study an electron gas in the presence of a cosine-wave potential. For this system we present results for the exchange-correlation hole and exchange-correlation energy density, and compare our findings with those from the local density approximation and generalized gradient approximation. 1. Introduction 2. The Adiabatic connection 3. Quantum Monte Carlo realization 3.1 Variational Monte Carlo 3.2 Fixed-density variance minimization 4. Cosine-wave jellium 4.1 Many-body wavefunction 5. Results and discussion 1 Introduction Density functional theory (DFT) (1, 2) is the main computational tool for the treatment of many-body effects in solid state electroni...
Abstract It is difficult to imagine that density functional theory (DFT) could be useful in improving quantum Monte Carlo methods, but it is easy to imagine that QMC could be useful in improving density functional methods. Since QMC can produce accurate wavefunctions and since wavefunctions lead to all other properties one should be able, in principle, to use QMC to predict the exchangecorrelation energy functional Exc[n(r)], the key (unknown) quantity required in DFT. In practice, that is not so easy, but this paper shows how QMC can take steps in that direction by providing insights gained from comparisons of results from QMC and DFT. The authors considered the case of the strongly inhomogeneous electron gas using variational Q:tvIC to obtain wavefunctions 1/J>. and energies E[7/J>.] for the ground state of the Hamiltonian iI>- = T +AVee+V >. associated with the coupling constant A in the range 0 to 1. Optimization of a Slater-Jastrow many-body wavefunction was made by minimizing a penalty function given by a weighted sum of the variance in local energies and a sum of squares of electron density differences. The variational parameters in 7/J>. and v>. were repeatedly varied to reach convergence followed by additional calculations to obtain the local values of the exchange-correlation energies e c and the exchange-correlation holes n c from which Exe and nxc could be determined.
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–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.