{"id":"1c0274b2-09f4-4abd-a60c-8b818e70f0db","arxiv_id":"2412.01896","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A time-dependent generalized Gibbs ensemble accurately describes weakly dissipative hardcore boson gases, becoming exact as dissipation slows.","lead":"This paper tests a shortcut description for quantum gases that slowly lose and gain atoms, comparing it against exact numerical simulations. It finds the shortcut becomes exact when losses and gains are slow, and stays accurate for weak dissipation, including when atoms flow between regions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exactness claim for γ→0 at fixed γt rests on an uncontrolled adiabatic limit: integrable relaxation is algebraic, and the paper provides no scaling extrapolation from its finite-γ data. The suggested check is a power-law error analysis at smaller γ and larger L.","rationale":"The reader's weakest assumption is exactly the adiabatic separation of timescales in the presence of algebraic relaxation, and I agree that this is the most load-bearing premise. The paper is honest: it explicitly flags the algebraic-relaxation issue in Sec. I and does not claim a proof, only numerical evidence. I found no internal inconsistency in the derivation of Eq. (19) or the analytic solution Eq. (21); the appendices are careful, and the approximate nature of the Renyi-2 non-Gaussianity measure is acknowledged in Sec. V B and Appendix B. The reason the exactness claim remains conditional is not a detected error but an unverified limit: the numerical data are consistent with convergence, but they do not establish the asymptotic scaling in γ. A systematic power-law error analysis at smaller γ and larger L would settle the issue. If the check shows a positive exponent and vanishing extrapolated error, this would justify promoting the verdict toward ACCEPT; without it, CONDITIONAL remains the appropriate assessment.","tokens_in":24524,"tokens_out":6964,"duration_ms":82083,"concrete_test":"Compute the integrated error Δρ(t)/n(t) defined in Eq. (33) at fixed τ=γt∈{0.5,1,2} for γ=0.02, 0.01, 0.005, 0.0025, using large enough chains (e.g. L=80 and L=160) with open boundary conditions, and fit Δρ/n to C γ^p. If p>0 and the fitted curve extrapolates to zero at γ=0, the exactness claim is supported; if p≈0, p<0, or the extrapolated offset is nonzero, Eq. (19) is not the γ→0 limit of the microscopic dynamics. The same runs can also report δ(ρ,ρ_GGE)/L at fixed τ to check that non-Gaussianity vanishes as γ→0.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (19) becomes exact in the limit γ→0 with τ=γt fixed, and accurate for weak finite dissipation. The argument requires that between dissipative events, which occur at intervals Δt~1/γ, the system locally returns to a GGE. The paper itself notes in Sec. I that local relaxation in integrable systems is algebraic rather than exponential, so the relaxation time is formally infinite. No bound or asymptotic estimate is provided for the residual non-GGE component after waiting a time 1/γ. The numerical evidence in Figs. 1, 3, and 4 reaches only γ=0.05 in the homogeneous case, and Fig. 5 uses only two points in the inhomogeneous case, with no extrapolation in γ. If the residual after time 1/γ decays as a power γ^α with α>0, exactness holds; if it decays only logarithmically or saturates, the limiting dynamics need not coincide with Eq. (19). This missing scaling check is the load-bearing gap. The non-Gaussianity proxy in Sec. V B is itself approximate, and the replacement tr(ρ log ρ_GGE)≈tr(ρ_GGE log ρ_GGE) is justified only under the assumption that the largest error component is non-Gaussianity, but this is secondary: the rapidity-distribution comparison already tests the central equation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper benchmarks the time-dependent GGE (t-GGE) approximation for a one-dimensional gas of hardcore bosons subject to weak one-body loss and gain, against tensor-network MPO simulations of the full Lindblad dynamics. It derives closed-form loss and gain functionals for the rapidity distribution, presents an analytic solution of the resulting evolution equation (Eq. (21)), and compares it with numerics in both homogeneous and inhomogeneous settings. In the homogeneous case the comparison uses the rapidity distribution and a Rényi-2 non-Gaussianity measure; in the inhomogeneous case the t-GGE assumption is combined with Generalized Hydrodynamics to describe a domain-wall quench. The paper's central claims are that the t-GGE equation (19) becomes exact as γ → 0 with γt fixed, that it is accurate for weak finite dissipation, and that dissipative GHD captures Euler-scale transport.","tokens_in":24941,"tokens_out":5404,"duration_ms":62741,"significance":"If the claims hold, this is a valuable and much-needed benchmark of an approximation that is widely used but rarely tested beyond small systems. The paper's strengths are the explicit analytic derivation of the dissipation functionals, the closed-form solution (21), the direct parameter-free comparison with MPO data, and the use of a non-Gaussianity measure to go beyond two-point functions. The homogeneous data (Fig. 1, Fig. 3) are visually convincing across several γ values, and the inhomogeneous GHD comparison (Fig. 5) is suggestive. The principal weakness is that the central exactness claim in the γ → 0 limit is not supported by any scaling or extrapolation analysis, a gap that is especially significant because the paper itself notes that local relaxation in integrable systems is algebraic rather than exponential.","major_comments":[{"comment":"The claim that Eq. (19) becomes exact as γ → 0 at fixed γt is the central result, but the paper provides no scaling analysis in γ. The authors themselves state in Sec. I that local relaxation towards a GGE is algebraic, so the adiabatic separation argument is not a proof; the numerical data in Fig. 3 reach only γ = 0.05. I request a systematic extrapolation: compute Δρ(t)/n(t) at fixed γt for a decreasing sequence of γ (e.g., 0.1, 0.05, 0.025, 0.0125) and fit to Δ ~ γ^α, reporting the exponent and its stability with L. Without this, the observed agreement cannot be distinguished from a small but nonzero residual error, and the word 'exact' is stronger than the evidence supports.","section":"§I, §V.A, Figs. 1 and 3"},{"comment":"The MPO simulations are described as 'numerically exact', but no convergence checks or error bars are reported for the bond dimension. The text states in Sec. IV that 'the convergence with respect to the bond dimension is immediate', yet no maximum χ, discarded weight, or χ-dependence is shown for any of the data in Figs. 1, 3, 4, or 5. Because the reported discrepancies are small at small γ, it is important to exclude the possibility that part of the agreement comes from truncation of the operator entanglement. Please provide convergence data for representative parameters (for instance ns = 0, γ = 0.05, γt = 2, and the inhomogeneous case ns = 0.75, γ = 0.05, t = L/2).","section":"§IV"},{"comment":"The non-Gaussianity measure used in Fig. 4 relies on the replacement tr(ρ log ρ_GGE) ≈ tr(ρ_GGE log ρ_GGE). The paper openly acknowledges in Appendix B that this is justified only under the assumption that the largest component of the error is non-Gaussianity. That assumption is plausible, but Fig. 4 is then used as evidence that the t-GGE is accurate beyond two-point functions, so the measure's own approximation weakens the conclusion. I ask the authors to quantify the error in this replacement for the studied system sizes, for instance by computing tr(ρ log ρ_GGE) directly for a few representative states, or to state more carefully that the non-Gaussianity plot is a proxy rather than a direct distance.","section":"§V.B, Eq. (35), Appendix B"},{"comment":"The inhomogeneous test uses only two pairs (L, γ) = (100, 0.1) and (200, 0.05), so the finite-size and finite-γ effects are varied simultaneously and cannot be separated. The collapse in Fig. 5 is visually convincing, but the claim that dissipative GHD becomes exact in the combined Euler and weak-dissipation limit would be considerably strengthened by adding a third point (or by fixing L and varying γ, and fixing γ and varying L, separately) and by reporting the quantitative deviation from the theoretical profile. This is a load-bearing issue because the inhomogeneous setting is one of the two main advertised results.","section":"§VI.B, Fig. 5"}],"minor_comments":[{"comment":"There is a typo: 'collpase' should be 'collapse'.","section":"Appendix C"},{"comment":"The notation is confusing because the analytic solution is written in terms of f(t) and g(t) inside an equation expressed in the rescaled time τ = γt. Please clarify whether t in f(t) and g(t) means the physical time or the rescaled time, and align the notation accordingly.","section":"Eq. (21) and surrounding text"},{"comment":"The label 'free dissipation' is used for the non-interacting fermion comparison; the caption should explain in one sentence that this refers to the model of Refs. [56,58] where gain/loss acts on free fermions and each mode relaxes independently, since the term is otherwise ambiguous.","section":"Fig. 5 caption"},{"comment":"The manuscript contains no data or code availability statement. Providing the numerical data for the curves in Figs. 1, 3, 4, and 5 would improve reproducibility and would allow independent verification of the claimed collapse.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the central idea is sound. The main risk is that the abstract and conclusion phrase the γ → 0 exactness more strongly than the numerical evidence supports without a scaling analysis. If the authors add the requested extrapolation in γ and the bond-dimension convergence checks, I would be willing to support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is the most substantial numerical test of the t-GGE program so far for a model with both gain and loss: the analytic solution in Eq. (21) genuinely extends the loss-only result of Ref. [34], and the MPO benchmarks reach system sizes and times that matter for cold-atom experiments. Second, the central claim about exactness in the gamma-to-0 limit at fixed gamma*t is not actually established. The paper says in Sec. I that local relaxation in integrable systems is algebraic, so the relaxation time is infinite; that means the adiabatic assumption needs a controlled argument, and the numerics stop at gamma=0.05 with no extrapolation in gamma. The stress-test note puts its finger on exactly this gap, and I think the gap is real.\n\nWhat is genuinely good: the derivation of the gain functional is transparent and the formulas are explicit; the rapidity-distribution comparisons in Figs. 1 and 3 are convincing at finite gamma, especially the collapse onto a function of gamma*t; the non-Gaussianity measure, while approximate, is a sensible way to go beyond two-point functions; and the inhomogeneous dissipative GHD benchmark in Fig. 5 is the first at this scale. The paper also earns credit for being honest: it flags the algebraic-relaxation problem and the approximation in Sec. V B itself.\n\nThe soft spots, in proportion. The main one is the uncontrolled adiabatic limit. To support \"exactness\" (even as numerical evidence), the authors should add a scaling analysis: fix tau=gamma*t, plot the error as a function of gamma for smaller gamma and larger L, and show a power-law decay with an estimate of the exponent. Without that, the honest reading is \"consistent with exactness for the studied range,\" not \"exact in the limit.\" Second, the non-Gaussianity measure relies on the replacement tr(rho log rho_GGE) approx tr(rho_GGE log rho_GGE); Appendix B shows the error is O(delta) under the assumption that the largest error is non-Gaussianity, but that assumption is not directly checked. This is secondary because the rapidity distribution already tests the central equation. Third, there are no bond-dimension convergence checks and no code/data artifacts; for a benchmarking paper, that is a moderate omission, not a fatal one.\n\nWho this is for: people working on dissipative integrable systems, GHD with losses/gains, and cold-atom implementations of open quantum gases. It deserves a serious referee. I would send it to review and ask for the gamma-scaling analysis and convergence data before acceptance; with those added, the paper would be a solid reference for the field.","headline":"A careful numerical benchmark of t-GGE for hardcore bosons with gain and loss, with real analytic extension and strong finite-gamma evidence, but the claimed gamma-to-0 exactness rests on an uncontrolled adiabatic limit that needs a scaling check.","tokens_in":25350,"tokens_out":3215,"would_cite":true,"duration_ms":34622,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Dissipation opens integrable models to a time-dependent GGE that becomes exact as the decay rate goes to zero.","keywords":["time-dependent Generalized Gibbs Ensemble","dissipative integrable systems","hardcore bosons","Lindblad master equation","Generalized Hydrodynamics","tensor networks","rapidity distribution","non-Gaussianity"],"falsifier":"A direct test would be to run the same MPO simulation at a much larger system size and longer time, and look for a deviation between the exact rapidity distribution and the t-GGE prediction as γt grows, at fixed small γ. If the error Δρ(t)/n(t) does not continue to shrink with decreasing γ, or if the collapse onto the γt scaling breaks down at any fixed γt, the exactness claim would be falsified. On the experimental side, measuring the rapidity distribution of a lossy 1D Bose gas prepared in a thermal state and comparing to Eq. (21) would be a decisive check.","tokens_in":24343,"feed_emoji":"⚛️","tokens_out":1632,"duration_ms":18906,"temperature":0.7,"pith_summary":"The paper argues that a simple description of a dissipative quantum gas—the time-dependent Generalized Gibbs Ensemble (t-GGE)—becomes exact in the limit of weak, adiabatic dissipation, and remains accurate for weak but finite dissipation. It tests this idea on a chain of hardcore bosons with one-body atom loss and gain, a model that is integrable when closed but not exactly solvable when open. Using tensor-network simulations for system sizes and times far beyond previous checks, the authors show that the rapidity distribution, which encodes the many-body state, collapses onto the t-GGE prediction as the dissipation rate γ goes to zero at fixed γt. They also verify that the state stays close to Gaussian, the key assumption behind the t-GGE, and that combining the approach with Generalized Hydrodynamics correctly captures transport in an inhomogeneous setup. The significance is that a closed, analytically tractable equation may describe the full open dynamics of an integrable many-body system in a regime relevant to cold-atom experiments.","feed_headline":"A closed equation predicts dissipative gas dynamics","feed_subtitle":"As atom loss and gain weaken, an integrable Bose gas relaxes to a time-dependent GGE that becomes exact, verified with tensor networks.","key_machinery":"The central object is the time-dependent Generalized Gibbs Ensemble (t-GGE), parameterized by the rapidity distribution ρ(t, k) of the Jordan-Wigner fermions that diagonalize the hardcore boson Hamiltonian. The key identity is Eq. (19), the closed nonlinear evolution equation for ρ(t, k) whose dissipation functional F[ρ] is expressed in terms of a Hilbert transform of the rapidity distribution. This functional captures the nonlocal-in-rapidity effect of the bosonic dissipators, which carry Jordan-Wigner strings that prevent Gaussianity from being preserved exactly. The t-GGE assumption—that between rare dissipative events the system locally relaxes to a GGE—turns the open dynamics into a closed equation for the GGE parameters.","core_discovery":"The central claim is that for a one-dimensional hardcore boson gas subject to weak one-body gain and loss, the t-GGE approximation—where the density matrix is taken to remain a Gaussian state of Jordan-Wigner fermions—becomes exact in the limit γ → 0 with γt fixed, and quantitatively accurate for weak but finite γ. This is established by deriving a closed evolution equation, Eq. (19), for the rapidity distribution, solving it analytically, and then benchmarking it against numerically exact MPO simulations. The same approximation, combined with Generalized Hydrodynamics, reproduces the emergent Euler-scale transport of a domain-wall initial state. The paper's evidence includes collapse of the rapidity distributions to a universal curve at fixed γt, agreement of the non-Gaussianity measure showing that the state becomes increasingly Gaussian as γ decreases, and matching density profiles in the inhomogeneous case.","pith_inferences":["The t-GGE framework likely extends to other integrable models with one-body dissipation where a closed equation for the rapidity distribution can be derived, such as the Lieb-Liniger gas with losses, where analogous formulas are already known.","The numerical evidence suggests that the algebraic relaxation problem does not invalidate the adiabatic assumption for this model at the studied scales; a rigorous proof might come from a modified ETH-like argument that accounts for the slow algebraic tails without requiring true exponential relaxation.","A testable extension would be to measure higher-order correlation functions (beyond the non-Gaussianity proxy) in a cold-atom experiment and compare their decay to the t-GGE prediction, which would probe the Gaussian approximation beyond two-point functions.","The observed cancellation at balanced gain and loss hints that there may be an exact hidden symmetry of the Lindbladian when ns = 1/2, which could be proven and used to simplify other open integrable models."],"forward_implications":["If the t-GGE equation is exact in the adiabatic limit, then the entire dissipative dynamics of the gas is captured by a single scalar function ρ(γt, k), computable in closed form for any initial state.","The analytic solution, Eq. (21), provides a fast and parameter-free prediction for rapidity distributions and two-point functions in cold-atom experiments with loss and gain, without needing full simulation.","Combining t-GGE with Generalized Hydrodynamics yields a dissipative GHD equation, Eq. (24), that correctly predicts large-scale density profiles in inhomogeneous settings.","The accuracy improves when loss and gain are balanced (ns = 1/2), because many terms in the master equation cancel, which may guide experimental choices to maximize predictability."],"supporting_citations":[{"why":"Provides the original derivation of the t-GGE equations for atom losses in the one-dimensional Bose gas, establishing the functional form that the present paper extends to gains.","marker":"[28]"},{"why":"Derives the loss-only version of the t-GGE equation for hardcore bosons; the present work adapts its analytical solution technique to include gains.","marker":"[34]"},{"why":"Introduces Generalized Hydrodynamics for integrable systems, the framework that the paper couples with t-GGE to describe inhomogeneous transport.","marker":"[50]"},{"why":"Co-introduces GHD, supplying the hydrodynamic equations that the paper extends with a dissipative term.","marker":"[51]"},{"why":"Provides earlier small-system quantum trajectory checks of t-GGE for lossy quantum gases, serving as the baseline that this paper improves upon with much larger system sizes.","marker":"[30]"},{"why":"Documents algebraic (non-exponential) local relaxation in integrable systems, the main theoretical objection to the t-GGE assumption that the paper addresses numerically.","marker":"[47]"}],"fun_headline_variants":["Time-dependent GGE exact for weakly dissipative bosons","Weak atom loss makes t-GGE exact in 1D Bose gas","t-GGE proven exact for bosons under adiabatic dissipation","Dissipative Bose gas relaxes to exact t-GGE at slow loss"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The adiabatic separation of timescales: after each rare dissipative event, the system locally relaxes to a GGE before the next event, which is nontrivial because local relaxation in integrable systems is algebraic rather than exponential, so the relaxation time is formally infinite.","fun_headline_variants_meta":{"raw":{"variants":["Time-dependent GGE exact for weakly dissipative bosons","Weak atom loss makes t-GGE exact in 1D Bose gas","t-GGE proven exact for bosons under adiabatic dissipation","Dissipative Bose gas relaxes to exact t-GGE at slow loss"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1462,"prompt_tokens":861,"completion_tokens":601,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":477,"completion_tokens_details":{"reasoning_tokens":526}},"tokens_in":477,"tokens_out":601,"duration_ms":6242,"temperature":1.0,"reasoning_tokens":526,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:51:19.503339+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to run the same MPO simulation at a much larger system size and longer time, and look for a deviation between the exact rapidity distribution and the t-GGE prediction as γt grows, at fixed small γ. If the error Δρ(t)/n(t) does not continue to shrink with decreasing γ, or if the collapse onto the γt scaling breaks down at any fixed γt, the exactness claim would be falsified. On the experimental side, measuring the rapidity distribution of a lossy 1D Bose gas prepared in a thermal state and comparing to Eq. (21) would be a decisive check.","supporting_citations":[{"cited_title":"Weimer, A","cited_arxiv_id":null,"evidence_quote":"Provides the original derivation of the t-GGE equations for atom losses in the one-dimensional Bose gas, establishing the functional form that the present paper extends to gains."},{"cited_title":"Rossini, A","cited_arxiv_id":null,"evidence_quote":"Derives the loss-only version of the t-GGE equation for hardcore bosons; the present work adapts its analytical solution technique to include gains."},{"cited_title":"Fagotti and F","cited_arxiv_id":null,"evidence_quote":"Co-introduces GHD, supplying the hydrodynamic equations that the paper extends with a dissipative term."},{"cited_title":"Lange, Z","cited_arxiv_id":null,"evidence_quote":"Provides earlier small-system quantum trajectory checks of t-GGE for lossy quantum gases, serving as the baseline that this paper improves upon with much larger system sizes."},{"cited_title":"Rigol, V","cited_arxiv_id":null,"evidence_quote":"Documents algebraic (non-exponential) local relaxation in integrable systems, the main theoretical objection to the t-GGE assumption that the paper addresses numerically."}],"review_version":1}