Exponential Spectral Scaling (ESS)

• Exponential Spectral Scaling is a concept defining exponential structure in spectral quantities with field-specific interpretations, such as subradiant decay rates in disordered waveguide QED.
• It captures the crossover from power-law to exponential decay induced by disorder, highlighting Anderson localization effects through finite-size scaling and critical exponents.
• ESS underpins practical methods in stellarator optimization and Mellin spectral theory, enabling efficient analysis of spectral behavior and multiscale dynamics.

Exponential Spectral Scaling (ESS) is a term used in several recent research contexts to denote exponential structure in spectral quantities, but the object being scaled differs by field. In disordered waveguide QED, ESS denotes the exponential system-size dependence of the typical subradiant decay rates, $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$, induced by positional disorder and identified with Anderson localization of subradiant eigenstates (Tian et al., 4 Apr 2026). In stellarator boundary optimization, ESS denotes a mode-dependent exponential rescaling of boundary Fourier coefficients that compresses mode-amplitude disparity and enables single-stage full-spectrum optimization (Jang et al., 19 Sep 2025). In Mellin spectral theory, the acronym is not introduced explicitly, but the explicit scale-invariant kernel $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$ realizes an exponential dependence in log scale whose Mellin spectrum is Lorentzian, thereby separating geometric from spectral scaling exponents (Jacobs et al., 1 Jun 2026). Related formulations appear in finite-size scaling with exponential basis sets for quantum criticality and in the exponential scaling limit of Anderson localization (Alharbi et al., 2013, Chulaevsky, 2015).

1. Range of meanings

The phrase has acquired a field-dependent meaning rather than a single universal definition. In each setting, ESS refers to exponential structure in a spectral representation, but the relevant variables, observables, and mechanisms are different.

Context	Quantity or variable	ESS meaning
Waveguide QED	Typical subradiant decay rates	$$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$
Mellin spectral theory	Ratio-dependent kernel shape / Mellin multiplier	$F(t)=c e^{-\sigma|\ln t|}$ gives Lorentzian $\widetilde F(\omega)$
Stellarator optimization	Boundary Fourier coefficients	$$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$

This divergence of usage matters technically. In waveguide QED, ESS is a statement about the disorder-averaged logarithm of decay rates rather than their arithmetic mean. In stellarator optimization, ESS is a variable-scaling strategy rather than a physical law. In Mellin spectral theory, ESS characterizes exponential dependence in $|\ln t|$, not an exponential law in the original coordinate itself. A related caution is that some adjacent papers do not use the acronym explicitly, even when their constructions are naturally interpreted as exhibiting exponential spectral scaling.

2. ESS in disordered waveguide QED

In a one-dimensional waveguide QED array, an effective non-Hermitian Hamiltonian governs collective radiative dynamics in the single-excitation sector,

$$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$

with decay rates $\Gamma_k=-\Im\omega_k$. For an ordered chain $x_m=md$, the subradiant eigenstates are extended standing waves $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$0, and the clean finite-size baseline is a power law: strong subradiant states near $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$1 obey $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$2, while weak subradiant states with $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$3 obey $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$4 (Tian et al., 4 Apr 2026).

Disorder is introduced through random positions $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$5 with $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$6. Because the resulting decay-rate distribution is broad and skewed, the relevant object is the typical decay rate

$F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$7

rather than the mean $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$8. For finite disorder $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$9, $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$0 exhibits a crossover: at small $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$1 the clean power law survives, whereas for $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$2 it crosses into

$$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$3

This exponential-in-$$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$4 regime is what the work identifies as ESS. Figures 2(a,b) show the crossover for both strong and weak subradiant sectors, while the mean decay rate remains algebraic, $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$5, even in disorder. One common misconception is therefore excluded directly by the formulation: ESS is not a statement about the mean spectrum, but about the typical spectrum.

The same work terms the crossover from clean power-law scaling to exponential scaling the subradiant scaling transition (SST). In the thermodynamic limit, any infinitesimal positional disorder drives the typical subradiant sector into ESS, whereas superradiant modes near $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$6 do not display the corresponding critical behavior.

3. Criticality, finite-size scaling, and the subradiant scaling transition

To characterize the SST quantitatively, the analysis introduces a finite-size characteristic scale

$$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$7

