Kubo-Thermalization Correspondence

• Kubo-Thermalization Correspondence is an exact relation linking short-time linear-response measurements to long-time thermalized magnetization in weakly driven quantum spin systems.
• It employs Boltzmann-weighted integrals and analytic continuation to connect measurable spectral rates to the zero crossing of steady-state magnetization.
• Experimental tests with ultracold Fermi polarons validate that spectral shifts and detailed-balance relations govern the dynamics of quantum equilibration.

The Kubo-Thermalization Correspondence denotes an exact relation between short-time linear-response spectroscopy and a long-time thermalized steady state for a weakly driven spin coupled to a thermal bath. In its explicit 2026 formulation, the correspondence states that the detuning $\Delta_0$ governing the zero crossing of the asymptotic magnetization is fixed by a Boltzmann-weighted integral transform of the short-time transition spectrum $R_\downarrow(\Delta)$, thereby linking a late-time thermalized observable to a quantity measured in the Kubo/Fermi’s Golden Rule regime (Huang et al., 7 May 2026). In a broader literature, the phrase also names looser, modified, or even failing links between response-theoretic structures, KMS analyticity, and quantum equilibration.

1. Exact statement and operative definition

In the formulation introduced for a weakly driven spin-$\tfrac12$ impurity, the system is described in the rotating frame by

$$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$

with a thermal bath Hamiltonian $\hat H^{\rm B}$ left unspecified and a spin-diagonal system-bath coupling

$$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$

The weak drive has detuning $\Delta$ and Rabi frequency $\Omega_0$; the main observable is the magnetization

$\mathcal M \equiv \langle \hat\sigma_z\rangle.$

The long-time side of the correspondence is the thermalized steady-state magnetization

$$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$

where $R_\downarrow(\Delta)$0, and $R_\downarrow(\Delta)$1 is defined by the zero crossing

$R_\downarrow(\Delta)$2

The short-time side is the linear-response transition spectrum $R_\downarrow(\Delta)$3, defined for a spin initially prepared in $R_\downarrow(\Delta)$4. After a short transient, the measured transition rate satisfies

$R_\downarrow(\Delta)$5

The exact correspondence is

$R_\downarrow(\Delta)$6

Equivalently, in units $R_\downarrow(\Delta)$7,

$R_\downarrow(\Delta)$8

This identifies the zero crossing of the long-time thermalized magnetization with a Boltzmann-weighted moment of the short-time linear-response spectrum. A related exact relation is

$R_\downarrow(\Delta)$9

which functions as a detailed-balance identity between forward and reverse spectra (Huang et al., 7 May 2026).

Two misconceptions are explicitly excluded by this formulation. First, the correspondence is not a statement that the long-time thermalized resonance is simply the spectral peak. Second, it is not a perturbative approximation around the final steady state; it remains valid even when the steady state differs substantially from the initial state, provided the stated assumptions hold.

2. Derivation from partition functions and analytic continuation

The derivation proceeds by showing that the same partition-function ratio controls both the thermalized steady state and the imaginary-time value of the short-time response function. In units $\tfrac12$0, the asymptotic magnetization is

$\tfrac12$1

Using the spin-diagonal form of the interaction, one obtains

$\tfrac12$2

Defining

$\tfrac12$3

this reduces exactly to the hyperbolic-tangent form above. The susceptibility at the zero crossing is then fixed: $\tfrac12$4

On the short-time side, the FGR spectrum is

$\tfrac12$5

and may be written as

$\tfrac12$6

with

$\tfrac12$7

It obeys the sum rule

$\tfrac12$8

The crucial analytic input is that $\tfrac12$9 is analytic in the strip

$$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$0

and continuous at the boundaries, provided the spectrum is bounded below. This allows evaluation at $$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$1: $$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$2 Using the Fourier representation,

$$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$3

which yields the correspondence.

The derivation is therefore neither a kinetic approximation nor an ETH-style asymptotic argument. Its content is exact within the stated assumptions because both the steady-state magnetization and the short-time spectrum are constrained by the same thermal trace ratio and the same analytic continuation structure (Huang et al., 7 May 2026).

3. Experimental realization in ultracold Fermi polarons

The correspondence was experimentally tested using a homogeneous ultracold gas of $$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$4 atoms in an optical box trap. The impurity spin-$$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$5 is encoded in two internal states, $$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$6 and $$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$7, while the bath is a third internal state $$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$8. The impurity fraction is

$$\hat H=\hat H^{\rm s}+\hat H^{\rm B}+\hat H^{\rm int}, \qquad \hat H^{\rm s}=-\frac12\hbar\Delta \hat\sigma_z+\frac12\hbar\Omega_0\hat\sigma_x,$$9

so spin-spin interactions are negligible and the bath is only weakly perturbed by the impurities. The bath parameters are

