Next Article in Journal
Non-Kolmogorovian Probabilities and Quantum Technologies
Next Article in Special Issue
Relevant Analytic Spontaneous Magnetization Relation for the Face-Centered-Cubic Ising Lattice
Previous Article in Journal
Iterative Joint Estimation Procedure of Channel and PDP for OFDM Systems
Previous Article in Special Issue
Generalized Solution of Inverse Problem for Ising Connection Matrix on d-Dimensional Hypercubic Lattice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

What Does It Take to Solve the 3D Ising Model? Minimal Necessary Conditions for a Valid Solution

by
Gandhimohan M. Viswanathan
1,2,*,
Marco Aurelio G. Portillo
1,3,
Ernesto P. Raposo
4 and
Marcos G. E. da Luz
3
1
Department of Physics, Federal University of Rio Grande do Norte, Natal 59078-970, RN, Brazil
2
National Institute of Science and Technology of Complex Systems, Federal University of Rio Grande do Norte, Natal 59078-970, RN, Brazil
3
Departmento de Física, Universidade Federal do Paraná, Curitiba 81531-980, PR, Brazil
4
Laboratório de Física Teórica e Computacional, Departamento de Física, Universidade Federal de Pernambuco, Recife 50670-901, PE, Brazil
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(11), 1665; https://doi.org/10.3390/e24111665
Submission received: 10 October 2022 / Revised: 3 November 2022 / Accepted: 8 November 2022 / Published: 15 November 2022
(This article belongs to the Special Issue Ising Model: Recent Developments and Exotic Applications II)

Abstract

:
An exact solution of the Ising model on the simple cubic lattice is one of the long-standing open problems in rigorous statistical mechanics. Indeed, it is generally believed that settling it would constitute a methodological breakthrough, fomenting great prospects for further application, similarly to what happened when Lars Onsager solved the two-dimensional model eighty years ago. Hence, there have been many attempts to find analytic expressions for the exact partition function Z, but all such attempts have failed due to unavoidable conceptual or mathematical obstructions. Given the importance of this simple yet paradigmatic model, here we set out clear-cut criteria for any claimed exact expression for Z to be minimally plausible. Specifically, we present six necessary—but not sufficient—conditions that Z must satisfy. These criteria will allow very quick plausibility checks of future claims. As illustrative examples, we discuss previous mistaken “solutions”, unveiling their shortcomings.

1. Introduction

Exact results are always welcome in science, even if they are simplified or idealized models of more realistic natural phenomena [1]. For example, an elucidating discussion about the general importance of analytical solutions in physics can be found in Ref. [2]. The Ising model of magnetism was originally proposed by Wilhelm Lenz in 1920 and exactly solved in one dimension by his graduate student Ernst Ising in 1924. A summary of the thesis was published in 1925 [3]. The Ising model long ago ceased to be a paradigm restricted only to magnetism models. Currently, it has found applications in diverse areas, from neuroscience to sea ice and voter models, to name a few [4]. Since Lars Onsager’s solution of the 2D Ising model in 1942 (which was published in 1944 [5]), statistical physics in general and equilibrium statistical mechanics in particular have experienced a great flourishing of powerful mathematical techniques [6,7], allowing considerable progress towards obtaining exact expressions for many relevant models [8,9,10].
In spite of all these advances, one of the most paradigmatic systems in statistical physics, namely, the ferromagnetic Ising model with interactions between nearest-neighbor two-state spins on the simple cubic lattice [11]—henceforth referred as the 3D Ising model—has withstood all attempts at exact solution. According to Rowlinson [12], the first claim of analytically cracking the 3D Ising model was presented at StatPhys 2, held in Paris in 1952, by John R. Maddox, who later became editor of Nature (see also Refs. [13,14,15]). Eliott W. Montroll, still during the conference, showed that the proposed expression could not be correct, by comparing it with the first few terms of the well established exact series expansions for the high and low temperature limit cases. It was identified later that the error was due to an incorrect application of the Jordan–Wigner transformation. Since then, new announcements of an exact solution have been made every few years, only to be systematically proved incorrect (e.g., Refs. [16,17,18,19]). What is seemingly lacking, in this historical context, is a set of clear-cut plausibility criteria that can be used to quickly verify whether or not a claimed solution is minimally worth considering. Therefore, we shall list a number of necessary—but not sufficient—conditions that a correct partition function Z must satisfy. Although some of these have already (in part or in full) appeared in the literature, here they are presented in an unified, rigorous, and comprehensive way. We emphasize that on the one hand these conditions can be used to refute unfounded claims—and thereby to identify any inappropriate protocols that are responsible for the incorrect features for Z. On the other hand, they might serve as a guide for the development of promising solution schemes. Finally, we briefly mention two reasons why the 3D Ising model is considered a major open problem, besides the obvious intellectually instigating fact that it still stands eighty years since Onsager solved the 2D version. First, despite the simple definitions of Ising-like models, their usefulness in studying a large number of diverse processes is overwhelming [11,20,21,22,23,24,25,26,27]. However, this should not be a surprise. Indeed, given the universality classes of classical spin models, of which the Ising is probably the most emblematic example, they all can be mapped to distinct instances of the logic problem of satisfiability (SAT) [28]. In this way, the analytic solution of the 3D case conceivably would represent a great boost for its already wide applicability. Second, the 3D Ising model near the critical point is very closely related to string and gauge field theories [29]. In particular, these associations can be analyzed through conformal bootstrap methods [30] (for a review, see, e.g., [31]). Hence, the eventual determination (or eventual disproof of existence) of an exact solution certainly will impact other important areas of physics.
All these potential perspectives involving a rigorous analytical expression for the partition function of the 3D Ising model certainly justify the establishment of a minimal set of clear-cut criteria for the validation of its exact solution, as presented below.

