REVIEW
BMN Spread Complexity Across Phase Transitions
High Energy Physics - Theory · Submitted 30 Aug 2026 · 72 references
Pith's verdict: The intermediate-phase result is a parametrization, not a prediction, because the spectral density is postulated with free parameters. The gapped and continuum phases are handled well.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
BMN spread complexity reveals three-phase scrambling crossover from bound states to maximal chaos.
The reading
What carries the argument
The key machinery is the charged Thermofield Double (cTFD) state of the BMN matrix model at leading order in the large mass deformation µ ≫ 1, combined with the Lanczos algorithm and the continuum limit n → x of the Krylov chain. The cTFD state incorporates both temperature β⁻¹ and chemical potential ν, and its return amplitude near the KMS analyticity boundary yields factorial growth of spectral moments, which translates to linear growth of Lanczos coefficients at high order. The continuum limit reduces the discrete Schrödinger equation to an advection partial differential equation $\partial_t \psi = -\partial_x (b(x) \psi)$ with velocity field $b(x) = \frac{\pi}{\beta} x + \gamma \sqrt{x}$. Solving this equation gives the wavepacket trajectory and the spread complexity across all three phases.
Extended reading notes
Core claim
The paper establishes that spread complexity in the BMN matrix model exhibits a three-phase structure governed by temperature relative to the mass gap µ. In the low-temperature gapped phase (βµ ≫ 1), the Lanczos coefficients $b_n$ are bounded (∼ µ), leading to quadratic growth $C(t) \sim \mu^2 t^2$. In the high-temperature continuum phase (βµ ≪ 1), $b_n \sim \frac{\pi}{\beta} n$, driving exponential growth $C(t) \sim \exp(\pi t/\beta)$ with Krylov Lyapunov exponent $\lambda_K = \pi/\beta$ that saturates the Maldacena-Shenker-Stanford chaos bound. In the intermediate regime (βµ ∼ 1, βν ∼ 1), the Lanczos coefficients acquire a sub-leading $\sqrt{n}$ correction: $b_n \sim \frac{\pi}{\beta} n + \gamma \sqrt{n}$, producing a sub-exponential crossover before asymptotic chaotic behavior sets in. The parameter $\gamma(\beta\mu,\beta\nu)$ encodes the competition between discrete mass-gap modes and the thermal continuum modulated by chemical potential ν.
Load-bearing premise
The derivation of the Lanczos coefficients in the intermediate regime ($b_n \sim \frac{\pi}{\beta} n + \gamma \sqrt{n}$) assumes a spectral density that artificially splits into a thermal continuum with a specific functional form (Eq. C.8) and discrete mass-gap poles with ad hoc amplitudes, without deriving these amplitudes from the BMN matrix model beyond leading order.
Editorial extensions
If this is right
- Spread complexity in the continuum phase of the BMN model serves as a direct microscopic probe of maximal chaos, with the Lyapunov exponent matching the universal bound $\lambda_L \leq 2\pi/\beta$ after accounting for the difference between operator and state complexity.
- The intermediate-temperature sub-exponential correction controlled by $\gamma(\beta\mu,\beta\nu)$ provides a controlled setting to study the onset of scrambling near a confinement-deconfinement transition in a holographic context.
- The cTFD construction offers a template for computing Krylov complexity in other massive matrix models with global $U(1)$ symmetries at leading order in the deformation.
- The analytical form of the Lanczos coefficients and the advection equation allow precise predictions for the time-dependence of spread complexity in each phase, which could be compared with lattice or Monte Carlo simulations of the BMN matrix model.
- The results extend the universal operator growth hypothesis to state complexity in a non-conformal theory with a mass gap, showing that the KMS analyticity strip width $\beta/2$ controls the asymptotic growth rate even when discrete states coexist with the continuum.
Reading between the lines
- The sub-leading $\gamma \sqrt{n}$ correction in the Lanczos coefficients might be testable in direct numerical simulations of the BMN matrix model at intermediate temperatures, where the spectral density computed via Hamiltonian truncation or lattice methods could be compared to the model in Eq. (C.8).
- The parameter $\gamma(\beta\mu,\beta\nu)$ could be reinterpreted as a renormalisation of the effective temperature due to the mass gap, potentially modifying the chaos bound at finite $\mu$; this would provide a concrete extension of the MSS bound to gapped systems.
- The three-phase spread complexity likely maps to distinct holographic duals: the LLM geometry (low $T$), a charged non-extremal D0-brane black hole (intermediate $T$), and an uncharged non-extremal D0-brane black hole (high $T$). The growth law of the wormhole volume or action in each geometry should match the corresponding $C(t)$ behaviour derived here.
- The transition from quadratic to sub-exponential to exponential growth may serve as a diagnostic for the confinement-deconfinement transition in the BMN matrix model, analogous to the thermal entropy jump seen in lattice studies.
Original abstract
We investigate the dynamical behavior and spread complexity of quantum states within the mass-deformed BMN matrix model equipped with an emergent $U(1)$ global charge at finite temperature $\beta^{-1}$ and chemical potential $\nu$. Focusing on the strong-coupling regime ($\mu \gg 1$), we model the dynamics via a charged Thermofield Double (cTFD) state and trace the evolution of spread complexity across three distinct thermodynamic regimes. In the low-temperature gapped phase ($\beta \mu \gg 1$), discrete mass-gap bound states dominate, trapping wavepacket dispersion and producing non-chaotic, oscillatory early-time growth. Conversely, in the high-temperature continuum phase ($\beta \mu \ll 1$), thermal excitations overwhelm the mass gap, driving a transition to a continuous advection field that exhibits maximal chaotic scrambling with a Krylov Lyapunov exponent $\lambda_K = \pi / \beta$ that saturates the universal bound. In the intermediate temperature regime ($\beta \mu \sim 1, \beta \nu \sim 1$), the interplay between mass-gap bound states and the continuous thermal background induces a sub-leading correction to the Lanczos coefficients $b_n \sim \frac{\pi}{\beta} n + \gamma \sqrt{n}$, governing a continuous sub-exponential crossover before full chaotic thermalization. Technical derivations regarding KMS boundary conditions, grand canonical spectral moments, residue analysis, and time-reversal symmetry breaking in Krylov space are detailed in four dedicated appendices.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We investigate the dynamical behavior and spread complexity of quantum states within the mass-deformed BMN matrix model equipped with an emergent $U(1)$ global charge at finite temperature $\beta^{-1}$ and chemical potential $\nu$. Focusing on the strong-coupling regime ($\mu \gg 1$), we model the dynamics via a charged Thermofield Double (cTFD) state and trace the evolution of spread complexity across three distinct thermodynamic regimes. In the low-temperature gapped phase ($\beta \mu \gg 1$), discrete mass-gap bound states dominate, trapping wavepacket dispersion and producing non-chaotic, oscillatory early-time growth. Conversely, in the high-temperature continuum phase ($\beta \mu \ll 1$), thermal excitations overwhelm the mass gap, driving a transition to a continuous advection field that exhibits maximal chaotic scrambling with a Krylov Lyapunov exponent $\lambda_K = \pi / \beta$ that saturates the universal bound. In the intermediate temperature regime ($\beta \mu \sim 1, \beta \nu \sim 1$), the interplay between mass-gap bound states and the continuous thermal background induces a sub-leading correction to the Lanczos coefficients $b_n \sim \frac{\pi}{\beta} n + \gamma \sqrt{n}$, governing a continuous sub-exponential crossover before full chaotic thermalization. Technical derivations regarding KMS boundary conditions, grand canonical spectral moments, residue analysis, and time-reversal symmetry breaking in Krylov space are detailed in four dedicated appendices.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We investigate the dynamical behavior and spread complexity of quantum states within the mass-deformed BMN matrix model equipped with an emergent $U(1)$ global charge at finite temperature $\beta^{-1}$ and chemical potential $\nu$. | contribution 0.97 | 0.86 | 0.04 | asserted 0.83 |
| s1 | Focusing on the strong-coupling regime ($\mu \gg 1$), we model the dynamics via a charged Thermofield Double (cTFD) state and trace the evolution of spread complexity across three distinct thermodynamic regimes. | method 0.81 | 0.93 | 0.07 | asserted 0.97 |
| s2 | In the low-temperature gapped phase ($\beta \mu \gg 1$), discrete mass-gap bound states dominate, trapping wavepacket dispersion and producing non-chaotic, oscillatory early-time growth. | result 1.00 | 0.84 | 0.10 | asserted 1.00 |
| s3 | Conversely, in the high-temperature continuum phase ($\beta \mu \ll 1$), thermal excitations overwhelm the mass gap, driving a transition to a continuous advection field that exhibits maximal chaotic scrambling with a Krylov Lyapunov exponent $\lambda_K = \pi / \beta$ that saturates the universal bound. | result 1.00 | 0.85 | 0.32 | asserted 1.00 |
| s4 | In the intermediate temperature regime ($\beta \mu \sim 1, \beta \nu \sim 1$), the interplay between mass-gap bound states and the continuous thermal background induces a sub-leading correction to the Lanczos coefficients $b_n \sim \frac{\pi}{\beta} n + \gamma \sqrt{n}$, governing a continuous sub-exponential crossover before full chaotic thermalization. | result 1.00 | 0.88 | 0.18 | asserted 1.00 |
| s5 | Technical derivations regarding KMS boundary conditions, grand canonical spectral moments, residue analysis, and time-reversal symmetry breaking in Krylov space are detailed in four dedicated appendices. | method 0.88 | 0.71 | 0.02 | asserted 0.99 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How the literature uses it
What this paper does with the arXiv papers it cites
What it cites · 56 sentences without a confident label
Builds on 5 papers
-
The purpose of this article is to extend the notion of spread complexity of states [ref] - [ref] in the context of BMN matrix models [ref] - [ref] at finite chemical potential $ν$ and temperature $β^-1$ [ref] .
Thermodynamics of the BMN matrix model at strong coupling · Introduction and General Idea · 0.86
-
The purpose of this article is to extend the notion of spread complexity of states [ref] - [ref] in the context of BMN matrix models [ref] - [ref] at finite chemical potential $ν$ and temperature $β^-1$ [ref] .
A Universal Operator Growth Hypothesis · Introduction and General Idea · 0.86
-
The purpose of this article is to extend the notion of spread complexity of states [ref] - [ref] in the context of BMN matrix models [ref] - [ref] at finite chemical potential $ν$ and temperature $β^-1$ [ref] .
Quantum chaos and the complexity of spread of states · Introduction and General Idea · 0.86
Uses the method or tool of 5 papers
-
To obtain the high-order spectral moments $ℳ_2n$, we begin with the spectral representation of the thermal return amplitude $ℛ(t) = ⟨TFD(0) e^-i ℋ t | TFD(0) ⟩$, which can be expressed as an integral of the following form [ref] [equation] where $f(ω)$ is the symmetric thermal spectral density.
A Universal Operator Growth Hypothesis · Moments and Lanczos coefficients · 0.93
-
We outline how to obtain the wave function (or the probability amplitude) $ψ_n(t)$ associated with the $n$th Krylov site $|K_n⟩$ [ref] .
Quantum chaos and the complexity of spread of states · Fuzzy sphere model and spread complexity: A review · 0.86
-
Given the above state of the art, we choose the Krylov basis at $t=0$ as [ref] [equation] where, $Ψ_0(x,y)$ is the Gaussian wavefunction of the ground state of the composite system and $|x,y⟩$ is the eigen basis in the position representation.
Krylov complexity for Lin-Maldacena geometries and their holographic duals · Fuzzy sphere model and spread complexity: A review · 0.97
Compares with 1 paper
-
Finally, substituting $x = 2μt$ leads directly to the early-time expansion of the spread complexity of states in the gapped phase of the matrix model [equation] which mimics previous results [ref] obtained at zero temperature [ref].
Krylov state complexity for BMN matrix model · Spread complexity in the gapped phase · 0.92
Cites as support 4 papers
-
In the continuum limit ($n →x$), this manifests as a continuous advection PDE, driving maximal quantum chaos and exponential scrambling $𝒞(t) ∼exp(πt / β)$ with Krylov Lyapunov exponent $λ_K = π/ β$ which saturates the universal bound [ref] . * Intermediate phase ($βμ∼1, βν∼1$): At intermediate scales, the interplay between discrete bound states and the continuous thermal background induces sub-exponential corrections governed by the parameter $γ$ that interpolates between discrete mass gap and the thermal continuum.
A bound on chaos · Summary and conclusions · 0.96
-
Instead, the Thermofield Double (TFD) state $|TFD⟩= 1/√(Z) ∑_n e^-βE_n/2 |n⟩_L ⊗|n⟩_R$ holographically maps to two decoupled, smooth LLM geometries [ref] - [ref] generated by background supergravity fluxes without black hole horizons.
The Supergravity Dual of the BMN Matrix Model · The complex oscillator model · 0.98
-
Instead, the Thermofield Double (TFD) state $|TFD⟩= 1/√(Z) ∑_n e^-βE_n/2 |n⟩_L ⊗|n⟩_R$ holographically maps to two decoupled, smooth LLM geometries [ref] - [ref] generated by background supergravity fluxes without black hole horizons.
Bubbling AdS space and 1/2 BPS geometries · The complex oscillator model · 0.98
Cites as background 12 papers
-
The BMN or Plane Wave Matrix Model [ref] is defined as[Throughout this paper, we choose to work with natural units where we set $ħ=c=l=1$.] [ref] [equation]
Chaos in the BMN matrix model · Fuzzy sphere model and spread complexity: A review · 0.98
-
The BMN or Plane Wave Matrix Model [ref] is defined as[Throughout this paper, we choose to work with natural units where we set $ħ=c=l=1$.] [ref] [equation]
Strings in flat space and pp waves from ${\cal N}=4$ Super Yang Mills · Fuzzy sphere model and spread complexity: A review · 0.98
-
Krylov and/or the spread complexity of states and operators [ref] - [ref] and its various applications including chaos [ref] - [ref] , holography [ref] - [ref] and black holes [ref] have garnered renewed attention in the recent years due to their several remarkable features[For a nice and comprehensive set of reviews on the subject see [ref] - [ref] .].
Random Matrix Theory for Complexity Growth and Black Hole Interiors · Introduction and General Idea · 1.00
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.