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

Physics Open

Analytical treatment of the Yukawa screened Coulomb interaction in a plane-wave basis

Priyanka Aggarwal, Ram Kuntal Hazra, Bharti Kapil, Shivalika Sharma, Igor Di Marco

Published Nov 26, 2024• 1 citation
Publisher

Abstract

The calculation of the many-electron (screened) Coulomb and exchange integrals is a very common task to perform in modern electronic structure theory. While an analytical treatment of the complete integrals is too complicated to be attempted directly, essential insight can be gained by focusing on the most significant contributions, as determined by a physical criterion. The monopole term of the screened Coulomb interaction is particularly important, since it defines the Hubbard interaction U among a set of localized electrons, which is an essential parameter for effective models of strongly correlated systems. Here, we derive an analytical solution for the matrix elements of the screened Coulomb interaction on a plane-wave basis, which is routinely used in the most common methods for electronic structure calculations. Screening is treated using a Yukawa potential, which is suitable to describe the valence electrons in the Fermi level region. For the solution of the integrals, the plane waves and the Yukawa potential are first expanded in Bessel functions of first and second kind. Then, by means of the lower and upper incomplete gamma functions, we are able to obtain closed-form integrals of a series that remains convergent for realistic parameters. Our exact solution for the radial integrals of the monopole term can find usage in plane-wave codes for electronic structure calculations, both as an output tool as well as within the computational cycle, as e.g., for many-body extensions of density-functional theory. Considering the importance of Bessel functions in solid state physics and electronic structure theory, it is also easy to foresee that our solution to the various integrals across this work may become useful to several other problems in the field. • Analytical method for matrix elements of the screened Coulomb interaction in plane wave. • Analysis of the monopole term, significant as the Hubbard U in strongly correlated systems. • Closed-form solutions achieved for both long-range and short-range components of integrals. • Method generalizable to evaluate higher-order terms. • Anticipated applicability to plane-wave codes for physical analysis and hybrid functional routines.

DOI
10.1016/j.physo.2024.100243

Related papers

Parametrization of the Coulomb interaction matrix with point-group symmetry

Coraline Letouzé, Guillaume Radtke, Benjamin Lenz +1 more

Coulomb integrals, i.e., matrix elements of bare or screened Coulomb interaction between one-electron orbitals, are fundamental objects in many approaches developed to tackle the challenging problem of calculating the electronic structure of strongly correlated materials. In this paper, Coulomb integrals are analyzed by considering both the point-group symmetry of the site occupied by the atom in the crystal or molecule and the permutation symmetries of the orbitals in the integrals. In particular, the case where one-electron orbitals form the basis of a general (i.e., a real, complex, or pseudoreal) irreducible representation is considered. Explicit formulas are provided to calculate all integrals of the interaction tensor in terms of a minimum set of independent ones. The effect of a symmetry breaking is also investigated by describing Coulomb integrals of a group in terms of those of one of its subgroups. We developed the specific example of $\mathrm{O}(3)$ as the larger group which can therefore be used to quantify the deviation of a specific system from the spherical symmetry. Possible applications of the presented framework include the calculation of solid-state and molecular spectroscopies via multiplet techniques, dynamical mean-field theory, or the $GW$ approximation.

Physical review. B./Physical review. B•Sep 18, 2023•3 citations

An ab-initio many-body method for electronic structure calculations of solids. I. Description of the method

I. Schnell, G. Czycholl, Robert C. Albers

We propose a new, alternative method for ab-initio calculations of the electronic structure of solids, which has been specifically adapted to treat many-body effects in a more rigorous way than many existing ab-initio methods. We start from a standard band-structure calculation for an effective one-particle Hamiltonian approximately describing the material of interest. This yields a suitable set of one-particle basis functions, from which well localized Wannier functions can be constructed using a method proposed by Marzari and Vanderbilt. Within this (minimal) basis of localized Wannier functions the matrix elements of the non-interacting (one-particle) Hamiltonian as well as the Coulomb matrix elements can be calculated. The result is a many-body Hamiltonian in second quantization with parameters determined from first principles calculations for the material of interest. The Hamiltonian is in the form of a multi-band Hamiltonian in second quantization (a kind of extended, multi-band Hubbard model) such that all the standard many-body methods can be applied. We explicitly show how this approach can be solved in the simplest many-body approximation, the mean-field Hartree-Fock approximation (HFA), which takes into account exact exchange and corrects for self-interaction effects.

arXiv:cond-mat/0304643•Apr 28, 2003•0 citations

The Yukawa potential: ground state energy and critical screening

James P. Edwards, Urs Gerber, Christian Schubert +2 more

