---
title: Complex-Temperature Ising Model
url: https://www.emergentmind.com/topics/complex-temperature-ising-model
type: topic
---

# Complex-Temperature Ising Model

The complex-temperature Ising model refers to the analytic continuation of the classical Ising model’s partition function, magnetization, and correlation functions to complex values of the inverse temperature $\beta = J/(k_B T)$ or, equivalently, to complex coupling parameters such as $\tanh\beta$. This extension uncovers rich mathematical and physical structures, including new “phases” in the complex plane, modulated correlations ("helimagnetism"), connections to quantum dynamics, and intricate relations among zeros of the partition function, analytic singularities, computational intractability, and the emergence of long-range order.

## 1. Analytic Continuation and Partition Function Structure

The Ising model on a lattice (or general graph) is defined by the Hamiltonian
\[
H[\{\sigma\}] = -\sum_{\langle i,j\rangle} J_{ij} \sigma_i \sigma_j - \sum_i h_i \sigma_i
\]
with $\sigma_i \in \{\pm 1\}$. The partition function at complex $\beta$ is
\[
Z(\beta) = \sum_{\{\sigma\}} \exp\left[-\beta H(\{\sigma\})\right], \quad \beta \in \mathbb{C}.
\]
On the square lattice, Onsager’s solution and its generalizations apply, with the free-energy density for complex $\beta$ in the isotropic case ($J_x = J_y = J$)
\[
-\beta f(\beta) = \ln 2 + \frac{1}{2\pi} \int_0^\pi d\theta \,\ln [C^2 - 2CS\cos\theta + 1]
\]
where $S = \sinh 2\beta$, $C = \cosh 2\beta$ [2001.05475]. For anisotropic couplings, further analytic continuation is required [2001.05475]. 

For general graphs of maximum degree $\Delta$, the complex-temperature partition function can be recast as
\[
Z_{\mathrm{Ising}}(G; \beta) = \sum_{\sigma} \prod_{\{u,v\} \in E} \beta^{\frac{1+\sigma(u)\sigma(v)}{2}}
\]
allowing deep combinatorial and algorithmic analyses [2105.00287].

## 2. Fisher Zeros and Phase Boundaries in the Complex Plane

The zeros of $Z(\beta)$ in the complex $\beta$-plane (“Fisher zeros”) are central to delineating phase boundaries, generalizing the role of critical points at real temperature. For the 2D isotropic Ising model, Fisher showed the zeros condense on two circles in the complex $t = \tanh\beta$-plane:
\[
|t - 1| = \sqrt{2} \quad\text{and}\quad |t + 1| = \sqrt{2}
\]
with intersections at $t = \pm i$ [2001.05475]. In the anisotropic model, these loci deform into arcs that define four regions: a ferromagnetic (FM) bubble, a paramagnetic (PM) region, and two non-ferromagnetic (NFM1, NFM2) regions. The Fisher zero structure defines bona fide phase boundaries in the thermodynamic limit [2001.05475, 1304.6314].

For finite-size 1D chains, zeros cluster differently: with open boundaries, all zeros stay at $y=-1$ in the $y = e^{2\beta J}$ plane, while for periodic boundary conditions, zeros condense towards the origin along rays as system size increases, directly impacting low-temperature expansions [1407.6919].

## 3. Correlation Functions and Spatially Modulated Order

In regions of the complex-temperature plane where the two leading eigenvalues of the transfer matrix become equimodular, the system exhibits spatially modulated (“helimagnetic”) correlations:
\[
\langle \sigma_0 \sigma_r \rangle \sim |\lambda_2/\lambda_1|^r e^{i r q}
\]
with $q = \arg(\lambda_2/\lambda_1)$ [1304.6314, 2001.05475]. In the 1D limit, $q = \arg(\tanh\beta)$. In the 2D anisotropic model, spatial modulation and long-range order can persist throughout distinct NFM regions of the plane; the FM phase maintains nonzero uniform magnetization only within a central “bubble.”

The magnetization itself generalizes to
\[
M(\beta_x, \beta_y) = 
\begin{cases}
\left[1 - (\sinh 2\beta_x \sinh 2\beta_y)^{-2}\right]^{1/8} & \text{if } |\sinh 2\beta_x \sinh 2\beta_y| > 1 \\
0 & \text{otherwise}
\end{cases}
\]
with the branch cut chosen for continuity with real-temperature results [2001.05475].

## 4. Dimensional Crossover: Ladders, Boundary Conditions, and Zeros

