---
title: Interacting Anderson Quantum Sun Model
url: https://www.emergentmind.com/topics/interacting-anderson-quantum-sun-model
type: topic
---

# Interacting Anderson Quantum Sun Model

The interacting Anderson Quantum Sun (AQS) model is a disordered many-body spin system in which a localized spin-\(1/2\) chain is coupled to a small, fully ergodic “sun” through exponentially decaying, \(U(1)\)-symmetric XY interactions. It was introduced to probe how an ergodic inclusion destabilizes Anderson localization and to test avalanche-based scenarios for delocalization. In addition to standard localized and ergodic phases, the model displays two nonstandard regimes: one with volume-law entanglement coexisting with intermediate spectral statistics, and another with Poisson level statistics together with sub-volume entanglement and rare-event-dominated systemwide correlations, interpreted as ergodic instabilities [2507.18286].

## 1. Model architecture and defining scales

The AQS Hamiltonian decomposes into chain, bath, and coupling terms,
\[
H = H_{\text{chain}} + H_{\text{bath}} + H_{\text{coupling}}.
\]
The localized chain is a disordered spin-\(1/2\) chain with open boundary conditions,
\[
H_{\text{chain}} = \sum_{j=1}^L h_j S_j^z + J \sum_{j=1}^{L-1} (S_j^x S_{j+1}^x + S_j^y S_{j+1}^y),
\]
where the fields are i.i.d. uniform random variables \(h_j \in [1-W,1+W]\) with \(W=0.5\), and \(J \ge 0\) controls intra-chain hopping. The model preserves the \(U(1)\) symmetry \(S^z_{\mathrm{tot}} = \sum_{\text{all spins}} S^z\) [2507.18286].

The ergodic bath, or “sun,” contains \(N=3\) spins and is represented by a Gaussian Orthogonal Ensemble random matrix restricted to conserve \(S^z_{\mathrm{tot}}\),
\[
H_{\text{bath}} = R_{\mathrm{QS}}, \qquad
R_{\mathrm{QS}} = (\beta/2)(B+B^T),
\]
with \(\beta = 0.3\), \(B_{ij}\) i.i.d. standard Gaussian, and all entries that violate \(U(1)\) symmetry set to zero. The chain and bath are connected through exponentially decaying XY couplings,
\[
H_{\text{coupling}} = \sum_{j=1}^L \alpha^{u_j}\big[S_{n_j}^x S_j^x + S_{n_j}^y S_j^y\big],
\]
where \(n_j \in \{1,\dots,N\}\) is a randomly chosen sun-spin index, and \(u_j \in [j-\zeta,j+\zeta]\) with \(\zeta=0.2\) sets the decay profile [2507.18286].

Two intrinsic length scales organize the physics. The first is the sun–chain coupling decay length,
\[
\Lambda_\alpha = \frac{1}{|\ln \alpha|},
\qquad \alpha^u = e^{-u/\Lambda_\alpha}.
\]
The second is the average Anderson localization length of the noninteracting chain,
\[
\xi_{\mathrm{AL}} = \frac{1}{\ln\!\left[1+\left(\frac{W}{W_0 J}\right)^2\right]},
\qquad W_0 \approx 1.2.
\]
Their competition controls whether localization survives, becomes unstable, or is replaced by thermal behavior. This suggests that the adjective “Anderson” refers to the localized-chain sector that seeds the model’s nonergodic behavior [2507.18286].

## 2. Numerical framework and observables

The main numerical analysis is performed in the largest \(U(1)\) symmetry sector, \(S^z_{\mathrm{tot}}=0\) for even \(N+L\) and \(S^z_{\mathrm{tot}}=-1/2\) for odd \(N+L\). Energies are taken from a mid-spectrum window \(\varepsilon \in [0.5-\varepsilon_0,\,0.5+\varepsilon_0]\) with \(\varepsilon_0=0.005\), where \(\varepsilon\) is rescaled to \([0,1]\), and at least \(30\) eigenvalues are retained per realization. Full exact diagonalization is used up to \(N+L \le 13\), and POLFED extends the calculations to \(N+L \le 20\). Disorder averaging ranges from \(10^4\) realizations for \(L=7\) to \(10^3\) for the largest sizes; finite-size critical points \(\alpha_L^c\) are extrapolated using
\[
\alpha_L^c = a/L + \alpha_\infty^c
\]
[2507.18286].

