Universal dynamics from a single-particle dark state
Pith reviewed 2026-05-19 21:29 UTC · model grok-4.3
The pith
A single-particle dark state at zero momentum in a dissipative spin chain leads to universal many-body scaling dynamics at long times.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
While the zero-momentum mode is dark at the single-particle level, it decays slowly as 1/log t as it becomes dressed by other modes through a dissipation-induced nonlinearity. The momentum distribution takes a universal scaling form in k sqrt(t), and the total density decays as 1/sqrt(t) log t. The same universal behavior persists for soft-core interactions.
What carries the argument
The dissipation-induced nonlinearity that dresses the single-particle dark state at zero momentum and generates the universal scaling in the many-body dynamics.
If this is right
- The momentum distribution adopts a universal scaling form when expressed in terms of k sqrt(t).
- The total particle density decays asymptotically according to 1/(sqrt(t) log t).
- The universal dynamics remain unchanged even when soft-core interactions replace hard-core ones.
- This framework accounts for the conflicting results reported in recent studies of similar dissipative systems.
Where Pith is reading between the lines
- Similar dressing mechanisms could appear in other open quantum platforms, such as arrays of Rydberg atoms or circuit QED systems, producing analogous scaling laws.
- Experiments with ultracold atoms could measure the predicted slow decay to confirm the role of the dark state.
- Extending the model to include driving fields might allow control over the effective nonlinearity and the resulting decay rates.
Load-bearing premise
The long-time behavior is controlled by the slow dressing of the zero-momentum dark state through the dissipation-induced nonlinearity rather than by other processes such as finite-size effects or initial condition details.
What would settle it
Numerical simulation or experimental measurement showing that the momentum distribution fails to collapse onto a single curve when scaled by k sqrt(t) at sufficiently long times would falsify the universal dynamics claim.
Figures
read the original abstract
Open quantum systems can host dark or subradiant states whose decay is highly suppressed. While these states have been extensively studied in the few-excitation regime, their impact on the many-body dynamics remains largely unexplored. Here, we study a spin chain subject to correlated dissipation on neighboring sites, which admits a single-particle dark state at zero momentum. We show that the single-particle dark state qualitatively alters the many-body dynamics at long times, and identify its distinct universal behavior. While the zero-momentum mode is dark at the single-particle level, it decays slowly as $1/\log t$ as it becomes dressed by other modes through a dissipation-induced nonlinearity. We demonstrate that the momentum distribution takes a universal scaling form in $k\sqrt{t}$, and the total density decays as $1/\sqrt{t}\log t$. Our results further elucidate the origin of the conflicting results in recent works. Finally, we corroborate the analytics with matrix product state simulations and show that the same universal behavior persists for soft-core interactions, underscoring the universality of the emergent dynamics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies a spin chain with correlated dissipation on neighboring sites that supports a single-particle dark state at zero momentum. It claims that this dark state is dressed at long times by other modes through a dissipation-induced nonlinearity arising from the Lindblad structure, producing a slow decay of the zero-momentum occupation as 1/log t. The momentum distribution is shown to collapse onto a universal scaling function of k sqrt(t), while the total density decays as 1/sqrt(t) log t. Analytical arguments are corroborated by matrix-product-state simulations, and the same scaling is reported to persist under soft-core interactions. The work also addresses the origin of conflicting results in earlier studies.
Significance. If the central scaling relations hold, the paper makes a useful contribution to open quantum many-body physics by demonstrating how a single-particle dark state can generate universal long-time dynamics through an effective nonlinearity. The explicit forms 1/log t, k sqrt(t), and 1/sqrt(t) log t, together with the numerical confirmation and the extension to soft-core interactions, provide concrete, testable predictions that could guide future experiments in dissipative spin systems.
major comments (2)
- The central claim that the dissipation-induced nonlinearity dominates the late-time dressing of the k=0 mode (producing the 1/log t decay and the k sqrt(t) scaling) rests on the assumption that higher-order or multi-particle processes remain subleading. The manuscript should supply a clearer scaling argument or bound showing why this effective term controls the asymptotics, rather than relying primarily on the existing MPS runs whose finite-time and finite-size reach may not yet probe the true t→∞ regime.
- The abstract states that the results elucidate conflicting findings in recent works, yet the manuscript does not appear to contain a direct, quantitative comparison (e.g., a table or figure overlay) that maps the present model onto those earlier setups. Without such a comparison, it remains unclear whether the reported universality resolves the discrepancies or merely reproduces them under different parameters.
minor comments (2)
- The scaling variable k sqrt(t) should be defined explicitly in the main text, including the precise normalization of the momentum k and any prefactors that may depend on the dissipation strength.
- A short paragraph summarizing the precise form of the neighboring-site jump operators and the single-particle dark-state wave function would help readers assess the generality of the nonlinearity without consulting the supplemental material.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and outline the revisions we will make to strengthen the presentation of our results.
read point-by-point responses
-
Referee: The central claim that the dissipation-induced nonlinearity dominates the late-time dressing of the k=0 mode (producing the 1/log t decay and the k sqrt(t) scaling) rests on the assumption that higher-order or multi-particle processes remain subleading. The manuscript should supply a clearer scaling argument or bound showing why this effective term controls the asymptotics, rather than relying primarily on the existing MPS runs whose finite-time and finite-size reach may not yet probe the true t→∞ regime.
Authors: We agree that an explicit scaling argument would make the dominance of the effective nonlinearity more transparent. In the revised manuscript we will add a dedicated paragraph deriving the scaling of higher-order multi-particle contributions relative to the leading dissipation-induced term, showing that they remain subleading for t ≫ 1 in the regime of interest. This analytical bound will be presented alongside the existing MPS data, which we maintain already capture the onset of the asymptotic regime but which we acknowledge are limited by finite simulation times. revision: yes
-
Referee: The abstract states that the results elucidate conflicting findings in recent works, yet the manuscript does not appear to contain a direct, quantitative comparison (e.g., a table or figure overlay) that maps the present model onto those earlier setups. Without such a comparison, it remains unclear whether the reported universality resolves the discrepancies or merely reproduces them under different parameters.
Authors: We thank the referee for highlighting this point. The manuscript discusses the origin of the conflicting results in the introduction and in the concluding section by contrasting the dissipation structure and initial conditions with those of prior studies. However, we agree that a direct, quantitative mapping would improve clarity. In the revision we will add a new figure that overlays the parameter regimes of the earlier works with our model and explicitly shows how the presence or absence of the single-particle dark state accounts for the differing long-time scalings. revision: yes
Circularity Check
No significant circularity; derivation is self-contained
full rationale
The paper starts from the Lindblad master equation for a spin chain with neighboring-site correlated dissipation, identifies the single-particle zero-momentum dark state by direct inspection of the jump operators, then constructs an effective nonlinear dressing term whose consequences (1/log t decay, k sqrt(t) scaling of the momentum distribution, and 1/sqrt(t) log t density decay) are derived analytically and cross-checked against independent matrix-product-state simulations. No equation reduces to a fitted parameter renamed as a prediction, no load-bearing premise rests solely on a self-citation, and no ansatz is imported without independent justification; the claimed universality follows from the structure of the dissipation operators and the long-time asymptotics of the resulting effective dynamics.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Correlated dissipation on neighboring sites admits a single-particle dark state at zero momentum.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
While the zero-momentum mode is dark at the single-particle level, it decays slowly as 1/log t as it becomes dressed by other modes through a dissipation-induced nonlinearity. ... momentum distribution takes a universal scaling form in k √t, and the total density decays as 1/√t log t
-
IndisputableMonolith/Foundation/BranchSelection.leanbranch_selection unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
1/Γ dρ_k/dt = −F_k({ρ_p}) with F_k containing the nonlocal integral terms |δ_p ρ_p / (δ_p − k)|^2 and 2m ∫ (ρ_k − ρ_p)/|δ_p − k|^2
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
- [2]
-
[3]
F. Verstraete, M. M. Wolf, and J. Ignacio Cirac, Na- ture Physics5, 633 (2009)
work page 2009
-
[4]
P. M. Harrington, E. J. Mueller, and K. W. Murch, Nat Rev Phys4, 660 (2022)
work page 2022
-
[5]
L. M. Sieberer, M. Buchhold, J. Marino, and S. Diehl, Rev. Mod. Phys.97, 025004 (2025)
work page 2025
- [6]
-
[7]
A. A. Svidzinsky, J.-T. Chang, and M. O. Scully, Phys. Rev. A81, 053821 (2010)
work page 2010
-
[8]
A. S. Sheremet, M. I. Petrov, I. V. Iorsh, A. V. Poshakinskiy, and A. N. Poddubny, Rev. Mod. Phys. 95, 015002 (2023)
work page 2023
-
[9]
A. Asenjo-Garcia, M. Moreno-Cardoner, A. Albrecht, H. J. Kimble, and D. E. Chang, Phys. Rev. X7, 031024 (2017)
work page 2017
- [10]
-
[11]
R. T. Sutherland and F. Robicheaux, Phys. Rev. A 94, 013847 (2016)
work page 2016
-
[12]
R. J. Bettles, S. A. Gardiner, and C. S. Adams, Phys. Rev. A94, 043844 (2016)
work page 2016
-
[13]
H. H. Jen, M.-S. Chang, and Y.-C. Chen, Phys. Rev. A94, 013803 (2016)
work page 2016
-
[14]
D. F. Kornovan, A. S. Sheremet, and M. I. Petrov, Phys. Rev. B94, 245416 (2016)
work page 2016
-
[15]
H. R. Haakh, S. Faez, and V. Sandoghdar, Phys. Rev. A94, 053840 (2016)
work page 2016
-
[16]
J. Ruostekoski and J. Javanainen, Phys. Rev. Lett. 117, 143602 (2016)
work page 2016
- [17]
-
[18]
T. S. do Espirito Santo, P. Weiss, A. Cipris, R. Kaiser, W. Guerin, R. Bachelard, and J. Schachen- mayer, Phys. Rev. A101, 013617 (2020)
work page 2020
-
[19]
A. Albrecht, L. Henriet, A. Asenjo-Garcia, P. B. Dieterle, O. Painter, and D. E. Chang, New J. Phys. 21, 025003 (2019)
work page 2019
-
[20]
S. J. Masson, I. Ferrier-Barbut, L. A. Orozco, A. Browaeys, and A. Asenjo-Garcia, Phys. Rev. Lett. 125, 263601 (2020)
work page 2020
- [21]
- [22]
- [23]
-
[24]
L. Henriet, J. S. Douglas, D. E. Chang, and A. Al- brecht, Phys. Rev. A99, 023802 (2019)
work page 2019
-
[25]
S. E. Begg and R. Hanai, Phys. Rev. Lett.132, 120401 (2024)
work page 2024
- [26]
-
[27]
A. March´ e, H. Yoshida, A. Nardin, H. Katsura, and L. Mazza, (2026), arXiv:2601.08606
-
[28]
T. Peyronel, O. Firstenberg, Q.-Y. Liang, S. Hof- ferberth, A. V. Gorshkov, T. Pohl, M. D. Lukin, and V. Vuleti´ c, Nature488, 57 (2012)
work page 2012
-
[29]
C. L. Baldwin, P. Bienias, A. V. Gorshkov, M. J. Gullans, and M. Maghrebi, Phys. Rev. A104, L041303 (2021)
work page 2021
- [30]
- [31]
-
[32]
See Supplemental Material for (I) the mapping of the dissipative XX spin chain to fermions; (II) the perturbative approach at weak dissipation based on time-dependent generalized Gibbs ensemble, includ- ing the derivation of the kinetic equation, its analyti- cal continuation, scaling solutions in different regimes, and the bosonic momentum distribution; ...
- [33]
- [34]
- [35]
- [36]
-
[37]
M. Ali, N. A. Kamar, A. Seif, and M. Maghrebi, Phys. Rev. Lett.136, 060404 (2026)
work page 2026
-
[38]
I. Bouchoule, L. Dubois, and L.-P. Barbier, Phys. Rev. A104, L031304 (2021)
work page 2021
- [39]
-
[40]
M. A. Cazalilla, R. Citro, T. Giamarchi, E. Orignac, and M. Rigol, Rev. Mod. Phys.83, 1405 (2011)
work page 2011
-
[41]
S. Paeckel, T. K¨ ohler, A. Swoboda, S. R. Manmana, U. Schollw¨ ock, and C. Hubig, Annals of Physics411, 167998 (2019)
work page 2019
-
[42]
A. J. Daley, Advances in Physics63, 77 (2014)
work page 2014
-
[44]
P. Calabrese, F. H. L. Essler, and M. Fagotti, Jour- nal of Statistical Mechanics: Theory and Experiment 2012, P07016 (2012)
work page 2012
-
[45]
Sachdev,Quantum Phase Transitions, 2nd ed
S. Sachdev,Quantum Phase Transitions, 2nd ed. (Cambridge University Press, Cambridge, 2011)
work page 2011
- [46]
-
[47]
N. Iorgov, V. Shadura, and Y. Tykhyy, Journal of Statistical Mechanics: Theory and Experiment2011, P02028 (2011) 6 Supplemental Material: Universal dynamics from a single-particle dark state In this Supplemental Material, we provide additional details on the results stated in the main text. In Section S.I, we provide the mapping to fermions for both the H...
work page 2011
-
[48]
Basis transformation 10
-
[49]
Analytic continuation 13
-
[50]
Dynamics from the GGE equation 14
Dissipative Tonks-Girardeau gas 14 C. Dynamics from the GGE equation 14
-
[51]
Analytical solution at long times 15
-
[52]
Full scaling solution 15
-
[53]
Evolution at short times 17
-
[54]
Universality at long times 19
-
[55]
Bosonic momentum distribution 19 S.III. MPS simulations 20 A. Technical implementation details 20 B. Further numerical results 21
-
[56]
Total density and trajectory-averaged entanglement entropy 22
-
[57]
Fermionic momentum space observables 22
-
[58]
Hard-core boson momentum observables 23
-
[59]
MAPPING TO FERMIONS In this section, we provide the mapping to fermions, and specifically derive Eq
Comparison with the GGE at short times 23 Supplemental References 25 S.I. MAPPING TO FERMIONS In this section, we provide the mapping to fermions, and specifically derive Eq. (4) of the main text. For convenience, we first apply a transformationσy,z i → −σ y,z i , aπrotation aroundx. Under this transformation, the Hamiltonian remains unchanged, while the ...
-
[60]
Basis transformation In the following analysis, we also need a change of basis states between the even and odd sectors. Since both ˆck fork∈Z ap and ˆcp forp∈Z p span the one-particle Hilbert space on a ring, they can be related to each other. Indeed, the corresponding momentum eigenstates are related via |k⟩= X p∈Zp ⟨p|k⟩|p⟩,(S.22) where ⟨p|k⟩= 1 L LX j=...
-
[61]
Wick’s theorem Here, we use Eq. (S.19) but substitute ˆL1 with ˆLL for convenience to find 1 Γ d dt ρk =⟨ ˆ˜c† ˆQkˆ˜c⟩ − ⟨ˆQkˆ˜c†ˆ˜c⟩.(S.27) where the operator ˆ˜c, defined in Eq. (S.11), is linear in fermionic operators. The second term on the rhs of the above equation does not involve a change of parity ( ˆ˜c†ˆ˜cconserves the particle number), and can b...
-
[62]
Analytic continuation As defined in Eq. (9) of the main text, one can make an analytic continuation to the unit disk via QA(z) = Z 2π 0 dk 2π Ak eik +z eik −z ⇒Q A(z) is analytic in|z|<1,and ReQ A(eik−0+ ) =A k .(S.40) A Kramers-Kronig relation follows as well: ImQ A(eik−0+ ) = Z ⧸ dλ 2π Aλ cot( k−λ 2 ),(S.41) thus reducing to the circular Hilbert transfo...
-
[63]
Dissipative Tonks-Girardeau gas Here, we briefly consider the dissipative Tonks-Girardeau gas. We start from the Lieb-Liniger Hamiltonian ˆH= Z x ˆΨ†(−∂2 x/2 +g ˆΨ† ˆΨ) ˆΨ,(S.50) where ˆΨ, ˆΨ† denote the bosonic field operators satisfying the commutation relation [ ˆΨ(x), ˆΨ†(y)] =δ(x−y); we have setℏ=m= 1. The hard-core limitg→ ∞recovers the Tonks-Girard...
-
[64]
Analytical solution at long times In this and the next subsection, we set Γ = 1 for convenience. We start with Eq. (6) of the main text at long wavelengths (1/|eip −1| 2 ≈p 2): d dt ρ0(t) =− 1 π2 Z dλρλ(t) 2 − 2m π Z ρ0(t)−ρ λ(t) λ2 ,(S.53) together with the scaling solutionρ k(t) =ρ 0(t)f(k √ t) [(Eq. (7) of the main text)] to obtain an analytical soluti...
-
[65]
(10) of the main text and expand it aroundz= 1
Full scaling solution To identify the full scaling solution, it is more convenient to work with the analytically continued function Q(t, z) in Eq. (10) of the main text and expand it aroundz= 1. Letz= 1 +wand, in an abuse of notation, Q(t, z)→Q(t, w) in the new coordinates. To the first nontrivial order derived in the previous section, we haveρ k(t)∝ 1 lo...
-
[66]
We thus require the coefficient of the singular term in Eq. (S.67) to vanish, which then yields A1 → √π 4 ⇒f 1(x) = π 2 e−x2 .(S.68) We thus obtain the universal coefficient as expected. 16 FIG. S.1. The functiong(x) defined in Eq. (S.69). We can then determine the corresponding momentum distribution at this order: f2(x) = ReF 2(0− +ix) =C 2e−x2 +g(x),(S....
-
[67]
We first consider a fully excited initial state
Evolution at short times In this subsection, we examine the dynamics of fermionic densities at short times, and argue that, for highly excited initial states, thek= 0 mode decays quickly. We first consider a fully excited initial state. Interestingly, thek= 0 (k=π) that is slow (fast) to decay at long times, is the fastest (slowest) to decay in this regim...
-
[68]
Universality at long times Here, we show that the asymptotic behavior is independent of the initial state. In Fig. S.4, we find that the inverseρ 0(t) exhibits the same logarithmic scaling down to the same coefficient, i.e., with the same slope in the log-linear plot. FIG. S.4. Universality of the late-time decay of the zero mode. Inverse zero-mode densit...
-
[69]
Bosonic momentum distribution The fermionic correlation function is given by ⟨ˆc† jc0⟩= Z π −π dk 2π ρk cos(kj) = ( a0, j= 0 1 2 a|j|, j̸= 0 (S.77) where we have usedρ k =ρ −k. The spin correlation function is given byC(j)≡ ⟨σ + j+lσ− l ⟩which can be written in terms of fermions as C(j)≡ * ˆc† j j−1Y l=1 (−1)ˆc† l ˆclˆc0 + .(S.78) 19 (a) (b) (c) FIG. S.5....
-
[70]
Total density and trajectory-averaged entanglement entropy Figure S.6(a) shows the unscaled total density for several system sizes, initial filling fractions and loss rates. The long-time dynamics appear independent of the specific initial conditions and loss rate, and the curves show the same trend at late times. For a bipartition into subsystemsAandB, t...
-
[71]
(7,12), on a finite open spin chain of lengthL= 110
Fermionic momentum space observables Figure S.7 tests the analytical predictions for the fermionic observables (after JW), Eqs. (7,12), on a finite open spin chain of lengthL= 110. The rescaled total density (n √ Γt)−1 is shown in Fig. S.7(a), where the straight line gives the logarithmic decay predicted by Eq. (12) of the main text vian= R dk 2π ρk. The ...
-
[72]
S.7(b,c) for the hard-core boson (spin) observables, i.e
Hard-core boson momentum observables Figure S.8 shows the analog of Fig. S.7(b,c) for the hard-core boson (spin) observables, i.e. without the Jordan-Wigner string. The trends mirror the fermionic case with similar fitted coefficients. The bosonic MPS data may also be compared with the GGE results of Fig. S.5; we again find agreement within at most a factor of 2
-
[73]
Comparison with the GGE at short times Figure S.9 compares the momentum distributionρ k from MPS against the GGE solution with the initial fillingn(0) = 0.8. Panel (a) shows that both methods reproduce theρ 0 ∼1/Γt→1/log(Γt) crossover around Γτ∼1, and panel (b) shows good quantitative agreement of the fullρ k at short times. This agreement is remarkable g...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.