2. The 3D Ising Model

Suppose a finite simple cubic lattice Z N 3 of size L × L × L , with L a finite positive integer. There are N = L 3 sites that can be tagged as l = ( l 1 , l 2 , l 3 ) for i = 1 , 2 , 3 indicating the ith spatial direction (with unit vector e ^ i ) and l i = 1 , , L . Consider also two disjoint lattices, Z N 3 ˜ and Z N 3 , of sites labeled, respectively, by l b = ( n 1 , n 2 , n 3 ) , with 2 n i L 1 , and l f = ( m 1 , m 2 , m 3 ) , where at least one m i is necessarily 1 or L. Note that Z N 3 = Z N 3 ˜ Z N 3 , for Z N 3 ˜ ( Z N 3 ) representing the “bulk” (“frontier” or boundary) sites of Z N 3 . In other words, Z N 3 is the border or surface, whereas Z N 3 ˜ is the inside region or interior of the lattice.
The spin variable σ l at the vertex l Z N 3 can assume only two possible values, namely, −1 and +1. The Hamiltonian of the 3D anisotropic Ising model on Z N 3 and with zero external magnetic field is given by
H N ( σ ) = i = 1 3 J i l b l b = e ^ i σ l b σ l b + l l f = ε ( l ) e ^ i σ l σ l f , = H N ( b ) ( σ ) + H N ( f ) ( σ ) ,
where ε ( l ) = ± 1 if l Z N 3 ˜ and ε ( l ) = + 1 if l Z N 3 . The quantities J 1 , 2 , 3 are the couplings constants in the three distinct spatial directions, i = 1 , 2 , 3 . The two terms in the last equality in Equation (1) can readily be identified as the energy contributions from the system bulk b and frontier f regions for a state σ of the system. In fact, each σ represents a possible distinct configuration of −1 s and +1 s along the sites of the whole lattice and characterizes a specific system state.
Different boundary conditions (BCs) can be imposed to the problem [32,33]. They essentially specify constraints on the spin configuration of the set { σ } Z N 3 (see below).
The canonical partition function is conventionally defined as Z N ( K 1 , K 2 , K 3 ) = σ exp [ β H N ( σ ) ] , where β = ( k B T ) 1 and K i = β J i is the ith reduced temperature. The sum is over all the possible spin configurations { σ } over Z N 3 (observing the specified BCs). The partition function per site in the thermodynamic limit, the object of our interest here, is defined as
Z ( K 1 , K 2 , K 3 ) = lim N Z N ( K 1 , K 2 , K 3 ) 1 / N .
The challenge of finding an exact analytic expression for the above Z has been called the “holy grail of statistical mechanics” [34]. To this day, it remains one of the most important unsolved problems in statistical physics.
The above seemingly simple model has a notorious history of errors and controversies. In 1925, soon after arriving at the exact solution of the 1D ferromagnetic Ising model (with an external field), Ising himself erroneously extrapolated to 2D one of his 1D findings [3]. Specifically, he generalized incorrectly to 2D and 3D his (correct) result that there can be no spontaneous magnetization in 1D for positive T. The historical details can be found in Ref. [35]. Another controversy related to the Ising model concerns the critical exponents for the 2D model. Up until Onsager’s exact solution of the model, it was assumed that the critical exponents were given by the mean field approximation. Onsager’s solution allowed the mistake to be caught [36]. A few years later, when Maddox in 1952 presented his claim of a solution to the 3D model, he applied the Jordan–Wigner transformation to effectively try to map 2D planes of spin operators into fermionic creation and annihilation operators [15]. The corresponding trick works for the 2D Ising model because linear chains of spin operators can indeed be mapped into fermionic creation and annihilation operators. However, this trick does not easily generalize to 2D planes. There are 2D generalizations of the Jordan–Wigner transforation, but thus far none of these methods have been shown to be useful for solving the 3D Ising model. All claims of an exact solution for the 3D model have been shown to contain errors (see Section 4).

3. Necessary Conditions for a Valid Solution

Our main goal in this contribution is to establish a set of six necessary conditions that must be satisfied by any prospective of an 3D Ising model exact analytic expression for Z ( K 1 , K 2 , K 3 ) . Although the list consists of necessary conditions, these are not also sufficient conditions. In other words, any claimed Z may be wrong even if all necessary conditions are satisfied. However, if even just one condition is violated, the claimed expression for Z is certainly wrong.
Next, we present the noted conditions in a order somehow going from the more basic and fundamental to the more technical and abstract.

3.1. Condition 1

In the thermodynamic limit the per-site partition function Z of the 3D Ising model must be independent of the boundary conditions [15,37,38].
Surprisingly, the very broad reach of this requirement seems not to be properly appreciated as one should expect. Some recent erroneous claims of an exact Z even have assumed that certain special BCs could violate Condition 1 (see Section 4). Given such misunderstandings, below we present a very general and rigorous (although concise) proof that this indeed must be the case.
The thermodynamic limit represented by Equation (2) is known to be equivalently stated in terms of sequences of subgraphs G k of Z 3 (e.g., see [39,40]). There is a large relative freedom in choosing the structure of the successive G k s, provided they satisfy three fundamental properties, known as van Hove’s assumptions [38].
Denoting the number of vertices of a finite lattice G as V ( G ) , these assumptions are as follows:
  • k G k = Z 3 ,
  • G k G k + 1 ,
  • lim k V ( G k ) / V ( G k ) = 0 , for G k representing the frontier of G k , namely, G k = { l G k | j G k , | l j | = 1 } .