Three classes of diagnostics are central. The first is spectral statistics, based on level spacings \(s_i = E_{i+1}-E_i\) and adjacent-gap ratios
\[
r_i = \frac{\min(s_i,s_{i+1})}{\max(s_i,s_{i+1})}.
\]
The corresponding benchmarks are \(\langle r\rangle \approx 0.386\) for Poisson and \(\langle r\rangle \approx 0.530\) for GOE statistics. Intermediate values of \(\langle r\rangle\), stable under finite-size extrapolation, are treated as evidence for a genuinely non-GOE spectral regime rather than crossover noise [2507.18286].

The second class is entanglement diagnostics. For a bipartition \(A\) containing the sun plus the \(\lfloor (N+L)/2\rfloor-N\) chain spins closest to the sun, the von Neumann entropy is
\[
S_A = -\mathrm{Tr}(\rho_A \ln \rho_A),
\qquad
\rho_A = \mathrm{Tr}_B |\psi\rangle\langle\psi|.
\]
It is normalized by the Page entropy \(S_P\), and its system-size dependence is parameterized as \(S(L)\propto L^b\), with \(b=0\) for area law, \(0<b<1\) for sub-volume scaling, and \(b=1\) for volume law. Complementary information comes from participation entropies
\[
P_q = \frac{1}{1-q}\ln \sum_{\alpha=1}^{\mathcal D} |\psi_\alpha|^{2q},
\qquad
P_1 = -\sum_\alpha |\psi_\alpha|^2 \ln |\psi_\alpha|^2,
\]
with scaling \(P_q = D_q \ln \mathcal D + b_q\). In the AQS, \(D_1=1\) in the volume-law phases, while \(0<D_q<1\) in the localized and sub-volume regimes [2507.18286].

The third class concerns spatial correlations and rare events. Connected correlators are defined by
\[
C_{ij}^{\gamma\gamma}
=
\big|\langle S_i^\gamma S_j^\gamma\rangle - \langle S_i^\gamma\rangle\langle S_j^\gamma\rangle\big|,
\qquad \gamma=x,z,
\]
with typical decay extracted from
\[
\overline{\ln C_{ij}^{\gamma\gamma}} = -|i-j|/\xi_\gamma + O(1).
\]
A rescaled observable,
\[
\ell_{ij}^{\gamma\gamma} = -\ln|4C_{ij}^{\gamma\gamma}|,
\]
is used to expose fat-tailed distributions at small \(\ell\), which signal rare-event-dominated correlations. The ratio \(\xi_x/\xi_z\) provides an additional instability indicator: in the AQS sub-volume regime it decreases with \(L\), unlike the \(\xi_x/\xi_z \to 2\) behavior reported for Anderson or strongly disordered XXZ chains [2507.18286].

The study is deliberately static. Thermalization is inferred from ETH-consistent eigenstate diagnostics rather than direct quench dynamics, although supplementary analyses involving Thouless time and spectral form factor are mentioned [2507.18286].

## 3. Phase structure and anomalous regimes

At fixed disorder \(W=0.5\), the phase diagram in the \((J,\alpha)\) plane contains four regimes, labeled \([A]\) through \([D]\) [2507.18286].

| Regime | Spectral statistics | Entanglement scaling |
|---|---|---|
| [A] | Poisson | Area law, \(b=0\) |
| [B] | Poisson | Sub-volume, \(0<b<1\) |
| [C] | Intermediate | Volume law, \(b=1\) |
| [D] | GOE | Volume law, \(b=1\) |

