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:cond-mat/9802051

Perturbation Expansion of the Partition Sum for any Temperature

Rolf Schumann

Published Feb 4, 1998• 0 citations
arXivHTMLPDF

Abstract

Based on the special properties of Liouville eigenoperators a perturbation theory for the partition sum is given. It is applicable for any temperature and includes the case of degenerate Hamiltonians. To demonstrate the reliability of the method, the second order correction to the atomic limit grand canonical potential of the Hubbard model is calculated and compared to results known from the literature.

DOI
10.48550/arxiv.cond-mat/9802051

Related papers

Finite-temperature many-body perturbation theory in the grand canonical ensemble

So Hirata, Punit K. Jha

A finite-temperature many-body perturbation theory is presented, which expands in power series the electronic grand potential, chemical potential, internal energy, and entropy on an equal footing. Sum-over-states and sum-over-orbitals analytical formulas for the second-order perturbation corrections to these thermodynamic properties are obtained in a time-independent, nondiagrammatic, algebraic derivation, relying on the sum rules of the Hirschfelder–Certain degenerate perturbation energies in a degenerate subspace as well as nine algebraic identities for the zeroth-order thermal averages of one- through four-indexed quantities and products thereof. They reproduce numerically exactly the benchmark data obtained as the numerical derivatives of the thermal-full-configuration-interaction results for a wide range of temperatures.

arXiv:2003.05540•The Journal of Chemical Physics•Jul 1, 2020•24 citations

Finite-temperature many-body perturbation theory for electrons: Algebraic recursive definitions, second-quantized derivation, linked-diagram theorem, general-order algorithms, and grand canonical and canonical ensembles

So Hirata

A comprehensive and detailed account is presented for the finite-temperature many-body perturbation theory for electrons that expands in power series all thermodynamic functions on an equal footing. Algebraic recursions in the style of the Rayleigh-Schrödinger perturbation theory are derived for the grand potential, chemical potential, internal energy, and entropy in the grand canonical ensemble and for the Helmholtz energy, internal energy, and entropy in the canonical ensemble, leading to their sum-over-states analytical formulas at any arbitrary order. For the grand canonical ensemble, these sum-over-states formulas are systematically transformed to sum-over-orbitals reduced analytical formulas by the quantum-field-theoretical techniques of normal-ordered second quantization and Feynman diagrams extended to finite temperature. It is found that the perturbation corrections to energies entering the recursions have to be treated as a nondiagonal matrix, whose off-diagonal elements are generally nonzero within a subspace spanned by degenerate Slater determinants. They give rise to a unique set of linked diagrams-renormalization diagrams-whose resolvent lines are displaced upward, which are distinct from the well-known anomalous diagrams of which one or more resolvent lines are erased. A linked-diagram theorem is introduced that proves the size-consistency of the finite-temperature many-body perturbation theory at any order. General-order algorithms implementing the recursions establish the convergence of the perturbation series toward the finite-temperature full-configuration-interaction limit unless the series diverges. The normal-ordered Hamiltonian at finite temperature sheds light on the relationship between the finite-temperature Hartree-Fock and first-order many-body perturbation theories.

arXiv:2106.10747•The Journal of Chemical Physics•Sep 7, 2021•17 citations

Finite-temperature many-body perturbation theory in the canonical ensemble

Punit K. Jha, So Hirata

Benchmark data are presented for the zeroth- through third-order many-body perturbation corrections to the electronic Helmholtz energy, internal energy, and entropy in the canonical ensemble in a wide range of temperature. They are determined as numerical λ-derivatives of the respective quantities computed by thermal full configuration interaction with a perturbation-scaled Hamiltonian, H[over ̂]=H[over ̂]_{0}+λV[over ̂]. Sum-over-states analytical formulas for up to the third-order corrections to these properties are also derived as analytical λ-derivatives. These formulas, which are verified by exact numerical agreement with the benchmark data, are given in terms of the Hirschfelder-Certain degenerate perturbation energies and should be valid for both degenerate and nondegenerate reference states at any temperature down to zero. The results in the canonical ensemble are compared with the same in the grand canonical ensemble.

arXiv:1910.07628•Physical review. E•Feb 6, 2020•25 citations

Low-temperature breakdown of many-body perturbation theory for thermodynamics

So Hirata

