Time-Dependent Glauber Rates

• Time-dependent Glauber rates are stochastic rules governing spin-flip dynamics by incorporating time-varying fields and spatial bias in systems like the Ising model.
• They bridge equilibrium relaxation with nonequilibrium dynamics by transitioning from detailed balance to global balance, influencing critical behavior and response functions.
• Their formulation facilitates the derivation of macroscopic phenomena, aiding research in phase transitions, driven quantum sources, and the development of linear dynamical field theories.

Time-dependent Glauber rates define the stochastic rules governing the evolution of spin systems, generalizing the classic Glauber dynamics to contexts where rates either depend on local time-dependent fields, spatial bias, or allow only approximate compliance with statistical physical constraints such as detailed balance. These rates play a central role in nonequilibrium statistical mechanics, enabling the study of driven Ising chains, single-electron quantum states, and the derivation of linear dynamical field theories. Their significance lies in their ability to interpolate between equilibrium relaxation governed by detailed balance and far-from-equilibrium scenarios controlled by global balance, providing a direct link between microscopic dynamics and macroscopic, often exactly solvable, phenomenology.

1. Definition and Functional Forms

Time-dependent Glauber rates originate from the extension of Markovian spin-flip dynamics, where the probability per unit time $w_n(\sigma,t)$ for a spin $\sigma_n = \pm 1$ to flip depends not only on its local neighborhood but also on explicit time-dependent parameters or external fields.

Generalized rate (Ising chain, one-parameter family):

$$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$

with $\alpha$ as time-scale, $\gamma = \tanh(2J/T)$ (enforcing consistency with the Ising coupling $J$ and temperature $T$), and $p\in [0,1]$ tuning spatial bias; $V = 2p-1$ provides a drift parameterization (Godreche, 2011).

With external time-dependent field $h_n(t)$:

$\sigma_n = \pm 1$0

where $\sigma_n = \pm 1$1. Here, explicit time-dependence arises solely from $\sigma_n = \pm 1$2, not from the bias itself (Godreche, 2011).

Optimal linearization (arbitrary graph):

$\sigma_n = \pm 1$3

with coefficients $\sigma_n = \pm 1$4 determined via Moore–Penrose regression to minimize detailed-balance violation (Sahoo et al., 2014). In isotropic cases, $\sigma_n = \pm 1$5 admits a closed form as a function of temperature ($\sigma_n = \pm 1$6), coupling ($\sigma_n = \pm 1$7), and coordination ($\sigma_n = \pm 1$8).

2. Conditions Imposed by Global and Detailed Balance

When the dynamics is symmetric ($\sigma_n = \pm 1$9), Glauber’s original rate satisfies detailed balance with respect to the Ising Gibbs measure:

$$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$0

This ensures reversibility on each bond (Godreche, 2011).

For generic $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$1 (“directed” or “asymmetric” cases), detailed balance fails. Global balance is imposed instead. The master equation, summed over all possible flips and weighted by the Ising-Gibbs stationary measure, requires:

$$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$2

where $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$3 denotes the configuration with spin $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$4 flipped, and $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$5 (Godreche, 2011). Detailed algebra constrains $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$6 (with $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$7) and yields the unique globally balanced rate for fully directed cases ($$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$8 or $$w_n(\sigma) = \frac{1}{2} \alpha [1 - \gamma\,\sigma_n (p\,\sigma_{n-1} + (1-p)\,\sigma_{n+1})]$$9):

$\alpha$0

Up to $\alpha$1, this is uniquely set by global balance alone (Godreche, 2011).

3. Manifestation of Time Dependence

Time dependence of Glauber rates emerges primarily through external, space- and time-dependent magnetic fields. In zero applied field, transition probabilities remain purely configuration-dependent (no explicit $\alpha$2). Under $\alpha$3, rates gain explicit $\alpha$4-dependence via $\alpha$5.

In kinetic Ising models, this means:

• Without field: Coarsening (aging) or stationary regime rates depend only on $\alpha$6.
• With $\alpha$7: $\alpha$8 encodes explicit temporal variation, relevant for linear response and nonequilibrium driving (Godreche, 2011, Sahoo et al., 2014).

For the linearized (optimal) forms relevant for analytic solution, the field-dependence causes the relaxation rates and susceptibilities to become explicit functions of $\alpha$9, with quadratic corrections in weak fields (Sahoo et al., 2014).

