Finite-temperature Many-body Perturbation Theory Research Articles

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm.

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose-Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh-Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born-Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick's theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn-Luttinger-type nonconvergence.

On the Chemical Potential of Many-Body Perturbation Theory in Extended Systems.

Finite-temperature many-body perturbation theory in the grand-canonical ensemble is fundamental to numerous methods for computing electronic properties at nonzero temperature, such as finite-temperature coupled-cluster. In most applications it is the average number of electrons that is known rather than the chemical potential. Expensive correlation calculations must be repeated iteratively in search for the interacting chemical potential that yields the desired average number of electrons. In extended systems with mobile charges the situation is particular, however. Long-ranged electrostatic forces drive the charges such that the average ratio of negative and positive charges is one for any finite chemical potential. All properties per electron are expected to be virtually independent of the chemical potential, as they are in an electric wire at different voltage potentials. This work shows that per electron, the exchange-correlation free energy and the exchange-correlation grand potential indeed agree in the infinite-size limit. Thus, only one expensive correlation calculation suffices for each system size, sparing the search for the interacting chemical potential. This work also demonstrates the importance of regularizing the Coulomb interaction such that each electron on average interacts only with as many electrons as there are electrons in the simulation, avoiding interactions with periodic images. Numerical calculations of the warm uniform electron gas have been conducted with the Spencer-Alavi regularization employing the finite-temperature Hartree approximation for the self-consistent field and linearized finite-temperature direct-ring coupled-cluster doubles for treating correlation.

General solution to the Kohn–Luttinger nonconvergence problem

A simple, but general solution is proposed for the Kohn–Luttinger problem, i.e., the nonconvergence of the finite-temperature many-body perturbation theory with its zero-temperature counterpart as temperature is lowered to zero under some circumstances. How this nonconvergence can be avoided by altering the reference wave function is illustrated numerically by using up to the fifth order of the perturbation theory.

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. We compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

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.

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.

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

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 grand potential and internal energy diverge as $T \to 0$. Contrary to earlier suggestions of renormalizability by the chemical potential $\mu$, this nonconvergence, first suspected by W. Kohn and J. M. Luttinger, is caused by the nonanalytic nature of the Boltzmann factor $e^{-E/k_\text{B}T}$ at $T=0$, also plaguing the canonical ensemble, which does not involve $\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.

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

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.

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

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.

Automatic computing of the grand potential in finite temperature many-body perturbation theory

A new program created in C/C[Formula: see text] language generates automatically the analytic expression of grand potential and prints it in Latex2e format and in textual data. Another code created in Mathematica language can translate the textual data into a mathematical expression and help any interested to evaluate the thermodynamic quantities in analytic or numeric forms.

On the Kohn–Luttinger conundrum

Kohn and Luttinger [Phys. Rev. 118, 41 (1960)] showed that the conventional finite-temperature extension of the second-order many-body perturbation theory had the incorrect zero-temperature limit in metals and, on this basis, argued that the theory was incorrect. We show that this inconsistency arises from the noninclusion of the temperature effect in the energies of the zeroth-order eigenstates of the perturbation theory, which causes not only the Kohn-Luttinger conundrum but also another inconsistency with the zero-temperature many-body perturbation theory, namely, the different rates of divergence of the correlation energy in a homogeneous electron gas (HEG). We propose a renormalized many-body perturbation theory derivable from the finite-temperature extension of the normal-ordered second quantization applied to the denominators of the energy expression, which involves the energies of the zeroth-order states, as well as to the numerators. The renormalized theory is shown to have the correct zero-temperature limit and the same rate of divergence in a HEG as the zero-temperature counterpart, and is, therefore, the correct finite-temperature many-body perturbation theory.

Strong-coupling expansions at finite temperatures: Application to quantum disordered and quantum critical phases

By combining conventional finite-temperature many-body perturbation theory with cluster expansions, we develop a systematic method to carry out high order arbitrary temperature perturbative calculations on a computer. The method is well suited to studying the thermodynamic properties of quantum disordered and quantum critical phases at finite temperatures. As an application, we calculate the magnetic susceptibility, internal energy, and specific heat of the bilayer Heisenberg model. It is shown that for a wide range of coupling constants these expansions show excellent convergence at all temperatures. Comparing the direct series (without extrapolations) for the bulk susceptibility to quantum Monte Carlo simulations we find an almost perfect agreement between the two methods even at the quantum critical coupling separating the dimerized and antiferromagnetic phases. The convergence fails only at very low temperatures, which are also difficult to reach by quantum Monte Carlo simulations.