---
title: Quench Spectroscopy Overview
url: https://www.emergentmind.com/topics/quench-spectroscopy
type: topic
---

# Quench Spectroscopy Overview

Quench spectroscopy denotes a class of nonequilibrium spectroscopic protocols in which a system is prepared in an initial state, subjected to a sudden quench at \(t=0\), and then analyzed through the real-time evolution of observables whose oscillatory structure encodes the excitation content of the post-quench problem. Across the literature, the quench may be a global parameter change, a local perturbation, a tunneling switch, a projective measurement, or a broadband THz pulse, while the readout may be equal-time correlation functions, local observables, nonlinear currents, optical absorption edges, Loschmidt echoes, or time-dependent spectral functions [2007.08381] [2112.12123] [1405.5865].

## 1. Scope and variants of the concept

The term is used for several closely related protocols rather than a single fixed experimental design. In lattice many-body systems, quench spectroscopy usually means preparing the ground state of an initial Hamiltonian \(H(\lambda_i)\), abruptly changing a parameter to \(\lambda_f\), and reconstructing excitation energies from the spatiotemporal Fourier transform of correlation functions or one-point observables [2007.08381]. In superconducting THz experiments, “quench-drive spectroscopy” uses a short, few-cycle quench pulse and a delayed multicycle drive pulse, and the nonlinear current is analyzed in both real time and quench–drive delay [2112.12123]. In topological and impurity settings, the quench can be a sudden local tunneling switch or a single-photon absorption event, with the spectral information encoded in the Loschmidt echo or in edge singularities of optical absorption [1405.5865] [1102.3982].

A distinct methodological extension is “measurement quench,” in which the Hamiltonian is left unchanged and the quench is produced solely by measurement-induced wavefunction collapse; the ensuing unitary dynamics under the original Hamiltonian then serves spectroscopic purposes [1806.09643]. Another extension is non-local quench spectroscopy, where non-local spin-string correlators are monitored because they coincide with fermionic correlators under the Jordan–Wigner mapping; this is required when the elementary excitations are non-local in the microscopic variables [2407.14802].

| Setting | Quench | Readout |
|---|---|---|
| Lattice many-body models | Global or local parameter change | Space-time Fourier transform of observables |
| THz superconducting spectroscopy | Few-cycle quench plus delayed multicycle drive | \(j(t,\tau)\), \(j(\omega,\tau)\), \(j(\omega,\Omega)\) |
| Majorana and impurity problems | Sudden local tunneling or optical absorption | \(\mathcal{L}(t)\), threshold absorption, \(A(\omega,t)\) |

A common misconception is that quench spectroscopy is simply pump–probe spectroscopy under another name. The superconducting literature explicitly distinguishes quench-drive spectroscopy from conventional pump–probe THz and from single-pulse third-harmonic generation: the quench is broadband and short, the drive is narrowband and multicycle, and the delay is itself used as a spectral knob [2112.12123]. Another misconception is that only local probes are relevant. In one-dimensional spin chains, local probes can access only continua, whereas non-local correlators are needed to recover sharp fermionic quasiparticle lines [2407.14802].

## 2. Formal structure and signal extraction

In its most standard form, the protocol begins from the ground state of an initial Hamiltonian \(H_i\), followed by unitary evolution under \(H_f\),
\[
|\psi(t)\rangle = e^{-i H_f t / \hbar} |\psi_0(\lambda_i)\rangle .
\]
Expanding the initial state in the eigenbasis of \(H_f\) yields, for an observable \(O\),
\[
O(t) = \sum_{n,m} c_n^* c_m e^{i (E_n^{(f)} - E_m^{(f)}) t / \hbar} \langle n_f |O| m_f\rangle ,
\]
so peaks in the Fourier transform occur at frequencies \(\omega_{nm} = (E_n^{(f)} - E_m^{(f)})/\hbar\) [2507.07699]. For weak quenches dominated by transitions from the ground state, the strongest features typically occur at \(\omega \simeq (E_n^{(f)} - E_0^{(f)})/\hbar\) [2507.07699].