For pure exponential scaling, $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$8, this $$\Gamma_k^{\rm typ}(N,W)\sim \exp[-N/\xi_\infty(W,k)]$$9 approaches $F(t)=c e^{-\sigma|\ln t|}$0 as $F(t)=c e^{-\sigma|\ln t|}$1. Numerically, $F(t)=c e^{-\sigma|\ln t|}$2 saturates rapidly in strong disorder, while in weak disorder it first grows $F(t)=c e^{-\sigma|\ln t|}$3 because of the residual power-law regime and then saturates once $F(t)=c e^{-\sigma|\ln t|}$4. In the clean limit, $F(t)=c e^{-\sigma|\ln t|}$5, showing that ESS disappears at $F(t)=c e^{-\sigma|\ln t|}$6 (Tian et al., 4 Apr 2026).

This behavior motivates the interpretation of $F(t)=c e^{-\sigma|\ln t|}$7 as a correlation-length-like quantity,

$F(t)=c e^{-\sigma|\ln t|}$8

and the associated finite-size scaling form

$F(t)=c e^{-\sigma|\ln t|}$9

Data collapse yields $\widetilde F(\omega)$0 for both strong and weak subradiant sectors, with $\widetilde F(\omega)$1 for strong subradiant states and $\widetilde F(\omega)$2 for weak subradiant states. The critical point is therefore at zero disorder strength,

$\widetilde F(\omega)$3

so the clean power-law scaling is infinitely fragile: any $\widetilde F(\omega)$4 ultimately drives ESS in the thermodynamic limit. Superradiant modes do not show the same finite-size scaling behavior; their cost function $\widetilde F(\omega)$5 remains large, indicating no SST.

The critical interpretation is significant because it places the onset of ESS within standard finite-size scaling methodology rather than presenting it as a purely numerical crossover. The transition is framed as a genuine disorder-driven critical phenomenon in the typical decay-rate sector.

4. Localization mechanism and boundary-radiation picture

The physical mechanism of ESS in the waveguide-QED problem is exposed by analyzing the inverse effective Hamiltonian,

$\widetilde F(\omega)$6

where $\widetilde F(\omega)$7 becomes a nearest-neighbor tight-binding chain with random parameters $\widetilde F(\omega)$8, and $\widetilde F(\omega)$9 places imaginary potentials only at the boundaries. This mapping converts the original infinite-range dipole interaction into a disordered one-dimensional chain with boundary radiation channels (Tian et al., 4 Apr 2026).

Using the spectral decomposition of $$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$0, the decay rate of a subradiant mode is related to the boundary population: $$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$1 for $$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$2. Positional disorder makes the effective tight-binding parameters random, and the subradiant eigenstates undergo Anderson localization,

$$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$3

Since emission occurs effectively through the boundaries in the inverse-Hamiltonian picture, a localized bulk excitation must tunnel to the chain ends. The boundary weight is therefore exponentially small in the system size, yielding

$$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$4

Averaging over the localization center $$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$5 gives the typical law

$$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$6

which the paper summarizes as the universal exponential scaling $$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$7 for both classes of subradiant states. This identifies the spectral scale extracted from ESS with the localization length up to a factor of two,

$$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$8

Finite-size scaling of $$x_{\text{scaled}(m,n)}=x_{\text{original}(m,n)}\exp[\alpha g(m,n)]$$9, extracted via the participation ratio $|\ln t|$0, yields $|\ln t|$1 for strong and $|\ln t|$2 for weak subradiant states, matching the exponents obtained from the spectral scale. After conformal mapping, the data for $|\ln t|$3 and $|\ln t|$4 collapse onto the same universal curve. In this sense, ESS is the spectral fingerprint of Anderson localization in the subradiant sector.

5. Mellin spectral theory and decoupled scaling exponents

A distinct use of the ESS idea appears in the spectral theory of scale-invariant operators on the multiplicative half-line $|\ln t|$5. There, a symmetric kernel satisfying

$|\ln t|$6

factorizes as

$|\ln t|$7

where $|\ln t|$8 is the geometric exponent and $|\ln t|$9 is a ratio-dependent shape function. The Mellin transform diagonalizes the associated operator, with generalized eigenfunctions

$$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$0

and eigenvalues given by the Mellin multiplier

$$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$1

The central message is the decoupling of the geometric exponent $$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$2 from an effective spectral exponent $$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$3 extracted from finite-dimensional truncations $$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$4 (Jacobs et al., 1 Jun 2026).

The explicit kernel

$$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$5

is exponential in log scale. Its Mellin multiplier is

$$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$6