From the above, it is also possible to define G k ˜ = { l G k | l G k } , which is the bulk or core graph associated to G k . As a trivial example of a G k satisfying all the above characteristics, we note the aforementioned limited cubic lattice Z N 3 , where k = N = L 3 .
Thus, Z ( K 1 , K 2 , K 3 ) can be written more generally as
Z ( K 1 , K 2 , K 3 ) = c lim k Z G k ( K 1 , K 2 , K 3 ) 1 V ( G k ) ,
where
Z G k ( K 1 , K 2 , K 3 ) = σ { σ } G k exp [ β H G k ( σ ) ]
and H G k ( σ ) is the natural extension of Equation (1) to G k . Further, c is a constant of a purely topological origin. It may differ from 1 depending on the characteristics of the chosen sequence G k . However, it should not alter the resulting physics associated to the obtained partition function. Therefore, c might be set equal to 1 for the sake of discussion.
An important aspect of the finite Ising model relates to the boundary conditions assumed for the G k s. A rather general formulation for typical BCs relies on the following construction. Let Ω ( G ) represent all the possible combinations of the spin configurations on the vertices of the finite G, i.e., an element of Ω is denoted by σ and is a map σ : G { 1 , 1 } . Consider then G k and a specific subset Γ B C k Ω ( G k ) . We say that Γ B C k determines the BCs on the Ising model if the allowed spin configurations σ B C belong to
Ω B C ( G k ) = Ω ( G k ˜ ) × Γ B C k .
For instance, for Γ B C k = Ω ( G k ) , we have the usual free BCs, namely, for any site in the frontier the spin value can assume both values, −1 and +1, without restrictions. On the other hand, for G k displaying torus (or periodic), cylindrical, Klein, twisted, screw, etc., topology, then the permitted configurations in Γ B C k are established by proper pairwise mappings in the form σ l f σ l f . Hence, the finite partition function with the BCs determined by Γ B C k reads
Z G k B C ( K 1 , K 2 , K 3 ) = σ Ω B C ( G k ) exp [ β ( H G k ( b ) ( σ ) + H G k ( f ) ( σ ) ) ] .
We can now easily show that in the proper limit the partition function is independent on the BCs. We first observe that, for the free BCs,
lim k ln [ Z G k free ] V ( G k )
exists (for a proof, see Ref. [38]). Let us denote the limit in Equation (5) as ln [ Z ( K 1 , K 2 , K 3 ) ] . Second, we determine a bound for H G k ( f ) ( σ ) in Equation (4). For any site l f in G k , the maximum number of l G k such that | l l f | = 1 is 5 for the simple cubic lattice, then we have that | H G k ( f ) | < 5 V ( G k ) . Third, by lemma 2.2.1 in [37]
| ln [ Z G k B C ] ln [ Z G k free ] | β | H G k ( f ) B C H G k ( f ) free | .
Finally, by dividing Equation (6) by V ( G k ) , considering the triangular inequality
| H G k ( f ) B C H G k ( f ) free | 2 × β 5 V ( G k )
and since the sequence G k is van Hove, then
lim k ln [ Z G k B C ] V ( G k ) = lim k ln [ Z G k free ] V ( G k ) .
Thus, we readily conclude that
lim k Z G k B C ( K 1 , K 2 , K 3 ) 1 / V ( G k ) = Z ( K 1 , K 2 , K 3 )
is well defined and independent on the BCs.

3.2. Condition 2

The claimed partition function per site of the 3D model Z ( K 1 , K 2 , K 3 ) must reduce to Onsager’s solution whenever one of the three reduced temperatures vanishes.
Indeed, suppose without loss of generality that J 3 = 0 , so that K 3 = 0 . Then, for l l indicating the sites which lie in the plane x 3 = l (whose set of spin configurations we represent by { σ l } ), the 3D Hamiltonian of size N = L × L × L = L 2 × L can be written as ( σ j 0 if j G k )
H L 2 × L ( σ ) = l = 1 L l l i = 1 2 J i σ l l σ l l + e ^ i = l = 1 L H L 2 ( σ l ) ,
with H L 2 ( σ l ) the energy associated to the plane x 3 = l for σ l a given distribution of spins in such plane. Note that the 3D Hamiltonian is now expressed as the sum of L independent and identical 2D Hamiltonians. In this case, the 3D partition function Z N = Z L 3 ( K 1 , K 2 , 0 ) = Z L 2 × L factors as
Z L 2 × L = σ 1 σ L l = 1 L exp [ β H L 2 ( σ l ) ] = ( Z L 2 ) L ,
where Z L 2 = Z L 2 ( K 1 , K 2 ) is the finite 2D partition function, whose per-site thermodynamic limit Z ( K 1 , K 2 ) is naturally defined as
Z ( K 1 , K 2 ) = lim L Z L 2 ( K 1 , K 2 ) 1 / L 2 .
Thus, from Equation (10), we get
Z ( K 1 , K 2 , 0 ) = lim L Z L 2 × L 1 / L 3 = lim L ( Z L 2 L ) 1 / L 3 = lim L Z L 2 1 / L 2 = Z ( K 1 , K 2 ) .
Thus, the 3D per-site partition function must reduce to Onsager’s solution when any of the three couplings is made to vanish.

3.3. Condition 3