We study the ground state energy and the critical screening parameter of the Yukawa potential in nonrelativistic quantum mechanics. After a short review of the existing literature on these quantities, we apply fifth-order perturbation theory to the calculation of the ground state energy, using the exact solutions of the Coulomb potential together with a cutoff on the principal number summations. We also perform a variational calculation of the ground state energy using a Coulomb-like radial wave function and the exact solution of the corresponding minimization condition. For not too large values of the screening parameter, close agreement is found between the perturbative and variational results. For the critical screening parameter, we devise a novel method that permits us to determine it to 10 digits. This is the most precise calculation of this quantity to date, and allows us to resolve some discrepancies between previous results.

arXiv:1706.09979•Progress of Theoretical and Experimental Physics•Jul 7, 2017•38 citations

Toward First Principles Electronic Structure Simulations of Excited States and Strong Correlations in Nano- and Materials Science

Anton Kozhevnikov, A. G. Eguiluz, T. C. Schulthess

Methods based on the many-body Green's function are generally accepted as the path forward beyond Kohn-Sham based density functional theory, in order to compute from first principles electronic structure of materials with strong correlations and excited-state properties in nano- and materials science. Here we present an efficient method to compute the screened Coulomb interactionW, the crucial and computationally most demanding ingredient in the GW method, within the framework of the all-electron Linearized Augmented Plane Wave method. We use the method to compute from first principles, within the constrained random phase approximation (c-RPA), the frequency-dependent screened Hubbard U-matrix defined for a Wannier basis in which we downfold the many-body Hamiltonian for La2CuO4, the canonical parent compound of several cuprate high-temperature superconductors. These results were computed at scale on the Cray XT5 at ORNL, sustaining 1.30 petaflop. We discuss the details of the algorithm and its implementation that allowed us to reach high efficiency and short time to solution on today's petaflop supercomputers.

International Conference for High Performance Computing, Networking, Storage and Analysis (SC)•Nov 1, 2010•14 citations

Bound states of the hydrogen-like atomic systems in plasma environments

Fatma Zohra Khaled, Mustafa Moumni, Mokhtar Falek

We conduct a non-relativistic study of plasma screening effects on hydrogen-like atomic systems using the Screened Coulomb Potential (SCP i.e. Yukawa potential). The radial Schrödinger equation is first reduced to a bi-confluent Heun (BCH) equation for the Killingbeck potential which is the truncated version of the SCP, and we write the exact BCH eigenfunctions and eigenenergies. We then study the limitations of these BCH solutions and obtain an analytic description valid for weak to moderate screening using the BCH functional form. The corrected eigenfunctions are constructed order by order up to O(k^3); they are expressed in terms of the Laguerre polynomials and reduces exactly to Coulomb eigenfunctions when k=0. Using these Functions, the energy spectrum is computed via two analytic methods: (i) direct evaluation of the full Yukawa Hamiltonian expectation value, and (ii) the Hellmann-Feynman theorem, yielding integral representation. All methods presented here provide explicit analytical formulas for both wavefunctions and eigenenergies valid for different ranges of the screening. These methods establish a powerful analytic framework for studying confined quantum systems. The thermodynamic properties are also derived using the BCH formulation.

quant-pharXiv:2607.19601•Jul 23, 2026

Many-body perturbation theory: The GW approximation

Christoph Friedrich, Arno Schindlmayr

In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.

