The ferromagnetic nearest-neighbor Ising model on the simple-cubic lattice in zero field is the paradigmatic exactly-formulated but unsolved problem of equilibrium statistical mechanics: since Onsager's 1944 exact free energy and Yang's 1952 spontaneous magnetization for the two-dimensional model, no mathematically justified expression for the free energy per spin or the spontaneous magnetization has been found in three dimensions. The canonical form of the problem, and the field's acceptance standards for a solution, were fixed by the 2008 comment of Wu, McCoy, Fisher and Chayes - written to refute Z.-D. Zhang's 2007 conjectured solution: an exact solution must determine every auxiliary temperature-dependent function in a mathematically justified way, must reduce to the exact lower-dimensional results under specialization, must reproduce the established high- and low-temperature series term by term within their radii of convergence, must be compatible with rigorous results on analyticity and crossover, and must yield the exact critical temperature, critical exponents, and correlation-range behavior. No mathematically justified exact solution has been accepted to date. The problem remains open.
Background
The Ising model is the minimal statistical-mechanical model of a magnet: each site of a lattice carries a spin , and pairs of nearest-neighbor sites interact with energy ( favors alignment). In zero external field the Hamiltonian on the simple-cubic lattice (each site has 6 nearest neighbors) is . All equilibrium thermodynamics follows from the partition function , where the sum runs over all spin configurations. The two central objects are the free energy per spin in the thermodynamic limit , and the spontaneous magnetization (the order parameter, obtained by letting a small field go to zero after the infinite-volume limit), which is nonzero below the critical temperature and vanishes at and above it. Writing , the model undergoes a continuous phase transition at a critical coupling known numerically as .
Ising solved the one-dimensional chain in 1925 (no phase transition). Onsager (1944) obtained the exact free energy of the two-dimensional square-lattice model in zero field - the first exact solution of a model exhibiting a phase transition - and Yang (1952) computed its spontaneous magnetization, which vanishes as , that is, with critical exponent (critical exponents are the powers governing the singular behavior of observables near ). In three dimensions, after more than eighty years of effort, no exact solution is known: there is no mathematically justified expression for or of the simple-cubic nearest-neighbor model. The critical data are known only numerically - by Monte Carlo sampling (stochastic simulation of spin configurations) and by the conformal bootstrap (a numerical method that constrains the scaling dimensions of the three-dimensional critical theory from the crossing symmetry of correlation functions) - and the model's universality class (the set of systems sharing its critical exponents because long-wavelength fluctuations forget microscopic details) covers uniaxial magnets, simple fluids, and binary mixtures, making it the reference point for three-dimensional critical phenomena.
Because the model attracts conjectured "solutions" more often than theorems, the field maintains sharp consistency tests, stated canonically in the 2008 comment of Wu, McCoy, Fisher and Chayes. First, series tests: the free energy and magnetization possess convergent expansions at high temperature in and at low temperature in , known to many orders for the simple-cubic lattice (the high-temperature series of Guttmann and Enting; the low-temperature magnetization series from Essam and Fisher); rigorous results guarantee that, within finite radii of convergence, these series represent the thermodynamic-limit functions in dimensions , so any exact expression must reproduce them term by term - the "cardinal, golden rule" of the comment. Second, dimensional reduction: setting one coupling to zero must recover the exact two-dimensional results, including the exponent , with physically correct crossover (smooth interpolation between the three- and two-dimensional regimes as the coupling is switched off). Third, auxiliary functions: any temperature-dependent function entering a representation must be determined by proven equations, not fixed by ansatz or by fitting to the very series it must reproduce.
The 2007-2008 episode that prompted these standards: Z.-D. Zhang (2007) advanced a conjectured solution (formulated for the orthorhombic generalization with three distinct couplings, of which the simple cubic is the isotropic case): a four-fold integral for the free energy containing undetermined weight functions that were set to zero by an "ansatz" (an assumption introduced without proof), and a magnetization formula giving the critical exponent . Wu, McCoy, Fisher and Chayes disproved the conjecture: the setting is incompatible with the high-temperature expansion it had been fitted to; the magnetization formula reduces in two dimensions to the exponent instead of Yang's ; and its low-temperature expansion disagrees with the exact series. They concluded that the three-dimensional model "remains to this date an outstanding unsolved problem", and that the conjectured values of , the critical exponents, and the correlation-range properties are false.
Problem Statement
For the ferromagnetic nearest-neighbor Ising model on the simple-cubic lattice in zero external field, with Hamiltonian , , , and : produce a mathematically justified exact solution in the thermodynamic limit, consisting of (1) an explicit expression for the free energy per spin , and (2) an explicit expression for the spontaneous magnetization below the critical point, together with the mathematically justified solution machinery that determines the spin correlations, from which the exact critical coupling , the critical exponents (e.g. ), and the true asymptotic behavior of the correlation range (the characteristic decay length of spin correlations, which diverges at ) follow exactly.
The expressions may be closed forms or rigorously convergent integral or series representations, provided convergence to the thermodynamic limit is proven; a derivation, not a fit, is required. As codified in the source comment, any valid solution must satisfy all of the following: (i) every auxiliary temperature-dependent function or parameter entering the expressions is determined in a mathematically justified way - not by ansatz, physical plausibility, or fitting; (ii) under specialization (e.g. one coupling tending to zero) the expressions reduce to the known exact lower-dimensional results - the Onsager free energy and the Yang spontaneous magnetization with two-dimensional exponent - with physically correct crossover; and (iii) the expressions reproduce, term by term, the established high-temperature expansions in and low-temperature expansions in within their rigorous radii of convergence, and are compatible with rigorous results on analyticity and with the standard understanding of dimensional crossover.
References
- Ernst Ising, 'Beitrag zur Theorie des Ferromagnetismus', Zeitschrift fur Physik 31, 253 (1925). DOI: 10.1007/BF02980577; https://doi.org/10.1007/BF02980577
- Lars Onsager, 'Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition', Physical Review 65, 117 (1944). DOI: 10.1103/PhysRev.65.117; https://doi.org/10.1103/PhysRev.65.117
- C. N. Yang, 'The Spontaneous Magnetization of a Two-Dimensional Ising Model', Physical Review 85, 808 (1952). DOI: 10.1103/PhysRev.85.808; https://doi.org/10.1103/PhysRev.85.808
- John W. Essam and Michael E. Fisher, 'Padé Approximant Studies of the Lattice Gas and Ising Ferromagnet below the Critical Point', Journal of Chemical Physics 38, 802-812 (1963). DOI: 10.1063/1.1733766; https://doi.org/10.1063/1.1733766
- A. J. Guttmann and I. G. Enting, 'Series studies of the Potts model. I. The simple cubic Ising model', Journal of Physics A: Mathematical and General 26, 807-821 (1993). DOI: 10.1088/0305-4470/26/4/010; https://doi.org/10.1088/0305-4470/26/4/010
- Z.-D. Zhang, 'Conjectures on the exact solution of three-dimensional (3D) simple orthorhombic Ising lattices', Philosophical Magazine 87, 5309-5419 (2007). DOI: 10.1080/14786430701646325; arXiv:0705.1045; https://doi.org/10.1080/14786430701646325
- F. Y. Wu, B. M. McCoy, M. E. Fisher and L. Chayes, 'Comment on a recent conjectured solution of the three-dimensional Ising model', Philosophical Magazine 88, 3093-3095 (2008). DOI: 10.1080/14786430802537738; arXiv:0811.3876; https://doi.org/10.1080/14786430802537738
- Z. D. Zhang, 'Response to "Comment on a recent conjectured solution of the three-dimensional Ising model"', Philosophical Magazine 88, 3097-3101 (2008). DOI: 10.1080/14786430802537720; https://doi.org/10.1080/14786430802537720
- F. Y. Wu, B. M. McCoy, M. E. Fisher and L. Chayes, 'Rejoinder to the Response to "Comment on a recent conjectured solution of the three-dimensional Ising model"', Philosophical Magazine 88, 3103 (2008). DOI: 10.1080/14786430802537779; https://doi.org/10.1080/14786430802537779
- Jacques H. H. Perk, 'Comment on "Conjectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices"', Philosophical Magazine 89, 761-764 (2009). DOI: 10.1080/14786430902776970; https://doi.org/10.1080/14786430902776970
- Z. D. Zhang, 'Response to the Comment on "Conjectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices"', Philosophical Magazine 89, 765-768 (2009). DOI: 10.1080/14786430902776988; https://doi.org/10.1080/14786430902776988
- Jacques H. H. Perk, 'Rejoinder to the Response to the Comment on "Conjectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices"', Philosophical Magazine 89, 769-770 (2009). DOI: 10.1080/14786430902776962; https://doi.org/10.1080/14786430902776962
- Jacques H. H. Perk, 'Erroneous solution of three-dimensional (3D) simple orthorhombic Ising lattices', arXiv:1209.0731 (2012). https://arxiv.org/abs/1209.0731
- Zhi-Dong Zhang, 'Mathematical structure of the three-dimensional (3D) Ising model', Chinese Physics B 22, 030513 (2013). DOI: 10.1088/1674-1056/22/3/030513; https://doi.org/10.1088/1674-1056/22/3/030513
- Jacques H. H. Perk, 'Comment on "Mathematical structure of the three-dimensional (3D) Ising model"', Chinese Physics B 22, 080508 (2013). DOI: 10.1088/1674-1056/22/8/080508; https://doi.org/10.1088/1674-1056/22/8/080508
- Michael Aizenman, Hugo Duminil-Copin and Vladas Sidoravicius, 'Random Currents and Continuity of Ising Model's Spontaneous Magnetization', Communications in Mathematical Physics 334, 719-742 (2015). DOI: 10.1007/s00220-014-2093-y; arXiv:1311.1937; https://doi.org/10.1007/s00220-014-2093-y
- Michael E. Fisher and Jacques H. H. Perk, 'Comments concerning the Ising model and two letters by N.H. March', Physics Letters A 380, 1339-1340 (2016). DOI: 10.1016/j.physleta.2015.09.055; https://doi.org/10.1016/j.physleta.2015.09.055
- Filip Kos, David Poland, David Simmons-Duffin and Alessandro Vichi, 'Precision Islands in the Ising and O(N) Models', Journal of High Energy Physics 2016, 36 (2016). DOI: 10.1007/jhep08(2016)036; arXiv:1603.04436; https://doi.org/10.1007/jhep08(2016)036
- Zhidong Zhang, Osamu Suzuki and Norman H. March, 'Clifford algebra approach of 3D Ising model', Advances in Applied Clifford Algebras 29 (2018). DOI: 10.1007/s00006-018-0923-2; https://doi.org/10.1007/s00006-018-0923-2
- Alan M. Ferrenberg, Jiahao Xu and David P. Landau, 'Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model', Physical Review E 97, 043301 (2018). DOI: 10.1103/PhysRevE.97.043301; arXiv:1806.03558; https://doi.org/10.1103/PhysRevE.97.043301
- Degang Zhang, 'Exact Solution for Three-Dimensional Ising Model', Symmetry 13, 1837 (2021). DOI: 10.3390/sym13101837; https://doi.org/10.3390/sym13101837
- Jacques H. H. Perk, 'Comment on "Exact Solution for Three-Dimensional Ising Model" by Degang Zhang followed by Response to the Reply by Degang Zhang', arXiv:2202.03136 (2021). https://arxiv.org/abs/2202.03136
- Gandhimohan M. Viswanathan, Marco Aurelio G. Portillo, Ernesto P. Raposo and Marcos G. E. da Luz, 'What Does It Take to Solve the 3D Ising Model? Minimal Necessary Conditions for a Valid Solution', Entropy 24, 1665 (2022). DOI: 10.3390/e24111665; https://doi.org/10.3390/e24111665
- Tuncer Kaya, 'Relevant spontaneous magnetization relations for the triangular and the cubic lattice Ising model', Chinese Journal of Physics 77, 2676-2683 (2022). DOI: 10.1016/j.cjph.2022.03.043; https://doi.org/10.1016/j.cjph.2022.03.043
- Jacques H. H. Perk, 'Comment on Zhang, D. Exact Solution for Three-Dimensional Ising Model. Symmetry 2021, 13, 1837', Symmetry 15, 374 (2023). DOI: 10.3390/sym15020374; https://doi.org/10.3390/sym15020374
- Degang Zhang, 'Reply to Perk, J.H.H. Comment on "Zhang, D. Exact Solution for Three-Dimensional Ising Model. Symmetry 2021, 13, 1837"', Symmetry 15, 375 (2023). DOI: 10.3390/sym15020375; https://doi.org/10.3390/sym15020375
- Zhidong Zhang, 'Exact solution of the three-dimensional (3D) Z2 lattice gauge theory', Open Physics 23, 20250215 (2025). DOI: 10.1515/phys-2025-0215; arXiv:2412.10412; https://doi.org/10.1515/phys-2025-0215
- D. Shapoval, Yu. Honchar, B. Delamotte, M. Dudka and Yu. Holovatch, 'Seven- and eight-loop critical exponents of the three-dimensional Ising model', arXiv:2607.07776 (2026). https://arxiv.org/abs/2607.07776