For the isotropic case, namely K = K i ( i = 1 , 2 , 3 ), any claimed per-site partition function must be analytic for 0 K < K c [41] and for K c < K < [42], where K c = 0.221 654 626 is the well known numerically estimated value of the critical temperature of the 3D Ising model.
Of course, when there are interactions besides nearest-neighbor or when there is an applied magnetic field, the range of analyticity in K for the partition function can be distinct from the above one (see, for instance, Refs. [38,41,42]).
Note also that an exact expression for Z should lead to an exact formula for K c . Thence, if one has derived a tentative exact Z, its analyticity should be relatively easy to test, e.g., for complex functions from the Cauchy–Riemann equations [43], and for real functions using standard techniques, as those described in [44]. We observe that the numerical value of the critical temperature is known with very high precision from Monte Carlo simulations and other numerical approaches [45,46,47,48].

3.4. Condition 4

For the isotropic case, the series expansion in the high (low) temperature limit—K small (large)—of the claimed solution must perfectly match the already established series to all known orders (see, e.g., Refs. [38,41]).
This is a direct consequence of the fact that in the domain of an analytic function, its series expansion around a fixed expansion point must be unique.
Let v i = tanh [ K i ] be the high temperature expansion variable. For the isotropic case, i.e., v = v i ( i = 1 , 2 , 3 ) , we define
Z ( K ) = 2 ( cosh [ K ] ) 3 Z high ( v ) .
For details, see, for instance, Equation (9) in Ref. [49]. This high temperature formulation is obtained by performing the sum over all states, and then using (for σ i = ± 1 )
exp [ K σ i σ j ] = cosh [ K ] + σ i σ j sinh [ K ] = cosh [ K ] ( 1 + σ i σ j tanh [ K ] ) .
This identity justifies the choice of the high temperature variable v = tanh [ K ] . It is well known (see Ref. [39] for the 2D analog) that Z high ( v ) is then given by a properly defined generating function for graphs of given length, all of whose nodes have even degrees and whose edges only connect nearest neighbors, in the limit of an infinite lattice. Such graphs have been variously called “closed graphs” or “multipolygons.”
The first few expansion terms of Z high ( v ) have been rigorously determined [38] via finite lattice methods [50]. From Condition 3, Z ( K ) is analytic for T > T c . Therefore, for large T, the function Z high ( v ) is also analytic. In this way, the series expansion of Z ( K ) / ( 2 ( cosh [ K ] ) 3 ) must coincide with the noted known terms.
Similarly, by setting u i = exp [ 4 K i ] as the low temperature variable, for the isotropic case we have u i = u ( i = 1 , 2 , 3 ). Then, one can write the partition function as
Z ( K ) = u 3 / 4 Z low ( u ) .
See, for instance, Equation (1.5) in Ref. [40]. This low temperature formulation simply ennumerates the number of configurations with a given domain wall area. Indeed, the quantity Z low is the generating function, in the limit of an infinite lattice, of closed surfaces of given area (see, e.g., Ref. [40]). The prefactor u 3 / 4 is due to the ground state energy not being zero. If each of the three bonds per site in the ground state contribute an energy β J = K , then the partition function per site for the ground state is exp [ 3 K ] = u 3 / 4 .
The first few exact terms of the low temperature expansion of Z low ( u ) have been calculated [50]. Therefore, similarly to the high temperature expansion, for T < T c low enough, any claimed exact solution should meet term-by-term the noted series.
The first few expansion terms of both Z high ( v ) and Z low ( u ) are known, thanks to the finite lattice method [50]. Ian Enting and Tom de Neef originally developed this innovative approach in the 1970s for generating series expansions, with applications to exact enumeration problems [51]. Since then, this powerful numerical method has led to significant advances. For example, Iwan Jensen used it to calculate the statistics of self-avoiding polygons on the square lattice [52]. Very recently, Nathan Clisby wrote an expository article concerning the method’s relevance to the enumerative combinatorics of lattice polymers [53]. It is arguably the most powerful algorithmic technique for obtaining exact series expansions to high order of models whose exact solution is not known, including, of course, the 3D Ising model.
Table 1 lists the first few non-zero terms of the high and low temperature expansions. Therefore, according to the present condition, any claimed Z that does not exactly comply with the terms in Table 1 cannot be the exact partition function for the 3D Ising model.

3.4.1. Condition 5

The claimed exact Z should display permutation symmetry and convexity on the reduced temperature variables K i = β J i .
The first property, namely permutation symmetry, as proposed in [15], implies that Z ( K 1 , K 2 , K 3 ) = Z ( K π ( 1 ) , K π ( 2 ) , K π ( 3 ) ) , where π ( i ) is a permutation for i = 1 , 2 , 3 . Since the final Z does not depend on the BCs or on the specific choice of the van Hove sequence, consider G k in Equation (3) as the simple cubic lattice Z N 3 with free BCs. Then, the proof is direct since the finite Z N ( K 1 , K 2 , K 3 ) trivially presents the aforementioned symmetry. Further (refer to Ref. [48]), if one can use Z to derive an equation for the critical β c , i.e., F ( β c J 1 , β c J 2 , β c J 3 ) = 0 , then also β c must be invariant under any permutation of J 1 , J 2 , and J 3 .
In Ref. [38], it has been shown that, for the isotropic case J i = J ( i = 1 , 2 , 3 ) , the partition function must be a convex function in β . For the anisotropic case, the proof follows exactly the same steps. Indeed, considering again G k the finite cubic lattice with free BCs, one finds that G k is convex (details omitted here). However, because the limit of a van Hove sequence of convex functions is also convex [54], the final result holds. As a consequence, one has the following. For α [ 0 , 1 ] , β a , β b ( 0 , ) and arbitrary J 1 , J 2 , J 3 , then for β d = α β a + ( 1 α ) β b
Z ( β d J 1 , β d J 2 , β d J 3 ) α Z ( β a J 1 , β a J 2 , β a J 3 ) + ( 1 α ) Z ( β b J 1 , β b J 2 , β b J 3 ) .

