Exact solution of the Glauber-Ising model on the finite-length semi-open chain

This paper derives the exact time-space correlation function for the finite-length semi-open Glauber-Ising model quenched to zero temperature, enabling the calculation of the dual coagulation-diffusion process's empty-interval probability and confirming consistency with dynamical finite-size scaling theory.

The Gist

Imagine a long, narrow hallway filled with people standing in a line. Each person can face either Left or Right. This is the "Ising model" in a nutshell: a simple line of spins (people) that can be in one of two states.

Now, imagine the lights go out, and everyone starts randomly changing their direction based on what their immediate neighbors are doing. This chaotic shuffling is the "Glauber dynamics." The paper by Malte Henkel asks a very specific question: How does this shuffling behave if the hallway has a specific length and a "hard wall" at one end?

Here is the breakdown of the paper's findings using everyday analogies:

1. The Setup: A Finite Hallway with a Wall

Usually, physicists study these lines as if they were infinite (going on forever) or circular (like a race track where the end connects to the start). But in the real world, things have ends.

• The Scenario: Imagine a hallway with $N$ spots.
• The Rules:
  • The person at the very left end is glued to the wall and cannot move (they are fixed).
  • The person at the very right end is told to leave (their influence must vanish).
  • Everyone in between is free to shuffle, but they are influenced by their neighbors.

The paper solves the math for exactly how the "agreement" (correlation) between the fixed person on the left and anyone else in the hallway changes over time.

2. The "Magic Mirror" Trick

The biggest headache in solving this math is the "hard wall" at the end. Standard math tools (like Fourier series) love circles or infinite lines, but they hate hard stops.

The author uses a clever trick called spatial symmetry.

• The Analogy: Imagine you have a mirror placed at the end of the hallway. Instead of trying to solve the problem with a wall, the author pretends the hallway continues through the mirror into a "ghost world."
• In this ghost world, the rules are different: if the real person is facing Left, the ghost person faces Right (or vice versa, depending on the math).
• By creating this "ghost hallway," the hard wall disappears, and the problem becomes a smooth, continuous wave that is much easier to solve. Once the math is done, the author folds the ghost world back into the real one to get the final answer.

3. The Result: How Fast Does Order Spread?

The paper calculates an exact formula for how the "mood" of the line settles down.

• The Finding: The way the people align depends on two things: how much time has passed and how long the hallway is.
• The Surprise: The author tested a "shortcut" method that physicists often use to guess these answers quickly. The shortcut assumes the hallway is so long that the walls don't matter yet.
  • The Verdict: The shortcut works great if the hallway is huge compared to how far the "mood" has spread. But if the hallway is short, the shortcut fails. The exact math shows that the "hard wall" changes the shape of the curve in a way the shortcut misses. It's like trying to predict traffic flow in a small cul-de-sac using a formula designed for a highway; the results look similar at first, but the details are wrong.

4. The Secret Connection: The "Empty Seat" Game

Here is the most fascinating part of the paper. The author reveals that this problem of people facing Left/Right is mathematically identical to a completely different game: The Empty Seat Game.

• The Game: Imagine a circular table (a ring) with seats. Some seats are empty; some have a person sitting in them.
• The Rule: If two people sit next to each other, they merge into one person (coagulation), leaving an empty seat behind. People also randomly hop to empty seats next to them (diffusion).
• The Connection: The paper proves that calculating the "agreement" between two people in the fixed-wall hallway is exactly the same as calculating the probability of finding a long stretch of empty seats on the circular table.
• Why it matters: This allows scientists to use the solution for the "people in a line" to instantly solve the "empty seats on a ring" problem, including how fast the number of people on the ring decreases over time.

Summary

In simple terms, this paper is a masterclass in solving a specific type of "shuffling" puzzle on a finite line with a wall.

1. It uses a mirror trick to turn a difficult "wall" problem into an easy "circle" problem.
2. It provides the exact recipe for how the system behaves, showing that simple shortcuts often fail when the system is small.
3. It reveals a secret twin: The behavior of this line of people is mathematically identical to how empty spaces behave in a game of merging particles on a ring.

The paper doesn't promise to cure diseases or build new engines; it simply provides a precise, mathematical map of how order and disorder evolve in a confined space, which serves as a perfect benchmark for future scientists to test their own theories against.

Technical Summary

Problem Statement
The paper addresses the exact solution of the time-space correlation function for the one-dimensional Glauber-Ising model quenched to temperature $T=0$ on a finite-length, semi-open lattice of size $N$. While the dynamics of the infinite system and periodic finite chains are well-understood, the treatment of non-periodic (specifically semi-open) boundary conditions presents a challenge. The standard equation of motion for the correlator $C_n(t)$ involves a constraint $C_0(t)=1$ (due to a fixed spin) and, in the semi-open case, a vanishing condition at the other end ($C_N(t)=0$). These constraints preclude the immediate application of standard Fourier series methods, which typically rely on periodicity. The authors aim to derive the exact correlation function for this geometry and analyze how finite-size effects modify the dynamical scaling behavior expected in the infinite-size limit. Additionally, the paper seeks to exploit the duality between the Glauber-Ising model and the 1D coagulation-diffusion process to determine exact empty-interval probabilities and particle concentrations on a periodic finite ring.

Methodology
The authors employ a spatial symmetry method to transform the boundary value problem into a form solvable by Fourier analysis. The core technique involves:

1. Analytic Continuation: Instead of solving the partial differential equation (PDE) directly with hard boundary conditions, the authors extend the physical correlation function $C(t; x)$ to negative spatial arguments using specific symmetry relations.
   • For the periodic case, the function is extended such that $C(t; -x) = 2 - C(t; x)$ and $C(t; x) = C(t; N-x)$, resulting in a periodicity of $2N$.
   • For the semi-open case (fixed spin at $x=0$, vanishing at $x=N$), the function is extended such that $C(t; -x) = 2 - C(t; x)$ and a linear transformation is applied to satisfy the vanishing boundary condition at $N$. This introduces an auxiliary function $B(t; x)$ which is shown to be anti-symmetric and periodic with period $2N$.
2. Fourier Series Representation: By embedding the physical constraints into the spatial symmetries of the analytically continued function, the authors construct a Fourier series representation. This allows the equation of motion (a diffusion equation in the continuum limit) to be solved exactly in Fourier space.
3. Duality Mapping: The resulting exact solution for the semi-open Glauber-Ising model is mapped to the empty-interval probability $E(t, x)$ of the dual 1D coagulation-diffusion process on a periodic ring. This mapping relies on the identity of the equations of motion (up to a time rescaling) and boundary conditions between the two systems.
4. Scaling Analysis: The exact solutions are analyzed in the scaling limit ($t \to \infty, N \to \infty, |x| \to \infty$) with finite scaling variables $u = |x|/\sqrt{t}$ and $v = N/\sqrt{t}$. The authors also test a simplified scaling ansatz $C_{scal}(t; x) = F(x/\sqrt{t})$ against the exact results to determine its validity range.

Key Contributions and Results

• Exact Correlation Functions: The paper derives the exact expression for the single-time correlator $C(t; x)$ on a semi-open chain. The solution is expressed in terms of Jacobi theta functions ($\vartheta_2, \vartheta_3$) and depends on the finite-size scaling variables $|x|/N$ and $t/N^2$.
  • For the semi-open case: $C_{semi}(t; x) = 1 - \int_0^{|x|/N} dv \, \vartheta_3(\frac{\pi}{2}v, e^{-\pi^2 t / 2N^2}) + \text{initial condition terms}$.
  • For the periodic case, a similar exact form is recovered, confirming consistency with previous literature.
• Finite-Size Scaling: The results demonstrate that the correlation functions satisfy dynamical finite-size scaling, $C(t; x; N) = F_C(|x|/\sqrt{t}, N/\sqrt{t})$. The universal scaling function $F_C$ is shown to be dependent on the boundary conditions.
  • The semi-open correlator vanishes exactly at $x=N$, whereas the periodic correlator decays exponentially to zero at $x=N/2$ (in the scaled variable) for large $N$.
  • The behavior of the semi-open correlator in the interval $0 \le |x|/N \le 1$ is analogous to the periodic correlator in $0 \le |x|/N \le 1/2$, though the boundary conditions induce distinct qualitative features (e.g., the "cusp" at $x=0$ due to hard Ising spins).
• Coagulation-Diffusion Application: Using the duality, the authors provide the exact empty-interval probability $E(t, x)$ for particles on a periodic ring and derive the time-dependent particle concentration $c(t)$.
  • The concentration is given by $c(t) = \frac{1}{N}\vartheta_3(0, e^{-2\pi^2 t/N^2}) + \text{initial condition terms}$.
  • The long-time decay $c(t) \sim t^{-1/2}$ is reproduced, consistent with the known behavior for $d=1$ coagulation processes.
• Validation of Scaling Ansatz: The paper tests a simplified scaling ansatz (leading to an error function form) against the exact solution.
  • The ansatz $C_{scal}(t; x) = 1 - \text{erf}(x/\sqrt{2t})/\text{erf}(N/\sqrt{2t})$ is found to be an excellent approximation only when the scaling variable $y = N/\sqrt{t} \gtrsim 2$.
  • For small $y$ (where the domain size $\ell(t) \sim \sqrt{t}$ becomes comparable to $N$), the simple ansatz fails to capture the crossover behavior, indicating that the single-variable scaling assumption is an over-simplification for finite geometries near the crossover regime.

Significance
The paper provides the first explicit exact results for dynamical non-equilibrium many-body systems with non-periodic (semi-open) boundary conditions. By successfully adapting the spatial symmetry method to open boundaries, the authors offer a rigorous benchmark for finite-size scaling theory in dynamics. The work clarifies how boundary conditions influence the form of universal scaling functions, a detail often obscured in infinite-system approximations. Furthermore, the exact solution for the coagulation-diffusion process on a finite ring offers a precise tool for studying particle concentration decay and empty-interval probabilities, which are relevant for understanding physical ageing and anomalous transport in one dimension. The authors suggest these techniques could serve as a background for more generic studies of finite-size effects and may be extendable to quantum dynamics and integrable chains, though such extensions are noted as potential future directions rather than immediate results of this work.