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

Advances in Pure Mathematics

Approximate Algebraic Solution of 3D Ising Model —Transfer-Matrix Method

Andras Szasz

Published Jan 1, 2026• 1 citation
Publisher

Abstract

The exact analytical solution of the three-dimensional (3D) Ising model remains a central unresolved problem in statistical mechanics. While exact solutions exist in one and two dimensions, the three-dimensional case has resisted closed-form treatment. While standard Kramers-Wannier duality fails to yield a self-dual spin model in 3D, instead mapping to a gauge theory, we introduce a quasi-duality transformation by perturbing the isotropic cubic lattice with a slight, localized spatial inhomogeneity. Specifically, this is introduced via the boundary correction term A n ( ? )= δ η,nN ( s η ? s η?N+1 ? ? s η ? s η+1 ? ) , which arises from flattening the 2D planes into 1D row-continuous chains. By enforcing a condition where the introduced inhomogeneity algebraically absorbs the leading-order topological mismatch of the dual gauge fields, we recover an effective self-dual relation. This algebraic solution yields a critical coupling of K c = 1 4 ln( 1+ 2 )≈0.22033 , which remarkably deviates by only ~0.6% from the widely accepted homogeneous Monte Carlo estimation of K c ( MC ) ?0.22165 . We discuss the derivation and outline future applications of this method to study phase transitions in complex, structurally disordered systems. The method provides insight into how weak inhomogeneity restores a form of dual symmetry that approximates the true critical manifold. Implications for critical phenomena, renormalization structure, and complex biological systems operating near criticality are discussed.

DOI
10.4236/apm.2026.166021

Related papers

Exact solution of the three-dimensional (3D) Z 2 lattice gauge theory

Zhidong Zhang

Abstract The gauge invariance is an important phenomenon widely occurring in various matters with many-body interactions. The gauge theories are very important for understanding the evolution of fundamental physical processes in particle physics, high-energy physics, and condensed matter physics. In this work, the origin of nonlocal effects is inspected, and the contributions of nontrivial topological structures to physical properties are investigated in detail for both the three-dimensional (3D) Ising model and the 3D Z 2 lattice gauge model. Then, the exact solution for the 3D Z 2 lattice gauge theory is derived by the duality between the two models. The partition function, the spontaneous magnetization, the true range of correlation, and the critical point are rigorously derived for the 3D Z 2 lattice gauge theory. The critical exponents of the 3D Z 2 lattice gauge theory are in the same universality class as the 3D Ising model, which are α = 0, β = 3/8, γ = 5/4, δ = 13/3, η = 1/8, and ν = 2/3. Several fundamental issues, such as dimensionality, duality, symmetry, manifold and degenerate states, are investigated for these many-body interacting spin systems. The connections with superfluid, superconductors, particle physics, etc. are evaluated. Furthermore, physical significances and mathematical aspects of the 3D Z 2 lattice gauge theory are discussed with respect to topology, geometry, and algebra. The present results for the 3D Z 2 lattice gauge theory are helpful for figuring out the basic features of 3D lattice gauge theories with other symmetries, such as U(1), SU(2), and SU(3). This work shines a light on the interdisciplinarity between many-body interacting systems, algebra, topology, and geometry.

arXiv:2412.10412•Open Physics•Jan 1, 2025•6 citations

Conjectures on the exact solution of three-dimensional (3D) simple orthorhombic Ising lattices

Z.D. Zhang