Regime \([A]\) is the conventional localized phase: Poisson spectrum, area-law entanglement, and no indication of extensive thermalization. Regime \([D]\) is the conventional ergodic phase: GOE level statistics and volume-law entanglement consistent with ETH. The AQS departs from the standard MBL/ETH dichotomy through regimes \([B]\) and \([C]\). In \([B]\), Poisson statistics persist while eigenstates already display sub-volume entanglement and rare-event-dominated systemwide correlations. In \([C]\), the eigenstates have volume-law entanglement and \(D_1=1\), yet \(\langle r\rangle\) remains strictly between Poisson and GOE [2507.18286].

For small \(J\), corresponding to strong disorder \(W/J \gtrsim 1.7\), the spectral and eigenstate transitions nearly coincide at \(\alpha^c \approx 0.75\), as in the \(U(1)\) Quantum Sun limit \(J=0\). For larger \(J\), a separation opens between the eigenstate transition \(\alpha_{S,\infty}^c\) and the spectral transition \(\alpha_{r,\infty}^c\): volume-law entanglement appears already at \(\alpha \gtrsim \alpha_{S,\infty}^c\), whereas GOE-like level repulsion only emerges at the larger scale \(\alpha \gtrsim \alpha_{r,\infty}^c\). This creates regime \([C]\), where eigenstate thermalization and spectral chaos are decoupled [2507.18286].

Inside the Poisson sector, the crossover between \([A]\) and \([B]\) is summarized by an empirical instability line,
\[
\Lambda_\alpha + z\,\xi_{\mathrm{AL}} = \zeta^\ast,
\qquad
z \approx 2,\quad \zeta^\ast \approx 3.7.
\]
This marks the onset of sub-volume entanglement and ergodic instabilities while the spectrum remains Poisson. A plausible implication is that the AQS supports an extended precursor regime in which the bath seeds long-range resonances before producing global spectral repulsion [2507.18286].

## 4. Rare events, avalanches, and the decoupling of ETH diagnostics

The physical interpretation centers on the competition between \(\Lambda_\alpha\), which measures how far the bath effectively reaches into the chain, and \(\xi_{\mathrm{AL}}\), which measures how readily the Anderson-localized chain can hybridize. Increasing either parameter promotes delocalization. The empirical line \(\Lambda_\alpha + z\xi_{\mathrm{AL}}=\zeta^\ast\) is consistent with the onset of rare resonant spots seeded by the bath, which correlate distant degrees of freedom without immediately generating full spectral level repulsion [2507.18286].

This mechanism is most visible in regime \([B]\). There, \(C^{zz}\) typically still decays exponentially, yet the histograms of \(\ell_{0,\lfloor L/2\rfloor}^{zz}\) develop pronounced power-law tails at small \(\ell\), involving a finite fraction of eigenstates. These tails are described as clear rare-event dominance, or Griffiths-like behavior, together with systemwide correlations. Simultaneously, the decrease of \(\xi_x/\xi_z\) with \(L\) is interpreted as an ergodic-like tendency inside a phase whose spectrum remains Poisson [2507.18286].

Regime \([C]\) sharpens the distinction between eigenstate and spectral criteria. Its eigenstates are ETH-like in the sense that \(S/S_P \to 1\) and \(D_1=1\), but the level statistics remain intermediate rather than GOE. The paper explicitly interprets this as a separation of eigenstate thermalization from spectral chaos, comparing the regime to mixed phase-space models and integrable-but-delocalized quadratic systems. In that sense, the AQS violates the usual “lockstep” expectation according to which Poisson spectra pair with area-law entanglement and GOE spectra with volume-law entanglement [2507.18286].

The resulting picture refines avalanche-based intuition. Rather than a single sharp conversion from localized to thermal behavior in all diagnostics at once, the AQS exhibits a staged destabilization: first rare-event-dominated correlations and sub-volume entanglement, then ETH-like volume-law eigenstates without full random-matrix spectral statistics, and only later complete GOE behavior. This suggests that in bath-coupled localized systems, different ETH indicators can separate over substantial parameter windows [2507.18286].

## 5. Relation to Quantum Sun models and localization theory

The AQS extends the earlier Quantum Sun construction by adding intra-chain XY hopping \(J\). In the \(U(1)\) Quantum Sun limit \(J=0\), a single critical coupling \(\alpha_c \approx 0.75\) separates localization and ergodicity, with relatively small finite-size effects. The AQS retains this scale at small \(J\) but develops a richer two-step thermalization structure once \(J\) is increased [2507.18286].