$\hat H^{\rm B}$0

and

$\hat H^{\rm B}$1

At $\hat H^{\rm B}$2, an rf field with detuning $\hat H^{\rm B}$3 and Rabi frequency $\hat H^{\rm B}$4 couples $\hat H^{\rm B}$5, while the bath is unaffected by the rf. The impurity-bath interactions are tuned via a Feshbach resonance with scattering lengths $\hat H^{\rm B}$6 and $\hat H^{\rm B}$7. The typical regime is one in which $\hat H^{\rm B}$8 is strongly interacting with the bath and $\hat H^{\rm B}$9 is weakly interacting.

The experiments measured both the short-time spectroscopy and the long-time steady-state magnetization. On the BCS side $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$0, the spectrum is narrow and almost symmetric, and the spectral peak $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$1 nearly coincides with $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$2. This is consistent with the limiting case

$$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$3

for which the correspondence gives $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$4. On the BEC side $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$5, the spectrum broadens strongly and $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$6 and $$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$7 differ clearly. The discrepancy grows with spectral width, and the experiments report

$$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$8

close to the Gaussian-model expectation

$$\hat H^{\rm int}=\hat n_\uparrow \hat O_\uparrow+\hat n_\downarrow \hat O_\downarrow, \qquad \hat n_{\uparrow,\downarrow}=\frac{1\pm \hat\sigma_z}{2}.$$9

This is the main empirical demonstration that the correspondence depends on the entire spectral lineshape rather than the peak position alone.

Direct use of the basic integral formula is experimentally difficult because the negative-detuning tail is exponentially amplified by $\Delta$0. For this reason, the measurements employed the symmetrized form involving both $\Delta$1 and $\Delta$2: $\Delta$3 with

$\Delta$4

where the $\Delta$5 sign is for $\Delta$6 and the $\Delta$7 sign is for $\Delta$8. Across the BCS-BEC crossover, the $\Delta$9 inferred from this symmetrized spectroscopy agreed closely with $\Omega_0$0 extracted from the long-time magnetization.

An independent thermalization check came from the slope at the zero crossing. The measured susceptibility

$\Omega_0$1

was found constant across interaction strength and equal to $\Omega_0$2, with $\Omega_0$3 independently obtained from time-of-flight thermometry. The correspondence was also tested on the metastable repulsive polaron branch at large positive $\Omega_0$4; when thermalization within that sector occurred faster than decay to lower states, the observed $\Omega_0$5 again agreed with the narrow-spectrum expectation (Huang et al., 7 May 2026).

4. KMS structure, detailed balance, and neighboring formulations

The correspondence is tightly connected to KMS analyticity but is not reducible to a generic statement that “thermalization implies KMS.” Its exact form depends on the analytic continuation of a very specific equilibrium response function,

$\Omega_0$6

to imaginary time $\Omega_0$7. The detailed-balance relation

$\Omega_0$8

is the immediate spectral counterpart of that structure.

In the wider literature, KMS analyticity has been used in several related but distinct ways. In relativistic QFT, the usual KMS condition has been localized into a local KMS condition that is equivalent, for analytic Hadamard states of the free Klein–Gordon field on Minkowski space, to the Buchholz–Ojima–Roos notion of local thermal equilibrium; the correspondence is pointwise and is expressed in terms of the local relative-variable two-point function (Gransee, 2016). In quantum quenches from regularized boundary states,

$\Omega_0$9

the initial state can satisfy only part of the KMS criterion, while local correlators at $\mathcal M \equiv \langle \hat\sigma_z\rangle.$0 were shown to satisfy the full KMS relation with effective inverse temperature

$\mathcal M \equiv \langle \hat\sigma_z\rangle.$1

in the cases analyzed (Guo, 2017). In SU(2)-symmetric chaotic many-body systems, a fine-grained KMS relation for individual energy eigenstates can be derived from non-Abelian ETH, with finite-size corrections that scale as $\mathcal M \equiv \langle \hat\sigma_z\rangle.$2 in some regimes and can become polynomially larger in others (Noh et al., 9 Jul 2025).

These neighboring constructions clarify the status of the exact spin-bath correspondence. It is neither a purely algebraic restatement of global KMS equilibrium nor an ETH-based asymptotic property of isolated eigenstates. It is an exact response-to-thermalization identity for a driven open quantum system under a specific set of assumptions.

5. Assumptions, limitations, and cases where analogous correspondences fail or change form

The exact 2026 correspondence holds under a sharply delimited set of hypotheses: a spin-$\mathcal M \equiv \langle \hat\sigma_z\rangle.$3 or more generally finite-level system, spin-diagonal system-bath coupling, a bath at temperature $\mathcal M \equiv \langle \hat\sigma_z\rangle.$4, weak driving, thermalization to the bath temperature, and analyticity of the response function in the strip $\mathcal M \equiv \langle \hat\sigma_z\rangle.$5. The Methods also give an $\mathcal M \equiv \langle \hat\sigma_z\rangle.$6-level generalization,

