Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Progress of Theoretical Physics

Cluster Transfer-Matrix Method for One-Dimensional Quantum Spin Systems

Hiroshi Betsuyaku

Published Apr 1, 1986• 19 citations
Publisher

Abstract

A formulation is given to map the one-dimensional S = 1/2 quantum spin systems onto the two-dimensional Ising model with many-spin interactions on the basis of the n-spin cluster decomposition (n CLD). The cluster transfer matrix is defined which characterizes the thermodynamical properties of the equivalent system. A prescription is presented to construct the cluster transfar matrix and to evaluate its maximum eigenvalue numerically. The method has been applied to the calculation of the energy and specific heat of the XY model for various values of the cluster size n and th Trotter size m. The convergence properties of the energy as a function of n and m have been studied in detail. The main results obtained are: (1) The low-temperature behavior is improved much more by considering larger spin clusters; (2) the energy depends on the spin cluster through the effective cluster size ne defined by ne =(n+1)/2; and (3) the energy converges linearly with 1/n2em2 except at very low temperatures.

DOI
10.1143/ptp.75.774

Related papers

Study of One-Dimensional Quantum Spin Systems by the Transfer-Matrix Method

Hiroshi Betsuyaku

The Suzuki-Trotter formula has been used to get the mth approximant to the classical representation of the partition function of the one-dimensional N-spin S=½ quantum spin systems. The equivalent two-dimensional (N×2m) Ising model with four-spin interactions has been studied in detail by using the numerically exact transfer-matrix method for T≥0.05 and m ≤8. The convergence properties have been examined in two different representations; checkerboard decomposition (CBD) and real-space decomposition (RSD). The spin correlation functions in RSD converge much faster than those in CBD. The limiting m→∞ behavior has been estimated from the extrapolation formula of the form: E(m)=E( ∞)+ α/ m2. The limiting values of the energy derived from the nearest-neighbor correlation agree with the correct values excellently.

Progress of Theoretical Physics•Feb 1, 1985•68 citations

Study of One-Dimensional Quantum Spin Systems by the Transger-Matrix Method. II: S = 1 Case

Hiroshi Betsuyaku, Terufumi Yokota

The one-dimensional S=1 quantum spin system has been mapped onto the equivalent two-dimensional Ising model on the basis of the equivalence theorem of Suzuki. The classical representation of the partition function has been expressed in terms of the maximum eigenvalue of the transfer matrix and evaluated by means of the numerically exact transfer-matrix method. The physical meaning of two kinds of correlation functions in RSD and CBD has been made clear in terms of the concept of the maximum eigenstate. Calculations have been made for both ferromagnetic and antiferromagnetic Heisenberg chains with several single-site anisotropy values, both positive and negative. Good agreements are found with some existing numerical results.

Progress of Theoretical Physics•Apr 1, 1986•30 citations

An Improved Transfer Matrix Method for Quantum Spin Systems

T. Tsuzuki

An improved transfer matrix method is proposed for calculating the partition function of quantum spin systems, which is expressed by the M-th approximant of the Suzuki-Trotter transformation with Trotter number M. The point is to incorporate the symmetry property of the Trotter subsystem in calculation. The size of our transfer matrix is reduced very much. Thermodynamics is described by the maximum eigenvalue of the transfer matrix. Our method is applied to the XY model with spin 1/2. We consider the cases of pair and tri spin decompositions with Trotter number up to 4 (M≤4). We can do our calculation analytically in those cases. In the pair spin decomposition we find results which have been obtained by the Monte Carlo simulation and by a transfer matrix method by computer. We find nice results in the tri spin decomposition. Our result in the case M=4 is as good as the result in M=8 by the previous Monte Carlo simulation on comparing them with the exact result of Katsura.

Progress of Theoretical Physics•Jun 1, 1985•23 citations

An Improved Transfer Matrix Method for Quantum Spin Systems. II: Isotropic Heisenberg Chain with Spin 1/2

T. Tsuzuki

