---
title: Pole-Skipping in Holographic Thermal Dynamics
url: https://www.emergentmind.com/topics/pole-skipping-phenomenon
type: topic
---

# Pole-Skipping in Holographic Thermal Dynamics

Pole-skipping is the phenomenon in which a retarded thermal Green’s function \(G^R(\omega,k)\) becomes intrinsically non-unique at special complex frequency-momentum pairs \((\omega_*,k_*)\) because a pole and a zero coincide, producing a local \(0/0\) ambiguity. In holographic realizations, this boundary indeterminacy is the image of a bulk near-horizon degeneracy: the linearized field equations admit an extra ingoing solution, so the usual regularity prescription fails to fix a unique response. In the energy-density channel, the leading upper-half-plane point encodes \(\lambda_L=2\pi T\) and the butterfly velocity, while more general probe sectors exhibit infinite towers at Matsubara frequencies in the lower half-plane. Recent work has further tied pole-skipping to shockwaves, entanglement-wedge geometry, and replica manifolds, enlarging its role from a diagnostic of holographic chaos to a structural feature of thermal correlators and black-hole response [1905.12014, 2504.08139].

## 1. Definition and near-horizon mechanism

At finite temperature, a standard diagnostic is the retarded energy-density correlator
\[
G^R(\omega,k)=\langle T_{00}T_{00}\rangle^{\mathrm{ret}}_\beta(\omega,k).
\]
Pole-skipping refers to special points where the correlator takes the schematic form
\[
G^R(\omega,k)\sim \frac{N(\omega,k)}{D(\omega,k)}\to \frac{0}{0},
\]
because the would-be pole and the zero of the numerator coincide. In holographic language, this occurs when the near-horizon recursion relations lose a constraint, so ingoing boundary conditions no longer select a unique bulk solution [2504.08139, 2209.04296].

For planar black holes in ingoing Eddington–Finkelstein coordinates,
\[
ds^2=-f(r)\,dv^2+2\,dv\,dr+r^2\,dx_i\,dx^i,
\]
metric perturbations of the form \(\delta g_{MN}(v,r,x)=e^{-i\omega v+ikx^1}\delta g_{MN}(r)\) admit a horizon expansion
\[
\delta g_{MN}(r)=\sum_{n=0}^{\infty}\delta g_{MN}^{(n)}(r-r_0)^n.
\]
In the energy-density channel, the leading nontrivial equation is \(E_{vv}\). At the first special point it decouples and reduces to the homogeneous transverse equation
\[
[-\partial_i^2+m^2]\,\delta g_{vv}^{(0)}(x)=0,
\]
whose nontrivial normalizable solution fixes \(k_*^2=-m^2\). This yields the first chaos-related point
\[
\omega_*=i\,2\pi/\beta,\qquad k_*=i\,m,
\]
and defines the pole-skipping butterfly velocity
\[
v_B^{\mathrm{PS}}\equiv \frac{\Im \omega_*}{\Im k_*}.
\]
In that setup, the pole-skipping points lie on a half-line \(\omega=i\,2\pi n/\beta\), \(\Re k=0\), and the lowest nontrivial point \(n=1\) carries the Lyapunov exponent \(\lambda_L=2\pi/\beta\) [2504.08139].

The same Frobenius logic applies more broadly. For generic bulk perturbations, the two indicial exponents near the horizon take the form \(\lambda_1=0\) and \(\lambda_2=i\omega/(2\pi T)\), so whenever \(i\omega/(2\pi T)\) is an integer one must test whether the logarithmic obstruction disappears. This produces a truncated near-horizon matrix \(M^{(n)}\), and pole-skipping occurs when
\[
\omega_n=-\,i\,2\pi T\,n,\qquad \det M^{(n)}(\omega_n,k_n)=0.
\]
This mechanism yields the standard lower-half-plane Matsubara tower for scalar, vector, and tensor channels [1909.09168, 1909.10223].

## 2. Regularity, non-uniqueness, and local analytic structure

The defining issue is not only degeneration of the recursion matrix but also regularity of the resulting horizon solutions. In the SAdS\(_4\) sound channel, the generic outgoing mode is singular at the future horizon in the upper-half \(\omega\)-plane: a gauge-invariant curvature invariant such as the linearized Kretschmann scalar remains finite for the ingoing branch but diverges for the would-be outgoing branch. At the special point, however, the regular singular point becomes an ordinary point and both independent solutions are regular. The ingoing mode therefore ceases to be uniquely defined precisely at the pole-skipping point [1905.12014].

From the boundary perspective, this produces a multi-valued retarded Green’s function. Approaching \((\omega_*,k_*)\) along different slopes in the complex \((\omega,k)\)-plane gives different limiting values. In the Rarita–Schwinger case, the local structure takes the standard form
\[
G^R(\omega_*+\delta\omega,k_*+\delta k)\sim \frac{\delta\omega-v_z\,\delta k}{\delta\omega-v_p\,\delta k},
\]
with \(v_{p,z}\) determined by the next-order horizon equations [2101.01490]. This slope dependence is the operational form of the \(0/0\) ambiguity.

A recurring structural feature is that the special frequencies are more rigid than the special momenta. Finite-coupling and higher-curvature analyses show that the Matsubara frequencies \(\omega_n=-i\,2\pi T\,n\) are not explicitly shifted by Gauss–Bonnet, \(R^2\), or \(R^4\) corrections, while the momenta \(q_n\) or \(k_n\) do shift. In some sectors, special points can even disappear at critical higher-derivative couplings; in the AdS\(_5\) tensor channel, this happens at \(\lambda_{\mathrm{GB}}=-1/8\) [1909.09168, 1909.10223]. This separation between protected \(\omega_n\) and theory-dependent \(k_n\) has become a central feature of the subject.

## 3. Relation to chaos, shockwaves, and entanglement

The original prominence of pole-skipping came from the energy-density correlator, where the leading upper-half-plane point is directly tied to many-body chaos:
\[
\omega_*=i\,\lambda_L,\qquad k_*=i\,\lambda_L/v_B.
\]
This identifies the Lyapunov exponent and butterfly velocity without computing an out-of-time-order correlator explicitly [1905.12014].

A deeper unification is furnished by the relation among pole-skipping, OTOC shockwaves, and entanglement-wedge reconstruction. In the late-time, large-region regime, the extremal HRT surface lies infinitesimally shifted off the bifurcation surface along the horizon. In Kruskal coordinates, the limiting surface can be written as
\[
U=0,\qquad V=F(x),\qquad F(x)\simeq e^{\nu |x|},
\]
and the extremal-surface equation reduces in the near-horizon limit to
\[
[-\partial^2+m^2]\,F(x)=0.
\]
The pole-skipping mode in Kruskal gauge has the form
\[
\delta ds^2_{\mathrm{PS}}=G(x)\left(\frac{dU^2}{U}+\cdots\right),\qquad G(x)=F(x),
\]
so the \(n\to 1\) replica geometry controlling late-time entanglement is precisely the pole-skipping mode. After regulating \(U\to U-i\epsilon\),
\[
\Im \frac{1}{U-i\epsilon}=\pi\,\delta(U),
\]
and the imaginary part of the pole-skipping perturbation becomes the standard gravitational shockwave computing the OTOC. Since the shockwave profile, the limiting HRT surface, and the pole-skipping mode obey the same transverse equation, one obtains the triple equality
\[
v_B^{\mathrm{PS}}=v_B^{\mathrm{OTOC}}=v_B^{\mathrm{EW}}
\]
in Einstein gravity and, more generally, in higher-derivative \(f(\mathrm{Riemann})\) theories with minimally coupled matter whenever the standard pole-skipping mode exists [2504.08139].

This coincidence is not confined to static isotropic branes. In the simply-spinning Myers-Perry–AdS\(_5\) plasma, there is a general proof of pole-skipping in the retarded energy-density Green’s function whenever the spatial profile of the perturbation satisfies the same shockwave equation that governs the OTOC. In the large-black-hole limit, the associated Lyapunov exponents and butterfly velocities can be extracted analytically, and the sound-mode dispersion relations numerically pass through the predicted pole-skipping locations [2211.00016].

## 4. Fields, symmetries, mass, and duality

Pole-skipping extends far beyond metric perturbations. For a minimally coupled Rarita–Schwinger field in AdS–Schwarzschild, the near-horizon \((\psi_v,\psi_r)\) subsystem becomes degenerate at
\[
\omega_*=i\,\pi T=i\,\lambda_L/2,\qquad k_*=-\,i\,m\,r_0.
\]
At that point, two linearly independent regular ingoing solutions exist. The same analysis suggests a spin-dependent pattern
\[
\omega_{*,s}\simeq 2\pi i T\,(s-1),\qquad s=0,1,3/2,2.
\]
For \(s=3/2\), the positive imaginary frequency is not a genuine instability, because \(k_*\) is also complex and the ingoing prescription is still imposed [2101.01490].

For massless \(U(1)\) \(p\)-forms in thermal AdS/CFT, the first-order pole-skipping points occur at
\[
\omega=-2\pi iT,\qquad
k^2=\pm\,\pi T\bigl[(d-2p)\,h'(r_0)+2h(r_0)\,Z'(r_0)/Z(r_0)\bigr],
\]
with the upper sign corresponding to the longitudinal mode and the lower sign to the transverse mode. The dependence on \(d-2p\) leads to a trans-mode equivalence under electromagnetic duality: the longitudinal sector of a \(p\)-form is mapped to the transverse sector of its \((d-p)\)-form dual, and the pole-skipping loci agree [2209.04296].

Mass terms enlarge the pattern. In massive \(p\)-form theories, extra first-order points appear, the massless limit separates into the massless longitudinal \(p\)-form branch and the massless transverse \((p-1)\)-form branch, and there is an additional zeroth-order point at
\[
\omega=0,\qquad k^2=-m^2 h(r_0).
\]
A complementary Stueckelberg analysis shows why the spectrum doubles: breaking gauge invariance promotes the would-be gauge sector to dynamical status, the number of skipped poles doubles, and the extra pole-skipping points are associated with the Stueckelberg fields. As the mass varies, some skipped-pole wave numbers move from a non-physical region with complex \(q\) to a physical region with real \(q\) [2209.04296, 2404.17354].

These examples establish that pole-skipping is sensitive both to the kinematics of horizon regularity and to the operator content of the bulk theory. The frequencies are often universal, but the momenta can depend on spin, mass, duality frame, and symmetry-breaking data.

## 5. Beyond static isotropic black branes

The phenomenon is not restricted to standard AdS black branes. In acoustic black holes, the analogue of pole-skipping appears for scalar perturbations at
\[
\omega_n=-\,i\,2\pi T\,n,\qquad n=1,2,\dots,
\]
with the first point at \(\omega_*=-i\,2\pi T\). In that setting, the lower-half-plane pole-skipping frequencies track the imaginary parts of quasinormal modes and encode dissipation and instability timescales of the fluid perturbations [2110.08074].

Pole-skipping also occurs without a horizon. In the AdS\(_5\) soliton, regularity at the tip replaces the ingoing horizon condition. For scalar, Maxwell, and gravitational perturbations with compact-circle momentum \(q_z\), the pole-skipping values are
\[
q_z=(s-1)-n,\qquad n=0,1,2,\dots,
\]
and at each such value one would-be normal mode is absent. The resulting “missing states” resolve otherwise puzzling features of the soliton spectrum and show that pole-skipping is not exclusively a black-hole phenomenon [2307.11178].

In holographic axion theories, the use of gauge-invariant variables reveals two distinct classes: regular pole-skipping points, where the near-horizon matrix remains finite and \(\det \mathcal M_n=0\), and singular pole-skipping points, where some entries diverge and a singular expansion is required. The lower-half-plane points are regular, whereas the leading spin-0 chaos point is singular. The same framework also produces a finite-\(\beta\) singular shear point at \((\omega_*,k_*)=(0,\pm i\bar\beta)\), and shows that if the pole-skipping momentum is purely real or purely imaginary at \(\mu=\beta=0\), it retains that character for \(\mu\neq 0\) and \(\beta\neq 0\) [2402.12951].

Spatial anisotropy and rotation deform, but do not erase, the structure. In anisotropic plasma, the sound-channel upper-half-plane point remains at \(\omega=+i\,2\pi T\), still yields \(\lambda_L=2\pi T\), and produces anisotropy-dependent corrections to \(v_B\); in the shear channel, the analytically continued momentum-diffusion branch passes through the first three pole-skipping points [2012.07710]. For rotating BTZ black holes, exact Green’s functions for spin-\(1/2\) and spin-1 fields give a full tower of pole-skipping points, and in the non-extremal case the leading frequency is
\[
\omega_{\mathrm{leading}}=\frac{2\pi i\,T_h\,(s-1+\nu\Omega)}{1-\Omega^2}.
\]
In the extremal case, a leading point survives only when \(\nu=s+1\), with
\[
\omega_{\mathrm{leading}}^{\mathrm{extremal}}=-\,2\pi i\,T_R\,(s+1)
\]
[2306.14805].

Pole-skipping has also been used as a probe of phase structure. In the 1RCBH model, the butterfly velocity probes critical behavior, pole-skipping occurs at the chaos point, and one finds \(v_B^2\ge c_s^2\) throughout the parameter space studied for both 1RCBH and AdS–RN [2305.00298]. In the Einstein–Maxwell–Chern–Simons model with a quantum critical point, the quantities
\[
O(n)=\sum_{i=1}^{2n}|\Re\,k_{n,i}|
\]
vanish in the disordered state and become positive in the ordered state, providing order parameters built directly from pole-skipping data [2307.16716].

## 6. Algebraic structures and spectral interpretations

A recent development is the recognition that pole-skipping points possess a hidden algebraic organization. For Klein–Gordon-type equations \((\nabla^2+V(r))\phi=0\), the determinant condition at level \(n\) is a degree-\(n\) polynomial in \(\mu\equiv k^2\), but only a subset of the \(\mu_{n,q}\) is independent: the rest satisfy universal homogeneous polynomial identities. In this framework, the full static planar black-hole metric can be reconstructed analytically from the infinite set of pole-skipping data \(\{(\omega_n,\mu_{n,q})\}\) by solving linear equations for the near-horizon coefficients, and even the near-horizon vacuum Einstein equations can be rewritten entirely in terms of the symmetric polynomials of the pole-skipping momenta [2507.13306].

This algebraic perspective clarifies a common misconception: pole-skipping is not merely a statement about poles. Zeros are equally essential. That point is especially clear in Kerr perturbation theory near algebraically special frequencies. There, anomalous quasinormal-mode behavior is explained by tracking poles and zeros of the Teukolsky Green-function building blocks across different Riemann sheets. The apparent disappearance of a Kerr quasinormal mode is due to pole-skipping: a quasinormal-mode pole is cancelled by a Matsubara-mode zero, producing a \(0/0\) structure. The accompanying bifurcation is resolved as an avoided crossing with resonant excitation rather than as a literal splitting of a single mode [2605.17840].

Related channel structures arise already in four-dimensional massive black holes with maximally symmetric horizons. There, pole-skipping points separate into algebraically special points and a second set common to the even and odd gravitational channels, with the Darboux transformation providing the relation between the two sectors. In that classification, the \(n=-1\) even-channel solution is the chaos point with \(\omega=i\lambda_L\) [2303.15921].

Across these developments, a consistent picture has emerged. Pole-skipping is determined by IR regularity data at the horizon or its analogue, survives substantial deformations of the bulk dynamics, constrains the analytic continuation of hydrodynamic and nonhydrodynamic response, and in the energy-density channel furnishes a direct bridge between retarded two-point functions, shockwaves, and entanglement structures [2209.04296, 2504.08139].

Source: https://www.emergentmind.com/topics/pole-skipping-phenomenon