arXiv:2402.13030 · The European Physical Journal Plus
Compressing the two-particle Green’s function using wavelets: Theory and application to the Hubbard atom
Emin Moghadas, Nikolaus Dräger, A. Toschi, Jiawei Zang, Matija Medvidović, Dominik Kiese, Andrew J. Millis, Anirvan M. Sengupta, Sabine Andergassen, Domenico Di Sante
Abstract Precise algorithms capable of providing controlled solutions in the presence of strong interactions are transforming the landscape of quantum many-body physics. Particularly, exciting breakthroughs are enabling the computation of non-zero temperature correlation functions. However, computational challenges arise due to constraints in resources and memory limitations, especially in scenarios involving complex Green’s functions and lattice effects. Leveraging the principles of signal processing and data compression, this paper explores the wavelet decomposition as a versatile and efficient method for obtaining compact and resource-efficient representations of the many-body theory of interacting systems. The effectiveness of the wavelet decomposition is illustrated through its application to the representation of generalized susceptibilities and self-energies in a prototypical interacting fermionic system, namely the Hubbard model at half-filling in its atomic limit. These results are the first proof-of-principle application of the wavelet compression within the realm of many-body physics and demonstrate the potential of this wavelet-based compression scheme for understanding the physics of correlated electron systems.
Abstract Precise algorithms capable of providing controlled solutions in the presence of strong interactions are transforming the landscape of quantum many-body physics. Particularly, exciting breakthroughs are enabling the computation of non-zero temperature correlation functions. However, computational challenges arise due to constraints in resources and memory limitations, especially in scenarios involving complex Green’s functions and lattice effects. Leveraging the principles of signal processing and data compression, this paper explores the wavelet decomposition as a versatile and efficient method for obtaining compact and resource-efficient representations of the many-body theory of interacting systems. The effectiveness of the wavelet decomposition is illustrated through its application to the representation of generalized susceptibilities and self-energies in a prototypical interacting fermionic system, namely the Hubbard model at half-filling in its atomic limit. These results are the first proof-of-principle application of the wavelet compression within the realm of many-body physics and demonstrate the potential of this wavelet-based compression scheme for understanding the physics of correlated electron systems.
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.
The Hubbard model is investigated in the framework of lattice density functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable of the many-body problem. A new approximation to the interaction-energy functional $W[\ensuremath{\gamma}]$ is proposed which is based on its scaling properties and which recovers exactly the limit of strong electron correlations at half-band filling. In this way, a more accurate description of W is obtained throughout the domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}},$ including the crossover from weak to strong correlations. As examples of applications results are given for the ground-state energy, charge-excitation gap, and charge susceptibility of the Hubbard model in one-, two-, and three-dimensional lattices. The performance of the method is demonstrated by comparison with available exact solutions, with numerical calculations, and with LDFT using a simpler dimer ansatz for W. Goals and limitations of the different approximations are discussed.
In this work we investigate some aspects of the physics of strongly correlated systems by taking into account both electron-electron and electron-phonon interactions as basic mechanisms for reproducing electronic correlations in real materials. The relevance of the electron-electron interactions is discussed in the first part of this thesis in the framework of a self-consistent theoretical approach, named Composite Operator Method (COM), which accounts for the relevant quasi-particle excitations in terms of a set of composite operators that appear as a result of the modification imposed by the interactions on the canonical electronic fields. We show that the COM allows the calculation of all the relevant Green’s and correlation functions in terms of a number of unknown internal parameters to be determined self-consistently. Therefore, depending on the balance between unknown parameters and self-consistent equations, exact and approximate solutions can be obtained. By way of example, we discuss the application of the COM to the extended t-U- J-h model in the atomic limit, and to the two-dimensional single-band Hubbard model. In the former case, we show that the COM provides the exact solution of the model in one dimension. We study the effects of electronic correlations as responsible for the formation of a plethora of different charge and/or spin orderings. We report the phase diagram of the model, as well as a detailed analysis of both zero and finite temperature single-particle and thermodynamic properties. As far as the single-band Hubbard model is concerned, we illustrate an approximated selfconsistent scheme based on the choice of a two-field basis. We report a detailed analysis of many unconventional features that arise in single-particle properties, thermodynamics and system’s response functions. We emphasize that the accuracy of the COM in describing the effects of electronic correlations strongly relies on the choice of the basis, paving the way for possible multi-pole extensions to the twofield theory. To this purpose, we also study a three-field approach to the single-band Hubbard model, showing a significant step forward in the agreements with numerical data with respect to the two-pole results. The role of the electron-phonon interaction in the physics of strongly correlated systems is discussed in the second part of this thesis. We show that in highly polarizable lattices the competition between unscreened Coulomb and Fröhlich interactions results in a short-range polaronic exchange term Jp that favours the formation of local and light pairs of bosonic nature, named bipolarons, which condense with a critical temperature well in excess of hundred kelvins. These findings, discussed in the framework of the so-called polaronic t-Jp model, are further investigated in the presence of a finite on-site potential ~U , coming from the competition between on-site Coulomb and Fröhlich interactions. We discuss the role of ~U as the driving parameter for a small-to-large bipolaron transition, providing a possible explanation of the BEC-BCS crossover in terms of the properties of the bipolaronic ground state. Finally, we show that a hard-core bipolarons gas, studied as a charged Bose-Fermi mixture, allows for the description of many non Fermi liquid behaviours, allowing also for a microscopic explanation of pseudogap features in terms of a thermal-induced recombination of polarons and bipolarons, without any assumption on preexisting order or broken symmetries. [edited by Author]
One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales -- for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections -- but strongly correlated electrons require a more complex description and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We derive nonequilibrium spectral moment sum rules for the first four moments of the single-particle Green's function. We then analytically confirm that this Green's function exactly satisfies the zeroth, first, and second moment sum rules, and provide numerical evidence that the third spectral moment sum rule also holds. Finally we calculate the two-particle Green's function in the atomic limit. These two functions are required inputs to create a self-consistent strong-coupling expansion for the lattice Green's functions of the Hubbard-Holstein model.
The self-consistent theory of the correlation effects in Highly Correlated Systems(HCS) is presented. The novel Irreducible Green's Functions(IGF) method is discused in detail for the Hubbard model and random Hubbard model. The interpolative solution for the quasiparticle spectrum, which is valid for both the atomic and band limit is obtained. The (IGF) method permits to calculate the quasiparticle spectra of many-particle systems with the complicated spectra and strong interaction in a very natural and compact way. The inelastic scattering corrections leads to the damping of the quasiparticles and are the main topic of the present consideration. The calculation of the damping has been done in a self-consistent way for both limits. For the random Hubbard model the weak coupling case has been considered and the self-energy operator has been calculated using the combination of the IGF method and Coherent Potential Approximation (CPA). The other applications of the method to s-f model, Anderson model, Heisenberg antiferromagnet, electron-phonon interaction models are discussed briefly.
We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.