For any purely fermionic system on a lattice in arbitrary space dimension we prove that there is a different but well-defined system on the same lattice, consisting of both bona fide fermions and bosons with an interaction depending on a parameter \ensuremath{\lambda}, which characterizes a same-site repulsion between particles. The energy spectrum and the scattering matrix of the former are identical to those in the finite energy sector of the latter, in the limit \ensuremath{\lambda}\ensuremath{\rightarrow}\ensuremath{\infty} (i.e., hard-core repulsion). This general theorem is applied to a resonance-boson model in a lattice fermion system, whose solution may have relevance to high-
Tc superconductivity.
We investigate the fermion spectrum within the Bose–Fermi–Hubbard model used for the description of boson-fermion mixtures of ultra-cold atoms in optical lattices. We used the method based on the Hubbard operator approach for an on-site basis. The equation for fermion Green’s function in the Bose–Fermi–Hubbard model is built; Green’s functions of higher orders are decoupled in the Hubbard-I approximation (the case of the strong on-site interaction). The corresponding spectral densities are calculated. In the case of hard-core bosons, the condition of appearance of additional bands in the fermion spectrum is investigated. It is shown that these bands exist only in the state with a Bose–Einstein condensate and appear because of the mixing of states with different numbers of bosons. These additional bands can be interpreted as a manifestation of composite excitations (when the appearance of a fermion on the site is accompanied by the simultaneous creation (or annihilation) of a boson).
A general system composed of many weakly interacting bosons and/or fermions is considered. The object is to develop a procedure for renormalizing the single-particle creation operators, so as to remove the interactions to successively higher orders of perturbation. The basic idea is that creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ must satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where
H is the Hamiltonian and ${\ensuremath{\omega}}_{i}$ are the particle energies. For the many-weakly-interacting-particle system, the zeroth-order boson and fermion creation operators satisfy these equations to zeroth order. It is shown that if the creation operators are renormalized so as to satisfy the Hamiltonian commutator equations to order
m, and also to satisfy the appropriate boson commutators and fermion anticommutators to order
m, then the problem is solved to order
m. In particular, the vectors formed by operating on the ground state with the renormalized creation operators, according to the usual boson and fermion occupation-number representation, are eigenvectors of
H and are orthonormal, all to order
m. A procedure is otained for finding the (
m+1)-order contributions to the particle creation energies, in terms of the
m-order operators. Explicit first-order calculations of these general results are provided for a system of bosons and a system of fermions, and these first-order results are shown to include similar results of Rayleigh-Schr\"odinger perturbation theory, the random-phase approximation, and the method of thermodynamic Green's functions. The problem of anharmonic lattice dynamics is studied in detail, and a method of undetermined coefficients is used to renormalize the phonon creation operators to first order. The phonon energies are calculated to second order, and this calculation shows that the interactions between renormalized phonons cannot be removed in second order. Statistical averages of the phonon energies give the energy shifts and lifetimes which have been calculated previously by various propagator techniques. In addition, the renormalized energies are used to calculate the temperature-dependent part of the Helmholtz free energy correct to second order. As a final example, electron-phonon interactions in a normal metal are studied. The electron and phonon creation operators are renormalized to first order, the particle energies are calculated to second order, and statistical averages of the particle energies recover the usual thermodynamic Green's-function results for energy shifts and lifetimes. These examples show the simplicity by which the renormalization procedure obtains a great amount of detailed information about the single-particle nature of a many-particle system.
Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green's function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear as promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green's function (GF2) method, where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in k-space are the key components of a computationally feasible algorithm. We apply this technique to the one-dimensional hydrogen lattice--a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.