The nearest-neighbor quantum antiferromagnetic (AF) Heisenberg model for spin-
1∕2 on a two-dimensional square lattice is studied in the auxiliary-fermion representation. Expressing spin operators by canonical fermionic particles requires a constraint on the fermion charge
Qi=1 on each lattice site
i, which is imposed approximately through the thermal average. The resulting interacting fermion system is first treated in mean-field theory (MFT), which yields an AF ordered ground state and spin waves in quantitative agreement with conventional spin-wave theory. At finite temperature a self-consistent approximation beyond mean field is required in order to fulfill the Mermin-Wagner theorem. We first discuss a fully self-consistent approximation, where fermions are renormalized due to fluctuations of their spin density, in close analogy to FLEX. While static properties like the correlation length, $\ensuremath{\xi}(T)\ensuremath{\propto}\mathrm{exp}(aJ∕T)$, come out correctly, the dynamical response lacks the magnon-like peaks which would reflect the appearance of short-range order at low
T. This drawback, which is caused by overdamping, is overcome in a ``minimal self-consistent approximation'' (MSCA), which we derive from the equations of motion. The MSCA features dynamical scaling at small energy and temperature and is qualitatively correct both in the regime of order-parameter relaxation at long wavelengths $\ensuremath{\lambda}>\ensuremath{\xi}$ and in the short-range-order regime at $\ensuremath{\lambda}<\ensuremath{\xi}$. We also discuss the impact of vertex corrections and the problem of pseudo-gap formation in the single-particle density of states due to long-range fluctuations. Finally we show that the (short-range) magnetic order in MFT and MSCA helps to fulfill the constraint on the local fermion occupancy.
A Wigner-Jordan (WJ) transformation is used to study the one-dimensional (1D) and two-dimensional (2D) quantum antiferromagnetic Heisenberg model. The advantage of using the Wigner-Jordan transformation is that it preserves all spin-commutation relations as well as the spin on-site exclusion principle. In the 1D case a nearest-neighbor covalent-bonding state of the WJ fermions is found to have a ground-state energy (-0.4351J per site) comparable with that from the Bethe-ansatz solution (-0.4431J per site), and a linear energy spectrum at low energies with velocity 1.6366J, in close agreement with the velocity obtained by Haldane for the quantum antiferromagnetic Heisenberg model with 1/
d2 interaction (1.5708J). The method used for studying the 1D model is then applied to the 2D Heisenberg model in a square lattice. The resulting state at finite temperature is the in-phase flux state, i.e., a flux state of the Wigner-Jordan spinless fermions with an in-phase fermion orbital current circulating around each elementary plaquette.The single-particle excitation spectrum, i.e., the energy dispersion of reversing orientation of a spin in the system, of the in-phase flux state is overall similar to that of spin waves, with significant difference near the edge point k=(\ensuremath{\pi},0), at which the excitation energy of the WJ fermion is zero, whereas that of the spin-wave excitation is 2J. The specific heat of the in-phase flux state predicts a correct temperature dependence over the entire temperature range, namely a
T2 dependence at low temperature, a peak near T/J=0.6, and a 1/
T2 decreasing at high temperature. It also gives excellent agreement with the specific heat calculated from numerical methods. In contrast, the spin-wave theory only correctly predicts a
T2 dependence at low temperature. The Raman spectrum of the in-phase flux state is calculated and shows significant improvement over that from spin-wave theories compared with the experimental spectrum of
La2CuO4. The exchange parameter, J, obtained from the comparison is 1060 ${\mathrm{cm}}^{\mathrm{\ensuremath{-}}1}$, in agreement with that obtained from analyzing the neutron-scattering data.
The one-dimensional Heisenberg model with S = 1/2 is treated with the use of the two-time Green's functions. The hierarchy of the equations of motion of the Green'f functions is decoupled at a stage one-step further than Tyablikov's decoupling. The thermal average of the spin component, ≪ Sz >, is set to zero, becaure the long-range order does not exist in one dimension. Instead, our Green's functions are expressed in terms of the correlation functions cn ≡4≪ S0zSnz >. The Green's function is essentially of the form representing undamped spin waves, whose spactrum depends on c1, c2 and one more parameter. They are determined by the requirement that c1 and c2 should be self-consistent and that c0 should be unity. The self-consistency equations have been solved analytically at high- and low-temperature limits, and also solved numerically in the whole range of the temperature. Thermodynamic quantities have been calculated using these solutions. It has turned out that the theory gives the correct high-temperature expansion for the thermodynamic quantities and the correlation functions. The latter is expressed by cn = (J / 4 kBT )n. In the case of ferromagnetic coupling, the correlation function cn at T = 0 is equal to 1/3 for all n's. This is what is expected from the correct ground state of the ferromagnetic Heisenberg system. The spin-wave spectrum at T = 0 also agrees with the correct one. At T ≪J/ kB, we find that the specific heat goes as T1/2 and the susceptibility goes as T-2. The gross feature of the temperature-dependence of the thermodynamic quantities agrees with Bonner and Fisher. In the case of the antiferromagnetic Heisenberg model, we find c1 = -0.55407 and c2 = 0.16100 at T = 0, which are fairly close to the exact values. The thermodynamic quantities are also in gross agreement with Bonner and Fisher.
The Hartree-Fock-RPA approach is applied to the 1D anti-ferromagnetic Heisenberg model in the Jordan-Wigner representation. Somewhat contrary to expectation, this leads to reasonable results for spectral functions and sum rules in the symmetry unbroken phase. In a preliminary application of Self-Consistent RPA to finite size chains strongly improved results are obtained. pacs: 75.40.Gb, 75.50.Ee, 75.10.Pq, 75.10.Jm The 1D Anti-Ferromagnetic Heisenberg Model (AFHM) [1] belongs to one of the most classic many body research fields. It was the first many-body model to be solved exactly by Bethe with his famous Bethe ansatz solution [2] and it has been a playground for testing and developing many body approaches ever since. In spite of tremendous progress in the understanding of many facets of the model, it continues to be a very active field of interest and research [3, 4]. In an attempt to apply extended RPA-like theories, as e.g. the recently developed Self-Consistent RPA (SCRPA) approach [5, 6], we found out, somewhat to our surprise, that even the standard RPA theory has so far not fully been developed in its application to the 1D-AFHM. This probably stems from the fact that the HF-RPA scheme is usually considered as unreliable in low dimensions [7]. However, it will be found that RPA shows interesting features even in this 1D-version of the model. On the other hand SCRPA shows very promising results in a number of cases [5] and in a first application for a small numbers of sites we here get very close agreement with results from exact diagonalization. We mostly consider the isotropic Heisenberg model with anti-ferromagnetic coupling, however some consideration will also be given to the anisotropic case, for instance to the
An approach for studying phase transitions in systems with quantum degrees of freedom is discussed. Starting from the microscopic Hamiltonian of a quantum model, we first derive a set of exact differential equations for the free energy and the correlation functions describing the effects of fluctuations on the thermodynamics of the system. These equations reproduce the full renormalization group structure in the neighborhood of a critical point keeping, at the same time, full information on the nonuniversal properties of the model. As a concrete application we investigate the phase diagram of a Heisenberg antiferromagnet in a staggered external magnetic field. At long wavelengths the known relationship to the quantum nonlinear $\ensuremath{\sigma}$ model naturally emerges from our approach. By representing the two point function in an approximate analytical form, we obtain a closed partial differential equation which is then solved numerically. The results in three dimensions are in good agreement with available quantum Monte Carlo simulations and series expansions. More refined approximations to the general framework presented here and few applications to other models are briefly discussed.
A simple real-space renormalization method yields the ground-state energy of the Heisenberg antiferromagnet. We find the ground-state energy per spin for
s=21 (-0.4438 in 1D, -0.6723 in 2D) and
s=1 (-1.388 in 1D and -1.907 in 2D) to three-figure accuracy, using properties of relatively small odd-numbered clusters. Our results provide reasonable proof for long-range order in the ground state of the
s=21 Heisenberg antiferromagnet in 2D.