We derive a variational cluster approximation for Heisenberg spin systems at finite temperature based on the ideas of the self-energy functional theory by Potthoff for fermionic and bosonic systems with local interactions. Partitioning the real system into a set of clusters, we find an analytical expression for the auxiliary free energy, depending on a set of variational parameters defined on the cluster, whose stationary points provide approximate solutions from which the thermodynamics of spin models can be obtained. We explicitly describe the technical details of how to evaluate the free energy for finite clusters and remark on specific problems and possible limitations of the method. To test the approximation we apply it to the antiferromagnetic spin chain and compare the results for varying cluster sizes and choices of variational parameters with the exact Bethe ansatz solution.
The motion of a single hole in a two-dimensional Heisenberg antiferromagnet (hAF) is studied in a representation where holes are described as spinless fermions (holons) and spins as normal bosons. Assuming long-range AF order the spin dynamics is treated in linear spin-wave theory. The formulation highlights the close relation with the conventional polaron problem. The holon Green's function is calculated self-consistently within the Born approximation using finite-cluster geometries for the numerical solution. As a remarkable result we find close agreement with the spectral function A(k,\ensuremath{\omega}) of a hole calculated by exact diagonalization methods. A(k,\ensuremath{\omega}) is characterized by a narrow quasiparticle (QP) peak at the low-energy side of the spectrum, which is well separated from the incoherent part for large enough J values. A complete characterization of our solution is given, including the spectral weight, the dispersion relation, and effective masses of the QP state. A finite-size-scaling study gives a nonvanishing spectral weight of the QP in the thermodynamic limit for values J/t typical for copper oxide superconductors. Our calculations indicate that the self-consistent Born approximation is a valuable scheme for characterizing the dynamics of a hole in the t-J model, even in the strong-coupling regime.
The thermodynamic effects of polaronic interactions in the Hubbard model are studied by high-temperature perturbation theory. Two electron-lattice interaction terms are added to the Hubbard Hamiltonian. The transfer term is treated as the perturbation, so the results are valid for narrow bandwidths and high temperatures. A one-particle Green's-function method is used to find an explicit expression for the secondorder change in the partition function, and a simple analytic expression is given for the weak-binding limit. The magnetic susceptibility and specific heat are computed numerically and given for several cases of polaron binding energies and molecular energy levels. The results are compared with the limit of no electron-lattice coupling and it is shown that they cannot be explained as a simple renormalization of the Hubbard parameters
t and
U.
We formulate a finite-temperature scheme of the variational cluster approximation (VCA) particularly suitable for an exact-diagonalization cluster solver. Based on the analytical properties of the single-particle Green's function matrices, we explicitly show the branch-cut structure of logarithm of the complex determinant functions appearing in the self-energy-functional theory (SFT) and whereby construct an efficient scheme for the finite-temperature VCA. We also derive the explicit formulas for entropy and specific heat within the framework of the SFT. We first apply the method to explore the antiferromagnetic order in a half-filled Hubbard model by calculating the entropy, specific heat, and single-particle excitation spectrum for different values of on-site Coulomb repulsion
U and temperature
T. We also calculate the
T dependence of the single-particle excitation spectrum in the strong coupling region, and discuss the overall similarities to and the fine differences from the spectrum obtained by the spin-density-wave mean-field theory at low temperatures and the Hubbard-I approximation at high temperatures. Moreover, we show a necessary and sufficient condition for the third law of thermodynamics in the SFT. On the basis of the thermodynamic properties, such as the entropy and the double occupancy, calculated via the
T and/or
U derivative of the grand potential, we obtain a crossover diagram in the
(U,T) plane, which separates a Slater-type insulator and a Mott-type insulator. Next, we demonstrate the finite-temperature scheme in the cluster-dynamical-impurity approximation (CDIA), i.e., the VCA with noninteracting bath orbitals attached to each cluster, and study the paramagnetic Mott metal-insulator transition in the half-filled Hubbard model. Formulating the finite-temperature CDIA, we first address a subtle issue regarding the treatment of the artificially introduced bath degrees of freedom which are absent in the originally considered Hubbard model. We then apply the finite-temperature CDIA to calculate the finite-temperature phase diagram in the
(U,T) plane. Metallic, insulating, coexistence, and crossover regions are distinguished from the bath-cluster hybridization-variational-parameter dependence of the grand-potential functional. We find that the Mott transition at low temperatures is discontinuous, and the coexistence region of the metallic and insulating states persists down to zero temperature. The result obtained here by the finite-temperature CDIA is complementary to the previously reported zero-temperature CDIA phase diagram.
The one-dimensional spin-1 Heisenberg antiferromagnetic chains are studied by a generalized two-time Green's-function method at finite temperature. 〈
Siz〉 is set to zero at each site to guarantee the lack of the long-range order in the one-dimensional Heisenberg system. The Green's functions are decoupled in terms of the correlation functions of spin operators. A set of self-consistent equations of the correlation functions are derived and solved numerically. Correlation functions of spin operators, thermodynamic properties such as the internal energy, and specific heat in the entire temperature region are obtained. The Haldane gap in the excitation spectrum appears naturally in the analytical result. When k is zero, there is a gap 2\ensuremath{\Delta}, which is 1.0J and close to the theoretical result of White and Huse. And there is a broad maximum in the temperature dependence of the specific heat. The properties of the thermodynamic quantities are consistent with the numerical results on the finite chains. \textcopyright{} 1996 The American Physical Society.