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:1809.04651

Feynman's solution of the quintessential problem in solid state physics

Kun Chen, Kristjan Haule

Published Sep 12, 2018• 4 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.1809.04651

Related papers

A combined variational and diagrammatic quantum Monte Carlo approach to the many-electron problem

Kun Chen, Kristjan Haule

Two of the most influential ideas developed by Richard Feynman are the Feynman diagram technique and his variational approach. Here we show that combining both, and introducing a diagrammatic quantum Monte Carlo method, 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 it to the quintessential problem of solid state, the uniform electron gas, which is at the heart of the density functional theory success in describing real materials, yet it has not been adequately solved for over 90 years. Our method allows us to calculate numerically exact momentum and frequency resolved spin and charge response functions. This method can be applied to a number of moderately interacting electron systems, including models of realistic metallic and semiconducting solids.

Nature Communications•Aug 19, 2019•70 citations

New Applications of Quantum Monte Carlo

Nicholas D. M. Hine

This thesis comprises a trio of separate yet thematically-linked investigations within the field of electronic structure that address challenging problems to which quantum Monte Carlo, the most accurate electronic structure method capable of treating large extended systems, has been applied only rarely or indeed never before. The first is the study of polarization and localization in many body systems. Density functional and variational quantum Monte Carlo methods are used to examine the behaviour of the many-electron localization length near band insulator to metal transitions in various oneand two-dimensional model systems. We examine when the localization length, known to be infinite in metals and finite in insulators and assumed to diverge at a metal-insulator transition, serves as a useful diagnostic of the transition. We also present a novel derivation of the formulae for polarization and localization in many body systems that expresses the connection of these two quantities to the concept of Many Body Wannier Functions. The second is a study of the surface energy of the electron gas. The electron gas, or ‘jellium’, is an important benchmark system for electronic structure methods and serves as a simple model for metals. We apply diffusion Monte Carlo (DMC), and a new method for determining surface energies without an unreliable comparison of bulk and slab results, to the calculation of the surface energy. Our method improves on previous calculations in a number of ways and removes many sources of error that plagued previous calculations. The final study applies DMC to the properties, specifically the formation energy, of point defects in Alumina. Alumina is an important technological material, and there are many important unanswered questions relating to its properties that stem from the behaviour and concentration of its point defects. We address many of the problems associated with the calculation of formation energies with electronic structure calculations, and present the first calculations of formation energies of charged point defects in DMC.

•Jan 1, 2007•0 citations

Diagrammatic Methods for Strongly-Correlated Electrons in Realistic Systems

Jia Li

This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.

•Jan 1, 2021•0 citations

A Quantum Monte Carlo Approach to the Adiabatic Connection Method

Maziar Nekovee, W. M. C. Foulkes, A. J. Williamson +3 more

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...

•Jan 1, 1999•8 citations

New variational approach to the many-body problem and application to the high-density electron gas

Yasutami Takada

A new variational approach is proposed for the study of the many-body problem with the purpose of combining the advantages of perturbation-theoretic and variational approaches, while avoiding their disadvantages. A trial wave function is constructed so that a general expansion for the expectation value of an arbitrary physical quantity can be derived. Each term in the expansion can be calculated with the use of the techniques of Feynman diagrams, but in contrast with usual perturbation-theoretic approaches, every term in the expansion always converges. The procedure for applying the method to real systems is illustrated by a calculation of the ground-state energy of the electron gas at high densities.

Physical Review A•Oct 1, 1983•20 citations

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel

J. P. F. LeBlanc, Kun Chen, Kristjan Haule +2 more

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

arXiv:2205.13595•Physical Review Letters•Dec 7, 2022•31 citations
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.
Two of the most influential ideas developed by Richard Feynman are the Feynman diagram technique and his variational approach. Here we show that combining both, and introducing a diagrammatic quantum Monte Carlo method, 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 it to the quintessential problem of solid state, the uniform electron gas, which is at the heart of the density functional theory success in describing real materials, yet it has not been adequately solved for over 90 years. Our method allows us to calculate numerically exact momentum and frequency resolved spin and charge response functions. This method can be applied to a number of moderately interacting electron systems, including models of realistic metallic and semiconducting solids.
This thesis comprises a trio of separate yet thematically-linked investigations within the field of electronic structure that address challenging problems to which quantum Monte Carlo, the most accurate electronic structure method capable of treating large extended systems, has been applied only rarely or indeed never before. The first is the study of polarization and localization in many body systems. Density functional and variational quantum Monte Carlo methods are used to examine the behaviour of the many-electron localization length near band insulator to metal transitions in various oneand two-dimensional model systems. We examine when the localization length, known to be infinite in metals and finite in insulators and assumed to diverge at a metal-insulator transition, serves as a useful diagnostic of the transition. We also present a novel derivation of the formulae for polarization and localization in many body systems that expresses the connection of these two quantities to the concept of Many Body Wannier Functions. The second is a study of the surface energy of the electron gas. The electron gas, or ‘jellium’, is an important benchmark system for electronic structure methods and serves as a simple model for metals. We apply diffusion Monte Carlo (DMC), and a new method for determining surface energies without an unreliable comparison of bulk and slab results, to the calculation of the surface energy. Our method improves on previous calculations in a number of ways and removes many sources of error that plagued previous calculations. The final study applies DMC to the properties, specifically the formation energy, of point defects in Alumina. Alumina is an important technological material, and there are many important unanswered questions relating to its properties that stem from the behaviour and concentration of its point defects. We address many of the problems associated with the calculation of formation energies with electronic structure calculations, and present the first calculations of formation energies of charged point defects in DMC.
This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.
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...
A new variational approach is proposed for the study of the many-body problem with the purpose of combining the advantages of perturbation-theoretic and variational approaches, while avoiding their disadvantages. A trial wave function is constructed so that a general expansion for the expectation value of an arbitrary physical quantity can be derived. Each term in the expansion can be calculated with the use of the techniques of Feynman diagrams, but in contrast with usual perturbation-theoretic approaches, every term in the expansion always converges. The procedure for applying the method to real systems is illustrated by a calculation of the ground-state energy of the electron gas at high densities.
Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.