pith. sign in

arxiv: 2605.16494 · v1 · pith:QZCX5E6Nnew · submitted 2026-05-15 · ❄️ cond-mat.quant-gas · cond-mat.stat-mech· quant-ph

Universal dynamics from a single-particle dark state

Pith reviewed 2026-05-19 21:29 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.stat-mechquant-ph
keywords dark statesopen quantum systemsdissipative spin chainsuniversal scalingmany-body dynamicssubradiant statescorrelated dissipationlong-time decay
0
0 comments X

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.

A spin chain with correlated dissipation on neighboring sites has a dark state for a single particle at zero momentum. This state remains protected at the few-body level but gets dressed by other modes at the many-body level through an effective nonlinearity created by the dissipation. The dressing causes the zero-momentum mode to decay slowly, resulting in a momentum distribution that scales universally with the variable k times the square root of time. The total density then falls off as one over the square root of t multiplied by the log of t. This picture matters because it reveals how seemingly protected states can still shape the collective long-time behavior in open quantum systems and clarifies earlier conflicting observations.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.16494 by Arghavan Safavi-Naini, Johannes Schachenmayer, Mohammad Maghrebi, Ruben Daraban.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Schematic representation of a spin chain (al [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Long-time dynamics of fermionic densities start [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Universal dynamics for the dissipative Bose [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. 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.
  2. 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)
  1. 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.
  2. 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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence of the zero-momentum dark state under the chosen dissipation and on the emergence of a dissipation-induced nonlinearity that dresses this mode at long times.

axioms (1)
  • domain assumption Correlated dissipation on neighboring sites admits a single-particle dark state at zero momentum.
    This is the foundational premise that allows the many-body analysis to proceed from the single-particle level.

pith-pipeline@v0.9.0 · 5726 in / 1284 out tokens · 48290 ms · 2026-05-19T21:29:38.351929+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation 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.lean branch_selection unclear
    ?
    unclear

    Relation 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

72 extracted references · 72 canonical work pages

  1. [1]

    Preskill, Quantum2, 79 (2018)

    J. Preskill, Quantum2, 79 (2018)

  2. [2]

    Diehl, A

    S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. P. B¨ uchler, and P. Zoller, Nature Physics4, 878 (2008)

  3. [3]

    Verstraete, M

    F. Verstraete, M. M. Wolf, and J. Ignacio Cirac, Na- ture Physics5, 633 (2009)

  4. [4]

    P. M. Harrington, E. J. Mueller, and K. W. Murch, Nat Rev Phys4, 660 (2022)

  5. [5]

    L. M. Sieberer, M. Buchhold, J. Marino, and S. Diehl, Rev. Mod. Phys.97, 025004 (2025)

  6. [6]

    Prasad and R

    S. Prasad and R. J. Glauber, Phys. Rev. A61, 063814 5 (2000)

  7. [7]

    A. A. Svidzinsky, J.-T. Chang, and M. O. Scully, Phys. Rev. A81, 053821 (2010)

  8. [8]

    A. S. Sheremet, M. I. Petrov, I. V. Iorsh, A. V. Poshakinskiy, and A. N. Poddubny, Rev. Mod. Phys. 95, 015002 (2023)

  9. [9]

    Asenjo-Garcia, M

    A. Asenjo-Garcia, M. Moreno-Cardoner, A. Albrecht, H. J. Kimble, and D. E. Chang, Phys. Rev. X7, 031024 (2017)

  10. [10]

    Guerin, M

    W. Guerin, M. O. Ara´ ujo, and R. Kaiser, Phys. Rev. Lett.116, 083601 (2016)

  11. [11]

    R. T. Sutherland and F. Robicheaux, Phys. Rev. A 94, 013847 (2016)

  12. [12]

    R. J. Bettles, S. A. Gardiner, and C. S. Adams, Phys. Rev. A94, 043844 (2016)

  13. [13]

    H. H. Jen, M.-S. Chang, and Y.-C. Chen, Phys. Rev. A94, 013803 (2016)

  14. [14]

    D. F. Kornovan, A. S. Sheremet, and M. I. Petrov, Phys. Rev. B94, 245416 (2016)

  15. [15]

    H. R. Haakh, S. Faez, and V. Sandoghdar, Phys. Rev. A94, 053840 (2016)

  16. [16]

    Ruostekoski and J

    J. Ruostekoski and J. Javanainen, Phys. Rev. Lett. 117, 143602 (2016)

  17. [17]

    Ruostekoski and J

    J. Ruostekoski and J. Javanainen, Phys. Rev. A96, 033857 (2017)

  18. [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)

  19. [19]

    Albrecht, L

    A. Albrecht, L. Henriet, A. Asenjo-Garcia, P. B. Dieterle, O. Painter, and D. E. Chang, New J. Phys. 21, 025003 (2019)

  20. [20]

    S. J. Masson, I. Ferrier-Barbut, L. A. Orozco, A. Browaeys, and A. Asenjo-Garcia, Phys. Rev. Lett. 125, 263601 (2020)

  21. [21]

    Zhang and K

    Y.-X. Zhang and K. Mølmer, Phys. Rev. Lett.122, 203605 (2019)

  22. [22]

    Zhang and K

    Y.-X. Zhang and K. Mølmer, Phys. Rev. Lett.125, 253601 (2020)

  23. [23]

    Zhang and K

    Y.-X. Zhang and K. Mølmer, Phys. Rev. Lett.128, 093602 (2022)

  24. [24]

    Henriet, J

    L. Henriet, J. S. Douglas, D. E. Chang, and A. Al- brecht, Phys. Rev. A99, 023802 (2019)

  25. [25]

    S. E. Begg and R. Hanai, Phys. Rev. Lett.132, 120401 (2024)

  26. [26]

    Pocklington and A

    A. Pocklington and A. A. Clerk, PRX Quantum6, 030349 (2025)

  27. [27]

    March´ e, H

    A. March´ e, H. Yoshida, A. Nardin, H. Katsura, and L. Mazza, (2026), arXiv:2601.08606

  28. [28]

    Peyronel, O

    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)

  29. [29]

    C. L. Baldwin, P. Bienias, A. V. Gorshkov, M. J. Gullans, and M. Maghrebi, Phys. Rev. A104, L041303 (2021)

  30. [30]

    Marino and S

    J. Marino and S. Diehl, Phys. Rev. Lett.116, 070407 (2016)

  31. [31]

    N. A. Kamar, M. Ali, and M. Maghrebi, (2025), arXiv:2505.22721

  32. [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. [33]

    Lange, Z

    F. Lange, Z. Lenarˇ ciˇ c, and A. Rosch, Nature Com- munications8, 15767 (2017)

  34. [34]

    Lange, Z

    F. Lange, Z. Lenarˇ ciˇ c, and A. Rosch, Phys. Rev. B 97, 165138 (2018)

  35. [35]

    Bouchoule, B

    I. Bouchoule, B. Doyon, and J. Dubail, SciPost Phys.9, 044 (2020)

  36. [36]

    Riggio, L

    F. Riggio, L. Rosso, D. Karevski, and J. Dubail, Phys. Rev. A109, 023311 (2024)

  37. [37]

    M. Ali, N. A. Kamar, A. Seif, and M. Maghrebi, Phys. Rev. Lett.136, 060404 (2026)

  38. [38]

    Bouchoule, L

    I. Bouchoule, L. Dubois, and L.-P. Barbier, Phys. Rev. A104, L031304 (2021)

  39. [39]

    Bouchoule and J

    I. Bouchoule and J. Dubail, Phys. Rev. Lett.126, 160603 (2021)

  40. [40]

    M. A. Cazalilla, R. Citro, T. Giamarchi, E. Orignac, and M. Rigol, Rev. Mod. Phys.83, 1405 (2011)

  41. [41]

    Paeckel, T

    S. Paeckel, T. K¨ ohler, A. Swoboda, S. R. Manmana, U. Schollw¨ ock, and C. Hubig, Annals of Physics411, 167998 (2019)

  42. [42]

    A. J. Daley, Advances in Physics63, 77 (2014)

  43. [44]

    Calabrese, F

    P. Calabrese, F. H. L. Essler, and M. Fagotti, Jour- nal of Statistical Mechanics: Theory and Experiment 2012, P07016 (2012)

  44. [45]

    Sachdev,Quantum Phase Transitions, 2nd ed

    S. Sachdev,Quantum Phase Transitions, 2nd ed. (Cambridge University Press, Cambridge, 2011)

  45. [46]

    Shirai and T

    T. Shirai and T. Mori, Phys. Rev. E101, 042116 (2020)

  46. [47]

    Iorgov, V

    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...

  47. [48]

    Basis transformation 10

  48. [49]

    Analytic continuation 13

  49. [50]

    Dynamics from the GGE equation 14

    Dissipative Tonks-Girardeau gas 14 C. Dynamics from the GGE equation 14

  50. [51]

    Analytical solution at long times 15

  51. [52]

    Full scaling solution 15

  52. [53]

    Evolution at short times 17

  53. [54]

    Universality at long times 19

  54. [55]

    MPS simulations 20 A

    Bosonic momentum distribution 19 S.III. MPS simulations 20 A. Technical implementation details 20 B. Further numerical results 21

  55. [56]

    Total density and trajectory-averaged entanglement entropy 22

  56. [57]

    Fermionic momentum space observables 22

  57. [58]

    Hard-core boson momentum observables 23

  58. [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 ...

  59. [60]

    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

    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=...

  60. [61]

    (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

    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...

  61. [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...

  62. [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...

  63. [64]

    We start with Eq

    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...

  64. [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...

  65. [66]

    (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

    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....

  66. [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...

  67. [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...

  68. [69]

    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

    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....

  69. [70]

    The long-time dynamics appear independent of the specific initial conditions and loss rate, and the curves show the same trend at late times

    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...

  70. [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 ...

  71. [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

  72. [73]

    Fishman, S

    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...