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...
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.
This article describes the variational and fixed-node diffusion quantum Monte Carlo methods and how they may be used to calculate the properties of many-electron systems. These stochastic wave-function-based approaches provide a very direct treatment of quantum many-body effects and serve as benchmarks against which other techniques may be compared. They complement the less demanding density-functional approach by providing more accurate results and a deeper understanding of the physics of electronic correlation in real materials. The algorithms are intrinsically parallel, and currently available high-performance computers allow applications to systems containing a thousand or more electrons. With these tools one can study complicated problems such as the properties of surfaces and defects, while including electron correlation effects with high precision. The authors provide a pedagogical overview of the techniques and describe a selection of applications to ground and excited states of solids and clusters.
Two of the most influential ideas developed by Richard Feynman are the Feynman diagram technique and his variational approach. The former provides a powerful tool to construct a systematic expansion for a generic interacting system, while the latter allows optimization of a perturbation theory using a variational principle. Here we show that combining a variational approach with a new diagrammatic quantum Monte Carlo method, both based on the Feynman's original ideas, results in a powerful and accurate solver to the generic solid state problem, in which a macroscopic number of electrons interact by the long range Coulomb repulsion. We apply the solver to the quintessential problem of solid state, the uniform electron gas (UEG), which is at the heart of the density functional theory (DFT) success in describing real materials, yet it has not been adequately solved for over 90 years. While some wave-function properties, like the ground state energy, have been very accurately calculated by the diffusion Monte Carlo method (DMC), the static and dynamic response functions, which are directly accessed by the experiment, and are crucial for development of non-local and non-adiabatic density functionals, remain poorly understood. Our method allows us to calculate the momentum-frequency resolved spin response functions for the first time, and to improve on the precision of the charge response function. The accuracy of both response functions is sufficiently high, so as to uncover previously missed fine structure in these responses. This method can be straightforwardly applied to a large number of moderately interacting electron systems in the thermodynamic limit, including realistic models of metallic and semiconducting solids.
to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.
The fixed-node diffusion quantum Monte Carlo (DMC) method is the most accurate method known for calculating the energies of large many-particle quantum systems. The key element of the method is the development of accurate trial many-body wavefunctions which control the statistical efficiency of the calculations and the accuracy obtained. Accurate wavefunctions can be obtained by building correlation effects on top of mean field descriptions such as density functional theory. The wavefunctions can be improved by introducing multi-determinants, pairing functions, and backflow transformations. The calculations are expensive, but the method scales well with system size and calculations on 1000 particles are possible. Some recent applications of the DMC method to atoms, molecules and solids will be presented.
Variational and diffusion Monte Carlo (VMC and DMC) are widely used techniques in the Monte Carlo study of many-body quantum systems. Through efficient use of ran- dom sampling to solve many-dimensional integrals, they can be used to reliably solve the Schrödinger equation. Within the wider ecosystems of electronic structure calculation, quantum chemistry, and quantum simulation, quantum Monte Carlo (QMC) methods are used in tandem with density functional theory (DFT) and other computational methods to deliver insight into a vast array of condensed matter systems. Of particular interest to this work are systems in the intermediate correlation regime, where methods such as simple band theory and the Hubbard model break down. This thesis proposes and tests an explicit three-body backflow term to improve the accuracy of VMC and DMC calculations by incorporating higher-order correlation in order to reduce fixed-node er- ror in DMC, finding statistically significant improvements in ground state energies and improved variance in a water molecule. Then, QMC is used to calculate accurate charge densities for weak and intermediately correlated materials, to test commonly used dens- ity functional approximations in a controlled environment, providing insight into where these approximations perform well, and where they fail. Finally, DMC is used to generate ground state energy data for a dipolar Fermi gas of dressed Rydberg atoms interacting in the intermediate correlation regime. This data is parametrised in terms of interaction strength, for use in generating exchange-correlation functionals, to enable the DFT study of inhomogeneous systems of dressed Rydberg atoms.