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

Study of linear response in Hubbard chains using Many-body Perturbation Theory

Irina Bostan

Published Jan 1, 2010• 0 citations
OpenAlex

Abstract

In this work the basic formalism of non-equilibrium Green’s functions is presented and then applied to study a Ward identity in linear response theory, namely the frequency sum-rule. It can be proven that the frequency sum-rule is satisfied when the quantities involved are calculated using perturbation theory within a conserving approximation for the self-energy. To illustrate this equality along with other properties of the response function, a numerical application that solves the Kadanoff-Baym equations for systems of Hubbard chains was used. The results showed that the frequency sum-rule was satisfied to the same extent by all the conserving approximations used as by the exact diagonalization numerical results. The density response function was analyzed diagrammatically for a series of conserving approximations for the self-energy and this demonstrated that even for a first order in perturbation theory approximation for the self-energy, the response function has a corresponding complex, third order in the perturbation diagrammatic structure.

Related papers

Exact self-energy of the many-body problem from conserving approximations

Yasutami Takada

A procedure is proposed to obtain the exact self-energy in the many-body problem. This algorithm is based on the formal analysis to reach the exact theory by repeated applications of an operator F to an arbitrarily chosen input self-energy represented as a functional of the dressed Green's function. The operator F is so defined that the microscopic conservation law for particle number is satisfied. The rigorous self-energy is obtained by the solution of an eigenfunction of F. Particular attention is paid to the relation between the present procedure and the Baym-Kadanoff framework of conserving approximations. By simplifying the procedure in F with use of the generalized Ward identity, we suggest a practical method to implement this algorithm rather easily in actual systems. In order to suggest future directions to improve on this practical method, the recently developed mean-field theory for the Hubbard model in the limit of high spatial dimensions is also discussed in the context of our theory.

Physical review. B, Condensed matter•Nov 1, 1995•24 citations

Artificial damping in the Kadanoff-Baym dynamics of small Hubbard chains

M. Puig von Friesen, Claudio Verdozzi, C.‐O. Almbladh

We perform a comparative study of exact and approximate time-evolved densities in small Hubbard chains. The approximate densities are obtained via many-body perturbation theory (Hartree-Fock, 2nd Born, GW and T-matrix approximations) within the framework of the time-dependent Kadanoff-Baym equations. Benchmarking approximate results against exact ones allows us to address two rather fundamental issues in the non equilibrium dynamics of strongly correlated systems. I) A characterisation of the performance of several standard MBAs in the non-equilibrium regime. Having a definite notion of how good a specific MBA can be is highly relevant to its application to cases (typically, infinite systems) where exact solutions are not available. Our results show that the T-matrix approximation is overall superior to the other MBAs, at all electron densities. II) A scrutiny of the whole idea of Many Body Perturbation Theory in the Kadanoff-Baym sense, when applied to finite systems. The surprising outcome of our study is that during the time evolution, the KBE develop an unphysical steady state solution. This is a genuinely novel feature of the time-dependent KBE, i.e. is not inherited from possible limitations/approximations in the calculation of the initial state. Our extensive numerical characterisation gives robust evidence that the problem occurs in general, whenever MBPT is applied to finite systems, and approximate self energies based upon infinite partial summations are used. We also offer some more conceptual and general consideration on the dependence of this behaviour on the number of particles and system size. This is followed by our conclusions and glimpses of future work.

Journal of Physics Conference Series•Apr 1, 2010•7 citations

Dynamical properties of theD=∞ Hubbard model

Dai S. Hirashima

Dynamical response functions of the Hubbard model on an infinite dimensional lattice (D=\ensuremath{\infty}) are investigated by the self-consistent second-order perturbation theory with respect to the Coulomb interaction U. The vertex corrections are properly considered so that the conservation laws are satisfied, i.e., the f-sum rule is satisfied. The explicit expression for the f-sum rule is obtained in the case of infinite dimensions. It is found that even in the second-order perturbation theory the electronic states show many-body features at finite temperature and/or frequency regions: the system crosses over from the Fermi-liquid state at the low temperature and low-frequency region to the incoherent regime at high temperatures.

Physical review. B, Condensed matter•Jun 15, 1993•8 citations