a Lorentzian of width $$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$7. The paper does not introduce a separate acronym “ESS,” but this example is precisely an exponential spectral scaling mechanism in the sense that exponential dependence in $$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$8 produces a non-power-law spectral line shape in Mellin space. The high-frequency tail is $$\hat{H}_{\rm eff} = -\frac{i\gamma}{2}\sum_{m,n=1}^{N} e^{i|x_m-x_n|k_0}\hat{\sigma}^{\dagger}_m\hat{\sigma}_n,$$9, but the spectral envelope is not scale-free; it carries the finite width $\Gamma_k=-\Im\omega_k$0.

This construction leads to a multicritical interpretation. The geometric exponent $\Gamma_k=-\Im\omega_k$1 controls the dilation law of the kernel, while the effective spectral exponent $\Gamma_k=-\Im\omega_k$2 depends on the Mellin line shape and the sampling of discrete frequencies in finite truncations. Hence $\Gamma_k=-\Im\omega_k$3 generically. Equality $\Gamma_k=-\Im\omega_k$4 corresponds to a simple critical fixed point of the renormalization group, whereas $\Gamma_k=-\Im\omega_k$5 signals multiple independent scaling dimensions. The same work also proves that exact discrete self-similarity on the lattice forces eigenvector collapse to rank one, which motivates the continuum Mellin formulation and clarifies why finite-size spectra must be interpreted as sampled continuum Mellin structure rather than literal lattice self-similarity.

6. ESS as mode-dependent scaling in stellarator optimization, and related formulations

In stellarator boundary optimization, ESS has a directly algorithmic meaning. Plasma boundaries are represented by double Fourier coefficients $\Gamma_k=-\Im\omega_k$6 and $\Gamma_k=-\Im\omega_k$7, whose magnitudes in optimized configurations exhibit near-exponential decay with increasing mode numbers. This creates a large mode amplitude disparity,

$\Gamma_k=-\Im\omega_k$8

typically $\Gamma_k=-\Im\omega_k$9, so low-order modes dominate nonlinear least-squares steps and high-order modes become numerically underweighted. ESS addresses this by rescaling each mode according to

$x_m=md$0

with the primary implementation using the $x_m=md$1 norm,

$x_m=md$2

which gives a square spectral decay profile in $x_m=md$3-space (Jang et al., 19 Sep 2025).

This rescaling compresses the dynamic range from $x_m=md$4 to $x_m=md$5, aligns with the natural spectral decay of physically meaningful configurations, and permits direct single-step optimization using the full Fourier spectrum. It replaces traditional Fourier continuation, which optimizes low modes first and introduces higher modes in multiple stages. Benchmark results on quasi-axisymmetric and quasi-helically symmetric configurations, using DESC and SIMSOPT, show that ESS eliminates arbitrary staging decisions, smooths the loss history, reduces sensitivity to initial conditions, avoids distorted or self-intersecting local minima, and reduces wall-clock time by a factor of $x_m=md$6 to $x_m=md$7. For a fixed QA benchmark with target objective $x_m=md$8, the reported times are $x_m=md$9 min for Fourier continuation and $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$00, $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$01, and $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$02 min for ESS with $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$03, $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$04, and $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$05, respectively. The recommended default is $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$06 with $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$07. A key conceptual clarification in this literature is that ESS is not a physics penalty and does not regularize by suppressing high modes; it is purely a variable scaling, so high-mode amplitudes may still grow if the optimizer finds them useful.

Two adjacent lines of work are closely related but terminologically distinct. In the finite-size scaling analysis of quantum criticality for the Hulthen potential, a meshfree spectral method uses exponential basis functions $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$08 with $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$09, $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$10 linearly spanning $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$11, so the basis covers eight orders of magnitude in decay rates. There the spectrum $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$12 is analyzed as a function of basis size $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$13, leading to estimates $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$14, $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$15, and $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$16; the paper does not use the ESS acronym, but it combines exponential basis structure with finite-size scaling of spectral observables (Alharbi et al., 2013). In the lattice Anderson model, a different but cognate notion appears as the exponential scaling limit, defined by

$F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$17

which encodes decay $F(t)=c\rho^{|\ln t|}=c e^{-\sigma|\ln t|}$18 for eigenfunction correlators and Green-function singularity probabilities. That work reformulates bootstrap multi-scale analysis as adaptive feedback scaling and proves such asymptotically exponential decay even for marginal disorder distributions weaker than any fixed Hölder regularity (Chulaevsky, 2015).

Taken together, these usages show that ESS is not a single canonical formalism but a family of exponential-scaling ideas attached to spectral objects: decay-rate spectra in open quantum systems, Mellin spectra of scale-invariant kernels, optimization variables in Fourier boundary representations, and closely related exponential-scaling limits or exponential-basis spectral analyses in localization and quantum criticality.