Transfer-matrix method (statistical mechanics)

Higher-dimensional models contain even more summations. For systems with more than a few particles, such expressions can quickly become too complex to work out directly, even by computer.

Instead, the partition function can be rewritten in an equivalent way. The basic idea is to write the partition function in the form

${\displaystyle {\mathcal {Z}}=\mathbf {v} _{0}\cdot \left\{\prod _{k=1}^{N}\mathbf {W} _{k}\right\}\cdot \mathbf {v} _{N+1}}$

where v₀ and v_N+1 are vectors of dimension p and the p × p matrices W_k are the so-called transfer matrices. In some cases, particularly for systems with periodic boundary conditions, the partition function may be written more simply as

${\displaystyle {\mathcal {Z}}=\operatorname {tr} \left\{\prod _{k=1}^{N}\mathbf {W} _{k}\right\}}$

where "tr" denotes the matrix trace. In either case, the partition function may be solved exactly using eigenanalysis. If the matrices are all the same matrix W, the partition function may be approximated as the N^th power of the largest eigenvalue of W, since the trace is the sum of the eigenvalues and the eigenvalues of the product of two diagonal matrices equals the product of their individual eigenvalues.

The transfer-matrix method is used when the total system can be broken into a sequence of subsystems that interact only with adjacent subsystems. For example, a three-dimensional cubical lattice of spins in an Ising model can be decomposed into a sequence of two-dimensional planar lattices of spins that interact only adjacently. The dimension p of the p × p transfer matrix equals the number of states the subsystem may have; the transfer matrix itself W_k encodes the statistical weight associated with a particular state of subsystem k − 1 being next to another state of subsystem k.

Importantly, transfer matrix methods allow to tackle probabilistic lattice models from an algebraic perspective, allowing for instance the use of results from representation theory.

As an example of observables that can be calculated from this method, the probability of a particular state ${\displaystyle m}$ occurring at position x is given by:

${\displaystyle \mathrm {Pr} _{m}(x)={\frac {\operatorname {tr} \left[\prod _{k=1}^{x}\mathbf {W} _{k}\mathbf {Pj} \prod _{k'=x+1}^{N}\mathbf {W} _{k'}\right]}{\operatorname {tr} \left[\prod _{k=1}^{N}\mathbf {W} _{k}\right]}}}$

Where ${\displaystyle Pj}$ is the projection matrix for state ${\displaystyle m}$, having elements ${\displaystyle Pj_{\mu \nu }=\delta _{\mu \nu }\delta _{\mu m}}$

Transfer-matrix methods have been critical for many exact solutions of problems in statistical mechanics, including the Zimm–Bragg and Lifson–Roig models of the helix-coil transition, transfer matrix models for protein-DNA binding, as well as the famous exact solution of the two-dimensional Ising model by Lars Onsager.

• Transfer operator

• Rodney J. Baxter (1982). Exactly solved models in statistical mechanics. Academic Press. ISBN 978-0-12-083182-1.
• Teif V.B. (2007). "General transfer matrix formalism to calculate DNA-protein-drug binding in gene regulation". Nucleic Acids Res. 35 (11): e80. doi:10.1093/nar/gkm268. PMC 1920246. PMID 17526526.
• Efremov AK, Winardhi RS, Yan J (2016). "Transfer-matrix calculations of DNA polymer micromechanics under tension and torque constraints". Phys. Rev. E. 94 (3) 032404. Bibcode:2016PhRvE..94c2404E. doi:10.1103/PhysRevE.94.032404. PMID 27739846.
• Efremov AK, Yan J (2018). "Transfer-matrix calculations of the effects of tension and torque constraints on DNA-protein interactions". Nucleic Acids Res. 46 (13): 6504–6527. doi:10.1093/nar/gky478. PMC 6061897. PMID 29878241.