Exact solution of a two-dimensional (2D) Ising model with the next nearest interactions

Zhidong Zhang

arxiv: 2601.10749 · v6 · pith:WJKK2WJGnew · submitted 2026-01-13 · ❄️ cond-mat.stat-mech

Exact solution of a two-dimensional (2D) Ising model with the next nearest interactions

Zhidong Zhang This is my paper

Pith reviewed 2026-05-16 15:25 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech

keywords 2D Ising modelnext-nearest-neighbor interactionsexact solutionpartition functionspontaneous magnetizationcritical temperaturetransfer matrixtriangular lattice

0 comments

The pith

The exact solution is derived for the two-dimensional Ising model with next-nearest-neighbor interactions at zero magnetic field, yielding closed-form partition function and spontaneous magnetization.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that a two-dimensional Ising model with next-nearest-neighbor interactions at zero magnetic field admits an exact solution. It reaches this result by first analyzing the transfer matrices in Clifford algebraic, transfer tensor, and schematic representations, then showing equivalence to a triangular Ising model supplemented by an interaction along a third axis. This equivalence permits the direct adaptation of solution methods previously used for three-dimensional Ising models. The resulting expressions quantify how added interactions and lattice topology raise the critical temperature relative to simpler Ising lattices.

Core claim

The two-dimensional Ising model with next-nearest interactions at zero field is equivalent to a triangular Ising model plus an interaction along the z axis. By modifying the approaches developed for the 3D Ising model, the partition function and the spontaneous magnetization are obtained exactly.

What carries the argument

Equivalence mapping of the 2D next-nearest-neighbor Ising lattice to a triangular Ising model coupled along a perpendicular axis, together with transfer-matrix analysis in three representations.

If this is right

The critical temperature rises when the number of interactions per unit cell increases.
Topological contributions in the lattice further elevate the critical point.
Exact thermodynamic quantities, including magnetization, become available without approximation for this interaction set.
Comparisons across different Ising lattices confirm that both interaction count and topology stabilize ordered phases.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same mapping strategy may extend to other two-dimensional lattices that include competing or longer-range couplings.
Real materials whose magnetic layers realize next-nearest interactions could display higher transition temperatures than nearest-neighbor models predict.
Finite-size numerical checks could directly verify the analytic magnetization formula in the vicinity of the critical point.

Load-bearing premise

The two-dimensional model with next-nearest interactions is equivalent to a triangular Ising model plus an interaction along the z axis, so that three-dimensional Ising solution methods apply directly.

What would settle it

A large-scale Monte Carlo simulation on the same lattice that produces a critical temperature or spontaneous-magnetization curve differing from the closed-form expressions derived in the paper.

Figures

Figures reproduced from arXiv: 2601.10749 by Zhidong Zhang.