It is shown analytically and numerically that the finite-temperature many-body perturbation theory in the grand canonical ensemble has zero radius of convergence at zero temperature when the energy ordering or degree of degeneracy for the ground state changes with the perturbation strength. When the degeneracy of the reference state is partially or fully lifted at the first-order Hirschfelder-Certain degenerate perturbation theory, the second-order grand potential and internal energy diverge as $T\ensuremath{\rightarrow}0$. Contrary to earlier suggestions of renormalizability by the chemical potential $\ensuremath{\mu}$, this nonconvergence, first suspected by Kohn and Luttinger, is caused by the nonanalytic nature of the Boltzmann factor ${e}^{\ensuremath{-}E/{k}_{\text{B}}T}$ at $T=0$, also plaguing the canonical ensemble, which does not involve $\ensuremath{\mu}$. The finding reveals a fundamental flaw in perturbation theory, which is deeply rooted in the mathematical limitation of power-series expansions and is unlikely to be removed within its framework.

arXiv:2006.00078•Physical Review A•Jan 26, 2021•10 citations

Analytical results of the one-dimensional Hubbard model in the high-temperature limit

I. C. Charret, E. V. Corrêa Silva, S. M. de Souza +3 more

We investigate the grand potential of the one-dimensional Hubbard model in the high-temperature limit, calculating the coefficients of the high-temperature expansion $(\ensuremath{\beta}$-expansion) of this function up to order ${\ensuremath{\beta}}^{4}$ by an alternative method. The results derived are analytical and do not involve any perturbation expansion in the hopping constant, being valid for arbitrary density of electrons in the one-dimensional model. In the half-filled case, we compare our analytical results for the specific heat and the magnetic susceptibility, in the high-temperature limit, with the ones obtained by Beni et al. [Phys. Rev. B 8, 3329 (1973)] and Takahashi's integral equations, showing that the latter result does not take into account the complete energy spectrum of the one-dimensional Hubbard model. The exact integral solution by J\"uttner et al. [Nucl. Phys. B 522, 471 (1998)] is applied to the determination of the range of validity of our expansion in $\ensuremath{\beta}$ in the half-filled case, for several different values of U.

arXiv:cond-mat/9910043•Physical review. B, Condensed matter•Oct 29, 2001•6 citations

Linked-cluster expansion around the atomic limit of the Hubbard model

Walter Metzner

We develop a perturbation expansion in the intersite hopping around the atomic limit of the Hubbard model. It is valid for arbitrary finite temperatures and interaction strengths. Diagrammatic rules that determine the grand-canonical potential and the Green's functions are derived. They reduce the calculation of any finite-order contribution to simple algebra. This opens the way for series extrapolations from computer-aided high-finite-order evaluations. Discrepancies in earlier expansions around the atomic limit are clarified. The present expansion scheme involves only connected diagrams with unrestricted lattice sums. This allows one to perform a vertex renormalization as for the linked-cluster expansion of the Ising model. The renormalized perturbation expansion can be used to construct self-consistent approximations which are automatically exact in the atomic limit. In the limit of high lattice dimensions, only fully two-particle reducible embeddings of diagrams on the lattice contribute. The single-particle properties of the infinite-dimensional Hubbard model reduce to those of independent tight-binding fermions hopping between dressed sites.