Partial self-consistency and analyticity in many-body perturbation theory: Particle number conservation and a generalized sum rule

Daniel Karlsson, Robert van Leeuwen

We consider a general class of approximations which guarantees the conservation of particle number in many-body perturbation theory. To do this we extend the concept of $\mathrm{\ensuremath{\Phi}}$ derivability for the self-energy $\mathrm{\ensuremath{\Sigma}}$ to a larger class of diagrammatic terms in which only some of the Green's function lines contain the fully dressed Green's function $G$. We call the corresponding approximations for $\mathrm{\ensuremath{\Sigma}}$ partially $\mathrm{\ensuremath{\Phi}}$ derivable. A special subclass of such approximations, which are gauge invariant, is obtained by dressing loops in the diagrammatic expansion of $\mathrm{\ensuremath{\Phi}}$ consistently with $G$. These approximations are number conserving but do not have to fulfill other conservation laws, such as the conservation of energy and momentum. From our formalism we can easily deduce whether commonly used approximations will fulfill the continuity equation, which implies particle number conservation. We further show how the concept of partial $\mathrm{\ensuremath{\Phi}}$ derivability plays an important role in the derivation of a generalized sum rule for the particle number, which reduces to the Luttinger-Ward theorem in the case of a homogeneous electron gas, and the Friedel sum rule in the case of the Anderson model. To do this we need to ensure that the Green's function has certain complex analytic properties, which can be guaranteed if the spectral function is positive-semidefinite. The latter property can be ensured for a subset of partially $\mathrm{\ensuremath{\Phi}}$-derivable approximations for the self-energy, namely those that can be constructed from squares of so-called half diagrams. For the case in which the analytic requirements are not fulfilled we highlight a number of subtle issues related to branch cuts, pole structure, and multivaluedness. We also show that various schemes of computing the particle number are consistent for particle number conserving approximations.

arXiv:1606.07486•Physical review. B./Physical review. B•Sep 13, 2016•14 citations

Kadanoff-Baym equations and approximate double occupancy in a Hubbard dimer

M. Puig von Friesen, Claudio Verdozzi, C.‐O. Almbladh

We use the Kadanoff-Baym equations within the framework of many-body perturbation theory to study the double occupancies in a Hubbard dimer. The double occupancies are obtained from the equations of motion of the single-particle Green's function. Our calculations show that the approximate double occupancies can become negative. This is a shortcoming of approximate, and yet conserving, many-body self-energies. Among the tested perturbation schemes, the T-matrix approximation is the only one providing double occupancies which are always positive.

arXiv:1009.2917•Sep 15, 2010•4 citations

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

Niko Säkkinen, Yang Peng, Heiko Appel +1 more

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

