This chapter, within the framework of classical statistical mechanics, discusses a family of models defined on one-dimensional lattices. It studies the simplest local examples: models that involve only interactions between nearest neighbours on the lattice. For such models, correlation functions can be calculated by a transfer matrix formalism. The chapter first describes some general properties of transfer matrices in one-dimensional models. This formalism is used to establish various properties of correlation functions, like the thermodynamic or infinite volume limit, the large-distance behaviour of the two-point correlation function, and introduces the very important concept of correlation length. Connected correlation functions, cumulants of the distribution, play a particularly important role. Indeed, these functions satisfy the cluster property, which characterizes their decay at large distance. The transfer matrix formalism is applied to the example of a Gaussian Boltzmann weight, which is studied in detail. The chapter calculates the partition function and correlation functions explicitly, and observes that , the correlation length diverges, making it possible to define a continuum limit. It shows that results of the continuum limit can be reproduced directly by solving a partial differential equation in which all traces of the initial lattice structure have disappeared. Finally, it exhibits a slightly more general class of models which share the same properties: divergent correlation length and continuum limit. Exercises are provided at the end of the chapter.
We develop a field theoretical approach to the classical two-dimensional (2D) models, particularly to 2D Ising model (2DIM) and
XYZ model, which is simple to apply for calculation of various correlation functions. We calculate the partition function of 2DIM and
XY model within the developed framework. Determinant representation of spin-spin correlation functions is derived using fermionic realization for the Boltzmann weights. The approach also allows formulation of the partition function of 2DIM in the presence of an external magnetic field.
The density-matrix renormalization group (DMRG) applied to transfer matrices allows it to calculate static as well as dynamical properties of one-dimensional quantum systems at finite temperature in the thermodynamic limit. To this end the quantum system is mapped onto a two-dimensional classical system by a Trotter-Suzuki decomposition. Here we discuss two different mappings: The standard mapping onto a two-dimensional lattice with checkerboard structure as well as an alternative mapping introduced by two of us. For the classical system an appropriate quantum transfer matrix is defined which is then treated using a DMRG scheme. As applications, the calculation of thermodynamic properties for a spin-1/2 Heisenberg chain in a staggered magnetic field and the calculation of boundary contributions for open spin chains are discussed. Finally, we show how to obtain real time dynamics from a classical system with complex Boltzmann weights and present results for the autocorrelation function of the XXZ-chain.
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.
Exact expressions for the partition function, spin pair correlation function, and susceptibility of the onedimensional isotropic classical Heisenberg model are obtained in zero external field with cyclic boundary conditions. It is shown that the methods used to derive these results enable the partition functions and susceptibilities of finite clusters of interacting classical spins to be evaluated in terms of the $3n\ensuremath{-}j$ symbols of Wigner. Exact results in one dimension are also obtained for the partition function and susceptibility of a "planar" classical Heisenberg model. In this model the spin vectors interact via a Heisenberg coupling but each spin vector is restricted to lie in a plane.The anisotropic classical Heisenberg model described by the Hamiltonian $\mathcal{H}=\ensuremath{-}\ensuremath{\Sigma}\stackrel{}{(\mathrm{ij})}2({{J}_{\mathrm{ij}}}^{x}{{s}_{i}}^{x}{{s}_{j}}^{x}+{{J}_{\mathrm{ij}}}^{y}{{s}_{i}}^{y}{{s}_{j}}^{y}+{{J}_{\mathrm{ij}}}^{z}{{s}_{i}}^{z}{{s}_{j}}^{z})\ensuremath{-}mH\ensuremath{\Sigma}\stackrel{N}{j=1}{{s}_{i}}^{z},$ where
six,
siy, and
siz are components of the unit vector
si, is also considered. A perturbation series for the zero-field free energy of the anisotropic model in one dimension with nearest-neighbor interactions
Jijx=Jijy=J and ${{J}_{\mathrm{ij}}}^{z}=\ensuremath{\gamma}J$ is developed in powers of $\ensuremath{\gamma}\ensuremath{-}1$ using the isotropic model as the unperturbed system. Detailed calculations are performed to third order in $\ensuremath{\gamma}\ensuremath{-}1$. It is found that the perturbation series for the energy per spin breaks down as $T\ensuremath{\rightarrow}0$. A high-temperature series expansion for the anisotropic model, which is valid for a general interaction potential and lattice, is derived by generalizing the methods developed by Horwitz and Callen for the Ising model. This series is rearranged to give a simplified diagram expansion. Finally, a practical technique for calculating the high-temperature series expansions of the zero-field free energy and susceptibility of the isotropic classical Heisenberg model is presented.
We show that there are novel generalizations of statistics for identical particles in one dimension. These arise due to possible boundary conditions on wave functions, or equivalently due to δ-function interactions, at points of coincidence of particle coordinates. Special choices of these boundary conditions describe bosons, fermions or paraparticles. The general solution for the boundary conditions involves vector-valued wave functions and statistics with non-Abelian features, even though the classical configuration space has an Abelian fundamental group. Physical models leading to such non-Abelian statistics, involving for example δ-function interactions of spins, are constructed. Properties under parity and time reversal of the new boundary conditions are studied. It is shown that the Bethe ansatz does not always give eigenstates of energy. When it does, the wave numbers in the ansatz for three or more particles must satisfy a weaker form of the Yang-Baxter equations and certain additional equations. The energy spectrum and wave functions for identical particles on R1 or S1 are discussed in some simple cases. Our work generalizes the previous work of Lieb and Liniger, Lieb, Yang and Leinaas and Myrheim.