Effective Hamiltonians are usually constructed by using canonical transformations or projection techniques. In contrast to this, we present a method for systems with arbitrary Hilbert space based on the introduction of cumulants. Cumulants guarantee size consistency, a property that is not always evident in other treatments. As a non-trivial example of use the derived method is applied to the strong-coupling limit of the half-filled Hubbard model on a general lattice in arbitrary spatial dimension for which the fourth-order expansion in t / U of the effective Hamiltonian is derived.
Cumulants represent a natural language for expressing macroscopic properties of many-body systems. The most important property of cumulants is that of size consistency, i.e. a cumulant expression for an extensive variable scales with the size of the system, independent of possible further approximations used in the evaluation procedure. Cumulants can be considered as a generalization of linked diagrams known from diagrammatic technique of many-body theory. In this paper we outline a recently introduced method based on cumulants in order to derive expressions for zero-temperature properties of many-particle systems, i.e. the ground-state energy, static expectation values and dynamical correlation functions. This cumulant formalism allows one to describe weakly and strongly correlated systems along the same lines. We show that the coupled-cluster method known from quantum chemistry can be derived from our cumulant approach. Finally, we demonstrate the usefulness of the cumulant method by applying it to examples from solid-state physics and quantum chemistry.
By applying projection techniques a unified description of weakly and strongly correlated electronic systems can be given. For that purpose the theory is formulated in Liouville instead of Hilbert space. By introducing an appropriate metric in terms of cumulants it is ensured that all energies are size consistent. An energy expression is derived which is a generalization of Goldstone’s linked cluster theorem for arbitrary Hamiltonians. It is suitable for applying the projection method of Zwanzig and Mori. The equations for the correlation energy within the Local Ansatz and within coupled electron pair approximation CEPA-2 follow straightforwardly from the present theory.
Introducing low-energy effective Hamiltonians is usual to grasp most correlations in quantum many-body problems. For instance, such effective Hamiltonians can be treated at the mean-field level to reproduce some physical properties of interest. Employing effective Hamiltonians that contain many-body correlations renders the use of perturbative many-body techniques difficult because of the overcounting of correlations. In this work, we develop a strategy to apply an extension of the many-body perturbation theory starting from an effective interaction that contains correlations beyond the mean-field level. The goal is to re-organize the many-body calculation to avoid the overcounting of correlations originating from the introduction of correlated effective Hamiltonians in the description. For this purpose, we generalize the formulation of the Rayleigh-Schrödinger perturbation theory by including free parameters adjusted to reproduce the appropriate limits. In particular, the expansion in the bare weak-coupling regime and the strong-coupling limit serves as a valuable input to fix the value of the free parameters appearing in the resulting expression. This method avoids double counting of correlations using beyond-mean-field strategies for the description of many-body systems. The ground state energy of various systems relevant for ultracold atomic, nuclear, and condensed matter physics is reproduced qualitatively beyond the domain of validity of the standard many-body perturbation theory. Finally, our method suggests interpreting the formal results obtained as an effective field theory using the proposed reorganization of the many-body calculation. The results, like ground state energies, are improved systematically by considering higher orders in the extended many-body perturbation theory while maintaining a straightforward polynomial expansion.
We extend and benchmark the recently developed Effective-Hamiltonian (EFFH) method [PRX Quantum 4, 020307 (2023)] as an approximation to the equilibrium state (``mean-force Gibbs state'') of a quantum system at strong coupling to a thermal bath. The EFFH method is an approximate framework. Through a combination of the reaction-coordinate mapping, a polaron transformation, and a controlled truncation, it imprints the system-bath coupling parameters into the system's Hamiltonian. First, we develop a variational EFFH technique. In this method, the system's parameters are renormalized by both the system-bath coupling parameters (as in the original EFFH approach) and bath temperature. Second, adopting the generalized spin-boson model, we benchmark the equilibrium state from the EFFH treatment against numerically exact simulations and demonstrate a good agreement for both polarization and coherences using the Brownian spectral function. Third, we contrast the (normal and variational) EFFH approach with the familiar (normal and variational) polaron treatment. We show that the two methods predict a similar structure for the equilibrium state, albeit the EFFH approach offers the advantage of simpler calculations and closed-form analytical results. Altogether, we argue that for temperatures comparable to the system's frequencies, the EFFH methodology provides a good approximation for the mean-force Gibbs state in the full range of system-bath coupling, from ultraweak to ultrastrong.
The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid
He3, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$.
F is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the
f(rij) embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements
Hmn,
Nmn of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the
Hmn are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (
Hmn) and (
Nmn), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (
Hmn) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (
Hmn). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (
Hmn). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any
F in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{},
Q-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid
He3 has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.