Zigzag ordering, defects, and anomalous relaxation in antiferromagnetic Kuramoto lattices

Technical Summary: Zigzag Ordering, Defects, and Anomalous Relaxation in Antiferromagnetic Kuramoto Lattices

Problem Statement
The paper investigates the nonequilibrium ordering dynamics of coupled Kuramoto oscillators with negative (antiferromagnetic) nearest-neighbor coupling. While positive coupling in the Kuramoto model is well-studied for its synchronization transitions, the dynamics of repulsive coupling—which favors phase differences of $\pi$ and induces antiphase ordering—remain less explored. Specifically, the authors address the lack of understanding regarding the coarsening dynamics, defect kinetics, and scaling properties of these systems. The central problem is to determine how deterministic nonlinear interactions and geometric frustration alone, without quenched disorder or stochastic noise, drive the system toward an ordered state and whether this process follows standard universality classes (such as the Ising model) or exhibits anomalous scaling.

Methodology
The study employs both numerical simulations and analytical continuum approximations:

• Numerical Simulations: The authors simulate 1D chains and 2D square lattices of $N$ oscillators governed by the Kuramoto equation with $K < 0$. The equations of motion are integrated using a fourth-order Runge–Kutta scheme.
  • Defect Definition: A defect is identified where the local slope sign ($s_i = \text{sign}(\phi_i - \phi_{i-1})$) fails to alternate, i.e., $|s_i + s_{i+1}| = 2$. The order parameter is the normalized defect density $D(t)$.
  • Persistence: Two types of persistence are measured: slope persistence ($P_s$), the probability a local gradient retains its initial sign, and field persistence ($P_\phi$), the probability a phase variable does not cross zero.
  • Parameters: Simulations vary system sizes ($N$), coupling strengths ($K$), and intrinsic frequencies ($\omega$) to test scaling independence.
• Analytical Approach: A continuum limit is derived via Taylor expansion of the discrete equations, resulting in a coarse-grained partial differential equation (PDE). Linear stability analysis is performed on this PDE to identify unstable modes and the dominant terms governing relaxation.

Key Results

1. One-Dimensional Dynamics (Anomalous Coarsening)

• Defect Decay: The defect density decays algebraically as $D(t) \sim t^{-1/4}$ at early times. This is significantly slower than the $t^{-1/2}$ decay typical of standard nonconserved phase ordering (e.g., Model A or Ising systems) with a dynamical exponent $z=2$.
• Dynamical Scaling: Finite-size scaling analysis reveals a dynamical exponent $z=2$, indicating that the characteristic relaxation time scales as $t_c \sim N^2$. This suggests diffusive domain growth ($L(t) \sim t^{1/2}$).
• Universality Class: The coexistence of $z=2$ and the defect decay exponent $\delta = 1/4$ defines a distinct universality class. The authors note that while $z=2$ is common, the combination with $\delta=1/4$ violates standard scaling relations where $\delta = 1/z$ for annihilating random walkers.
• Persistence: Both slope and field persistence follow a stretched-exponential form, $P(t) \sim \exp(-c t^{1/4})$. The stretching exponent $\alpha = 1/4$ exactly matches the defect decay exponent $\delta$.
• Saturation: In finite systems, the defect density and persistence saturate at a time $t_c \sim N^2$, leaving a residual, nearly equidistant defect configuration. In the thermodynamic limit, long-range zigzag order is recovered.
• Role of Frequency: The scaling exponents are independent of the intrinsic frequency $\omega$ and coupling magnitude $K$ (provided $K<0$), which only rescale the time axis. However, a non-zero $\omega$ prevents the saturation of persistence (memory is eventually lost), whereas $\omega=0$ leads to a frozen metastable state.

2. Two-Dimensional Dynamics (Geometric Frustration)

• In 2D, geometric frustration inhibits the formation of global antiferromagnetic order.
• While local checkerboard patterns emerge, the system becomes trapped in metastable states characterized by extended domain walls (defect lines) rather than point defects.
• The defect density exhibits an initial transient decay but eventually saturates at a non-zero value, preventing complete ordering within the simulated time scales. The elimination of these extended interfaces requires collective rearrangements that are kinetically suppressed.

Theoretical Explanation
The authors derive a continuum equation for the phase field $\phi(x,t)$:
$$\partial_t \phi = -|K| \partial_x^2 \phi - \frac{|K|}{12} \partial_x^4 \phi$$

• Instability: The negative diffusion term ($\partial_x^2$) drives an instability, while the biharmonic term ($\partial_x^4$) stabilizes short wavelengths. This competition selects a zigzag wavelength.
• Scaling Origin: The $t^{-1/4}$ decay is attributed to the dominance of the biharmonic term (surface-diffusion-like behavior) at intermediate scales, which governs defect smoothing. However, the long-wavelength dynamics are eventually dominated by the diffusive $\partial_x^2$ term, leading to the $z=2$ scaling for the global equilibration time. The authors argue that the defect motion is "slaved" to the slower smoothing of the zigzag background, resulting in the anomalous $\delta = 1/4$ exponent despite the diffusive $z=2$ background.
• Persistence Derivation: A kinetic argument links the persistence exponent to the defect density. Assuming the encounter rate of defects scales as $\Gamma(t) \sim t^{-3/4}$ (derived from $D(t) \sim t^{-1/4}$ and domain size $\xi(t) \sim t^{1/2}$), the cumulative hazard integrates to $t^{1/4}$, yielding the stretched-exponential persistence.

Significance and Claims
The paper claims that deterministic nonlinear dynamics and geometric frustration are sufficient to generate slow relaxation and anomalous scaling in many-body systems, without the need for quenched disorder or stochastic noise.

• Breakdown of Standard Scaling: The work demonstrates a breakdown of standard dynamic scaling relations in antiferromagnetic Kuramoto lattices, revealing a unique universality class characterized by $\delta = 1/4$ and $z=2$.
• Mechanism of Slow Relaxation: It establishes that the slow relaxation arises from the constrained mobility of defects within an evolving zigzag background, rather than from external randomness.
• Dimensionality Dependence: The study highlights how dimensionality fundamentally alters the outcome: 1D systems eventually reach long-range order with anomalous kinetics, while 2D systems are trapped in frustrated, metastable states with persistent domain walls.
• Theoretical Insight: The continuum approximation successfully explains the coexistence of quartic amplitude decay and quadratic time scaling, providing a theoretical framework for the observed exponents.

The authors conclude that these results offer a new perspective on coarsening dynamics in continuous phase systems, showing that frustration alone can lead to complex, slow relaxation phenomena distinct from standard spin models.