**Figure 1.** Figure 1: Illustration of a rectangular Ising model with the next nearest interactions on a 66 lattice. The nearest neighboring interactions 𝐾1 and 𝐾2 are represented by black lines, while the next nearest interactions 𝐾3 and 𝐾4 are represented by blue and red lines respectively [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗

**Figure 2.** Figure 2: Illustration of a triangular Ising model with interactions 𝐾1, 𝐾2 and 𝐾3 on a 66 lattice plus an interaction 𝐾4 along the z direction at each lattice point. The [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

**Figure 3.** Figure 3: Temperature dependence of the spontaneous magnetization of the 2D Ising model with the next nearest interactions when 𝐾1 = 𝐾2 and 𝐾3 = 𝐾4 . The curves from left to right correspond to the next nearest interaction 𝐾3 = 0, 0.5𝐾1 , 𝐾1 , 1.0001𝐾1, 1.25𝐾1 and 1.5𝐾1, respectively [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗

read the original abstract

The exact solution of a two-dimensional (2D) Ising model with the next nearest interactions at zero magnetic field is derived. At first, the transfer matrices are analyzed in three representations, i.e., Clifford algebraic representation, transfer tensor representation and schematic representation, to inspect nontrivial topological structures in this system. The system is equivalent to a triangular Ising model plus an interaction along the z axis, so that the approaches developed for the 3D Ising model are modified to be appropriable for solving the exact solution of the 2D Ising model with the next nearest interactions. The partition function and the spontaneous magnetization are obtained. The comparison with the exact solutions of other Ising lattices reveals that either the increase of the number of interactions in a unit cell or the presence/increase of topological contributions enhances the critical point of the Ising lattices. The results obtained in this work are helpful for understanding the physical properties of the 2D magnetic materials.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The exact-solution claim for this 2D NNN Ising model hinges on a mapping to a triangular-plus-z model whose solvability is not obvious, since standard 3D Ising methods do not yield closed forms.

read the letter

The paper's main move is to treat the 2D Ising lattice with next-nearest couplings as equivalent to a triangular Ising model plus an interaction along the z direction, then adapt transfer-matrix techniques from 3D Ising work. It analyzes the transfer matrices in Clifford-algebra, tensor, and schematic representations to extract topological features, and it reports closed-form expressions for the partition function and spontaneous magnetization at zero field. It also compares critical points across several lattices and concludes that more interactions or added topology push the transition temperature higher. Those comparisons are straightforward and could serve as a small benchmark for people modeling 2D magnets. The topological inspection of the transfer matrices is a sensible step that might help classify similar systems. The central difficulty is the mapping itself. The 3D Ising model has no known exact solution, so any adaptation only produces a true closed form if the equivalence introduces a special cancellation or reduction that the abstract does not spell out. Without the algebra showing how the z-term is handled and why the result remains exact rather than approximate, the claim rests on an unverified step. The abstract presents the equivalence as given, which leaves open the possibility that the final expressions either contain hidden assumptions or are not fully rigorous. This work is aimed at researchers who care about exactly solvable cases in statistical mechanics and at those studying how lattice geometry affects criticality in 2D materials. A reader who wants concrete expressions for this specific lattice might extract the comparison data, but anyone needing a verified benchmark should treat the main result as provisional until the mapping is checked line by line. The paper deserves peer review because the lattice is concrete and the claim is testable; a referee can focus on whether the modified 3D methods actually close the solution or whether the equivalence merely restates the problem.

Referee Report

2 major / 2 minor

Summary. The manuscript claims an exact solution for the 2D Ising model with next-nearest-neighbor interactions at zero field. Transfer matrices are analyzed in Clifford algebraic, transfer-tensor, and schematic representations to reveal topological structures. The model is mapped to a triangular Ising model plus a z-axis interaction, allowing modification of 3D Ising methods to derive the partition function and spontaneous magnetization. Comparisons with other Ising lattices show that more interactions or topological contributions raise the critical point.

Significance. An exact closed-form solution for this 2D Ising variant would be significant for statistical mechanics, providing a benchmark for 2D magnetic materials and clarifying how interaction range and topology affect criticality. The use of multiple representations and the mapping approach, if rigorously justified, could offer reusable techniques for related models.

major comments (2)

The load-bearing step is the asserted equivalence of the 2D NNN model to a triangular Ising model plus z-axis interaction, followed by direct modification of 3D Ising techniques (Clifford, transfer-tensor, schematic). Since the 3D Ising model has no known closed-form solution, the manuscript must explicitly demonstrate what property of this mapping renders the problem exactly solvable rather than inheriting the intractability; without that justification the derived partition function cannot be exact. (Abstract; section on equivalence and transfer-matrix analysis)
The abstract states that the partition function and spontaneous magnetization are obtained exactly, yet the provided description supplies no explicit final expressions, verification against known limits (e.g., reduction to Onsager solution when NNN coupling vanishes), or numerical checks. These must be supplied with derivation steps to substantiate the claim.

minor comments (2)

Notation for the three representations (Clifford, transfer tensor, schematic) should be introduced with explicit definitions and a table comparing their structures for clarity.
The comparison of critical points across lattices would benefit from a table listing the models, interaction counts, topological features, and reported Tc values.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback. We address the major comments point by point below, with plans to revise the manuscript accordingly to strengthen the justification and presentation of the exact solution.

read point-by-point responses

Referee: The load-bearing step is the asserted equivalence of the 2D NNN model to a triangular Ising model plus z-axis interaction, followed by direct modification of 3D Ising techniques (Clifford, transfer-tensor, schematic). Since the 3D Ising model has no known closed-form solution, the manuscript must explicitly demonstrate what property of this mapping renders the problem exactly solvable rather than inheriting the intractability; without that justification the derived partition function cannot be exact. (Abstract; section on equivalence and transfer-matrix analysis)

Authors: We agree that the mapping requires more explicit justification to clarify why it permits an exact solution. In the revised manuscript we will expand the equivalence section to detail how the specific topological structures identified via the Clifford algebraic, transfer-tensor, and schematic representations of the transfer matrices enable a factorization. This factorization, arising from the z-axis interaction in the mapped model, allows the modified 3D techniques to yield closed-form results without inheriting the general intractability of the 3D Ising model. Step-by-step reasoning for this property will be added. revision: yes
Referee: The abstract states that the partition function and spontaneous magnetization are obtained exactly, yet the provided description supplies no explicit final expressions, verification against known limits (e.g., reduction to Onsager solution when NNN coupling vanishes), or numerical checks. These must be supplied with derivation steps to substantiate the claim.

Authors: We accept that the explicit forms and verifications should be more prominent. The derivations are present in the full text, but we will revise to include the closed-form expressions for the partition function and spontaneous magnetization in the main body, together with a dedicated verification subsection. This will demonstrate the reduction to the Onsager solution when the next-nearest-neighbor coupling is set to zero, include the key derivation steps, and add numerical checks against Monte Carlo simulations for representative parameter values. revision: yes

Circularity Check

1 steps flagged

Exact solution rests on asserted equivalence to triangular-plus-z model, allowing modification of 3D Ising methods

specific steps

self definitional [Abstract]
"The system is equivalent to a triangular Ising model plus an interaction along the z axis, so that the approaches developed for the 3D Ising model are modified to be appropriable for solving the exact solution of the 2D Ising model with the next nearest interactions."

The equivalence is taken as an un-derived premise. The partition function and spontaneous magnetization are then obtained by modifying 3D Ising techniques on the basis of this premise. The final expressions are therefore equivalent to re-expressing the assumed mapping rather than deriving new results from the 2D model equations, satisfying the self-definitional pattern.

full rationale

The paper asserts without derivation that the 2D NNN Ising model is equivalent to a triangular Ising model plus z-axis interaction, then modifies approaches from the 3D Ising model to obtain the partition function and magnetization. Because the 3D Ising model has no known closed-form solution, the claimed exact results reduce to consequences of this input equivalence rather than an independent derivation from the original 2D Hamiltonian. The visible text provides no equations showing how the mapping is constructed or verified, making the central claim dependent on the premise by construction. This is partial circularity (score 6) rather than total, as the paper does contain analysis of transfer matrices in multiple representations.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the unproven equivalence between the 2D next-nearest model and a triangular lattice plus z-interaction, plus the assumption that 3D Ising solution techniques remain valid after this mapping.

axioms (1)

domain assumption The 2D Ising model with next-nearest interactions is equivalent to a triangular Ising model plus an interaction along the z axis.
Explicitly invoked in the abstract as the justification for adapting 3D Ising methods.

pith-pipeline@v0.9.0 · 5457 in / 1377 out tokens · 36519 ms · 2026-05-16T15:25:36.168007+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith.Constants phi_golden_ratio / phi_sq_eq echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

The critical point for the ferromagnetic-paramagnetic phase transition is located at the golden ratio xc = e−2Kc = (√5−1)/2 , 1/Kc = 4.15617384.. for the cubic Ising lattice
IndisputableMonolith.Foundation.AlexanderDuality alexander_duality_circle_linking echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

The system is equivalent to a triangular Ising model plus an interaction along the z axis... nontrivial topological structures

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.