An improved transfer matrix method which has been proposed by the author is applied to the isotropic Heisenberg chain with spin 1/2. Both ferro- and antiferromagnetic couplings are considered. We first explain how to find the size of transfer matrix in consideration of the symmetry of Trotter subsystem with Trotter number M. We then compute the internal energy, specific heat, magnetization and susceptibility in the presence as well as absence of magnetic field in the cases of the pair and tri spin decompositions with M=2, 4 and 6. Our results are qualitatively explained on the basis of the level schemes of the clusters of decomposition. We find good and rapid convergence to the infinite chain limit of Bonner and Fisher with increasing Trotter number in the tri spin decomposition except at very low temperatures, when we examine the zero field limit. Convergence with Trotter number in the case of finite field is also rapid in either decomposition. However there are interesting discrepancies, though small, between the results of Bonner and Fisher and ours in the tri spin decomposition, which remain as one of future problems to be investigated. Finally we present brief discussion about strategy and tactics for getting improved results by our method.

Progress of Theoretical Physics•Feb 1, 1986•11 citations

Modified effective-field approach to low-dimensional spin-1/2 systems

G. Kamieniarz, Raf Dekeyser, Grzegorz Musiał +2 more

Two different generalization schemes for the molecular field approximation are described. Both contain the possibility to incorporate generalized effective fields that induce correlations between the degrees of freedom at the boundary of a chosen cluster. Applications are presented for low-dimensional spin-1/2 systems. For the Ising models the critical couplings and critical exponents as well as the spontaneous magnetization curve are determined. An application of the method to the quantum $S$=$\frac{1}{2}$ one-dimensional Heisenberg model is presented and reliable low-temperature estimates of the specific heat are evaluated. The data are combined with the quantum transfer matrix predictions for large system sizes and this leads to the following prediction for the low-temperature specific heat: ${C=AT}^{\ensuremath{\alpha}}$ with $\ensuremath{\alpha}=0.35\ifmmode\pm\else\textpm\fi{}0.07$.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jul 1, 1997•7 citations

Investigation of Magnetic Properties of Spin 5/2 Ising Chain by Using Transfer Matrix Method

Mehmet Gökhan Şensoy, Mehmet Batı

A magnetic property of the one dimensional spin 5/2 Ising model under the magnetic field has been investigated by means of transfer matrix method. Thermodynamic response functions are also obtained for varying values of scaling temperature and scaling magnetic field. Entropy and heat capacity of the system were calculated by benefiting from the temperature dependencies of Helmholtz free energy. Our simulation results demonstrate that as the strength of the magnetic field is increased, heat capacity tend to shift to the relatively higher temperature regions, and these findings are consistent with previous results for the low spin values in one dimensional Ising systems.

Sakarya University Journal of Science•Oct 23, 2018•1 citation
A formulation is given to map the one-dimensional S = 1/2 quantum spin systems onto the two-dimensional Ising model with many-spin interactions on the basis of the n-spin cluster decomposition (n CLD). The cluster transfer matrix is defined which characterizes the thermodynamical properties of the equivalent system. A prescription is presented to construct the cluster transfar matrix and to evaluate its maximum eigenvalue numerically. The method has been applied to the calculation of the energy and specific heat of the XY model for various values of the cluster size n and th Trotter size m. The convergence properties of the energy as a function of n and m have been studied in detail. The main results obtained are: (1) The low-temperature behavior is improved much more by considering larger spin clusters; (2) the energy depends on the spin cluster through the effective cluster size ne defined by ne =(n+1)/2; and (3) the energy converges linearly with 1/n2em2 except at very low temperatures.
The Suzuki-Trotter formula has been used to get the mth approximant to the classical representation of the partition function of the one-dimensional N-spin S=½ quantum spin systems. The equivalent two-dimensional (N×2m) Ising model with four-spin interactions has been studied in detail by using the numerically exact transfer-matrix method for T≥0.05 and m ≤8. The convergence properties have been examined in two different representations; checkerboard decomposition (CBD) and real-space decomposition (RSD). The spin correlation functions in RSD converge much faster than those in CBD. The limiting m→∞ behavior has been estimated from the extrapolation formula of the form: E(m)=E( ∞)+ α/ m2. The limiting values of the energy derived from the nearest-neighbor correlation agree with the correct values excellently.
The one-dimensional S=1 quantum spin system has been mapped onto the equivalent two-dimensional Ising model on the basis of the equivalence theorem of Suzuki. The classical representation of the partition function has been expressed in terms of the maximum eigenvalue of the transfer matrix and evaluated by means of the numerically exact transfer-matrix method. The physical meaning of two kinds of correlation functions in RSD and CBD has been made clear in terms of the concept of the maximum eigenstate. Calculations have been made for both ferromagnetic and antiferromagnetic Heisenberg chains with several single-site anisotropy values, both positive and negative. Good agreements are found with some existing numerical results.
An improved transfer matrix method is proposed for calculating the partition function of quantum spin systems, which is expressed by the M-th approximant of the Suzuki-Trotter transformation with Trotter number M. The point is to incorporate the symmetry property of the Trotter subsystem in calculation. The size of our transfer matrix is reduced very much. Thermodynamics is described by the maximum eigenvalue of the transfer matrix. Our method is applied to the XY model with spin 1/2. We consider the cases of pair and tri spin decompositions with Trotter number up to 4 (M≤4). We can do our calculation analytically in those cases. In the pair spin decomposition we find results which have been obtained by the Monte Carlo simulation and by a transfer matrix method by computer. We find nice results in the tri spin decomposition. Our result in the case M=4 is as good as the result in M=8 by the previous Monte Carlo simulation on comparing them with the exact result of Katsura.
An improved transfer matrix method which has been proposed by the author is applied to the isotropic Heisenberg chain with spin 1/2. Both ferro- and antiferromagnetic couplings are considered. We first explain how to find the size of transfer matrix in consideration of the symmetry of Trotter subsystem with Trotter number M. We then compute the internal energy, specific heat, magnetization and susceptibility in the presence as well as absence of magnetic field in the cases of the pair and tri spin decompositions with M=2, 4 and 6. Our results are qualitatively explained on the basis of the level schemes of the clusters of decomposition. We find good and rapid convergence to the infinite chain limit of Bonner and Fisher with increasing Trotter number in the tri spin decomposition except at very low temperatures, when we examine the zero field limit. Convergence with Trotter number in the case of finite field is also rapid in either decomposition. However there are interesting discrepancies, though small, between the results of Bonner and Fisher and ours in the tri spin decomposition, which remain as one of future problems to be investigated. Finally we present brief discussion about strategy and tactics for getting improved results by our method.
Two different generalization schemes for the molecular field approximation are described. Both contain the possibility to incorporate generalized effective fields that induce correlations between the degrees of freedom at the boundary of a chosen cluster. Applications are presented for low-dimensional spin-1/2 systems. For the Ising models the critical couplings and critical exponents as well as the spontaneous magnetization curve are determined. An application of the method to the quantum SSS=12\frac{1}{2}21​ one-dimensional Heisenberg model is presented and reliable low-temperature estimates of the specific heat are evaluated. The data are combined with the quantum transfer matrix predictions for large system sizes and this leads to the following prediction for the low-temperature specific heat: ${C=AT}^{\ensuremath{\alpha}}$ with $\ensuremath{\alpha}=0.35\ifmmode\pm\else\textpm\fi{}0.07$.
A magnetic property of the one dimensional spin 5/2 Ising model under the magnetic field has been investigated by means of transfer matrix method. Thermodynamic response functions are also obtained for varying values of scaling temperature and scaling magnetic field. Entropy and heat capacity of the system were calculated by benefiting from the temperature dependencies of Helmholtz free energy. Our simulation results demonstrate that as the strength of the magnetic field is increased, heat capacity tend to shift to the relatively higher temperature regions, and these findings are consistent with previous results for the low spin values in one dimensional Ising systems.