Earlier work on Quantum Sun models established that these star-graph systems can support a sharp ergodicity-breaking phase transition and an energy-dependent many-body mobility edge. In the original and \(U(1)\)-preserving Quantum Sun models, the critical coupling at \(\varepsilon=1/2\) is near \(\alpha_c \approx 0.75\), and the mobility-edge line \(\alpha_c(\varepsilon)\) rises away from the spectral center because the many-body density of states is exponentially suppressed in the tails [2308.01073]. The AQS does not replace that picture; rather, it introduces a localized-chain sector whose Anderson length \(\xi_{\mathrm{AL}}\) competes explicitly with the bath reach \(\Lambda_\alpha\), thereby exposing unconventional regimes between fully localized and fully ergodic behavior [2507.18286].

A complementary rigorous result concerns the small-\(\alpha\) limit of the Quantum Sun model. For sufficiently small \(\alpha\), one can prove Fock-space many-body localization and Poisson spectral statistics, with eigenstates overwhelmingly aligned with product \(Z\)-configurations beyond a random finite core. In that analysis the physically expected transition scale is \(\alpha_c = 1/\sqrt{2}\), although the proof applies only for \(\alpha \le \alpha_0 < 1/2\) [2506.13511]. This rigorous localized regime is consistent with the AQS phase \([A]\), but the AQS findings indicate that adding an Anderson-localized chain and finite \(J\) creates intermediate structures not captured by the simpler lockstep localization picture [2507.18286].

The comparison to XXZ chains is also diagnostic. At comparable parameters, the AQS sub-volume regime \([B]\)—Poisson statistics together with sub-volume entanglement and rare-event-dominated correlations—is reported to be absent in XXZ. The paper attributes this difference to the special role of the hub-and-spoke bath coupling and the noninteracting localization seed provided by the chain [2507.18286].

## 6. Finite-size constraints, experimental relevance, and open problems

The numerical evidence is based on relatively small total sizes, \(N+L \le 20\), with a bath as small as \(N=3\). The authors nevertheless emphasize that finite-size effects are modest enough to support extrapolations of \(\alpha_{r,\infty}^c\) and \(\alpha_{S,\infty}^c\), especially because the underlying Quantum Sun geometry is known to exhibit limited finite-size drift. Even so, the ultimate fate of regime \([B]\) in the thermodynamic limit remains open: the data show robust ergodic instabilities, but whether an eventual avalanche fully thermalizes that region is not settled [2507.18286].

Another limitation is diagnostic rather than numerical. The analysis is static, so the dynamical properties and timescales associated with the unconventional phases are not directly computed. Open questions include avalanche growth, Thouless time across the distinct regimes, and a microscopic theory for why intermediate spectral statistics in \([C]\) coexist with ETH-like eigenstates [2507.18286].

The proposed experimental setting is a platform with tunable hub-and-spoke couplings and controllable disorder, such as trapped-ion chains with ancilla spins acting as the bath. The relevant observables are entanglement scaling \(S(L)\), spectral \(r\)-statistics obtained from spectroscopy or Floquet quasienergies, and long-range connected spin correlations \(C_{ij}^{\gamma\gamma}\) together with their fat-tailed distributions. These observables map directly onto the distinctions among regimes \([A]\)–\([D]\) [2507.18286].

In summary, the interacting Anderson Quantum Sun model is significant because it exposes two distinct failures of the conventional MBL/ETH lockstep. It realizes a Poisson-spectral phase with sub-volume entanglement and rare-event-dominated correlations, and a volume-law phase with intermediate rather than GOE spectral statistics. This combination makes it a useful benchmark for studying how ergodic inclusions destabilize Anderson-localized structures, how avalanche scenarios unfold in finite systems, and how different diagnostics of thermalization can decouple in interacting disordered matter [2507.18286].

Source: https://www.emergentmind.com/topics/interacting-anderson-quantum-sun-model