3.4.2. Condition 6

The claimed exact solution must bring clarity to the conundrum related to the behavior of the partially resummed high temperature expansion of the anisotropic partition function [49,55].
Based on the anomalous behavior of the partially resummed series solution of the 3D Ising model, it is now believed by many (but not all) that the 3D Ising model might “not be solvable” [55]. Here, by not solvable we mean that Z may not be a differentiably finite function (known in the mathematical jargon as a D-finite or holonomic function). Recall that F ( z 1 , z 2 , , z N ) is said holonomic if F is analytic in all the variables and satisfies a linear differential equation, whose coefficients are polynomials in z 1 , z 2 , , z N (see, e.g., [56]). Hence, any claimed exact holonomic Z must be able to explain how and why the anisotropic 3D Ising model has a high-temperature series that, upon partial resummation, seems to indicate non-D-finiteness.

4. Some Previously Claimed Solutions in the Literature

We briefly review how previous advanced solutions have failed to satisfy the above conditions, thus not representing the correct exact Z ( K 1 , K 2 , K 3 ) .
As already noted, the proposal by Maddox in 1952 violated Condition 4 for the series expansions, as did the ones by Das [16], Lou and Wu [17], and Z.-D. Zhang [18]. The serious errors in the latter were also extensively addressed in Refs. [15,57]. Moreover, except for the solution proposed by Zhang, all others also violate Condition 3. In fact, Zhang’s solution only seems to satisfy it because the numerical value of the critical temperature is imposed as an ansatz, built into his construction. Still, Zhang’s critical temperature of K c 0.221 658 637 208 698 [18] (taken from a conjecture of Rosengren [58]) differs from the best known numerical estimate K c 0.221 654 626. It may seem that the discrepancy is small, but an exact solution should give the exact critical temperature to arbitrary precision, i.e., to all decimal places. For the exact critical temperature, it is not acceptable to tolerate a discrepancy in the 10th or even 10 10 th decimal place.
Crucially, none of the above claimed solutions minimally attend Condition 6. As explained above, the resummed high temperature series of the anisotropic 3D model seems to show an anomalous behavior, strongly suggesting non-D-finiteness. However, the claimed solutions all behave normally under resummation of the anisotropic high temperature series.
In 2021, D. Zhang [19] (not to be confused with Z.-D. Zhang) made another claim, promptly criticized in [48]. It is easy to check that the assertion in [19] violates Conditions 1 and 4. The claimed solution also fails to bring new insights regarding Condition 6. For Condition 3, the predicted critical temperature disagrees drastically with the known numerically estimated value. Finally, regarding Condition 2, Zhang writes (inaccurately) that:
When the interaction energy in the third dimension vanishes, Onsager’s exact solution of the 2D Ising model is recovered immediately. This guarantees the correctness of the exact solution of the 3D Ising model [emphasis added].
In fact, there is no such guarantee. Condition 2 is a necessary but not a sufficient condition for a solution to be correct. For example, the expression (for v i as defined before)
ln [ Z ] = ln [ 2 cosh [ K 1 ] cosh [ K 2 ] cosh [ K 3 ] ] + 1 2 ( 2 π ) 3 π π d k 1 d k 2 d k 3 × ln [ ( 1 + v 1 2 ) ( 1 + v 2 2 ) ( 1 + v 3 2 ) 2 v 1 ( 1 v 2 2 ) ( 1 v 3 2 ) cos [ k 1 ] 2 v 2 ( 1 v 1 2 ) ( 1 v 3 2 ) cos [ k 2 ] 2 v 3 ( 1 v 2 2 ) ( 1 v 1 2 ) cos [ k 3 ] ] ,
correctly reduces to Onsager’s solution of the 2D model if any one of the three K i are made to vanish. However, this expression clearly is not correct because it violates Conditions 1, 3, 4, and 6 above. See also the famous (but wrong) conjecture of Mark Kac, e.g., in Ref. [39].

5. Final Remarks and Conclusions