Finite-width Ising ladders interpolate between 1D chains and 2D planes. The ladder partition function $Z_M$ (for $M$ legs) is determined by the eigenvalue spectrum of a $2^M \times 2^M$ transfer matrix. Zeros of $Z_M$ arise where two eigenvalues become equimodular $|\lambda_1| = |\lambda_2|$, marking phases with persistent spatial modulation [2001.05475, 1304.6314]. As $M \to \infty$, the loci of zeros fill extended areas between Fisher circles, resulting in a phase diagram with regions of spatially modulated order.

Boundary conditions dramatically affect the analytic structure of zeros. For 1D periodic chains, zeros approach the real axis as $1/N$ and strongly alter low-temperature behavior compared to open chains, where zeros remain static [1407.6919]. Thus, both boundary conditions and finite-size scaling encode subtle analytic information even in the absence of finite-temperature transitions.

## 5. Algorithmic and Quantum Information Aspects

Complex-temperature Ising partition functions connect directly to computational complexity and quantum measurement protocols. For bounded-degree graphs, there exists a large zero-free region in the complex $\beta$-plane,
\[
\left|\frac{\beta-1}{\beta+1}\right| < \tan\left(\frac{\pi}{4\Delta-4}\right)
\]
within which a deterministic fully polynomial-time approximation scheme (FPTAS) is provably available; outside a strictly larger circle, approximation becomes $\#P$-hard [2105.00287]. The existence of zeros directly implies algorithmic intractability.

Quantum circuits of constant depth enable the *measurement* of complex-temperature Ising partition functions, exploiting commuting entangling gates and analytic continuation from imaginary time protocols. Wick rotations relate measured amplitudes $A(\alpha)$ for imaginary “Hamiltonian time” to $Z(i\alpha)$, with precise statistical error bounds achievable via central-limit theorems [1208.3918]. This approach allows practical probing of partition functions and the associated Jones polynomial evaluations in moderate-sized systems.

Table: Zones of Algorithmic Tractability for $Z_{\mathrm{Ising}}(G;\beta)$ [2105.00287]

| Region in $\beta$-plane                                           | Approximation Complexity     | Zeros Present?      |
|-------------------------------------------------------------------|-----------------------------|---------------------|
| $\left|\frac{\beta-1}{\beta+1}\right| < \tan\left(\frac{\pi}{4\Delta-4}\right)$ | FPTAS                      | Zero-free           |
| $\left|\frac{\beta-1}{\beta+1}\right| > 1/\sqrt{\Delta-1}$        | $\#P$-hard                  | Zeros occur         |
| Intermediate annulus                                              | No FPTAS, evidence of zeros | Zeros numerically   |

Within the zero-free disc, Taylor expansion and Barvinok–Patel–Regts cluster expansion strategies guarantee efficient computation, while zeros and associated Fisher phase boundaries in the complex plane demarcate the computational transition.

## 6. Special Unitary Points and Quantum Dynamical Connections

Points in the complex $\beta$-plane where $\tanh\beta = \pm i$ are “special unitary points.” Here, the transfer matrix $T(\beta)$ is (after normalization) unitary, and all eigenvalues reside on the unit circle. Physically, these points admit a quantum-dynamics interpretation: $T(\beta)$ realizes a unitary quantum evolution, with the partition function representing a perfectly coherent quantum circuit [2001.05475]. In tensor-network language, each Boltzmann weight becomes a unitary gate and the classical partition sum generalizes to a quantum circuit amplitude.

These unitary points coincide with multicritical junctions of the phase diagram (Fisher circles' intersection points) and are crucial for mapping between equilibrium statistical mechanics at complex temperature and protocols in quantum information science.

## 7. Thermodynamic and Physical Consequences

Analytic continuation endows the Ising model with a more intricate phase diagram, defined not by thermodynamic phases but by regions of holomorphy, loci of degenerate transfer-matrix eigenvalues, and modulated correlation functions (“helimagnetism”). The partition function's zeros correspond to both analytic singularities and computational hardness boundaries. Thermodynamic observables such as free energy, internal energy, and specific heat become multivalued and may possess complex branch cuts and discontinuities across these boundaries [1304.6314, 2001.05475].

In the 2D limit, only the Fisher circles (specific loci in the complex $\tanh\beta$-plane) are associated with singularities in a given topological sector; for finite-width ladders and 1D chains with periodic boundary conditions, zeros fill extended lines or regions, fundamentally altering physical and analytic behavior even in the absence of true transitions at positive real temperature [1407.6919].

The complex-temperature Ising model thus serves as a rich theoretical testbed, connecting statistical mechanics, complex analysis, computational complexity theory, and quantum information science.

Source: https://www.emergentmind.com/topics/complex-temperature-ising-model