•Jan 1, 2006•19 citations
The calculation of the many-electron (screened) Coulomb and exchange integrals is a very common task to perform in modern electronic structure theory. While an analytical treatment of the complete integrals is too complicated to be attempted directly, essential insight can be gained by focusing on the most significant contributions, as determined by a physical criterion. The monopole term of the screened Coulomb interaction is particularly important, since it defines the Hubbard interaction U among a set of localized electrons, which is an essential parameter for effective models of strongly correlated systems. Here, we derive an analytical solution for the matrix elements of the screened Coulomb interaction on a plane-wave basis, which is routinely used in the most common methods for electronic structure calculations. Screening is treated using a Yukawa potential, which is suitable to describe the valence electrons in the Fermi level region. For the solution of the integrals, the plane waves and the Yukawa potential are first expanded in Bessel functions of first and second kind. Then, by means of the lower and upper incomplete gamma functions, we are able to obtain closed-form integrals of a series that remains convergent for realistic parameters. Our exact solution for the radial integrals of the monopole term can find usage in plane-wave codes for electronic structure calculations, both as an output tool as well as within the computational cycle, as e.g., for many-body extensions of density-functional theory. Considering the importance of Bessel functions in solid state physics and electronic structure theory, it is also easy to foresee that our solution to the various integrals across this work may become useful to several other problems in the field. • Analytical method for matrix elements of the screened Coulomb interaction in plane wave. • Analysis of the monopole term, significant as the Hubbard U in strongly correlated systems. • Closed-form solutions achieved for both long-range and short-range components of integrals. • Method generalizable to evaluate higher-order terms. • Anticipated applicability to plane-wave codes for physical analysis and hybrid functional routines.
Coulomb integrals, i.e., matrix elements of bare or screened Coulomb interaction between one-electron orbitals, are fundamental objects in many approaches developed to tackle the challenging problem of calculating the electronic structure of strongly correlated materials. In this paper, Coulomb integrals are analyzed by considering both the point-group symmetry of the site occupied by the atom in the crystal or molecule and the permutation symmetries of the orbitals in the integrals. In particular, the case where one-electron orbitals form the basis of a general (i.e., a real, complex, or pseudoreal) irreducible representation is considered. Explicit formulas are provided to calculate all integrals of the interaction tensor in terms of a minimum set of independent ones. The effect of a symmetry breaking is also investigated by describing Coulomb integrals of a group in terms of those of one of its subgroups. We developed the specific example of O(3)\mathrm{O}(3)O(3) as the larger group which can therefore be used to quantify the deviation of a specific system from the spherical symmetry. Possible applications of the presented framework include the calculation of solid-state and molecular spectroscopies via multiplet techniques, dynamical mean-field theory, or the GWGWGW approximation.
We propose a new, alternative method for ab-initio calculations of the electronic structure of solids, which has been specifically adapted to treat many-body effects in a more rigorous way than many existing ab-initio methods. We start from a standard band-structure calculation for an effective one-particle Hamiltonian approximately describing the material of interest. This yields a suitable set of one-particle basis functions, from which well localized Wannier functions can be constructed using a method proposed by Marzari and Vanderbilt. Within this (minimal) basis of localized Wannier functions the matrix elements of the non-interacting (one-particle) Hamiltonian as well as the Coulomb matrix elements can be calculated. The result is a many-body Hamiltonian in second quantization with parameters determined from first principles calculations for the material of interest. The Hamiltonian is in the form of a multi-band Hamiltonian in second quantization (a kind of extended, multi-band Hubbard model) such that all the standard many-body methods can be applied. We explicitly show how this approach can be solved in the simplest many-body approximation, the mean-field Hartree-Fock approximation (HFA), which takes into account exact exchange and corrects for self-interaction effects.
We study the ground state energy and the critical screening parameter of the Yukawa potential in nonrelativistic quantum mechanics. After a short review of the existing literature on these quantities, we apply fifth-order perturbation theory to the calculation of the ground state energy, using the exact solutions of the Coulomb potential together with a cutoff on the principal number summations. We also perform a variational calculation of the ground state energy using a Coulomb-like radial wave function and the exact solution of the corresponding minimization condition. For not too large values of the screening parameter, close agreement is found between the perturbative and variational results. For the critical screening parameter, we devise a novel method that permits us to determine it to 10 digits. This is the most precise calculation of this quantity to date, and allows us to resolve some discrepancies between previous results.
Methods based on the many-body Green's function are generally accepted as the path forward beyond Kohn-Sham based density functional theory, in order to compute from first principles electronic structure of materials with strong correlations and excited-state properties in nano- and materials science. Here we present an efficient method to compute the screened Coulomb interactionW, the crucial and computationally most demanding ingredient in the GW method, within the framework of the all-electron Linearized Augmented Plane Wave method. We use the method to compute from first principles, within the constrained random phase approximation (c-RPA), the frequency-dependent screened Hubbard U-matrix defined for a Wannier basis in which we downfold the many-body Hamiltonian for La2CuO4, the canonical parent compound of several cuprate high-temperature superconductors. These results were computed at scale on the Cray XT5 at ORNL, sustaining 1.30 petaflop. We discuss the details of the algorithm and its implementation that allowed us to reach high efficiency and short time to solution on today's petaflop supercomputers.
We conduct a non-relativistic study of plasma screening effects on hydrogen-like atomic systems using the Screened Coulomb Potential (SCP i.e. Yukawa potential). The radial Schrödinger equation is first reduced to a bi-confluent Heun (BCH) equation for the Killingbeck potential which is the truncated version of the SCP, and we write the exact BCH eigenfunctions and eigenenergies. We then study the limitations of these BCH solutions and obtain an analytic description valid for weak to moderate screening using the BCH functional form. The corrected eigenfunctions are constructed order by order up to O(k^3); they are expressed in terms of the Laguerre polynomials and reduces exactly to Coulomb eigenfunctions when k=0. Using these Functions, the energy spectrum is computed via two analytic methods: (i) direct evaluation of the full Yukawa Hamiltonian expectation value, and (ii) the Hellmann-Feynman theorem, yielding integral representation. All methods presented here provide explicit analytical formulas for both wavefunctions and eigenenergies valid for different ranges of the screening. These methods establish a powerful analytic framework for studying confined quantum systems. The thermodynamic properties are also derived using the BCH formulation.
In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.