The last few decades have witnessed significant developments [59,60,61] aimed at obtaining an exact expression for the 3D Ising model. In the absence of strong theoretical results pointing otherwise, such steady progress should dispel the false myth regarding the (non)solvability of the 3D model (see below).
First, we emphasize that the ferromagnetic 3D Ising model with nearest neighbor interactions is not an NP-complete problem. Recall that a problem is said to be NP-complete when it can be used to simulate any problem classified as NP. A problem belongs to the NP complexity class when a proposed solution can be checked in polynomial time by a deterministic Turing machine. It is an open question to determine if NP problems can also be not just checked but also solved in polynomial time by a deterministic Turing machine (P complexity class). It is true that there is a theorem concerning NP-completeness of the Ising model due to Sorin Istrail [59]. Nonetheless, it refers to the 3D Ising spin glass with arbitrary interactions, not to the ferromagnetic model. Moreover, the problem addressed in [59] relates to finding the ground state. For the ferromagnetic case, the ground state is trivial, viz., with all spins aligned (so doubly degenerate).
Second, although there is strong numerical evidence of non-solvability of the 3D Ising model in terms of D-finite functions (see Condition 6), mathematical proofs for this supposition are still lacking. If indeed this would be the case, an exact analytic solution based on nonholonomic functions could still be possible. Actually, many experts have been careful to make clear that the aforementioned non-solvability of the 3D Ising model is conjectural. Barry Cipra, writing in Science [34], has stated that
It might still be possible to find exact answers for some special cases of the Ising model, Istrail notes. In particular, the ferromagnetic case of the 3D Ising model may turn out to be simple enough to solve.
And third, it is also not true that the progress is too slow or that the problem is hopelessly too difficult. Nor is it a waste of time—quite the contrary. In the preface to Polygons, Polyominoes and Polycubes, Anthony J. Guttmann writes [62],
This is indeed a golden age for studying such problems. With powerful computers and new algorithms, unimaginable numerical precision in our estimates of properties of many of these models is now possible. On the mathematical side, we are developing tools for solving increasingly complex functional equations, while the theory of conformal invariance, and the developments around stochastic Löwner evolution have given us powerful tools to predict, and in some cases to prove, new results. The scientific community in this field is divided into those who think we will never solve the problem, of say the perimeter or area generating function of self-avoiding polygons in two dimensions, and those who think that we will. I am firmly in the latter camp …
Finally, we briefly discuss other three dimensional lattice systems. Condition 1 is general and valid for all lattice systems, so long as the interactions in the Hamiltonian are nearest-neighbor. In contrast, conditions 2, 4, and 5 as formulated are specific to the simple cubic lattice. However, we can expect that there should be analogs of these conditions for each lattice system. The same should be true for Condition 3 and the numerical value of K c . Condition 6, however, is the most difficult to generalize. Very little is known about the analog of Condition 6 for other lattice systems.
Summarizing, we have reviewed, systematized, and enlarged a set of necessary conditions characterizing a potentially exact Z for the 3D Ising model. We have arranged this set into a single framework. Obviously, this set does not per se establish a concrete protocol that can solve the Ising system. Nevertheless, if even a single criterion is violated, one can be 100% certain that the methodology followed is fatally flawed. In this sense, the advance reported here has the potential to guide the maturing of future attempts to obtain the true Z.
We emphasize that the discussion about misguided attempts in the literature presented here by no means has the intention of criticizing these authors. Our purpose is solely to illustrate the subtleties and intricacies of the problem, which has deceived even some of the most respectable researchers. Our discussion thus makes clear the real necessity of clear-cut tests to check the plausibility of future claims.
Finally, we note an eventual (although improbable) curious consequence of our results. An exact analytic expression for Z should observe all the previously addressed requirements. However, it could be the case that such a function—observing the full set of conditions 1–6—cannot exist. A proof, of course, would settle negatively the possibility of an analytical Z. However, we conjecture that the six conditions are not mutually inconsistent.

Author Contributions

All authors contributed equally overall. Conceptualization, G.M.V. and M.G.E.d.L.; formal analysis, G.M.V., M.A.G.P., E.P.R., M.G.E.d.L.; investigation, G.M.V., M.A.G.P., E.P.R., M.G.E.d.L.; writing—original draft preparation, G.M.V., M.A.G.P., E.P.R., M.G.E.d.L. All authors have read and agreed to the published version of the manuscript.

Funding

