We extend the formalism of the thermodynamic two-time Green's functions to nonextensive quantum statistical mechanics. Working in the optimal Lagrangian multiplier representation, the q -spectral properties and the methods for a direct calculation of the two-time q Green's functions and the related q -spectral density ( q measures the nonextensivity degree) for two generic operators are presented in strict analogy with the extensive (q=1) counterpart. Some emphasis is devoted to the nonextensive version of the less known spectral density method whose effectiveness in exploring equilibrium and transport properties of a wide variety of systems has been well established in conventional classical and quantum many-body physics. To check how both the equations of motion and the spectral density methods work to study the q -induced nonextensivity effects in nontrivial many-body problems, we focus on the equilibrium properties of a second-quantized model for a high-density Bose gas with strong attraction between particles for which exact results exist in extensive conditions. Remarkably, the contributions to several thermodynamic quantities of the q -induced nonextensivity close to the extensive regime are explicitly calculated in the low-temperature regime by overcoming the calculation of the q grand-partition function.
We extend the formalism of the thermodynamic two-time Green's functions to nonextensive quantum statistical mechanics. Working in the optimal Lagrangian multiplier representation, the q -spectral properties and the methods for a direct calculation of the two-time q Green's functions and the related q -spectral density ( q measures the nonextensivity degree) for two generic operators are presented in strict analogy with the extensive (q=1) counterpart. Some emphasis is devoted to the nonextensive version of the less known spectral density method whose effectiveness in exploring equilibrium and transport properties of a wide variety of systems has been well established in conventional classical and quantum many-body physics. To check how both the equations of motion and the spectral density methods work to study the q -induced nonextensivity effects in nontrivial many-body problems, we focus on the equilibrium properties of a second-quantized model for a high-density Bose gas with strong attraction between particles for which exact results exist in extensive conditions. Remarkably, the contributions to several thermodynamic quantities of the q -induced nonextensivity close to the extensive regime are explicitly calculated in the low-temperature regime by overcoming the calculation of the q grand-partition function.
The traditional basis of description of many-particle systems in terms of Green functions is here generalized to the case when the system is nonextensive, by incorporating the Tsallis form of the density matrix indexed by a nonextensive parameter q. This is accomplished by expressing the many-particle q Green function in terms of a parametric contour integral over a kernel multiplied by the usual grand canonical Green function which now depends on this parameter. We study one- and two-particle Green functions in detail. From the one-particle Green function, we deduce some experimentally observable quantities such as the one-particle momentum distribution function and the one-particle energy distribution function. Special forms of the two-particle Green functions are related to physical dynamical structure factors, some of which are studied here. We deduce different forms of sum rules in the q formalism. A diagrammatic representation of the q Green functions similar to the traditional ones follows because the equations of motion for both of these are formally similar. Approximation schemes for one-particle q Green functions such as Hartree and Hartree-Fock schemes are given as examples. This extension enables us to predict possible experimental tests for the validity of this framework by expressing some observable quantities in terms of the q averages.
Abstract Green-function techniques are among the most versatile tools for investigating quantum many-body systems at finite temperature. Yet, standard textbook treatments often fragment the formalism—separating bosons and fermions, systems with massless and massive particles, or restricting proofs to special Hamiltonians—thereby obscuring the structural unity of the theory. This work develops a general, Hamiltonian-independent framework for both real- and imaginary-time Green-functions, valid for any combination of bosonic and fermionic degrees of freedom, including mixed systems such as interacting electron–phonon models. The key conceptual advance is the nature index η , which complements the conventional statistical sign ζ . While ζ encodes particle statistics, η distinguishes canonical quantization of systems with massless particles ( η = 0) from second quantization of systems with massive particles ( η = 1). This classification removes ambiguities in chemical-potential terms and unifies the treatment of photons, phonons, electrons, and their mixtures under a single notation. From this basis, we formulate generalized ensemble constructions, generalized Lehmann-representations, and a fluctuation-dissipation theorem valid for arbitrary Hamiltonians. Several results—such as the monotonicity of particle number in the chemical potential, an explicit formula for the quadratic contribution to the internal energy from the retarded Green-function, and the Hartree-contribution to the self-energy in electron–phonon systems—are either absent from, or only sparsely mentioned in, existing literature. We further introduce Green’s Chart, a visual and algebraic map that systematically connects all Green-function types through analytic continuation, Fourier-transforms, and Wick-rotations, enabling all to be generated from the Matsubara Green-function. Applications to both non-interacting and interacting models demonstrate the framework’s ability to compute densities of states, thermodynamic potentials, and spectral functions with a consistent formalism. The results provide a complete and reproducible toolkit for equilibrium many-body theory, consolidating scattered results into a unified structure suitable for both pedagogical and computational purposes.
Schrödinger suggested that thermodynamical functions cannot be based on the gratuitous allegation that quantum-mechanical levels (typically the orthogonal eigenstates of the Hamiltonian operator) are the only allowed states for a quantum system [E. Schrödinger, Statistical Thermodynamics (Courier Dover, Mineola, 1967)]. Different authors have interpreted this statement by introducing density distributions on the space of quantum pure states with weights obtained as functions of the expectation value of the Hamiltonian of the system. In this work we focus on one of the best known of these distributions and prove that, when considered in composite quantum systems, it defines partition functions that do not factorize as products of partition functions of the noninteracting subsystems, even in the thermodynamical regime. This implies that it is not possible to define extensive thermodynamical magnitudes such as the free energy, the internal energy, or the thermodynamic entropy by using these models. Therefore, we conclude that this distribution inspired by Schrödinger's idea cannot be used to construct an appropriate quantum equilibrium thermodynamics.
A brief but clear and complete account is given for two analytical methods of calculating correlation functions in quantum statistical physics from first principles—the widely used method of two-time temperature Green's functions (the GF method) and a new, "direct algebraic" (DA) method. The mathematical and technical clarity and simplicity of the DA method and its resulting practical value are demonstrated for the five most widely used models in quantum statistical physics. Since the DA method is an exactly self-consistent method (in the sense that the expansion coefficients in the equations of motion are chosen from the requirement that the Jacobi operator identity be satisfied exactly), it in principle affords the possibility of an internal check, which is not possible in the GF method. Like the GF method, the DA method permits calculation of the spectra of possible elementary excitations and, hence, of the density of single-particle energy states corresponding to them.
Description of electron correlation is crucial to chemical accuracy in quantum chemical calculations. However, the interaction of electrons in a system is an insoluble many-body problem. Constructing approximations to describe electron correlation is a challenging task, but there has been success in the quantum chemistry community using wavefunction methods and density functional theory. However, there are still challenges to be overcome in areas such as theoretical solid state chemistry, which requires description of large systems and the use of finite temperature. Large system sizes and finite temperature can be difficult to treat solely with the currently available methods. Therefore, a new class of methods, based on the temperature-dependent Green's function, is implemented and explored. This thesis is toward investigating the use of temperature-dependent Green's functions for ab-initio quantum chemistry. That is, we are treating a quantum chemical Hamiltonian with realistic electron interactions. While this formalism has been applied to model systems in the condensed matter community, its has been used much less by the quantum chemistry community. Therefore, the numerical behavior and accuracy of Green's function methods for quantum chemical calculations is relatively unknown. This work investigates the ability of the temperature-dependent Green's function, which is an ensemble formalism, to give access to temperature dependent thermodynamic quantities such as free energy and entropy when calculated in a second-order and perturbative manner (GF2). We find that this method is able to give good accuracy for lower temperatures and excellent accuracy for higher temperatures for a molecular case and is able to qualitatively describe a simple model of a solid. The results of this work are presented in chapter 3. Although Green's functions have a clear connection with spectra at zero temperature, it is not straight forward to obtain spectra from the finite temperature Green's function, which is calculated on the imaginary frequency axis. Therefore, we must investigate methods to obtain spectral quantities in a consistent and reliable manner from the imaginary axis. Chapter 4 investigates several methods to do so and we compare our results with experimental and highly accurate benchmark data. We find that it is possible to obtain spectral quantities that can differ by several electron volts, even if the same level of theory is used to obtain the Green's function. This reiterates that finding a spectrum from the imaginary axis is nontrivial and that one must exercise caution when comparing spectral quantities that were calculated using different techniques, even if they are treated with the same theoretical accuracy. Accessing larger systems with the Green's function requires the use of quantum embedding. Quantum embedding describes nontrivial electron interaction between a highly accurate ``active space'' or ``impurity'' and a larger, lower level ``environment''. It is challenging to construct an impurity solver that is reliable at low temperatures. In chapter 5 of this work, we implement and test a temperature dependent configuration interaction impurity solver for quantum embedding. This solver can be used in quantum embedding schemes such as dynamical mean field theory and self-energy embedding theory for larger systems. Overall, this thesis has made progress toward using Green's functions for ab-initio quantum chemistry at finite temperature. Groundwork has been laid for using this formalism to calculate thermodynamics and spectra using a Green's function with realistic electron interactions and to explore quantum embedding using an impurity solver at low temperatures.
The traditional basis of description of the many-particle systems in terms of the Green functions is here generalized when the system is nonextensive by incorporating the Tsallis form of the density matrix, indexed by a nonextensive parameter, q. This extension enables us to predict possible experimental tests for the validity of this framework by expressing some observable quantities in terms of the q averages.