We report conjectures on the three-dimensional (3D) Ising model of simple orthorhombic lattices, together with details of calculations for a putative exact solution. Two conjectures, an additional rotation in the fourth curled-up dimension and weight factors on the eigenvectors, are proposed to serve as a boundary condition to deal with the topologic problem of the 3D Ising model. The partition function of the 3D simple orthorhombic Ising model is evaluated by spinor analysis, employing these conjectures. Based on the validity of the conjectures, the critical temperature of the simple orthorhombic Ising lattices could be determined by the relation of KK* = KK′ + KK″ + K′K″ or sinh 2K · sinh 2(K′ + K″ + (K′K″/K)) = 1. For a simple cubic Ising lattice, the critical point is putatively determined to locate exactly at the golden ratio , as derived from K* = 3K or sinh 2K · sinh 6K = 1. If the conjectures would be true, the specific heat of the simple orthorhombic Ising system would show a logarithmic singularity at the critical point of the phase transition. The spontaneous magnetization of the simple orthorhombic Ising ferromagnet is derived explicitly by the perturbation procedure, following the conjectures. The spin correlation functions are discussed on the terms of the Pfaffians, by defining the effective skew-symmetric matrix Aeff. The true range κx of the correlation and the susceptibility of the simple orthorhombic Ising system are determined by procedures similar to those used for the two-dimensional Ising system. The putative critical exponents derived explicitly for the simple orthorhombic Ising lattices are α = 0, β = 3/8, γ = 5/4, δ = 13/3, η = 1/8 and ν = 2/3, showing the universality behaviour and satisfying the scaling laws. The cooperative phenomena near the critical point are studied and the results based on the conjectures are compared with those of approximation methods and experimental findings. The putative solutions have been judged by several criteria. The deviations of the approximation results and the experimental data from the solutions are interpreted. Based on the solution, it is found that the 3D-to-2D crossover phenomenon differs with the 2D-to-1D crossover phenomenon and there is a gradual crossover of the exponents from the 3D to the 2D values. Special attention is also paid to the extra energy caused by the introduction of the fourth curled-up dimension and the states at/near infinite temperature revealed by the weight factors of the eigenvectors. The physics beyond the conjectures and the existence of the extra dimension are discussed. The present work is not only significant for statistical and condensed matter physics, but also fill the gap between the quantum field theory, cosmology theory, high-energy particle physics, graph theory and computer science.

arXiv:0705.1045•The Philosophical Magazine A Journal of Theoretical Experimental and Applied Physics•Oct 30, 2007•127 citations

Transfer matrix and Monte Carlo tests of critical exponents in lattice models

J. Kaupužs

The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems of square geometry with periodic boundaries oriented either along <10> or along <11> direction have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.

arXiv:cond-mat/0201221•Jan 14, 2002•0 citations

Mathematical structure of the three-dimensional (3D) Ising model

Zhang Zhi-dong

An overview of the mathematical structure of the three-dimensional (3D) Ising model is given, from the viewpoints of topologic, algebraic and geometric aspects. By analyzing the relations among transfer matrices of the 3D Ising model, Reidemeister moves in the knot theory, Yang-Baxter and tetrahedron equations, the following facts are illustrated for the 3D Ising model: 1) The complexified quaternion basis constructed for the 3D Ising model represents naturally the rotation in a (3 + 1) - dimensional space-time, as a relativistic quantum statistical mechanics model, which is consistent with the 4-fold integrand of the partition function by taking the time average. 2) A unitary transformation with a matrix being a spin representation in 2^(nlo)-space corresponds to a rotation in 2nlo-space, which serves to smooth all the crossings in the transfer matrices and contributes as the non-trivial topologic part of the partition function of the 3D Ising model. 3) A tetrahedron relation would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model, and its existence is guaranteed also by the Jordan algebra and the Jordan-von Neumann-Wigner procedures. 4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases Φx, Φy, and Φz. The relation with quantum field and gauge theories, physical significance of weight factors are discussed in details. The conjectured exact solution is compared with numerical results, and singularities at/near infinite temperature are inspected. The analyticity in β = 1/(kB T) of both the hard-core and Ising models has been proved for β > 0, not for β = 0. Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.

arXiv:1305.2956•Chinese Physics B•Mar 1, 2013•62 citations

The bond-algebraic approach to dualities

Emilio Cobanera, Gerardo Ortíz, Zohar Nussinov

An algebraic theory of dualities is developed based on the notion of bond algebras. It deals with classical and quantum dualities in a unified fashion explaining the precise connection between quantum dualities and the low temperature (strong-coupling)/high temperature (weak-coupling) dualities of classical statistical mechanics (or (Euclidean) path integrals). Its range of applications includes discrete lattice, continuum field and gauge theories. Dualities are revealed to be local, structure-preserving mappings between model-specific bond algebras that can be implemented as unitary transformations, or partial isometries if gauge symmetries are involved. This characterization permits us to search systematically for dualities and self-dualities in quantum models of arbitrary system size, dimensionality and complexity, and any classical model admitting a transfer matrix or operator representation. In particular, special dualities such as exact dimensional reduction, emergent and gauge-reducing dualities that solve gauge constraints can be easily understood in terms of mappings of bond algebras. As a new example, we show that the ℤ2 Higgs model is dual to the extended toric code model in any number of dimensions. Non-local transformations such as dual variables and Jordan–Wigner dictionaries are algorithmically derived from the local mappings of bond algebras. This permits us to establish a precise connection between quantum dual and classical disorder variables. Our bond-algebraic approach goes beyond the standard approach to classical dualities, and could help resolve the long-standing problem of obtaining duality transformations for lattice non-Abelian models. As an illustration, we present new dualities in any spatial dimension for the quantum Heisenberg model. Finally, we discuss various applications including location of phase boundaries, spectral behavior and, notably, we show how bond-algebraic dualities help constrain and realize fermionization in an arbitrary number of spatial dimensions.

arXiv:1103.2776•Advances In Physics•Sep 29, 2011•135 citations

Mathematical structure of the three-dimensional(3D) Ising model

张志东

An overview of the mathematical structure of the three-dimensional(3D) Ising model is given from the points of view of topology,algebra,and geometry.By analyzing the relationships among transfer matrices of the 3D Ising model,Reidemeister moves in the knot theory,Yang-Baxter and tetrahedron equations,the following facts are illustrated for the 3D Ising model.1) The complex quaternion basis constructed for the 3D Ising model naturally represents the rotation in a(3+1)-dimensional space-time as a relativistic quantum statistical mechanics model,which is consistent with the 4-fold integrand of the partition function obtained by taking the time average.2) A unitary transformation with a matrix that is a spin representation in 2 n·l·o-space corresponds to a rotation in 2n·l·o-space,which serves to smooth all the crossings in the transfer matrices and contributes the non-trivial topological part of the partition function of the 3D Ising model.3) A tetrahedron relationship would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model,and its existence is guaranteed by the Jordan algebra and the Jordan-von Neumann-Wigner procedures.4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases φx,φy,and φz.The relationship with quantum field and gauge theories and the physical significance of the weight factors are discussed in detail.The conjectured exact solution is compared with numerical results,and the singularities at/near infinite temperature are inspected.The analyticity in β=1/(kBT) of both the hard-core and the Ising models has been proved only for β>0,not for β=0.Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.