We thank CNPq (grants 304532/2019-3, 305062/2017-4 and 302051/2018-0), CAPES, and FACEPE for funding. MAGP thanks the Brazilian National Institute of Science and Technology of Complex Systems for funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Jacques H. H. Perk for feedback and for very helpful discussions about analyticity at K = 0 . G.M.V. thanks V.M. Kenkre for discussions about boundary conditions during his stay at the Consortium of the Americas for Interdisciplinary Science, at the University of New Mexico, in 2008–2009.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gershenfeld, N. The Nature of Mathematical Modeling; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
  2. Borwein, J.M.; Crandall, R.E. Closed forms: What they are and why we care. Not. Am. Math. Soc. 2013, 60, 50–65. [Google Scholar] [CrossRef] [Green Version]
  3. Ising, E. Beitrag zur Theorie des Ferromagnetismus. Z. Physik 1925, 31, 253. [Google Scholar] [CrossRef]
  4. Tranquillo, J. An Introduction to Complex Systems: Making Sense of a Changing World; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  5. Onsager, L. Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition. Phys. Rev. 1944, 65, 117. [Google Scholar] [CrossRef]
  6. Kotecky, R. (Ed.) Methods of Contemporary Mathematical Statistical Physics; Springer: Berlin, Germany, 2009. [Google Scholar]
  7. Thompson, C.J. Mathematical Statistical Mechanics; Princeton University Press: Princeton, NJ, USA, 2015. [Google Scholar]
  8. Sutherland, B. Beautiful Models: 70 Years of Exactly Solved Quantum Many-Body Problems; World Scientific: Singapore, 2005. [Google Scholar]
  9. Mattis, D.C. The Many-Body Problem: An Encyclopedia of Exactly Solved Models in One Dimension; World Scientific: Singapore, 1993. [Google Scholar]
  10. Lieb, E.H.; Mattis, D.C. Mathematical Physics in One Dimension; Academic Press Inc.: New York, NY, USA, 1966. [Google Scholar]
  11. Singh, S.P. The Ising Model: Brief Introduction and Its Application. In Solid State Physics—Metastable, Spintronics Materials and Mechanics of Deformable Bodies—Recent Progress; Sivasankaran, S., Nayak, P.K., Günay, E., Eds.; IntechOpen: London, UK, 2020. [Google Scholar] [CrossRef] [Green Version]
  12. Van der Waals, J.D. On the Continuity of the Gaseous and Liquid States; Rowlinson, J.S., Ed.; Dover: New York, NY, USA, 2004. [Google Scholar]
  13. Maddox, J.R. Talk Presented at StatPhys 2, Paris, France. 1952. [Google Scholar]
  14. Maddox, J.R. Changement de Phases Societé de Chimie Physique; Presses Universitaires de France: Paris, France, 1952. [Google Scholar]
  15. Perk, J.H.H. Comment on ‘Conjectures on exact solution of three-dimensional (3D) simple orthorhombic Ising lattices’. Philos. Mag. 2009, 89, 761–764. [Google Scholar] [CrossRef] [Green Version]
  16. Das, D.D. Exact Partition Function of Ising Model in Magnetism in One, Two and Three Dimensions in Nonzero Field. Indian J. Phys. 1970, 44, 244. [Google Scholar]
  17. Lou, S.L.; Wu, S.H. Three-dimensional Ising model and transfer matrices. Chin. J. Phys. 2000, 38, 841–854. [Google Scholar]
  18. Zhang, Z.-D. Conjectures on the exact solution of three-dimensional (3D) simple orthorhombic Ising lattices. Philos. Mag. 2007, 87, 5309–5419. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, D. Exact solution for three-Dimensional Ising model. Symmetry 2021, 13, 1837. [Google Scholar] [CrossRef]
  20. Majewski, J.; Li, H.; Ott, J. The Ising model in physics and statistical genetics. Am. J. Hum. Gen. 2001, 69, 853–862. [Google Scholar] [CrossRef] [Green Version]
  21. Binek, C. Ising-Type Antiferromagnets; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  22. Suzuki, S.; Inoue, J.-I.; Chakrabarti, B.K. Quantum Ising Phases and Transitions in Transverse Ising Models, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  23. McCoy, B.M.; Wu, T.T. The Two-Dimensional Ising Model, 2nd ed.; Dover: New York, NY, USA, 2014. [Google Scholar]
  24. Adler, M. Monte Carlo Simulations of the Ising Model; Anchor Academic Publishing: Hamburg, Germany, 2016. [Google Scholar]
  25. Fadil, Z. Semi-Infinite Ising Model by the Renormalization Group: Applicable in Nanotechnology and Spintronics; Lap Lambert Academic Publishing: Beau Bassin, Mauritius, 2020. [Google Scholar]
  26. Farah, A. The applications of the Ising model in statistical thermodynamics and quantum mechanics. Eur. Acad. Res. 2020, 8, 2229–2237. [Google Scholar]
  27. Lipowski, A. (Ed.) Special Issue “Ising Model: Recent Developments and Exotic Applications”. Entropy 2022. Available online: https://www.mdpi.com/journal/entropy/special_issues/Ising_Model (accessed on 9 October 2022).
  28. Las Cuevas, G.D.; Cubitt, T.S. Simple universal models capture all classical spin physics. Science 2016, 351, 1180–1183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Polyakov, A.M. Gauge Fields and Strings; Harwood Academic Publishers: London, UK, 1987. [Google Scholar]
  30. Showk, S.-E.; Paulos, M.F.; Poland, D.; Rychkov, S.; Duffin, D.-S.; Vichi, A. Solving the 3D Ising model with the conformal bootstrap. Phys. Rev. D 2012, 86, 025022. [Google Scholar] [CrossRef]
  31. Rychkov, S. 3D Ising model: A view from the conformal bootstrap island. Comptes Rendus. Phys. 2020, 21, 185–198. [Google Scholar] [CrossRef]
  32. Wu, M.-C.; Hu, C.-K. Exact partition functions of the Ising model on M × N planar lattices with periodic—Aperiodic boundary conditions. J. Phys. A 2002, 35, 5189–5206. [Google Scholar] [CrossRef]
  33. Izmailian, N.S.; Hu, C.-K. Finite-size effects for the Ising model on helical tori. Phys. Rev. E 2007, 76, 041118. [Google Scholar] [CrossRef] [Green Version]
  34. Cipra, B. Statistical physicists phase out a dream. Science 2000, 288, 1561–1562. [Google Scholar] [CrossRef] [Green Version]
  35. Brush, S.G. History of the Lenz-Ising model. Rev. Mod. Phys. 1967, 39, 883–893. [Google Scholar] [CrossRef]
  36. Longuet-Higgins, H.C.; Fisher, M.E. Lars Onsager. 27 November 1903–5 October 1976. Biogr. Mem. Fellows R. Soc. 1978, 24, 443–471. [Google Scholar]
  37. Ruelle, D. Statistical Mechanics: Rigorous Results; World Scientific: Singapore, 1999. [Google Scholar]
  38. Friedli, S.; Velenik, Y. Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  39. Feynman, R.P. Statistical Mechanics: A Set of Lectures; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  40. Gattringer, C.R.; Jaimungal, S.; Semenoff, G.W. Loops, surfaces and Grassmann representation in two- and three-dimensional Ising models. Int. J. Mod. Phys. A 1999, 14, 4549–4574. [Google Scholar] [CrossRef] [Green Version]
  41. Ott, S. Weak mixing and analyticity of the pressure in the Ising model. Commun. Math. Phys. 2020, 377, 675–696. [Google Scholar] [CrossRef] [Green Version]
  42. Ott, S. Weak mixing and analyticity in Random Cluster and low temperature Ising models. arXiv 2020, arXiv:2003.05879. [Google Scholar]
  43. Marsden, J.E.; Hoffman, M.J. Basic Complex Analysis, 3rd ed.; W. H. Freeman: New York, NY, USA, 1998. [Google Scholar]
  44. Krantz, S.G.; Parks, H.R. A Primer of Real Analytic Functions, 2nd ed.; Birkhäuser: Boston, MA, USA, 2002. [Google Scholar]
  45. Talapov, A.L.; Blöte, H.W.J. The magnetization of the 3D Ising model. J. Phys. A 1996, 29, 5727–5733. [Google Scholar] [CrossRef] [Green Version]
  46. Blöte, H.W.J.; Heringa, J.R.; Hoogland, A.; Meyer, E.W.; Smit, T.S. Monte Carlo renormalization of the 3D Ising model: Analyticity and convergence. Phys. Rev. Lett. 1996, 76, 2613–2616. [Google Scholar] [CrossRef] [Green Version]
  47. Ferrenberg, A.M.; Xu, J.; Landau, D.P. Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model. Phys. Rev. E 2018, 97, 043301. [Google Scholar] [CrossRef] [PubMed]
  48. Perk, J.H.H. Comment on “Exact Solution for Three-Dimensional Ising model” by Degang Zhang. arXiv 2021, arXiv:2202.03136. [Google Scholar]
  49. Hansel, D.; Maillard, J.M.; Oitmaa, J.; Velgakis, M.J. Analytical properties of the anisotropic cubic Ising model. J. Stat. Phys. 1987, 48, 69–80. [Google Scholar] [CrossRef]
  50. Guttmann, A.J.; Enting, I.G. Series studies of the Potts model. I: The simple cubic Ising model. J. Phys. A 1993, 26, 807–822. [Google Scholar] [CrossRef] [Green Version]
  51. De Neef, T.; Enting, I.G. Series expansions from the finite lattice method. J. Phys. A 1977, 10, 801–806. [Google Scholar] [CrossRef]
  52. Jensen, I. Size and area of square lattice polygons. J. Phys. A 2000, 33, 3533–3544. [Google Scholar] [CrossRef] [Green Version]
  53. Clisby, N. Enumerative combinatorics of lattice polymers. Not. Am. Math. Soc. 2021, 68, 504–515. [Google Scholar] [CrossRef]
  54. Niculescu, C.P.; Persson, L.-E. Convex Functions and Their Applications, 2nd ed.; Springer: Cham, Switzerland, 2018. [Google Scholar]
  55. Guttmann, A.J.; Enting, I.G. Solvability of some statistical mechanical systems. Phys. Rev. Lett. 1996, 76, 344. [Google Scholar] [CrossRef] [PubMed]
  56. Van Der Hoven, J. Fast evaluation of holonomic functions near and in regular singularities. J. Symb. Comput. 2001, 31, 717–743. [Google Scholar] [CrossRef] [Green Version]
  57. Wu, F.Y.; McCoy, B.M.; Fisher, M.E.; Chayes, L. Comment on a recent conjectured solution of the three-dimensional Ising model. Phil. Mag. 2008, 88, 3093–3095. [Google Scholar] [CrossRef]
  58. Rosengren, A. On the combinatorial solution of the Ising model. J. Phys. A. 1986, 19, 1709–1714. [Google Scholar] [CrossRef]
  59. Istrail, S. Statistical mechanics, three-dimensionality and NP-completeness: I. Universality of intractability for the partition function of the Ising model across non-planar lattices. In STOC ’00: Proceedings of the Thirty-Second Annual ACM Symposium on Theory of Computing; Yao, F., Luks, E., Eds.; ACM Press: New York, NY, USA, 2000; pp. 87–96. [Google Scholar]
  60. Regge, T.; Zecchina, R. Combinatorial and topological approach to the 3D Ising model. J. Phys. A 2000, 33, 741–761. [Google Scholar] [CrossRef]
  61. Cimasoni, D. A generalized Kac-Ward formula. J. Stat. Mech. 2010, 2010, P07023. [Google Scholar] [CrossRef]
  62. Guttman, A.J. (Ed.) Polygons, Polyominoes and Polycubes; Springer: Dordrecht, The Netherlands, 2009. [Google Scholar]