For spatially resolved settings, the central object is often a quench spectral function obtained by Fourier transforming in both space and time. In the local-quench formalism,
\[
G(k,\omega) = (2\pi)^{D+1} \sum_{n,n'} \rho_i^{n'n} \langle n| O |n'\rangle \delta(E_{n'} - E_n - \omega)\delta(P_{n'} - P_n - k),
\]
which makes explicit that the peaks arise from energy and momentum differences between post-quench eigenstates populated by the quench and connected by the chosen observable [2007.08381]. In this sense, quench spectroscopy is selection-rule dependent: the measured operator determines which branches are visible.

Several concrete readout formulas recur across implementations. In the dipolar XY Rydberg simulator, linear spin-wave theory gives
\[
S(k,t) \approx C_k + A_k \cos[2\omega_k t],
\]
so fitting early-time oscillations of the time-dependent structure factor directly yields \(\omega_k\) [2311.11726]. In superconducting quench-drive spectroscopy, the measured current depends on both real time \(t\) and delay \(\tau\), and the two-dimensional spectrum is
\[
S(\omega,\Omega) = \int dt \int d\tau \, e^{i(\omega t + \Omega \tau)} j(t,\tau) ,
\]
with \(\omega\) associated with harmonics of the drive and \(\Omega\) the delay frequency conjugate to \(\tau\) [2112.12123]. In boundary and topological versions, the relevant quantity can instead be the Loschmidt echo,
\[
\mathcal{L}(t) = \left|\langle \psi(0)|\psi(t)\rangle\right|^2,
\]
whose long-time asymptotics diagnose the post-quench boundary condition and low-energy fixed point [1405.5865].

These formulations share a structural feature emphasized in the broader quench-dynamics literature: the quench itself acts as the pump, while the time-resolved observable is the spectroscopic probe. This viewpoint naturally connects oscillation frequencies, light-cone propagation, dephasing, and prethermal plateaus to the excitation structure and to the stability or breakdown of generalized Gibbs descriptions [1703.09740].

## 3. Lattice many-body implementations and quantum simulators

Global and local quenches in one-dimensional lattice models provided an early systematic formulation of quench spectroscopy as a space-time Fourier analysis of nonequilibrium observables. In translationally invariant models, local quenches are sufficient to reconstruct momentum-resolved spectra from one-point observables, and the resulting quench spectral function can recover single-quasiparticle dispersions as well as energy-difference branches. For cosine-like bands \(E_k = \Delta - v \cos k\), the strongest conserving-quasiparticle features follow \(\omega(k) = \pm 2 v \sin(k/2)\) [2007.08381].

A direct experimental realization was demonstrated in a two-dimensional dipolar XY Rydberg simulator. There the system is initialized in coherent spin states aligned along \(y\), the quench is the sudden onset of dipolar XY exchange dynamics, and the equal-time \(zz\) correlations are Fourier transformed to obtain a time-dependent structure factor [2311.11726]. The ferromagnetic case exhibits elementary excitations behaving as linear spin waves, long-lived oscillations, super-ballistic correlation fronts \(d \sim t^2\), and a nonanalytic small-\(k\) dispersion \(\omega_k \propto \sqrt{|k|}\). The antiferromagnetic case instead displays \(\omega_k \propto k\), a clear light cone with \(v_g \approx J a/\hbar \approx 24 \,\mathrm{m/s}\), and strong damping attributed in the Supplemental Material to one-to-three magnon decay channels allowed by quartic Holstein–Primakoff terms [2311.11726].

Quench spectroscopy has also been used as a phase diagnostic in disordered bosonic systems. In the one-dimensional disordered Bose–Hubbard model, the momentum-resolved quench spectral function distinguishes the Mott insulator, Bose glass, and superfluid by the presence or absence of low-frequency weight and of a sound-like linear mode. Operationally, the Mott insulator is characterized by \(|G(0,0)| = 0\), whereas the Bose glass and superfluid have \(|G(0,0)| \neq 0\); the superfluid alone retains a well-defined linear lower edge \(\omega(k) \approx c_s |k|\) at small \(k\) [2105.05774]. Spatially resolved spectroscopy via local density fluctuations supplies a second diagnostic through the typical gapped-region size \(\xi\), with \(\xi=L\) in the Mott phase, \(\xi \approx 0\) in the superfluid, and \(0<\xi/L<1\) in the Bose glass [2105.05774].

The trapped-boson implementation shows how confinement modifies otherwise homogeneous quench spectra. In one-dimensional bosons in an optical lattice with harmonic trapping, the Mott-insulating spectrum retains a gap near \(\hbar\omega \approx U\) but acquires a broadened band with width \(w_{\mathrm{DMRG}} \approx 0.164\,U\) in DMRG and \(w_{\mathrm{exp}} \approx 0.180 \pm 0.019\,U\) experimentally, together with a suppression beyond \(k a \approx \pi/2\) [2507.07699]. The superfluid response remains gapless near \(k a = 0\), and the inverse quench with larger amplitude yields the clearest trapped Mott spectrum [2507.07699].

One-dimensional spin chains introduce a sharper distinction between local and non-local probes. In Jordan–Wigner-mappable models, local spin observables typically excite multiparticle continua, whereas the non-local fermionic quench spectral function
\[
Q_F(k,\omega) = 2\pi F(k)\left[\delta(\omega - 2\epsilon_k) + \delta(\omega + 2\epsilon_k)\right]
\]
yields sharp peaks at twice the fermionic quasiparticle energy [2407.14802]. This is the central result of non-local quench spectroscopy of fermionic excitations in 1D quantum spin chains. The XX-chain case gives peaks at \(\omega = \pm 4|J \cos k|\), while local probes instead display broad continua with edges characteristic of two-particle processes [2407.14802].

The same dynamical logic extends to modes that are not isolated poles in linear response. In the gapless phase of the \((1+1)\)D XXZ spin chain, nonequilibrium quench spectroscopy identifies a symmetry-preserving underdamped amplitude mode on top of a Luttinger-liquid ground state, visible in the oscillations of the \(U(1)\)-symmetric observable \(\mathcal{O}^{zz}\). At \(\Delta_f=0\), the dominant oscillation has period \(T \approx 2.5/J\), corresponding to \(\omega_0 \approx 2.5J\), and the amplitude scales with the quench strength \(|\Delta_f-\Delta_i|\) [2601.20926]. The Bethe-ansatz analysis traces the mode to specific families of string excitations, especially two 2-strings at \(k=0\) [2601.20926].

A conceptually different but experimentally economical route is the measurement quench. Here a local projective measurement changes the state while leaving the Hamiltonian unchanged; spectroscopy then proceeds from the post-measurement time trace. In Ising chains this yields oscillations governed by matrix elements \(|\langle E_n|\sigma_j^x|E_0\rangle|^2\), while in a Kondo spin-chain emulation it provides nonequilibrium scaling collapse at fixed \(N/\xi_K\) and an exponential extraction of the Kondo cloud through \(\Delta \bar m_j(J') \sim e^{-j/\xi_K}\), with \(\xi_K \sim e^{A/J'}\) and \(A \approx 0.18\) [1806.09643].

## 4. Superconducting and quench-drive spectroscopy

In superconductors, quench spectroscopy takes a particularly developed form in THz nonlinear optics. The single-band BCS formulation of quench-drive spectroscopy combines a short, broadband quench \(E_q(t)\) with a delayed, multicycle narrowband drive \(E_d(t-\tau)\), and records the nonlinear current \(j(t,\tau)\) under
\[
A(t;\tau)=A_q(t+\tau)+A_d(t).
\]
The dynamics are computed in the Anderson-pseudospin framework,
\[
\partial_t s_k = 2 b_k(t) \times s_k(t), \qquad
\Delta(t)=V\sum_k (s_k^x+i s_k^y),
\]
with the electromagnetic field entering through minimal coupling and the dominant third-order response summarized by
\[
j^{(3)}(t)=A(t)[\chi_{\rho\rho}*A^2](t)
\]
in inversion-symmetric \(s\)-wave superconductors [2112.12123].

The distinctive object is the two-dimensional spectrum \(j(\omega,\Omega)\), obtained by Fourier transforming with respect to both real time and quench–drive delay. Pure third-harmonic generation produces vertical features at \(\Omega=0\) and \(\omega=\pm\omega_d,\pm 3\omega_d\), whereas quench–drive mixing yields diagonal lines \(\omega=\pm n\omega_d-\Omega\) with \(n=0,1,2\) [2112.12123]. Along these diagonals, the signal directly tracks the nonlinear susceptibility \(\chi_{\rho\rho}(\Omega)\), allowing extraction of resonances near \(\Omega \approx 2\Delta/\hbar\) without scanning the drive frequency. The third harmonic shows \(\tau\)-dependent modulation, visible as sidebands at \(-\omega_d\), \(-2\omega_d\), and \(-3\omega_d\) in \(j(\Omega,\omega=3\omega_d)\), and polarization control can suppress or preserve specific pathways [2112.12123].

The cuprate version adapts the same logic to a clean \(d\)-wave BCS superconductor with \(\Delta(k)=\Delta_0(\cos k_x-\cos k_y)/2\). In that setting, the Higgs mode is broadened by nodal quasiparticles, but the delay-dependent third-harmonic signal still carries a peak in the modulation spectrum at \(\Omega_{\Delta t}\approx 2\Delta/\hbar\), together with sidebands near \(\omega \approx 2\Delta/\hbar \pm \omega_d\) [2203.01242]. The theory also models a phase-fluctuating pseudogap state by attaching random momentum-dependent phases to the Cooper pairs; moderate incoherence strongly suppresses THG and removes the sharp \(\Omega_{\Delta t}\approx 2\Delta/\hbar\) signature, while stronger pulses can transiently increase \(\Delta\) and partially restore phase coherence [2203.01242].

A symmetry-resolved generalization makes the role of irreducible representations explicit. In the phase-fluctuating \(d\)-wave state above \(T_c\), the equilibrium order parameter vanishes because global phase coherence is lost, yet a short quench can transiently align phases of pre-existing pairs. The resulting third-order current factorizes into a susceptibility and Raman-like vertices,
\[
\chi^{(3)}_{ijkl}=\chi^{(3)}_{\rho\rho}\gamma_{ij}\gamma_{kl},
\]
which can be decomposed into \(A_{1g}\), \(B_{1g}\), and \(B_{2g}\) channels by polarization geometry [2402.02120]. A central result is that fully incoherent pairs with an equilibrium vanishing order parameter can still produce a third-harmonic signal after a quench. In the resulting two-dimensional spectrum \(J(\omega,\nu)\), the \(\nu=0\) vertical features are drive-only, diagonal lines encode quench–drive mixing, \(B_{1g}\) dominates THG, and \(B_{2g}\) controls quench-induced gap enhancement [2402.02120].

These superconducting variants make especially clear that quench spectroscopy is not merely the time-domain counterpart of linear response. The delay axis in quench-drive spectroscopy is itself a conjugate spectral variable, and in cuprates it has been proposed as a route to detect induced or increased phase coherence in the pseudogap regime through the reappearance of a delay-locked Higgs modulation peak and revived THG [2203.01242].

## 5. Boundary, impurity, and topological realizations

A different branch of quench spectroscopy uses local boundary quenches to probe topological and impurity excitations. In the Majorana protocol, a metallic lead is suddenly coupled to the end of a topological superconducting wire. The post-quench low-energy Hamiltonian is a chiral fermion coupled to an end Majorana zero mode, and the central observable is the Loschmidt echo \(\mathcal{L}(t)\) [1405.5865]. For \(t \gg 1/T^\star\),
\[
\mathcal{L}(t) \sim t^{-\alpha}, \qquad \alpha = \frac14,
\]
with the exponent fixed by a boundary-condition change from normal reflection to Andreev reflection and by a boundary condition changing operator of dimension \(h_{\mathrm{BCC}}=1/16\) [1405.5865]. This exponent remains robust for interacting Luttinger-liquid leads with \(g>1/2\), for spinful and multichannel leads in the weak-tunneling regime, and even in the presence of an additional bound level when the Majorana coupling is nonzero [1405.5865]. In a trivial gapped phase, by contrast, the exponent is non-universal and small, \(\alpha_{\mathrm{triv}} \sim \mathrm{const}/\Delta^2\), so the quench distinguishes a Majorana end mode from a trivial subgap state in the same device geometry [1405.5865].

Optical quench spectroscopy of Kondo correlations realizes a related idea in semiconductor quantum dots tunnel-coupled to a fermionic reservoir. Absorption of a single photon on the \(X^{-}\) transition suddenly changes the local Hamiltonian by creating a doubly occupied dot and effectively turning off the spin-exchange interaction with the reservoir [1102.3982]. The absorption rate near threshold obeys
\[
I(\omega)\propto (\omega-\omega_{\mathrm{th}})^{-\eta},
\]
where \(\eta\) is set by the change of scattering phase shifts across the quench. In the Kondo-singlet to local-singlet quench, the dominant value is \(\eta \simeq 0.5\) at \(B=0\), while at \(B_{\mathrm{ext}}=1\,\mathrm{T}\) the extracted exponents are \(\eta \approx 0.66\) for the blue transition and \(\eta \approx 0.31\) for the red transition [1102.3982]. The associated edge singularity is an optical manifestation of Anderson orthogonality catastrophe, and the scaling collapse of \(A(\nu)\) versus \(\nu/T_K\) demonstrates Kondo universality in an optical spectrum [1102.3982].

Time-dependent impurity spectral functions supply another readout layer. For the Anderson impurity model after a sudden quench, the time-dependent numerical renormalization group yields retarded, advanced, lesser, and greater Green functions for several time-reference choices, and defines
\[
A(\omega,t) = -\frac{1}{\pi}\operatorname{Im} G^r(\omega,t), \qquad
N(\omega,t)=\frac{G^<(\omega,t)}{2\pi i},
\]
with \(N(\omega,t)\) directly determining time-resolved photoemission for a Gaussian probe pulse [1912.08474]. In the symmetric-to-symmetric Kondo quench considered there, the high-energy satellite peaks rearrange on the charge timescale \(1/\Gamma\), whereas the Kondo resonance develops on the much longer spin timescale \(1/T_K\) [1912.08474]. The different choices of time reference shift how the transient buildup appears, but the long-time spectral function tends to the final equilibrium one, while the occupied density of states shows imperfect thermalization characteristic of single-quench TDNRG [1912.08474].

These impurity and boundary cases clarify that quench spectroscopy need not rely on spatial Fourier transforms. The same general logic can be implemented through threshold power laws, occupied density of states, or overlap dynamics, provided the sudden local change prepares a coherent superposition that is diagnostic of the post-quench fixed point.

## 6. Extensions, limitations, and outlook

The scope of quench spectroscopy has expanded beyond isolated Hermitian systems. In a resonant Bose gas quenched to strong interactions near a Feshbach resonance, post-quench dynamics shows many-body Rabi-like oscillations between the condensate and finite-momentum quasiparticle pairs, followed by a pre-thermalized nonequilibrium steady state with a broad stationary momentum distribution, a \(k^{-4}\)-type tail through Tan’s contact, and characteristic signatures in the structure function and RF spectroscopy [1601.01397]. In this setting, quench spectroscopy is less about resolving a single dispersion branch than about diagnosing interaction strength, prethermalization times \(t_{\mathrm{pth}}(k)\approx 1/E_k\), and short-range correlations through their temporal fingerprints [1601.01397].

The method has also been generalized explicitly to dissipative and non-Hermitian lattice models. For the open Bose–Hubbard chain confined in the superfluid phase, a sudden global quench of dissipations and repulsive interactions produces equal-time correlators whose quench spectral function still reconstructs the isolated post-quench Bogoliubov dispersion, supplemented by a \(k=0\) background from the loss quench [2412.00637]. For the non-Hermitian transverse-field Ising chain in the paramagnetic phase, the quench spectral function reveals peaks at \(\omega=\pm 2\,\mathrm{Re}\,\mathcal{A}_k\), so the spectroscopy tracks twice the real part of the complex single-particle energy, while the imaginary part governs the exponential envelope [2412.00637]. This suggests that the notion of quench spectroscopy is broader than unitary Hamiltonian evolution, although the interpretation of linewidths and envelopes becomes model dependent.

Several limitations recur across implementations. Finite observation windows impose a nominal resolution \(\Delta\omega \approx 2\pi/T\), finite size broadens or splits peaks, and boundary reflections limit the useful time window in real-space reconstructions [2507.07699] [2407.14802]. In superconducting quench-drive spectroscopy, moderate fields, short or narrowband pulses, and slow relaxation underpin the quasiequilibrium interpretation, while strong-field, strongly nonequilibrium regimes may require Keldysh nonequilibrium Green’s functions [2112.12123]. In trapped bosonic systems, harmonic confinement broadens the Mott band and mixes momentum components [2507.07699]. In Majorana boundary quenches, the quench must be effectively instantaneous relative to \(1/T^\star\), the superconducting gap must remain clean, and finite Majorana splitting \(\delta_M\) restricts the available power-law window [1405.5865]. In superconducting pseudogap settings, phase incoherence can strongly suppress THG and remove the sharp delay-frequency peak [2203.01242].

Another persistent issue is operator selectivity. Local probes can be sufficient for bosonic quasiparticles or for one-point local quench spectroscopy in translation-invariant models, but they can completely miss non-local quasiparticles, as shown by the 1D fermionic spin-chain case [2007.08381] [2407.14802]. Conversely, non-local observables can be experimentally costly because they require parity strings or site-resolved readout. This tension between selectivity and accessibility is intrinsic to the method rather than incidental.

The recent literature points toward several concrete extensions already articulated in the underlying works: multi-band and unconventional superconductors with Leggett or Bardasis–Schrieffer modes, different power-law exponents and frustrated geometries in quantum simulators, higher dimensions, dissipative and non-Hermitian settings, and full pump-pump-probe or more general multidimensional coherent protocols [2112.12123] [2311.11726] [2412.00637]. A plausible implication is that quench spectroscopy is evolving from a niche post-quench diagnostic into a general framework for extracting excitation structure from controlled nonequilibrium dynamics, with the precise observable and transform chosen to match the conserved quantities, symmetry sector, and experimental platform.

Source: https://www.emergentmind.com/topics/quench-spectroscopy