$\mathcal M \equiv \langle \hat\sigma_z\rangle.$7

showing that the underlying logic is not restricted to two-level impurities. The exactness does not remove practical limitations: direct use of the unsymmetrized integral can be experimentally unstable because the negative-detuning tail is exponentially amplified, and in metastable situations the correspondence applies only on timescales where thermalization within the relevant sector outpaces decay out of it (Huang et al., 7 May 2026).

Broader work on the subject also shows that “Kubo–thermalization correspondence” is not a universal principle with a single form. In a chaotic trimer-monomer Bose-Hubbard model, the standard Kubo/LRT estimate for the diffusion coefficient $\mathcal M \equiv \langle \hat\sigma_z\rangle.$8 can fail badly even though the dynamics remains diffusive at coarse-grained level; the actual coefficient is a smaller quantum value $\mathcal M \equiv \langle \hat\sigma_z\rangle.$9 obtained from a resistor-network treatment of a sparse transition graph, and the suppression factor

$$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$0

encodes an anomalous $$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$1-dependence (Khripkov et al., 2014). For finite systems, the physically meaningful Kubo formula for transport requires explicit coupling to baths or leads in the unperturbed Liouvillian; otherwise one encounters singular degenerate-state contributions that are artifacts of applying infinite-system logic to a finite open device (Wu et al., 2010). In engineered quantum-circuit reservoirs with non-orthogonal reservoir-qubit eigenstates, a modified KMS relation

$$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$2

can support complex-balanced thermalization toward a non-Gibbs nonequilibrium steady state rather than ordinary detailed-balance relaxation (Mao et al., 8 Jan 2026).

Setting	Response–thermalization link	Status
Weakly driven spin coupled to thermal bath	$$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$3 fixed by Boltzmann transform of $$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$4	Exact under stated assumptions
Sparse chaotic Bose-Hubbard network	Thermalization rate not set by $$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$5	Standard Kubo prediction fails
Finite-size open transport system	Kubo response requires explicit baths/leads	Open-system reformulation
Engineered non-Hermitian reservoir circuit	Modified KMS yields complex-balanced stationary state	Generalized nonequilibrium version

A persistent misconception is therefore that any thermalizing dynamics should automatically admit a standard Kubo prediction for its rate. The literature cited above shows that exact relations, breakdowns, and modified correspondences all occur, depending on microscopic connectivity, openness, and the character of the reservoir.

6. Broader significance and relation to holographic and hydrodynamic thermalization

The exact spin-bath correspondence sits within a larger research landscape in which short-time response structure and long-time equilibration are often linked only indirectly. In the D1D5 CFT, large gravitational redshift is explicitly distinguished from genuine thermalization: coherent infall with large redshift is not yet thermalization, whereas resonant multi-particle processes with secular $$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$6 growth are proposed as microscopic precursors of spreading over accessible states (Hampton et al., 2019). In semiclassical large-$$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$7 $$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$8 CFT, vacuum-block thermalization is diagnosed by monodromy of a Fuchs equation, and an effective temperature is extracted from

$$\mathcal M_\infty(\Delta)=\tanh\left(\frac{\beta\hbar(\Delta-\Delta_0)}{2}\right),$$9

but no Kubo formula is derived (Vos, 2018). In time-dependent AdS/CFT more broadly, far-from-equilibrium boundary dynamics is dual to bulk collapse, while late-time relaxation is governed by black-hole quasinormal modes, which are the poles of retarded correlators and thus the natural bridge to linear-response physics (Lindgren, 2019).

At the level of effective hydrodynamics, equilibrium data can determine a restricted class of Kubo formulas without addressing thermalization dynamics at all. For relativistic normal fluids and superfluids, the equilibrium partition function or Goldstone effective action generates zero-frequency Euclidean correlators, which coincide with the retarded correlators relevant for non-dissipative Kubo formulas at $R_\downarrow(\Delta)$00; dissipative transport, by contrast, requires genuine time-dependent analysis (Chapman et al., 2013). This establishes a strong equilibrium–response correspondence, but not a theory of approach to equilibrium.

A plausible implication is that exact links between short-time spectra and long-time thermalized observables are exceptional rather than generic. The 2026 Kubo-Thermalization Correspondence is notable precisely because it delivers an exact statement in a domain where most other connections are asymptotic, approximate, model-dependent, or deliberately modified. In that restricted but rigorous sense, it provides a rare benchmark for how linear-response data can encode a nontrivial property of the eventual thermalized state (Huang et al., 7 May 2026).