🔢 mathematics

Approximate thermodynamics of the two-dimensional Ising model in an external magnetic field

This paper introduces an approximate analytical framework for the thermodynamics of the two-dimensional Ising model in an external magnetic field by utilizing a quaternion representation and small-angle approximation to derive closed-form expressions for key properties like free energy and magnetization, which are validated against BKL algorithm simulations.

Original authors: Mathieu Gaudreault

Published 2026-07-21

📖 3 min read🧠 Deep dive

CC BY 4.0

Original authors: Mathieu Gaudreault

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world made of tiny, invisible magnets called "spins." In some materials, these spins are like a row of soldiers: at low temperatures, they all stand at attention, pointing in the same direction, creating a strong magnetic field. But as the material heats up, the soldiers get jittery and start pointing in random directions, losing their collective power. This tug-of-war between order and chaos is the heart of a famous puzzle in physics called the Ising model. Scientists have been solving this puzzle for decades, but there's a tricky version that has stumped them: what happens when you add an outside force, like a giant magnet, pushing on these spins? It's like trying to predict how a crowd of people will move when a strong wind is blowing through them. The math gets incredibly messy because the wind breaks the perfect symmetry of the crowd, making the usual tricks for solving the puzzle fail.

This is where a researcher named Mathieu Gaudreault steps in with a clever workaround. Instead of trying to solve the messy, windy equation exactly, he proposes a "small angle" shortcut. Think of it like trying to describe a wobbly, spinning top. If the top is spinning wildly, it's hard to describe. But if you look at it from a slightly tilted angle where the wobble looks tiny, you can use simple math to guess how it behaves. Gaudreault uses a special kind of math called "quaternions" (which are like 4D numbers used to describe rotations) to tilt his view of the magnetic spins. By assuming the wind isn't too strong and the temperature is high enough, he turns the impossible, messy equation into a simpler one that looks like the ones we already know how to solve.

The paper doesn't claim to have found the perfect, exact answer for every possible situation. Instead, it offers a set of approximate formulas that act like a very good map for a specific territory: when the temperature is high and the external magnetic field isn't overwhelming. Gaudreault tested his new map by running massive computer simulations using an algorithm called BKL, which acts like a super-fast digital laboratory. The results were promising: in the regions where his "small angle" assumption held true, his formulas matched the computer simulations almost perfectly. He calculated how much the material would magnetize, how sensitive it would be to changes in the magnetic field, and how much heat it would absorb.

However, the paper is careful to point out where the map breaks down. Near the "critical temperature"—the exact moment the material switches from ordered to disordered—the approximation gets fuzzy. Just as a flat map of the Earth distorts the poles, Gaudreault's method smooths out the sharp, dramatic changes that happen right at the tipping point. He also notes that while his method works for a flat, square grid of spins, it might need tweaking for other shapes or for the even more complex 3D world. Ultimately, this work doesn't solve the entire mystery of magnets in a wind, but it provides a powerful new tool to understand them when the conditions are just right, bridging the gap between the messy reality of external fields and the elegant math of the past.

Technical Summary: Approximate Thermodynamics of the Two-Dimensional Ising Model in an External Magnetic Field

Problem Statement
The two-dimensional (2D) Ising model is a cornerstone for understanding phase transitions and critical phenomena. While the model is exactly solvable in one dimension and in two dimensions with zero external magnetic field ( $H=0$ ), the 2D case with a non-zero external field ( $H \neq 0$ ) remains unsolved analytically. The introduction of an external field breaks the system's symmetry, rendering the transfer matrix non-centrosymmetric. This loss of symmetry prevents the application of standard exact solution techniques (such as those by Onsager and Kaufman) used for the $H=0$ case, leaving a gap in closed-form analytical expressions for the system's thermodynamics, particularly regarding susceptibility near the critical temperature.

Methodology
The authors propose an approximate analytical framework based on a quaternion representation of the 2D Ising model, building upon the geometric interpretations introduced by Onsager and Kaufman. The core of the methodology involves the following steps:

Quaternion Representation: The transfer matrix $T$ for the 2D lattice is expressed using a basis of Pauli matrices ( $X, Y, Z$ ). The interaction energy and external field terms are decomposed into matrices $V_1$ and $V_2$ , where $T = V_2 V_1$ .
Rotated Frame and Small-Angle Approximation: The authors introduce a rotation in the quaternion space defined by an angle $\phi$ . By applying a small-angle approximation ( $\phi \to 0$ ), the rotated frame aligns with the reference frame. This approximation is valid under conditions of high temperature ( $T \gg 2(H+J)$ ) and weak fields, specifically where $2\beta J \ll 1$ and $2\beta H \ll 1$ .
Rescaling: To ensure the system responds correctly to the full physical field $H$ in the limit $H \to 0$ , the field parameter $\psi = \beta H$ is rescaled by a factor of 4 ( $\psi \to 4\psi$ ) within the approximation.
Restoration of Centrosymmetry: Under this approximation, the complex transfer matrix simplifies into a centrosymmetric form. This structural change allows the authors to utilize the known eigenvalue solutions developed by Kaufman for the $H=0$ case.
Derivation of Thermodynamics: The largest eigenvalue of the simplified transfer matrix is used to derive the free energy. From this, analytical expressions for magnetization ( $M$ ), magnetic susceptibility ( $\chi$ ), internal energy ( $U$ ), and specific heat ( $C$ ) are derived. These expressions involve elliptic integrals and depend on the external field strength.

Key Results
The paper presents closed-form approximate expressions for the thermodynamic quantities of the 2D Ising model in an external field.

Validation: The analytical results were compared against numerical simulations performed using the BKL (Bortz-Kalos-Lebowitz) algorithm on square lattices ( $N=20\times20$ and $50\times50$ ).
Agreement Regions: The approximate expressions show excellent agreement with numerical simulations in regions where the small-angle approximation holds (high temperatures and small fields). Specifically, the absolute difference between analytical and numerical values was found to be less than 0.22 for $H \in [-0.1, 0.1]$ at $T=0.4$ (where $J=0.01$ ), and the agreement interval widened as temperature increased.
Limitations: The approximation fails near the critical temperature ( $T_c$ ) and for strong fields where the symmetry breaking is significant. Consequently, the approximate susceptibility at $H=0$ remains smooth around $\beta_c J$ rather than exhibiting the divergence seen in exact solutions and numerical simulations. The model predicts a sharp transition singularity at $\kappa' = 0$ , which does not reflect the smooth crossover observed in non-zero fields in numerical data.

Significance and Claims
The authors claim this work provides the first comprehensive set of analytical expressions for the thermodynamics of the 2D Ising model in an external field that explicitly include the dependence on field strength and reduce to the known exact solutions when $H=0$ .

Geometric Insight: The study highlights the utility of the quaternion representation and geometric rotations in simplifying complex statistical mechanical problems.
Practical Utility: The derived expressions offer a computationally efficient alternative to numerical simulations for high-temperature and weak-field regimes.
Future Directions: The authors suggest that the framework could be extended to anisotropic cases, different lattice types (triangular, hexagonal), and potentially the three-dimensional Ising model, provided a quaternion representation exists for that dimension. They also note that while the approximation fails near $T_c$ , the geometric approach might inspire new methods for tackling the exact solution of the non-centrosymmetric transfer matrix.

The paper concludes that while the approximation has clear boundaries of validity, it successfully bridges the gap between the exactly solvable $H=0$ case and the unsolved $H \neq 0$ case, offering a tractable analytical tool for specific thermodynamic regimes.

More like this