{"id":"95219f29-5258-410c-af79-9e7e2288d13b","arxiv_id":"2509.10835","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Tensor-network simulations of the massive Schwinger model show Bjorken-like hydrodynamics at small m/g and sharp dynamical order parameters marking the parity-breaking phase transition near m/g=0.33 at θ=π.","lead":"This paper simulates the real-time evolution of the massive Schwinger model using tensor networks, starting from a localized energy burst. It finds signs of Bjorken-like fluid flow at small fermion mass and dynamical markers of the parity-breaking transition in the electric field.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Bjorken-flow claim is currently supported only by qualitative v_z≈z/t and Π→0 plots; the decisive Milne-space analysis is deferred to an unavailable companion, and the fixed-volume convergence tests cannot detect boundary reflections.","rationale":"I read the paper in good faith. The tensor-network methods are standard, the convergence checks are extensive, and the EoS benchmark against exact diagonalization is a real strength. The phase-transition section is suggestive and the dynamical order parameters are plausible, though the fixed external critical point and limited time window leave room for protocol dependence. The most load-bearing unresolved issue is the hydrodynamic identification, because it is the paper's headline claim and its quantitative support is deferred to a companion that is not available in this submission. The reader's weakest assumption identified the finite-size/boundary aspect; I agree with that partially and add that the Milne-space observable needed to define Bjorken flow is not shown. This does not require changing the reader's conditional verdict: the appropriate outcome is still conditional acceptance, with the missing Milne/volume analysis as the condition. I am not recommending rejection because the underlying data could plausibly support the claim, and the qualitative plots are consistent with it.","tokens_in":21241,"tokens_out":14306,"duration_ms":139275,"concrete_test":"Re-analyze the stored TEBD wavefunctions for m/g = 0.1, θ = 0 to compute T^{ττ} in Milne coordinates (τ = sqrt(t^2 - z^2), η = atanh(z/t)) and plot ε(τ, η) for η ∈ {-1, -0.4, 0.4, 1} over t ∈ [0, 25/g]. If the curves are independent of η and follow the ideal-hydro evolution with the computed EoS, Bjorken flow is supported; if they spread or deviate, the claim is weakened. Separately, rerun with N = 200 at a = 0.5/g (doubling physical volume while keeping lattice spacing fixed) and compare the central-region ε(τ) and Π(t) up to t = 25/g; any difference from N = 100 indicates boundary contamination.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Bjorken flow is defined by boost invariance: in Milne coordinates (τ, η), all dimensionless observables should be independent of η, and ε should evolve only with τ. The present paper does not show this. The evidence for hydrodynamics in Sec. IV.B and Fig. 4 is qualitative: the velocity profile is visually close to v_z≈z/t and the bulk pressure Π tends to zero in a central window. But v≈z/t is a necessary, not sufficient, signature of Bjorken flow; it is shared by self-similar 1+1D expansions such as simple rarefaction waves from a localized source. The quantitative Milne-space analysis is explicitly deferred to the companion letter [53], whose arXiv identifier is a placeholder (2509.xxxxx) and therefore not checkable in this submission. A second, independent gap: the lattice has physical length Na = 50/g and light-crossing time t = 25/g, which is exactly the end of the simulation window. The lattice-size test in Appendix A (N = 80–140) keeps the total volume Na fixed, so it tests discretization but cannot detect boundary reflections. Thus the central hydrodynamic claim is underdetermined by the evidence in this paper as it stands.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports real-time tensor network simulations (TEBD) of the massive Schwinger model on a lattice of N=100 sites with spacing a=0.5/g, starting from a gauge-invariant, charge-neutral localized excitation obtained by applying U(pi,nu_50)U(pi,nu_51) to the ground state. It computes local expectation values of the energy-momentum tensor, extracts Landau-frame energy density, flow velocity, and effective pressure, and constructs the thermal equation of state by Thermofield-Double purification. The main physical claims are: (i) at small m/g the system displays Bjorken-like collective expansion, characterized by v_z about z/t and decaying bulk pressure, with the behavior degrading as m/g increases; and (ii) at theta=pi the electric-field and charge-density dynamics show a sharp change near m/g approximately 0.33, interpreted as a dynamical signal of the spontaneous parity-breaking transition, with string breaking and screening on one side and stable charge separation on the other. Numerical convergence with respect to Lmax, bond dimension, time step, and lattice spacing is documented in Appendix A, and the equation of state is benchmarked against exact diagonalization.","tokens_in":21498,"tokens_out":9187,"duration_ms":83069,"significance":"The numerics are careful: Appendix A reports convergence over electric-field cutoff, time step, bond dimension, and lattice spacing; the equation of state is benchmarked against exact diagonalization up to N=16; and SVD truncation errors and energy conservation are quantified. If the hydrodynamic interpretation is correct, the paper offers a concrete (1+1)-dimensional gauge-theory realization of hydrodynamization and a dynamical probe of the parity phase transition. However, the central Bjorken-flow claim is currently supported only by qualitative features, with the quantitative Milne-space analysis deferred to a companion letter that is not checkable in this submission. As it stands, the manuscript is therefore a valuable numerical study of the Schwinger model, but the headline hydrodynamic conclusion is underdetermined by the evidence presented.","major_comments":[{"comment":"The identification of Bjorken flow rests on visual inspection of v_z approximately z/t and Pi going to zero. Bjorken flow in (1+1)D is defined by boost invariance, meaning all dimensionless observables should be independent of Milne rapidity eta and depend only on tau; v_z approximately z/t is a necessary but not sufficient condition, since self-similar non-hydrodynamic expansions can share it. The manuscript explicitly defers the Milne-space analysis to the companion Letter [53], whose arXiv identifier (2509.xxxxx) is a placeholder and cannot be checked. Please include the fixed-eta analysis in this manuscript, for example plots of epsilon(tau,eta), P+Pi(tau,eta), or T^tau tau and T^eta eta at fixed eta over a wide rapidity window, or explicitly weaken the claim from 'exhibits Bjorken flow' to 'shows Bjorken-flow-like signatures'.","section":"Section IV.B, Fig. 4"},{"comment":"The test with N=80 to 140 keeps Na=50/g fixed, so it changes the lattice spacing while holding the total volume constant. It therefore cannot detect finite-volume artifacts or boundary reflections. With open boundary conditions, the light-crossing time from the central excitation to either boundary is t=25/g, exactly the final time of the simulation, so the late-time portions of the full profiles in Fig. 4 may be affected by boundary physics, and no test at a larger physical volume is provided. Please add a fixed-a, larger-Na comparison for m/g=0.1, or a boundary-contamination analysis, or restrict the hydrodynamic statement to the central region and to times before reflected signals return.","section":"Appendix A, lattice-independence test"},{"comment":"The power-law fit uses x_c=0.33 as a fixed input taken from known static results, so the fit is a consistency check rather than an independent dynamical determination of the critical point; the exponent, amplitude, and averaging window are chosen or fitted. Moreover, the two 'order parameters' are both derived from the same electric-field data and are not statistically independent, so their 'close agreement' should not be presented as independent confirmation. Please state these limitations explicitly and present the scaling result as evidence consistent with, rather than a determination of, the known transition.","section":"Section IV.C, Fig. 5"}],"minor_comments":[{"comment":"The quantity tau_ini used for the constant-proper-time lines is not defined in the text; please define it and state the reference time from which tau is measured.","section":"Fig. 4 caption"},{"comment":"Table I lists D=500 for all real-time evolution, while the Fig. 8 caption states D=250 for m/g=0.1 and 2.0 and D=500 for m/g=0.5; please clarify the production bond dimensions used for each parameter set.","section":"Table I and Fig. 8 caption"},{"comment":"The statement that 'the medium spreads out within the lightcone' is difficult to verify from the color-scale panels of Fig. 4; consider overlaying the light-cone boundaries on the energy-density panels.","section":"Section IV.B"},{"comment":"The term 'order parameter' is used loosely for the time-averaged central electric field and the growth rate of total electric field; neither is a true order parameter for spontaneous parity breaking in a finite system, so consider using 'dynamical indicators' and identifying one as the primary observable.","section":"Section IV.C"}],"recommendation":"major_revision","confidential_remarks":"The strongest and most novel claim, quantitative Milne-space hydrodynamics, is outsourced to a companion with a placeholder identifier. If the split-publication strategy is intentional, the current submission does not stand alone for its main claim, and I would ask the editor to require either integration of that analysis or a correspondingly softened abstract and conclusion before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things worth knowing up front. The numerics are careful and the phase-transition part of the paper is the stronger half. The Bjorken-flow claim, as it stands, is plausible but not yet demonstrated.\n\nWhat is new and good: the systematic mass scan at theta=pi with two specific order parameters (time-averaged central electric field and growth rate of the total electric field) rising sharply near m/g ~ 0.33 is a genuinely useful dynamical probe of the parity transition. The tensor-network work is methodical: gauge-invariant MPS, convergence tests covering Lmax, bond dimension, time steps, EoS benchmarked against exact diagonalization, and reported SVD errors and energy conservation. That is real evidence and the paper deserves credit for it. The authors are also honest about limits, explicitly saying the critical exponent is not determined rigorously and deferring Milne-space analysis to the companion.\n\nThe soft spots are real but mostly concentrated in the hydrodynamic interpretation. The evidence for Bjorken flow is qualitative: v_z close to z/t and bulk pressure tending to zero in a central window. Those are necessary signatures of boost-invariant expansion, but v ~ z/t is also shared by simple self-similar 1+1D rarefaction waves; the Milne-space check (observables independent of spacetime rapidity) is exactly what would discriminate, and it is deferred to a companion Letter whose arXiv ID is a placeholder. On top of that, the lattice length is 50/g and the light-crossing time from the center to the boundary is 25/g, which is the end of the simulation window. The lattice-size test in Appendix A varies N while keeping total volume fixed, so it checks discretization but cannot rule out boundary reflections contaminating late-time data. Those two gaps make the central hydrodynamic claim underdetermined in this submission.\n\nThe phase-transition part is more robust. The critical point is taken from prior literature and the power-law exponent is fitted, so the scaling claim is a consistency check rather than an independent determination, but the sharp simultaneous rise of two distinct observables is a clean signal. I do not see a circularity problem; using the known critical point as input is acceptable when the goal is to show a dynamical probe sees the transition.\n\nOverall: the paper is worth engaging with. The hydro claim needs either the Milne-space analysis included or a separate finite-volume test with varying total length but fixed lattice spacing. The phase-transition content will probably survive scrutiny. Send it to peer review, but the referee should insist on seeing the Milne-space analysis or an equivalent boost-invariance check before the Bjorken-flow claim is accepted.","headline":"Careful tensor-network study of the Schwinger model whose phase-transition results are solid but whose Bjorken-flow claim is under-supported until the deferred Milne-space analysis is made available.","tokens_in":22019,"tokens_out":2513,"would_cite":false,"duration_ms":25213,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Light fermions turn the massive Schwinger model into a fluid and expose its parity-breaking transition","keywords":["massive Schwinger model","tensor networks","TEBD","real-time dynamics","hydrodynamization","boost-invariant flow","parity symmetry breaking","string breaking"],"falsifier":"Evolve $m/g=0.1$ at $\\theta=0$ on a larger lattice (say $N=200$) with the same energy and check at $t=25/g$ whether the bulk viscous pressure still decays to zero and the velocity still follows $v^z\\approx z/t$ before reflections return from the boundaries; if $\\Pi$ stays at its early magnitude or boundary echoes dominate the profile, the boost-invariant flow claim is falsified. For the parity transition, scan $m/g$ in small steps across 0.33 at $\\theta=\\pi$: the claim predicts both electric-field order parameters stay consistent with zero below the threshold and rise sharply above it, so a smooth monotone rise or a jump at a different mass would falsify the dynamical-transition claim.","tokens_in":21042,"feed_emoji":"🌊","tokens_out":10118,"duration_ms":81157,"temperature":0.7,"pith_summary":"This paper follows, with a tensor-network simulation, what happens after a localized blob of energy is placed on top of the vacuum of the massive Schwinger model—a one-space, one-time-dimensional gauge theory that has confinement and chiral dynamics in common with quantum chromodynamics. It reports two findings. When the fermion mass is small compared with the gauge coupling ($m/g = 0.1$), the energy spreads inside a light cone with a flow velocity close to $v^z \\approx z/t$ and the bulk viscous pressure decays toward zero, so the system hydrodynamizes in the manner assumed for the central region of heavy-ion collisions. This fluid-like behavior weakens as $m/g$ increases and is replaced by oscillations. At the topological angle $\\theta = \\pi$, the electric field and charge density expose a dynamical counterpart of the spontaneous parity-breaking transition near $m/g \\approx 0.33$: below the threshold, string breaking screens charges and the central electric field decays; above it, the primary charges propagate as nearly free particles connected by a persistent electric field.","feed_headline":"Light fermions turn a 1D gauge theory into a fluid","feed_subtitle":"At small mass the expansion is fluid-like; the electric field then marks the m/g≈0.33 parity transition.","key_machinery":"The calculations are driven by a gauge-invariant matrix product state that encodes the $U(1)$ Gauss law blockwise in the tensor network, evolved in real time with the time-evolving block decimation (TEBD) algorithm; the electric field is truncated with a cutoff large enough for convergence. The hydrodynamic analysis uses the symmetric, gauge-invariant stress tensor, decomposed in the energy-comoving frame into energy density, flow velocity, and effective pressure $P + \\Pi$, with the equilibrium equation of state $P(\\varepsilon)$ computed separately by a doubled-space (Thermofield Double) purification, so that the bulk viscous pressure $\\Pi$ isolates the off-equilibrium response. The phase-transition analysis uses two electric-field observables, the time-averaged central field and the growth rate of the total field, which serve as dynamical order parameters whose sharp rise near $m/g \\approx 0.33$ at $\\theta=\\pi$ signals the spontaneous breaking of parity.","core_discovery":"On its own terms, the paper claims that a localized energy excitation in the massive Schwinger model evolves like a relativistic fluid whenever the coupling dominates the mass: the stress tensor, decomposed in the energy-comoving frame, gives a velocity profile $v^z \\approx z/t$ and a bulk viscous pressure $\\Pi$ that decreases in time and tends to zero, matching the standard boost-invariant expansion scenario. As $m/g$ grows, $\\Pi$ stays sizable and the velocity becomes oscillatory, marking the breakdown of the fluid description. Independently, at $\\theta = \\pi$ the evolution of the electric field and charge density separates into two regimes across $m/g \\approx 0.33$: in the parity-restored phase the electric string breaks repeatedly, producing secondary pairs that screen the primary charges, so the central field oscillates and decays; in the parity-broken phase the string is stabilized, primary charges separate almost freely, and a persistent electric field connects them. The paper quantifies this with two order parameters—the time-averaged central electric field $\\langle E(0,t)\\rangle$ and the time-averaged growth rate of the total electric field $dE_{\\mathrm{total}}/dt$—both of which stay near zero below the threshold and rise sharply above it, with a power-law fit consistent with a critical transition.","pith_inferences":["Extension: the same two electric-field observables could be measured on a cold-atom quantum simulator with site-resolved charge and electric field, turning the sharp rise near $m/g\\approx 0.33$ into a laboratory witness of the parity transition.","Extension: scanning the topological angle $\\theta$ between 0 and $\\pi$ in this setup could test whether the stiffening of the equation of state at the parity transition shifts the mass threshold for the onset of fluidity.","Extension: the paper quotes a power-law fit as evidence of critical scaling, but extracting a reliable critical exponent would require a finite-size scaling study across several lattice sizes rather than one lattice."],"forward_implications":["The light-mass Schwinger model becomes a first-principles testbed for hydrodynamization: the full approach to equilibrium can be watched through the decay of the bulk viscous pressure.","The threshold in $m/g$ where the bulk pressure stops vanishing and the velocity becomes oscillatory also marks where the fluid description loses validity.","At $\\theta=\\pi$ and $m/g>0.33$, the persistent electric field and nearly free charge separation give a real-time, finite-volume analogue of deconfinement at zero temperature.","The two electric-field order parameters supply a dynamical definition of the parity-breaking transition that does not rely on equilibrium order parameters and could be measured in a simulator."],"supporting_citations":[{"why":"Defines the massive Schwinger model whose real-time dynamics are simulated.","marker":"[3]"},{"why":"Supplies the time-evolving block decimation algorithm used for ground-state preparation and real-time evolution.","marker":"[22]"},{"why":"Establishes the zero-temperature phase diagram and the parity-breaking transition near $m/g\\approx 0.33$ at $\\theta=\\pi$.","marker":"[23]"},{"why":"Provides the doubled-space purification method used to compute the equilibrium equation of state.","marker":"[49]"},{"why":"Supports the background-field and string-tension picture of the parity-broken phase with stable electric strings.","marker":"[50]"},{"why":"Adds independent lattice evidence for the parity-breaking transition at $\\theta=\\pi$.","marker":"[51]"},{"why":"Defines the reference boost-invariant expansion profile used to assess the hydrodynamic behavior.","marker":"[52]"},{"why":"Gives the lattice-continuum mass shift that sets the physical fermion mass in the simulations.","marker":"[45]"}],"fun_headline_variants":["Massive Schwinger model shows fluid flow and a parity transition","1D gauge theory mimics fluid when coupling dominates mass","Tensor network reveals fluid dynamics and string breaking in Schwinger model","Heavy fermions break fluidity; parity flips at m/g≈0.33","From fluid to string-breaking: Schwinger model phase transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The hydrodynamic conclusion rests on treating a 100-site open-boundary lattice with a two-site phase rotation in the middle as a fair stand-in for an infinite, effectively boost-invariant expanding medium over the simulated time window; if finite-size effects or boundary reflections shape the early expansion, the fluid interpretation weakens.","fun_headline_variants_meta":{"raw":{"variants":["Massive Schwinger model shows fluid flow and a parity transition","1D gauge theory mimics fluid when coupling dominates mass","Tensor network reveals fluid dynamics and string breaking in Schwinger model","Heavy fermions break fluidity; parity flips at m/g≈0.33","From fluid to string-breaking: Schwinger model phase transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000749,"raw_usage":{"total_tokens":3344,"prompt_tokens":960,"completion_tokens":2384,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":2295}},"tokens_in":576,"tokens_out":2384,"duration_ms":14854,"temperature":1.0,"reasoning_tokens":2295,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:52:55.812990+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evolve $m/g=0.1$ at $\\theta=0$ on a larger lattice (say $N=200$) with the same energy and check at $t=25/g$ whether the bulk viscous pressure still decays to zero and the velocity still follows $v^z\\approx z/t$ before reflections return from the boundaries; if $\\Pi$ stays at its early magnitude or boundary echoes dominate the profile, the boost-invariant flow claim is falsified. For the parity transition, scan $m/g$ in small steps across 0.33 at $\\theta=\\pi$: the claim predicts both electric-field order parameters stay consistent with zero below the threshold and rise sharply above it, so a smooth monotone rise or a jump at a different mass would falsify the dynamical-transition claim.","supporting_citations":[{"cited_title":"Gauge Invariance and Mass. 2","cited_arxiv_id":null,"evidence_quote":"Defines the massive Schwinger model whose real-time dynamics are simulated."},{"cited_title":"Ma- trix Product Density Operators: Simulation of Finite- Temperature and Dissipative Systems,","cited_arxiv_id":null,"evidence_quote":"Provides the doubled-space purification method used to compute the equilibrium equation of state."},{"cited_title":"The Massive Schwinger Model on the Lattice Studied via a Local Hamiltonian Monte Carlo Method,","cited_arxiv_id":null,"evidence_quote":"Adds independent lattice evidence for the parity-breaking transition at $\\theta=\\pi$."},{"cited_title":"The Massive Schwinger Model on a Lattice: Background Field, Chiral Symmetry and the String Ten- sion,","cited_arxiv_id":null,"evidence_quote":"Supports the background-field and string-tension picture of the parity-broken phase with stable electric strings."},{"cited_title":"Highly Relativistic Nucleus-Nucleus Col- lisions: The Central Rapidity Region,","cited_arxiv_id":null,"evidence_quote":"Defines the reference boost-invariant expansion profile used to assess the hydrodynamic behavior."}],"review_version":1}