Physical review. B, Condensed matter•Apr 1, 1991•239 citations
Based on the special properties of Liouville eigenoperators a perturbation theory for the partition sum is given. It is applicable for any temperature and includes the case of degenerate Hamiltonians. To demonstrate the reliability of the method, the second order correction to the atomic limit grand canonical potential of the Hubbard model is calculated and compared to results known from the literature.
A finite-temperature many-body perturbation theory is presented, which expands in power series the electronic grand potential, chemical potential, internal energy, and entropy on an equal footing. Sum-over-states and sum-over-orbitals analytical formulas for the second-order perturbation corrections to these thermodynamic properties are obtained in a time-independent, nondiagrammatic, algebraic derivation, relying on the sum rules of the Hirschfelder–Certain degenerate perturbation energies in a degenerate subspace as well as nine algebraic identities for the zeroth-order thermal averages of one- through four-indexed quantities and products thereof. They reproduce numerically exactly the benchmark data obtained as the numerical derivatives of the thermal-full-configuration-interaction results for a wide range of temperatures.
A comprehensive and detailed account is presented for the finite-temperature many-body perturbation theory for electrons that expands in power series all thermodynamic functions on an equal footing. Algebraic recursions in the style of the Rayleigh-Schrödinger perturbation theory are derived for the grand potential, chemical potential, internal energy, and entropy in the grand canonical ensemble and for the Helmholtz energy, internal energy, and entropy in the canonical ensemble, leading to their sum-over-states analytical formulas at any arbitrary order. For the grand canonical ensemble, these sum-over-states formulas are systematically transformed to sum-over-orbitals reduced analytical formulas by the quantum-field-theoretical techniques of normal-ordered second quantization and Feynman diagrams extended to finite temperature. It is found that the perturbation corrections to energies entering the recursions have to be treated as a nondiagonal matrix, whose off-diagonal elements are generally nonzero within a subspace spanned by degenerate Slater determinants. They give rise to a unique set of linked diagrams-renormalization diagrams-whose resolvent lines are displaced upward, which are distinct from the well-known anomalous diagrams of which one or more resolvent lines are erased. A linked-diagram theorem is introduced that proves the size-consistency of the finite-temperature many-body perturbation theory at any order. General-order algorithms implementing the recursions establish the convergence of the perturbation series toward the finite-temperature full-configuration-interaction limit unless the series diverges. The normal-ordered Hamiltonian at finite temperature sheds light on the relationship between the finite-temperature Hartree-Fock and first-order many-body perturbation theories.
Benchmark data are presented for the zeroth- through third-order many-body perturbation corrections to the electronic Helmholtz energy, internal energy, and entropy in the canonical ensemble in a wide range of temperature. They are determined as numerical λ-derivatives of the respective quantities computed by thermal full configuration interaction with a perturbation-scaled Hamiltonian, H[over ̂]=H[over ̂]_{0}+λV[over ̂]. Sum-over-states analytical formulas for up to the third-order corrections to these properties are also derived as analytical λ-derivatives. These formulas, which are verified by exact numerical agreement with the benchmark data, are given in terms of the Hirschfelder-Certain degenerate perturbation energies and should be valid for both degenerate and nondegenerate reference states at any temperature down to zero. The results in the canonical ensemble are compared with the same in the grand canonical ensemble.
It is shown analytically and numerically that the finite-temperature many-body perturbation theory in the grand canonical ensemble has zero radius of convergence at zero temperature when the energy ordering or degree of degeneracy for the ground state changes with the perturbation strength. When the degeneracy of the reference state is partially or fully lifted at the first-order Hirschfelder-Certain degenerate perturbation theory, the second-order grand potential and internal energy diverge as $T\ensuremath{\rightarrow}0$. Contrary to earlier suggestions of renormalizability by the chemical potential $\ensuremath{\mu}$, this nonconvergence, first suspected by Kohn and Luttinger, is caused by the nonanalytic nature of the Boltzmann factor ${e}^{\ensuremath{-}E/{k}_{\text{B}}T}$ at T=0T=0T=0, also plaguing the canonical ensemble, which does not involve $\ensuremath{\mu}$. The finding reveals a fundamental flaw in perturbation theory, which is deeply rooted in the mathematical limitation of power-series expansions and is unlikely to be removed within its framework.
We investigate the grand potential of the one-dimensional Hubbard model in the high-temperature limit, calculating the coefficients of the high-temperature expansion $(\ensuremath{\beta}$-expansion) of this function up to order ${\ensuremath{\beta}}^{4}$ by an alternative method. The results derived are analytical and do not involve any perturbation expansion in the hopping constant, being valid for arbitrary density of electrons in the one-dimensional model. In the half-filled case, we compare our analytical results for the specific heat and the magnetic susceptibility, in the high-temperature limit, with the ones obtained by Beni et al. [Phys. Rev. B 8, 3329 (1973)] and Takahashi's integral equations, showing that the latter result does not take into account the complete energy spectrum of the one-dimensional Hubbard model. The exact integral solution by J\"uttner et al. [Nucl. Phys. B 522, 471 (1998)] is applied to the determination of the range of validity of our expansion in $\ensuremath{\beta}$ in the half-filled case, for several different values of U.
We develop a perturbation expansion in the intersite hopping around the atomic limit of the Hubbard model. It is valid for arbitrary finite temperatures and interaction strengths. Diagrammatic rules that determine the grand-canonical potential and the Green's functions are derived. They reduce the calculation of any finite-order contribution to simple algebra. This opens the way for series extrapolations from computer-aided high-finite-order evaluations. Discrepancies in earlier expansions around the atomic limit are clarified. The present expansion scheme involves only connected diagrams with unrestricted lattice sums. This allows one to perform a vertex renormalization as for the linked-cluster expansion of the Ising model. The renormalized perturbation expansion can be used to construct self-consistent approximations which are automatically exact in the atomic limit. In the limit of high lattice dimensions, only fully two-particle reducible embeddings of diagrams on the lattice contribute. The single-particle properties of the infinite-dimensional Hubbard model reduce to those of independent tight-binding fermions hopping between dressed sites.