arXiv:1507.04726•The Journal of Chemical Physics•Dec 15, 2015•36 citations
In this work the basic formalism of non-equilibrium Green’s functions is presented and then applied to study a Ward identity in linear response theory, namely the frequency sum-rule. It can be proven that the frequency sum-rule is satisfied when the quantities involved are calculated using perturbation theory within a conserving approximation for the self-energy. To illustrate this equality along with other properties of the response function, a numerical application that solves the Kadanoff-Baym equations for systems of Hubbard chains was used. The results showed that the frequency sum-rule was satisfied to the same extent by all the conserving approximations used as by the exact diagonalization numerical results. The density response function was analyzed diagrammatically for a series of conserving approximations for the self-energy and this demonstrated that even for a first order in perturbation theory approximation for the self-energy, the response function has a corresponding complex, third order in the perturbation diagrammatic structure.
A procedure is proposed to obtain the exact self-energy in the many-body problem. This algorithm is based on the formal analysis to reach the exact theory by repeated applications of an operator F to an arbitrarily chosen input self-energy represented as a functional of the dressed Green's function. The operator F is so defined that the microscopic conservation law for particle number is satisfied. The rigorous self-energy is obtained by the solution of an eigenfunction of F. Particular attention is paid to the relation between the present procedure and the Baym-Kadanoff framework of conserving approximations. By simplifying the procedure in F with use of the generalized Ward identity, we suggest a practical method to implement this algorithm rather easily in actual systems. In order to suggest future directions to improve on this practical method, the recently developed mean-field theory for the Hubbard model in the limit of high spatial dimensions is also discussed in the context of our theory.
We perform a comparative study of exact and approximate time-evolved densities in small Hubbard chains. The approximate densities are obtained via many-body perturbation theory (Hartree-Fock, 2nd Born, GW and T-matrix approximations) within the framework of the time-dependent Kadanoff-Baym equations. Benchmarking approximate results against exact ones allows us to address two rather fundamental issues in the non equilibrium dynamics of strongly correlated systems. I) A characterisation of the performance of several standard MBAs in the non-equilibrium regime. Having a definite notion of how good a specific MBA can be is highly relevant to its application to cases (typically, infinite systems) where exact solutions are not available. Our results show that the T-matrix approximation is overall superior to the other MBAs, at all electron densities. II) A scrutiny of the whole idea of Many Body Perturbation Theory in the Kadanoff-Baym sense, when applied to finite systems. The surprising outcome of our study is that during the time evolution, the KBE develop an unphysical steady state solution. This is a genuinely novel feature of the time-dependent KBE, i.e. is not inherited from possible limitations/approximations in the calculation of the initial state. Our extensive numerical characterisation gives robust evidence that the problem occurs in general, whenever MBPT is applied to finite systems, and approximate self energies based upon infinite partial summations are used. We also offer some more conceptual and general consideration on the dependence of this behaviour on the number of particles and system size. This is followed by our conclusions and glimpses of future work.
Dynamical response functions of the Hubbard model on an infinite dimensional lattice (D=\ensuremath{\infty}) are investigated by the self-consistent second-order perturbation theory with respect to the Coulomb interaction U. The vertex corrections are properly considered so that the conservation laws are satisfied, i.e., the f-sum rule is satisfied. The explicit expression for the f-sum rule is obtained in the case of infinite dimensions. It is found that even in the second-order perturbation theory the electronic states show many-body features at finite temperature and/or frequency regions: the system crosses over from the Fermi-liquid state at the low temperature and low-frequency region to the incoherent regime at high temperatures.
We consider a general class of approximations which guarantees the conservation of particle number in many-body perturbation theory. To do this we extend the concept of $\mathrm{\ensuremath{\Phi}}$ derivability for the self-energy $\mathrm{\ensuremath{\Sigma}}$ to a larger class of diagrammatic terms in which only some of the Green's function lines contain the fully dressed Green's function GGG. We call the corresponding approximations for $\mathrm{\ensuremath{\Sigma}}$ partially $\mathrm{\ensuremath{\Phi}}$ derivable. A special subclass of such approximations, which are gauge invariant, is obtained by dressing loops in the diagrammatic expansion of $\mathrm{\ensuremath{\Phi}}$ consistently with GGG. These approximations are number conserving but do not have to fulfill other conservation laws, such as the conservation of energy and momentum. From our formalism we can easily deduce whether commonly used approximations will fulfill the continuity equation, which implies particle number conservation. We further show how the concept of partial $\mathrm{\ensuremath{\Phi}}$ derivability plays an important role in the derivation of a generalized sum rule for the particle number, which reduces to the Luttinger-Ward theorem in the case of a homogeneous electron gas, and the Friedel sum rule in the case of the Anderson model. To do this we need to ensure that the Green's function has certain complex analytic properties, which can be guaranteed if the spectral function is positive-semidefinite. The latter property can be ensured for a subset of partially $\mathrm{\ensuremath{\Phi}}$-derivable approximations for the self-energy, namely those that can be constructed from squares of so-called half diagrams. For the case in which the analytic requirements are not fulfilled we highlight a number of subtle issues related to branch cuts, pole structure, and multivaluedness. We also show that various schemes of computing the particle number are consistent for particle number conserving approximations.
We use the Kadanoff-Baym equations within the framework of many-body perturbation theory to study the double occupancies in a Hubbard dimer. The double occupancies are obtained from the equations of motion of the single-particle Green's function. Our calculations show that the approximate double occupancies can become negative. This is a shortcoming of approximate, and yet conserving, many-body self-energies. Among the tested perturbation schemes, the T-matrix approximation is the only one providing double occupancies which are always positive.
We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.