Transfer-matrix construction in 1D

Rewrite 1D Boltzmann weights as powers of a matrix to compute partition functions, free energies, and correlations.
Transfer-matrix construction in 1D

For one-dimensional, nearest-neighbor lattice models, the transfer matrix turns the sum over all into linear algebra. It provides an exact route to the , the , and decay of in 1D.

Nearest-neighbor form and the transfer matrix

Consider a 1D chain of length NN with spins σiS\sigma_i\in S and a nearest-neighbor energy of the form

H(σ1,,σN)=i=1NU(σi,σi+1),σN+1=σ1 H(\sigma_1,\dots,\sigma_N)=\sum_{i=1}^{N} U(\sigma_i,\sigma_{i+1}), \qquad \sigma_{N+1}=\sigma_1

(periodic boundary conditions). Fix .

Define the transfer matrix TT (a indexed by SS) by

Tab=exp ⁣(βU(a,b)). T_{ab}=\exp\!\big(-\beta\,U(a,b)\big).

More generally, if the Hamiltonian includes on-site terms V(σi)V(\sigma_i), one often splits them symmetrically into neighboring factors so that Tab=exp(β[U(a,b)+12V(a)+12V(b)])T_{ab}=\exp(-\beta[U(a,b)+\tfrac12 V(a)+\tfrac12 V(b)]); this keeps the representation exact and symmetric.

Partition function as a trace

The periodic-chain partition function is

ZN=σ1,,σNexp ⁣(βH(σ1,,σN))=Tr(TN), Z_N=\sum_{\sigma_1,\dots,\sigma_N}\exp\!\big(-\beta H(\sigma_1,\dots,\sigma_N)\big) =\mathrm{Tr}\big(T^N\big),

where Tr\mathrm{Tr} is the . (For open boundary conditions, one obtains a bilinear form ZN=,TN1rZ_N=\langle \ell, T^{N-1} r\rangle with boundary vectors encoding the end weights.)

This representation is a concrete instance of “compute thermodynamics from logZ\log Z,” as in .

Thermodynamic limit and dominant eigenvalue

Let the eigenvalues of TT satisfy λ0λ1|\lambda_0|\ge |\lambda_1|\ge\cdots. Under mild positivity assumptions (typical for physical transfer matrices), λ0>0\lambda_0>0 is the unique largest eigenvalue. Then, in the NN\to\infty,

1NlogZNlogλ0,f=1βlogλ0, \frac{1}{N}\log Z_N\longrightarrow \log \lambda_0, \qquad f = -\frac{1}{\beta}\log \lambda_0,

where ff is the free energy per site (a 1D analogue of the thermodynamic density encoded by ).

Correlations and correlation length

Transfer matrices also control spatial correlations. For suitable local observables AA and BB at separation rr, one typically finds asymptotic decay

A0BrA0Br    (λ1λ0)ras r, \langle A_0 B_r\rangle - \langle A_0\rangle\langle B_r\rangle \;\propto\; \left(\frac{\lambda_1}{\lambda_0}\right)^r \quad\text{as } r\to\infty,

so the associated decays exponentially. The corresponding is

ξ1=logλ1λ0. \xi^{-1} = -\log\left|\frac{\lambda_1}{\lambda_0}\right|.

Physical interpretation

The transfer matrix is a “propagator” that advances the system by one lattice step. Because ZNZ_N is governed by a largest eigenvalue, thermodynamic quantities in 1D with finite-range interactions are typically analytic in parameters for TT with strictly positive entries—one reason phase transitions are absent in many 1D short-range models.