Table 1. The first ten nonzero coefficients in the high temperature expansion of Z high ( v ) = n = 0 a n v n and in the low temperature expansion of Z low ( u ) = n = 0 b n u n (see main text). The coefficients were obtained by Guttmann and Enting using the finite lattice method [50].
Table 1. The first ten nonzero coefficients in the high temperature expansion of Z high ( v ) = n = 0 a n v n and in the low temperature expansion of Z low ( u ) = n = 0 b n u n (see main text). The coefficients were obtained by Guttmann and Enting using the finite lattice method [50].
n (High) a n n (Low) b n
0101
4331
62253
81926−3
102046715
1224,8538−30
14329,3349101
164,649,60110−261
1868,884,35611807
201,059,830,11212−2308
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Viswanathan, G.M.; Portillo, M.A.G.; Raposo, E.P.; da Luz, M.G.E. What Does It Take to Solve the 3D Ising Model? Minimal Necessary Conditions for a Valid Solution. Entropy 2022, 24, 1665. https://doi.org/10.3390/e24111665

AMA Style

Viswanathan GM, Portillo MAG, Raposo EP, da Luz MGE. What Does It Take to Solve the 3D Ising Model? Minimal Necessary Conditions for a Valid Solution. Entropy. 2022; 24(11):1665. https://doi.org/10.3390/e24111665

Chicago/Turabian Style

Viswanathan, Gandhimohan M., Marco Aurelio G. Portillo, Ernesto P. Raposo, and Marcos G. E. da Luz. 2022. "What Does It Take to Solve the 3D Ising Model? Minimal Necessary Conditions for a Valid Solution" Entropy 24, no. 11: 1665. https://doi.org/10.3390/e24111665

APA Style

Viswanathan, G. M., Portillo, M. A. G., Raposo, E. P., & da Luz, M. G. E. (2022). What Does It Take to Solve the 3D Ising Model? Minimal Necessary Conditions for a Valid Solution. Entropy, 24(11), 1665. https://doi.org/10.3390/e24111665

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop