A Study of the Transfer Matrix Method for the Classical Statistical Mechanics of One Dimensional Systems

E. T. Y. Hsieh

Wenlong Wang, Rogelio Díaz-Méndez, Raudys R. Capdevila

Abstract The aim of this work is to present a formulation to solve the one-dimensional Ising model using the elementary technique of mathematical induction. This formulation is physically clear and leads to the same partition function form as the transfer matrix method, which is a common subject in the introductory courses of statistical mechanics. In this way our formulation is a useful tool to complement the traditional more abstract transfer matrix method. The method can be straightforwardly generalised to other short-range chains, coupled chains and is also computationally friendly. These two approaches provide a more complete understanding of the system, and therefore our work can be of broad interest for undergraduate teaching in statistical mechanics.

V. I. Skripnik

Chunhoa Chhenh

The concept of the corner transfer matrix (CTM) was first discovered by Baxter in 1968, when he derived a set of variational matrix equations, which for finite matrix sizes, could be numerically solved to obtain a sequence of approximations for the statistical mechanical properties of a system of monomers and dimers on a rectangular lattice [1]. It was not until 1978, however, in a seminal paper entitled "Variational Approximations for Square Lattice Models in Statistical Mechanics" [4], that Baxter outlined his CTM variational method, which brought to light the potential power of the former objects to obtain numerics and series expansions for unsolved models and to calculate the order parameter of solved ones. Subsequent numerical work led to the realisation that the method, though general, was not very efficient; and increasing efficiency required making model specific modifications, which restricted the transferability of the resulting algorithm to other models. The CTM variational method was thus not widely adopted, despite holding much promise. More recently, however, Nishino and Okunishi discovered that White's density matrix renormalisation group (DMRG) algorithm [55, 56] could be efficiently extended to study two-dimensional classical lattice models, if the density matrix was approximated by Baxters CTMs. Numerical tests of their algorithm, the corner transfer matrix renormalisation group (CTMRG) method, were met with much success. Notable among these is its implementation within the infinite projected entangled-pair states (iPEPS) algorithm of Orus and Vidal to simulate the ground state of infinite two-dimensional quantum systems. In this thesis we review CTM derived variational methods: Baxter's original CTM iterative method, and developments of the CTM concept within the DMRG algorithm by Nishino and Okunishi, which led to the CTMRG method to calculate critical phenomena of classical and quantum systems. This will begin with elucidating the theoretical formalisms of both methods in two dimensions and their generalisations to three dimensions; followed by review of important application, namely, within the iPEPS algorithm of Orus and Vidal to numerically study infinite planar quantum systems.

V. I. Skripnik

Hiroshi Ueda, Andrej Gendiar, Tomotoshi Nishino

We report a physical background of the wave function prediction in the infinite system density matrix renormalization group (DMRG) method, from the view point of two-dimensional vertex model, a typical lattice model in statistical mechanics.Singular value decomposition applied to rectangular corner transfer matrices naturally draws matrix product representation for the maximal eigenvector of the row-to-row transfer matrix.The wave function prediction can be expressed as the insertion of an approximate half-column transfer matrix.This insertion process is in accordance with the scheme proposed by McCulloch recently.