{"id":"9affc3c9-56f6-4fa0-87f6-9156c177bc89","arxiv_id":"2509.10855","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In the 1+1D massive Schwinger model, tensor network simulation of a localized excitation reveals Bjorken-like hydrodynamic flow for small fermion mass, but not for large mass.","lead":"A quantum simulation of the massive Schwinger model shows that a localized energy burst evolves into a flow matching Bjorken hydrodynamics when the coupling is strong. The result is a step toward deriving fluid behavior of quark-gluon plasma from first-principle quantum dynamics rather than classical assumptions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Missing D and a convergence checks make the MPS evolution unreliable; truncation-induced dissipation could mimic hydrodynamization, so the central claim is not yet established.","rationale":"The paper's strongest claim is that a pure-state strong-coupling Schwinger model shows Bjorken-type hydrodynamization. The observable evidence is entirely numerical: T^{mu nu} profiles from MPS evolution. Therefore, the reliability of that evolution is the gate-keeping assumption. The authors fix D=500 and a=0.5/g based on previous studies but provide no convergence data for this particular initial condition and observable. MPS truncation is an approximation that can introduce artificial dissipation: discarding small Schmidt values removes entanglement correlations and can drive the state toward a lower-entanglement, more 'thermal-looking' regime. This is especially dangerous in a real-time simulation of a quench from a localized UV excitation, where the initial state generates extensive entanglement. If the entropy across typical cuts approaches log2(D), the simulation loses fidelity, and the decline of Pi/(epsilon+P) may be spurious. The secondary issues—fitted Bjorken normalization, lack of event averaging, and the companion-paper reference (arXiv:2509.xxxxx) for numerical details—are real but would only affect the strength of the interpretation; they are not as fundamental as the numerical convergence question. The manuscript itself acknowledges the single-state nature of the analysis, but that is a feature of the pure-state claim, not a fatal flaw. Since the central concern can be settled by additional runs, the condition attached to the paper is appropriate. The reader's CONDITIONAL verdict is retained; no adjustment is needed.","tokens_in":11313,"tokens_out":9132,"duration_ms":83519,"concrete_test":"Repeat the quench with D=1000 and D=1500 at a=0.5/g, and with a=0.25/g (N=200) at D=1000. Compute T^{tt}, T^{tz}, T^{zz} at z=0 for tau=5,10,20 g^-1; also track the time tau_hydro when |Pi|/(epsilon+P) first remains below 0.1 for at least 2 g^-1. If the profiles shift by more than 10% or tau_hydro changes by more than 2 g^-1, the hydrodynamization seen in Figs. 3-5 is not converged. Additionally, report the largest discarded Schmidt weight and the central-cut entanglement entropy S(t); if S approaches log2(D) or the discarded weight exceeds 10^-4, D=500 is saturated and the observed relaxation may be truncation-induced.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that a strongly coupled Schwinger model hydrodynamizes and locally thermalizes—rests entirely on the fidelity of the MPS real-time evolution with bond dimension D=500 and lattice spacing a=0.5/g. No convergence checks in D or a are reported. This is not a minor technical gap: fixed-bond-dimension TEBD is a dissipative approximation. Each truncation discards small Schmidt weights, effectively performing an environment-induced projection that removes entanglement and can generate apparent entropy growth and relaxation. In an isolated unitary system, such truncation-induced decoherence can mimic the approach to local thermal equilibrium even when the true state is far from it. The worry is sharpened because the initial excitation (Eq. 5) has width a=0.5/g, i.e., exactly one lattice spacing; the initial UV structure couples strongly to high-momentum modes whose entanglement must be resolved. If D=500 becomes saturated, the late-time 'hydrodynamic' profiles (Figs. 3–5) could be an artifact of the retained subspace rather than a property of the continuum gauge theory. Since the paper offers a first-principles statement about quantum dynamics, this numerical provenance issue is the load-bearing assumption. The fitting of the Bjorken normalization and the single-state fluctuations are secondary; they would weaken the interpretation but not invalidate it if the data are converged.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses matrix product state (MPS) time-evolving block decimation to simulate the real-time quantum evolution of the massive Schwinger model on a lattice with N=100 sites, lattice spacing a=0.5/g, and bond dimension D=500. Starting from a local excitation placed on the vacuum, it computes the full stress-energy tensor T^{\\mu\\nu}, defines the fluid velocity in the Landau-Lifshitz frame, and compares the resulting energy density, velocity, and bulk pressure with the generalized Bjorken flow with finite collision time. For the strong-coupling case m=g/10, the authors report that v_z approaches z/t, the rapidity profile develops a plateau, the mid-rapidity energy density follows \\varepsilon \\propto \\tau^{-1-c_s^2} with c_s^2=0.5 from the independently computed equation of state, and the bulk pressure decays within a \\tau^{-2} envelope; they interpret the timescale \\tau~10 g^{-1} as the onset of hydrodynamics and find that this behavior is absent for m=2g.","tokens_in":11563,"tokens_out":4211,"duration_ms":41909,"significance":"If the numerical results are converged, this would be a notable first-principles demonstration that an isolated, strongly coupled gauge theory can hydrodynamize without external sources or semiclassical assumptions, and the m=g/10 versus m=2g contrast provides a falsifiable, qualitative prediction about when Bjorken-like behavior emerges. The paper's strengths include computing the equation of state and c_s^2 independently on the same lattice rather than fitting the exponent, checking energy-momentum conservation at the operator level, and directly using the stress tensor rather than proxy observables. The main weakness is numerical provenance: the central claim rests on a single bond dimension D=500 and a single lattice spacing a=0.5/g with no reported convergence checks, and the comparison to the Bjorken power law uses an unspecified amplitude normalization. These issues are fixable, so the manuscript is suitable for major revision rather than rejection.","major_comments":[{"comment":"The central claim rests on real-time evolution with D=500 and a=0.5/g, but the manuscript reports no convergence checks in either the bond dimension or the lattice spacing. Fixed-bond-dimension TEBD is a dissipative approximation: each truncation discards Schmidt weights and can induce apparent relaxation and entropy growth, which could mimic hydrodynamization even if the true unitary evolution is far from local equilibrium. The concern is sharpened because the initial excitation in Eq. (5) has width a, exactly one lattice spacing, so the initial state contains UV structure whose entanglement must be resolved. Please add convergence tests (for example D=250, 500, and 1000 at fixed volume, and a=0.25/g, 0.5/g, and 0.75/g at fixed physical volume) and report the truncation error or discarded weight over the full simulated time; without these, the late-time profiles in Figs. 3-5 cannot be taken as properties of the continuum model.","section":"Model and Numerical Setup"},{"comment":"The comparison between the simulated energy density and the Bjorken power law \\varepsilon \\propto \\tau^{-1-c_s^2} requires an amplitude normalization that is not specified in the text. If the normalization is simply matched to the simulation at one time, the agreement is partly a fit and the power-law claim is weakened. Please state the normalization procedure explicitly, or show the unnormalized product \\varepsilon(t) \\tau^{1+c_s^2} so that the reader can judge whether the exponent is actually reproduced without a free amplitude.","section":"Hydrodynamic Evolution, Fig. 5(b)"},{"comment":"All results are obtained from a single pure initial state, and the text itself notes O(0.1) quantum fluctuations. Because the stated motivation is to mimic event-by-event heavy-ion collisions, a single state cannot distinguish state-specific fluctuations from generic hydrodynamic behavior. Please provide either an average over several initial excitations or a quantitative sensitivity study varying the excitation width and profile, so that the claimed onset time \\tau~10 g^{-1} is not tied to one specially chosen initial condition.","section":"Summary and Discussion"}],"minor_comments":[{"comment":"The initial state is defined by a phase profile \\phi(z) of width a, but the text does not state the resulting deposited energy or its relation to the scales in the problem; adding this would make the analogy to heavy-ion energy deposition more concrete.","section":"Model, Eq. (5)"},{"comment":"Reference [98] is listed as arXiv:2509.xxxxx; please update it to the actual companion paper identifier before publication so that the numerical details can be verified.","section":"References"},{"comment":"The caption of Fig. 3 does not define all curves or panels clearly; please label the three rows as (a), (b), and (c) and specify the plotted quantity in each panel.","section":"Hydrodynamic Evolution, Fig. 3"},{"comment":"The rapidity-distribution panel would benefit from a legend or color bar identifying the values of \\tau-\\tau_{\\rm ini}, since the caption lists several times but does not indicate which curve corresponds to which time.","section":"Hydrodynamic Evolution, Fig. 5(a)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the hep-ph community if the convergence issue is resolved. I would request the convergence checks and normalization details before further consideration; the single-state issue is important but could be addressed with a modest additional computation. The companion paper reference should be available to the referee, since the methods are partly delegated to it."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Interesting thing here is that the paper computes the full T^mu nu in real time and directly checks the hallmark signatures of Bjorken flow: the v_z = z/t velocity profile, the rapidity plateau, the energy density power law with the independently computed c_s^2, and the decay of bulk pressure. That is more than previous tensor-network studies of jets and collective motion in the Schwinger model, and the mass dependence (strong coupling / small mass shows hydro-like behavior, large mass does not) is a clean qualitative contrast. The EoS is computed separately on the same lattice, so the Bjorken exponent is not fitted; that is a genuine strength. The authors are also honest that they run a single pure state and that late-time fluctuations become visible, which they attribute to event-by-event averaging in nuclear collisions. Where I agree with the skeptic: the refusal to show convergence checks is the weak point, and it is not a minor omission. Fixed-bond-dimension TEBD is a dissipative approximation. Each truncation removes small Schmidt weights; that can mimic entropy growth and relaxation even if the true unitary dynamics would not. The initial excitation has width equal to one lattice spacing, so the early-time evolution has to resolve short-distance entanglement, and D=500 might saturate. Without varying D and a, the apparent hydrodynamic onset in Figs. 3-5 could in principle be an artifact of the retained subspace. That said, the observed signatures are more structured than simple relaxation to local equilibrium: the velocity field tracks the Bjorken form over many sites, and the bulk pressure oscillates and then falls as tau^{-2}. These are nontrivial and hard to fake, but the numerical provenance still needs to be demonstrated. The other soft spots are milder. The Bjorken curve in Fig. 5(b) is presumably normalized to the simulation amplitude, but the text does not say so; that is a fitting detail that should be stated. The current results are for a single initial condition, so all the quantum fluctuations in the figures are not statistical error bars. And the comparison with the related tensor-network paper by Janik et al. (Ref. [86]) is too terse; what exactly is the new observable or interpretation here deserves a sentence or two. Bottom line: this is a plausible and potentially important proof-of-principle calculation, but the central claim rests on numerical fidelity that the paper does not yet establish. I would send it to a serious referee, asking specifically for convergence checks in bond dimension and lattice spacing and for a clear statement of the normalization in Fig. 5(b). If those checks survive, this becomes a nice contribution. If not, the result could be fortuitous. The paper is worth engaging; it just needs another round of evidence.","headline":"A serious tensor-network study that plausibly sees Bjorken-like hydrodynamic onset in the massive Schwinger model from a local excitation, but missing convergence checks in bond dimension and lattice spacing leave the central claim not yet fully established.","tokens_in":791,"tokens_out":2111,"would_cite":true,"duration_ms":30589,"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":"A strongly coupled quantum gauge theory can develop Bjorken-like flow on its own, the paper claims.","keywords":["massive Schwinger model","Bjorken flow","hydrodynamization","tensor networks","matrix product states","real-time lattice gauge theory","thermalization","quark-gluon plasma"],"falsifier":"Repeat the same evolution with bond dimension $D=1000$ or $D=2000$ and with lattice spacing $a=0.25/g$; if the energy-density scaling $\\varepsilon\\propto\\tau^{-1-c_s^2}$ or the bulk-pressure decay $\\Pi\\propto\\tau^{-2}$ changes materially, the reported hydrodynamic onset is a truncation artifact. Averaging over many localized-excitation initial states would also test whether the large-rapidity quantum fluctuations are event-by-event noise that disappears in the ensemble.","tokens_in":11111,"feed_emoji":"🌊","tokens_out":9130,"duration_ms":75653,"temperature":0.7,"pith_summary":"Using a real-time tensor-network simulation of the massive Schwinger model, this paper asks whether a pure gauge-theory state can turn into a fluid by itself. It prepares the vacuum with one localized energy bump, mimicking energy deposited by colliding nuclei, and follows the full unitary evolution. For strong coupling (fermion mass $m=g/10$), the system shows exactly the signatures of generalized Bjorken flow: the flow velocity approaches $v_z=z/t$, the energy density forms a widening rapidity plateau and decays with proper time as $\\varepsilon\\propto\\tau^{-1-c_s^2}$, and the bulk pressure falls inside the Navier-Stokes bound $\\Pi\\propto\\tau^{-2}$. For weak coupling and heavy mass ($m=2g$), none of this appears. The point is that local thermal equilibrium and hydrodynamization would then emerge from first-principle quantum dynamics, not from external sources or classical approximations.","feed_headline":"Strong coupling turns a quantum gauge theory into a flowing fluid","feed_subtitle":"A first-principle simulation of a QCD analog finds fluid behavior emerging from a localized excitation.","key_machinery":"The load-bearing object is the full energy-momentum tensor $T^{\\mu\\nu}(t,z)$, computed as the expectation value of lattice stress-tensor operators minus vacuum contributions and decomposed in the Landau-Lifshitz frame into energy density $\\varepsilon$, fluid velocity $u^\\mu$, and bulk pressure $\\Pi$. The test is whether these components satisfy the generalized Bjorken-flow solution: $u^z=z/t$, $\\varepsilon$ depending only on proper time $\\tau$, and $\\Pi$ bounded by the Navier-Stokes formula $\\zeta s\\theta\\propto\\tau^{-2}$. Numerically, this is carried by staggered (Kogut-Susskind) fermions on a 100-site lattice and by evolving a matrix product state with bond dimension $D=500$ using the time-evolving block decimation algorithm; the equation of state $P(\\varepsilon)$ and the speed of sound $c_s^2$ are computed separately from finite-temperature traces.","core_discovery":"The central claim is that in the strongly coupled massive Schwinger model at $m=g/10$, a pure state prepared as a localized excitation on the vacuum evolves, without external sources, into a state whose energy-momentum tensor matches the generalized Bjorken solution with a finite collision time. The simulation finds $v_z=z/t$ inside the light cone, a mid-rapidity plateau in energy density that widens with time, mid-rapidity $\\varepsilon$ scaling as $\\tau^{-1-c_s^2}$ with $c_s^2\\approx 0.5$ taken from the equation of state, and a bulk pressure that starts large, then oscillates and decays within the Navier-Stokes band $\\Pi\\propto\\tau^{-2}$ after $\\tau-\\tau_{\\mathrm{ini}}\\gtrsim 1/g$. Local thermal equilibrium, judged by the effective pressure matching the thermal equation of state, is reached on a timescale $\\tau\\sim 10/g$, which matches the thermalization time of the quantum momentum distribution. The same protocol at $m=2g$ shows neither Bjorken energy scaling nor decaying bulk-pressure oscillations, supporting the conclusion that strong coupling, equivalently small fermion mass, drives the onset of hydrodynamic behavior.","pith_inferences":["The paper leaves implicit that averaging over many localized-excitation initial states, analogous to event averaging in heavy-ion collisions, should suppress the $O(0.1)$ quantum fluctuations seen in the large-rapidity tails and sharpen the Bjorken comparison; this is directly testable with the same simulation setup.","Because 1+1 dimensions have no shear tensor, this test probes only bulk viscosity; whether the same rapid onset survives in 2+1 or 3+1 non-Abelian theories, where shear modes dominate, remains an open question the paper does not address.","A natural falsifying check is to repeat the simulation at larger bond dimension and smaller lattice spacing; the absence of published convergence checks means the present claim is only as strong as the truncation control.","One could probe local thermalization more stringently by extracting local entropy density from entanglement or particle-number fluctuations, going beyond the stress-tensor criteria used here."],"forward_implications":["If the central claim holds, a pure, isolated, strongly coupled gauge theory can hydrodynamize without any external source, semi-classical input, or assumption of local equilibrium, meaning the fluid description is emergent.","Hydrodynamization, judged by small viscous corrections and a Bjorken velocity profile, sets in before full global thermalization, consistent with the picture of early fluid behavior in far-from-equilibrium heavy-ion collisions.","The contrast between $m=g/10$ and $m=2g$ identifies strong coupling and light fermions as the microscopic ingredient responsible: enhanced scattering and pair production accelerate thermalization, while heavy fermions do not thermalize on the accessible timescale.","The onset time $\\tau\\sim 10/g$ coincides with the thermalization time of the quantum momentum distribution, suggesting that a single dynamical mechanism controls both hydrodynamization and thermalization.","The calculation provides a benchmark set of observables, $\\varepsilon$, $u^\\mu$, and $\\Pi$, for future real-time non-perturbative simulations of hydrodynamization in higher-dimensional gauge theories."],"supporting_citations":[{"why":"Defines Bjorken flow and its scaling laws ($v_z=z/t$, $\\varepsilon\\propto\\tau^{-1-c_s^2}$) that the simulation is matched against.","marker":"[100]"},{"why":"Generalizes the Bjorken solution to finite collision time and predicts the rapidity-plateau structure that the simulation observes.","marker":"[101]"},{"why":"Supplies the thermalization time of the quantum momentum distribution in the Schwinger model, which the hydrodynamic-onset timescale is compared with.","marker":"[87]"},{"why":"Supplies the Landau-Lifshitz stress-tensor decomposition and the Navier-Stokes relation $\\Pi=\\zeta s\\theta$ used to define and bound the bulk pressure.","marker":"[8]"},{"why":"Establishes the lattice spacing $a=0.5/g$ as a regime where thermalization has been observed in the Schwinger model, fixing the lattice cutoff.","marker":"[96]"},{"why":"Supplies the time-evolving block decimation algorithm used to evolve the matrix product state in real time.","marker":"[97]"},{"why":"Establishes matrix product state representations for gauge field theories, underpinning the tensor-network simulation of the Schwinger model.","marker":"[92]"},{"why":"Companion paper defining the numerical discretization, operator representation, initial-state preparation, real-time evolution, and finite-temperature equation-of-state calculations on which the stress-tensor profiles rely.","marker":"[98]"}],"fun_headline_variants":["Strong coupling drives Bjorken flow onset in Schwinger model","Quantum evolution of QCD analog shows hydrodynamic onset","Thermalization and Bjorken flow in strongly coupled Schwinger model","From localized quantum state to hydrodynamic flow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on the assumption that the simulation's truncation limit, bond dimension $D=500$, and lattice spacing $a=0.5/g$ faithfully represent the true real-time quantum evolution; no convergence checks with $D$ or $a$ are shown, so if the truncation error is substantial, the apparent hydrodynamization could be a numerical artifact rather than a physical property of the continuum model.","fun_headline_variants_meta":{"raw":{"variants":["Strong coupling drives Bjorken flow onset in Schwinger model","Quantum evolution of QCD analog shows hydrodynamic onset","Thermalization and Bjorken flow in strongly coupled Schwinger model","From localized quantum state to hydrodynamic flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00051,"raw_usage":{"total_tokens":2482,"prompt_tokens":948,"completion_tokens":1534,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":1469}},"tokens_in":564,"tokens_out":1534,"duration_ms":10091,"temperature":1.0,"reasoning_tokens":1469,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:52:09.234847+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same evolution with bond dimension $D=1000$ or $D=2000$ and with lattice spacing $a=0.25/g$; if the energy-density scaling $\\varepsilon\\propto\\tau^{-1-c_s^2}$ or the bulk-pressure decay $\\Pi\\propto\\tau^{-2}$ changes materially, the reported hydrodynamic onset is a truncation artifact. Averaging over many localized-excitation initial states would also test whether the large-rapidity quantum fluctuations are event-by-event noise that disappears in the ensemble.","supporting_citations":[],"review_version":2}