4. Symmetric, Asymmetric, and Continuum Limits

Time-dependent Glauber rates bridge between symmetric, detailed-balance-preserving, and asymmetric, globally-balanced forms:

Case	Explicit Rate Formula	Balance Condition
$\gamma = \tanh(2J/T)$0 (symmetric)	$\gamma = \tanh(2J/T)$1	Detailed balance
$\gamma = \tanh(2J/T)$2 or $\gamma = \tanh(2J/T)$3 (directed)	$\gamma = \tanh(2J/T)$4	Global balance only
$\gamma = \tanh(2J/T)$5 (intermediate)	Interpolates via $\gamma = \tanh(2J/T)$6 in the general rate	Interpolates between

The continuous interpolation affects key dynamical observables:

• For $\gamma = \tanh(2J/T)$7, correlation and response decay rates remain identical; the limit fluctuation-dissipation ratio ($\gamma = \tanh(2J/T)$8) decreases from $\gamma = \tanh(2J/T)$9 (equilibrium) to $J$0 at the critical bias (Godreche, 2011).
• Beyond $J$1, the response decays faster, and $J$2.

In the continuum (linearized) limit, time-dependent Glauber rates give rise to the time-dependent Ginzburg–Landau (TDGL) equation:

$J$3

with $J$4 and $J$5, directly connecting microscopic rates with macroscopic relaxation and pattern formation (Sahoo et al., 2014).

5. Observable Consequences: Magnetization, Correlation, Response, and Fluctuation-Dissipation

Magnetization: The evolution equation in the directed Ising chain takes the form:

$J$6

with solutions obtainable via Fourier-Laplace analysis. The bias parameter $J$7 modifies propagation, introducing a drift term (Godreche, 2011).

Equal-time correlations: The equation

$J$8

is independent of $J$9 and $T$0; thus, the domain growth law and scaling at $T$1 coincide with the symmetric case. This suggests bias affects dynamic but not static spatial structures (Godreche, 2011).

Two-time correlations and response: Asymmetric rates generate nontrivial temporal correlations with decay characteristics crossing over at the critical bias $T$2. The linear response function satisfies a generalized fluctuation-dissipation relation:

$T$3

Reducing to the equilibrium FDT at $T$4, this highlights the central dynamical alteration induced by directedness (Godreche, 2011).

Relaxation and fluctuation-dissipation theorem (LGM): In the linear Glauber model, the exact relaxation time in uniform field $T$5 is:

$T$6

showing quadratic $T$7-dependence (Sahoo et al., 2014). The fluctuation-dissipation relation is preserved despite approximate balance conditions:

$T$8

6. Quantum and Single-Particle Generalizations

Time-dependent Glauber rates also underpin the theory of single-electron quantum sources, where the relevant quantities are the first-order and second-order (Glauber) correlation functions:

$T$9

For periodically driven sources, this is formulated using Floquet theory. In the adiabatic case, $p\in [0,1]$0 is a symmetric Lorentzian in time, with lifetime $p\in [0,1]$1 and coherence time $p\in [0,1]$2—the Fourier-transform-limited regime. In non-adiabatic pumping, $p\in [0,1]$3 becomes asymmetric and decays exponentially. The emission probability rate $p\in [0,1]$4 connects directly to experimentally measurable current pulses (Haack et al., 2012).

A pure on-demand single-electron state ensures $p\in [0,1]$5 vanishes, realizing full anti-bunching and substantiating the quantumness of the source (Haack et al., 2012).

7. Implications and Applications

Time-dependent Glauber rates establish a rigorous framework for analyzing nonequilibrium processes in spin systems, stochastic processes on networks, and quantum sources. Their flexibility enables the study of non-reversible dynamics (directed Ising chains (Godreche, 2011)), analytic derivation of mean-field and continuum kinetic equations (linear and optimal LGM (Sahoo et al., 2014)), and fully quantum-coherent time-dependent transport (Floquet-driven sources (Haack et al., 2012)).

A notable implication is that nontrivial time dependence (in the absence of external field) arises not from spatial bias, but from the explicit introduction of a time-varying field or drive. This distinction clarifies the terminology: “time-dependent Glauber rates” conventionally refers to rates gaining $p\in [0,1]$6-dependence through an external field, not merely through asymmetric neighbor coupling (Godreche, 2011).