中国物理B:英文版•Jan 1, 2013•0 citations
The exact analytical solution of the three-dimensional (3D) Ising model remains a central unresolved problem in statistical mechanics. While exact solutions exist in one and two dimensions, the three-dimensional case has resisted closed-form treatment. While standard Kramers-Wannier duality fails to yield a self-dual spin model in 3D, instead mapping to a gauge theory, we introduce a quasi-duality transformation by perturbing the isotropic cubic lattice with a slight, localized spatial inhomogeneity. Specifically, this is introduced via the boundary correction term A n ( ? )= δ η,nN ( s η ? s η?N+1 ? ? s η ? s η+1 ? ) , which arises from flattening the 2D planes into 1D row-continuous chains. By enforcing a condition where the introduced inhomogeneity algebraically absorbs the leading-order topological mismatch of the dual gauge fields, we recover an effective self-dual relation. This algebraic solution yields a critical coupling of K c = 1 4 ln( 1+ 2 )≈0.22033 , which remarkably deviates by only ~0.6% from the widely accepted homogeneous Monte Carlo estimation of K c ( MC ) ?0.22165 . We discuss the derivation and outline future applications of this method to study phase transitions in complex, structurally disordered systems. The method provides insight into how weak inhomogeneity restores a form of dual symmetry that approximates the true critical manifold. Implications for critical phenomena, renormalization structure, and complex biological systems operating near criticality are discussed.
Abstract The gauge invariance is an important phenomenon widely occurring in various matters with many-body interactions. The gauge theories are very important for understanding the evolution of fundamental physical processes in particle physics, high-energy physics, and condensed matter physics. In this work, the origin of nonlocal effects is inspected, and the contributions of nontrivial topological structures to physical properties are investigated in detail for both the three-dimensional (3D) Ising model and the 3D Z 2 lattice gauge model. Then, the exact solution for the 3D Z 2 lattice gauge theory is derived by the duality between the two models. The partition function, the spontaneous magnetization, the true range of correlation, and the critical point are rigorously derived for the 3D Z 2 lattice gauge theory. The critical exponents of the 3D Z 2 lattice gauge theory are in the same universality class as the 3D Ising model, which are α = 0, β = 3/8, γ = 5/4, δ = 13/3, η = 1/8, and ν = 2/3. Several fundamental issues, such as dimensionality, duality, symmetry, manifold and degenerate states, are investigated for these many-body interacting spin systems. The connections with superfluid, superconductors, particle physics, etc. are evaluated. Furthermore, physical significances and mathematical aspects of the 3D Z 2 lattice gauge theory are discussed with respect to topology, geometry, and algebra. The present results for the 3D Z 2 lattice gauge theory are helpful for figuring out the basic features of 3D lattice gauge theories with other symmetries, such as U(1), SU(2), and SU(3). This work shines a light on the interdisciplinarity between many-body interacting systems, algebra, topology, and geometry.
We report conjectures on the three-dimensional (3D) Ising model of simple orthorhombic lattices, together with details of calculations for a putative exact solution. Two conjectures, an additional rotation in the fourth curled-up dimension and weight factors on the eigenvectors, are proposed to serve as a boundary condition to deal with the topologic problem of the 3D Ising model. The partition function of the 3D simple orthorhombic Ising model is evaluated by spinor analysis, employing these conjectures. Based on the validity of the conjectures, the critical temperature of the simple orthorhombic Ising lattices could be determined by the relation of KK* = KK′ + KK″ + K′K″ or sinh 2K · sinh 2(K′ + K″ + (K′K″/K)) = 1. For a simple cubic Ising lattice, the critical point is putatively determined to locate exactly at the golden ratio , as derived from K* = 3K or sinh 2K · sinh 6K = 1. If the conjectures would be true, the specific heat of the simple orthorhombic Ising system would show a logarithmic singularity at the critical point of the phase transition. The spontaneous magnetization of the simple orthorhombic Ising ferromagnet is derived explicitly by the perturbation procedure, following the conjectures. The spin correlation functions are discussed on the terms of the Pfaffians, by defining the effective skew-symmetric matrix Aeff. The true range κx of the correlation and the susceptibility of the simple orthorhombic Ising system are determined by procedures similar to those used for the two-dimensional Ising system. The putative critical exponents derived explicitly for the simple orthorhombic Ising lattices are α = 0, β = 3/8, γ = 5/4, δ = 13/3, η = 1/8 and ν = 2/3, showing the universality behaviour and satisfying the scaling laws. The cooperative phenomena near the critical point are studied and the results based on the conjectures are compared with those of approximation methods and experimental findings. The putative solutions have been judged by several criteria. The deviations of the approximation results and the experimental data from the solutions are interpreted. Based on the solution, it is found that the 3D-to-2D crossover phenomenon differs with the 2D-to-1D crossover phenomenon and there is a gradual crossover of the exponents from the 3D to the 2D values. Special attention is also paid to the extra energy caused by the introduction of the fourth curled-up dimension and the states at/near infinite temperature revealed by the weight factors of the eigenvectors. The physics beyond the conjectures and the existence of the extra dimension are discussed. The present work is not only significant for statistical and condensed matter physics, but also fill the gap between the quantum field theory, cosmology theory, high-energy particle physics, graph theory and computer science.
The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems of square geometry with periodic boundaries oriented either along <10> or along <11> direction have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.
An overview of the mathematical structure of the three-dimensional (3D) Ising model is given, from the viewpoints of topologic, algebraic and geometric aspects. By analyzing the relations among transfer matrices of the 3D Ising model, Reidemeister moves in the knot theory, Yang-Baxter and tetrahedron equations, the following facts are illustrated for the 3D Ising model: 1) The complexified quaternion basis constructed for the 3D Ising model represents naturally the rotation in a (3 + 1) - dimensional space-time, as a relativistic quantum statistical mechanics model, which is consistent with the 4-fold integrand of the partition function by taking the time average. 2) A unitary transformation with a matrix being a spin representation in 2^(nlo)-space corresponds to a rotation in 2nlo-space, which serves to smooth all the crossings in the transfer matrices and contributes as the non-trivial topologic part of the partition function of the 3D Ising model. 3) A tetrahedron relation would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model, and its existence is guaranteed also by the Jordan algebra and the Jordan-von Neumann-Wigner procedures. 4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases Φx, Φy, and Φz. The relation with quantum field and gauge theories, physical significance of weight factors are discussed in details. The conjectured exact solution is compared with numerical results, and singularities at/near infinite temperature are inspected. The analyticity in β = 1/(kB T) of both the hard-core and Ising models has been proved for β > 0, not for β = 0. Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.
An algebraic theory of dualities is developed based on the notion of bond algebras. It deals with classical and quantum dualities in a unified fashion explaining the precise connection between quantum dualities and the low temperature (strong-coupling)/high temperature (weak-coupling) dualities of classical statistical mechanics (or (Euclidean) path integrals). Its range of applications includes discrete lattice, continuum field and gauge theories. Dualities are revealed to be local, structure-preserving mappings between model-specific bond algebras that can be implemented as unitary transformations, or partial isometries if gauge symmetries are involved. This characterization permits us to search systematically for dualities and self-dualities in quantum models of arbitrary system size, dimensionality and complexity, and any classical model admitting a transfer matrix or operator representation. In particular, special dualities such as exact dimensional reduction, emergent and gauge-reducing dualities that solve gauge constraints can be easily understood in terms of mappings of bond algebras. As a new example, we show that the ℤ2 Higgs model is dual to the extended toric code model in any number of dimensions. Non-local transformations such as dual variables and Jordan–Wigner dictionaries are algorithmically derived from the local mappings of bond algebras. This permits us to establish a precise connection between quantum dual and classical disorder variables. Our bond-algebraic approach goes beyond the standard approach to classical dualities, and could help resolve the long-standing problem of obtaining duality transformations for lattice non-Abelian models. As an illustration, we present new dualities in any spatial dimension for the quantum Heisenberg model. Finally, we discuss various applications including location of phase boundaries, spectral behavior and, notably, we show how bond-algebraic dualities help constrain and realize fermionization in an arbitrary number of spatial dimensions.
An overview of the mathematical structure of the three-dimensional(3D) Ising model is given from the points of view of topology,algebra,and geometry.By analyzing the relationships among transfer matrices of the 3D Ising model,Reidemeister moves in the knot theory,Yang-Baxter and tetrahedron equations,the following facts are illustrated for the 3D Ising model.1) The complex quaternion basis constructed for the 3D Ising model naturally represents the rotation in a(3+1)-dimensional space-time as a relativistic quantum statistical mechanics model,which is consistent with the 4-fold integrand of the partition function obtained by taking the time average.2) A unitary transformation with a matrix that is a spin representation in 2 n·l·o-space corresponds to a rotation in 2n·l·o-space,which serves to smooth all the crossings in the transfer matrices and contributes the non-trivial topological part of the partition function of the 3D Ising model.3) A tetrahedron relationship would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model,and its existence is guaranteed by the Jordan algebra and the Jordan-von Neumann-Wigner procedures.4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases φx,φy,and φz.The relationship with quantum field and gauge theories and the physical significance of the weight factors are discussed in detail.The conjectured exact solution is compared with numerical results,and the singularities at/near infinite temperature are inspected.The analyticity in β=1/(kBT) of both the hard-core and the Ising models has been proved only for β>0,not for β=0.Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.