{"id":"b2585229-92b5-43f2-86c3-4f688d99fb4a","arxiv_id":"2510.05346","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A quench from a TPQ state in the XX chain produces a double-plateau entanglement profile, exactly derived by CFT on a Klein bottle, Gaussian numerics, and a quasiparticle picture.","lead":"The paper computes the exact real-time entanglement entropy after a quench from a thermal pure quantum (TPQ) state in a free-fermion XX chain, finding a universal double-plateau evolution. A generalist might care because it provides an exactly solvable example of thermalization in an integrable system that proceeds by coherent dephasing rather than chaos.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'exact' claim is not yet established for the spin TPQ state: the derivation and numerics use a single Gaussian fermionic crosscap state, while the Jordan-Wigner image of Eq. (1) is a superposition of two Gaussian states.","rationale":"The reader identified the Gaussian-subspace simplification as the weakest assumption. My independent reading of the manuscript confirms this is the most load-bearing gap. The main text explicitly limits the numerical benchmark to a Gaussian version of the crosscap state, and the supplemental derivation of the BCS representation (Eqs. S106–S133) starts from that Gaussian lattice state, not from the Jordan-Wigner image of the spin crosscap in Eq. (1). Since the actual spin state is a superposition of two Gaussian states, the exact formula Eq. (4) is not demonstrated for the spin TPQ state. A concrete numerical test on small accessible sizes can settle whether the two-component interference is negligible in the scaling regime. The reader's CONDITIONAL verdict is appropriate: if the simplification proves negligible, the paper would merit acceptance; currently the exactness claim is conditional on an unverified assumption. I therefore recommend no change to the reader's verdict.","tokens_in":23118,"tokens_out":26389,"duration_ms":211691,"concrete_test":"Perform exact evolution of the strict spin crosscap state of Eq. (1) for L=16, 24, 32, 48 with l=L/4 and β chosen so that β_CFT = 2π v_F β/L matches the Fig. 1 values (e.g., β_CFT ≈ 0.25). Express Eq. (1) as the two-component fermionic Gaussian superposition identified in Ref. [15], evolve each component with the XX Hamiltonian, and compute S_A(t) from the resulting two-Gaussian state (or use exact diagonalization for L ≤ 32). Compare against the single-Gaussian simulation and Eq. (4). If the maximum deviation over the two plateaus does not decay with L or exceeds the plotted agreement, the exactness claim for the spin TPQ state fails; if it decays as O(1/L) or vanishes, the simplification is benign.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (4) is the exact entanglement entropy after a quench from the spin-1/2 XX-chain TPQ state defined by Eq. (1). However, the paper's own numerical benchmark and the BCS/quasiparticle derivation use the simplified Gaussian crosscap state |C> = ∏(1 + c†_j c†_{j+L/2})/√2 |0>, not the Jordan-Wigner transform of the spin crosscap of Eq. (1). As the paper notes, that Jordan-Wigner transform is a superposition of two distinct fermionic Gaussian states (Ref. [15]). The CFT calculation is performed for a single crosscap state, and no argument is given that the two-component superposition has the same Klein-bottle two-point function or the same reduced density matrix after the quench. Entanglement entropy is not linear in the state, so the entropy of the superposition cannot be inferred from the entropy of one component; interference terms can substantially alter the reduced density matrix. The paper explicitly defers a quantitative treatment of the full superposition, meaning the exactness claim for the actual spin TPQ state lacks support. This is not a disagreement with the calculation for the Gaussian state, but a gap between what is computed and what is claimed. The sign typo in Eq. (7) is secondary; the missing superposition treatment is the load-bearing issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the quench dynamics of the spin-1/2 XX chain starting from a thermal pure quantum (TPQ) state defined by imaginary-time evolution of a crosscap state, Eqs. (1)–(2). It claims an exact formula Eq. (4) for the time-dependent entanglement entropy, exhibiting a double-plateau structure, and supports it by three complementary routes: a CFT calculation of a vertex-operator two-point function on a Klein bottle, an exact numerical solution of the matrix Riccati equation for the fermionic covariance matrix, and a quasiparticle formula Eq. (10). The paper further argues that local observables dephase to the canonical Gibbs ensemble and interprets the non-monotonic entanglement as the macroscopic signature of dephasing of anomalous pairing correlations. The appendices contain a detailed derivation of the CFT two-point function, including a proof of the key operator identity (S12), the Riccati evolution, and the BCS representation of the TPQ state.","tokens_in":23362,"tokens_out":10205,"duration_ms":85959,"significance":"If the result were established for the stated spin TPQ state, it would be a valuable exact example of coherent, non-chaotic equilibration in an integrable lattice model. The paper has real strengths: the CFT computation is detailed and self-contained; the Riccati solution is derived in the appendix; the quasiparticle entropy density s(k) follows from the BCS form rather than being fitted; and the three independent approaches agree in Figs. 1–2. However, the exactness currently demonstrated is for a single Gaussian fermionic crosscap state, not for the spin crosscap state advertised in Eq. (1). Because the authors explicitly defer the full Jordan–Wigner superposition, the central claim is wider than the evidence presented.","major_comments":[{"comment":"The load-bearing gap is stated by the authors themselves: the JW transform of the strict spin crosscap of Eq. (1) is a superposition of two fermionic Gaussian states, while the numerical simulation and the quasiparticle derivation use the single Gaussian state |C> = ∏(1+c†_j c†_{j+L/2})/√2 |0>. Entanglement entropy is not a linear functional of the state, so agreement for one component cannot establish Eq. (4) for the superposition; the two components can interfere in the reduced density matrix. In addition, the JW transformation is a nonlocal unitary, so the fermionic subsystem entropy computed for the Gaussian component is not automatically equal to the spin subsystem entropy of Eq. (1). The statement that the simplification 'does not affect the qualitative features' is not sufficient for the quantitative exactness claim. The authors should either compute the full two-Gaussian superpos","section":"Numerical Benchmark / Eq. (2)"},{"comment":"The CFT derivation computes the vertex-operator two-point function on a Klein bottle for the single crosscap boundary state defined by the constraint (S6). No mapping is given between this CFT crosscap state and the JW image of the lattice spin crosscap of Eq. (1). The numerical benchmark uses the same simplified Gaussian state, so the three-way agreement in Figs. 1–2 validates the Gaussian model, not necessarily the spin TPQ state. For the claimed exactness, the authors need to identify which lattice object the CFT crosscap corresponds to, or treat the two components explicitly.","section":"CFT approach, Eq. (4) and Supplemental S16–S27"}],"minor_comments":[{"comment":"There is a sign inconsistency between the main-text solution and the appendix solution of the Riccati equation. Eq. (7) reads Γ(τ) = (cos(Hτ)Γ_0 + sin(Hτ))(cos(Hτ) − Γ_0 sin(Hτ))^{-1}, while the derived result in Eq. (S91) is (cos(Hτ)Γ_0 − sin(Hτ))(sin(Hτ)Γ_0 + cos(Hτ))^{-1}. These differ; please correct the main-text formula or explain the discrepancy.","section":"Eq. (7) and Eq. (S91)"},{"comment":"The abstract says the result is for a 'free-fermion system', while the introduction and the central formula are framed for the 'spin-1/2 XX chain'. After clarifying the Gaussian-state caveat, the model statement should be made precise so that the exactness claim is unambiguous.","section":"Abstract/Introduction"},{"comment":"The periodic time variable τ_k is defined only verbally as t modulo L/|v(k)|. Please give an explicit definition, including the treatment of velocities v(k)=0 at k=0,π and the branch choices near revivals, to make the quasiparticle formula reproducible.","section":"Eq. (10)"},{"comment":"There is a typo in the paragraph after Eq. (S10): 'anihation' should be 'annihilation'. Please proofread the supplement.","section":"Supplemental Material"}],"recommendation":"major_revision","confidential_remarks":"This is a potentially strong paper, but the advertised exactness for the spin TPQ state is not supported by the calculations, which are performed for a single Gaussian fermionic crosscap state. I would ask the editor to require either a quantitative treatment of the two-Gaussian superposition or an explicit narrowing of the title/abstract to the fermionic Gaussian model. The sign inconsistency in Eq. (7) should also be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things up front. The CFT work is genuine: the two-point function on the Klein bottle is derived with an explicit proof of the key identity (S12), and the theta-function formula (4) agrees well with the Gaussian numerics in the tested cases. But the paper's central exactness claim does not cover the TPQ state it defines. Eq. (1) is the spin crosscap state, while every calculation in the paper uses a single Gaussian fermionic crosscap state. The paper admits the Jordan-Wigner image of Eq. (1) is a superposition of two Gaussian states (Ref. [15]) and defers a quantitative treatment. Entanglement entropy is not linear in the state, so this is not pedantic; the reduced density matrix of the spin state could differ from the single-Gaussian result. As written, the paper solves a well-defined but different problem: the quench from the fermionic Gaussian TPQ state.\n\nWhat is new and good: the finite-temperature beta > 0 treatment and the double-plateau structure. The quasiparticle formula (10) extends Ref. [16], as the authors acknowledge, but the BCS pair-structure derivation in the appendix cleanly justifies the Fermi-Dirac occupations and the dephasing-to-Gibbs mechanism. The three approaches agree with each other, which is the strongest evidence in the paper. The Riccati imaginary-time evolution is standard but used correctly. The explicit proof of the Klein-bottle identity (S12) is careful formal work.\n\nSoft spots. The Gaussian-superposition gap is the main one, and it is substantial rather than cosmetic. The sign discrepancy between Eq. (7) and Eq. (S91) is minor but needs fixing: an exact main-text formula that disagrees with its own derivation invites distrust. Releasing code would help, but the numerics look plausible. The 'third paradigm' language oversells it; the stationary state is a GGE that happens to coincide with the Gibbs ensemble here. The novelty is the entanglement route and the double-plateau signature, not a new equilibration paradigm.\n\nBottom line: this deserves a serious referee. The likely outcome is major revision—either show the two Gaussian components of the spin crosscap yield the same block reduced density matrices, or restate the claims for the fermionic model. The core calculation is worth keeping in the literature. I'd bring it to a reading group for the Jordan-Wigner discussion.","headline":"Solid CFT and numerics for a fermionic Gaussian TPQ quench with a new double-plateau entanglement signature; the exactness claim for the spin-chain TPQ state outruns the calculation.","tokens_in":23924,"tokens_out":12715,"would_cite":true,"duration_ms":97051,"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":"The paper derives an exact theta-function formula for entanglement entropy after a quench from a thermal pure quantum state in the XX chain, predicting a universal double-plateau profile.","keywords":["thermal pure quantum state","entanglement entropy","quench dynamics","XX chain","crosscap state","conformal field theory","quasiparticle picture","anomalous pairing"],"falsifier":"Run the exact Gaussian-state evolution for the full superposition of the spin crosscap state at L=500, β=20, l=80 and compare S_A(t) to Fig. 1; a deviation beyond the symbol size would falsify the quantitative CFT formula. A complementary test is a cold-atom or superconducting-qubit realization of the antipodal-pair state and a direct measurement of the two plateaus.","tokens_in":22922,"feed_emoji":"⚛️","tokens_out":5239,"duration_ms":87188,"temperature":0.7,"pith_summary":"Thermal pure quantum (TPQ) states are deterministic pure states that look locally thermal. The paper asks what happens when such a state in the integrable spin-1/2 XX chain is evolved with the same Hamiltonian that prepared it. The answer is an exact, closed-form entanglement profile with a double plateau, not the usual linear growth and saturation. The derivation runs through a conformal field theory two-point function on a Klein bottle, exact Gaussian-state numerics from the matrix Riccati equation, and a quasiparticle picture of antipodally entangled pairs. Together these show that the stationary state is the canonical Gibbs ensemble, reached by coherent dephasing of anomalous pairing correlations.","feed_headline":"Quench from thermal pure state shows exact double-plateau entropy","feed_subtitle":"Theta functions, Gaussian numerics, and quasiparticles agree on a non-monotonic route to Gibbs equilibrium.","key_machinery":"The machinery is the crosscap state |C>, a product of maximally entangled antipodal pairs, and its imaginary-time evolved version |Ψβ>. The load-bearing identity is the normalized two-point function of vertex twist fields on the Klein bottle, Eq. (S26), which yields the theta-function formula; numerically the matrix Riccati equation dΓ/dτ = -H - ΓHΓ evolves the Gaussian covariance matrix exactly; and conceptually the quasiparticle picture with s(k) as the binary entropy of a Fermi-Dirac mode occupation explains the subtraction in Eq. (10).","core_discovery":"The central claim is that S_A(t,σ), the von Neumann entropy of an interval of length σ at time t after the quench, is exactly given by Eq. (4), a ratio of Jacobi theta functions with modulus iβ/2π, and that this function displays two plateaus: an initial plateau before antipodal pairs enter the interval, a decrease as pairs become fully contained, and a rise to a second plateau as they exit. The same profile is obtained exactly from the covariance-matrix evolution and is quantitatively captured by the quasiparticle formula Eq. (10). The paper also proves that after dephasing the reduced state equals the Gibbs ensemble, because the crosscap initial state's symmetries force all higher conserve","pith_inferences":["The single-Gaussian simplification of the crosscap state is the paper's stated gap; evaluating the full two-component superposition would quantify whether the plateau heights shift for the actual spin chain.","Non-monotonic entanglement could serve as an experimentally accessible witness of pairing-coherence dephasing in quantum simulators, since it requires only measuring S_A(t) rather than anomalous correlators.","The exact structure may persist approximately in interacting integrable models, providing a way to test how integrability governs equilibration beyond free fermions; this remains an extension the paper does not make."],"forward_implications":["Universal double-plateau: the theta-function profile is independent of microscopic details within the free-boson/Dirac-fermion universality class.","Gibbs without GGE: because the conserved charges beyond the Hamiltonian vanish for this initial state, the stationary GGE reduces to the canonical Gibbs ensemble, a rare exact lattice example.","Dephasing mechanism: the anomalous pairing ⟨c_k c_{−k}⟩ rotates at frequency 2E(k) and dephases; observables sensitive to pairing, such as ⟨c_j c_{j+1}⟩, relax exactly to thermal values.","Direct transfer: the same CFT and Gaussian methods extend to mutual information and entanglement negativity, with the quasiparticle picture adapted by minor changes."],"fun_headline_variants":["Thermal quench entropy: exact double-plateau profile","Free-fermion quench yields exact double-plateau entropy","Double-plateau entropy from thermal pure quench","Entropy after thermal quench: exact twin plateaus"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative predictions rely on treating the crosscap state as a single fermionic Gaussian state; if the two Gaussian components of the exact spin-state superposition interfere, the plateau heights and timescales could differ from Eq. (4).","fun_headline_variants_meta":{"raw":{"variants":["Thermal quench entropy: exact double-plateau profile","Free-fermion quench yields exact double-plateau entropy","Double-plateau entropy from thermal pure quench","Entropy after thermal quench: exact twin plateaus"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000712,"raw_usage":{"total_tokens":2972,"prompt_tokens":606,"completion_tokens":2366,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":350,"completion_tokens_details":{"reasoning_tokens":2297}},"tokens_in":350,"tokens_out":2366,"duration_ms":13396,"temperature":1.0,"reasoning_tokens":2297,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T11:19:07.838654+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the exact Gaussian-state evolution for the full superposition of the spin crosscap state at L=500, β=20, l=80 and compare S_A(t) to Fig. 1; a deviation beyond the symbol size would falsify the quantitative CFT formula. A complementary test is a cold-atom or superconducting-qubit realization of the antipodal-pair state and a direct measurement of the two plateaus.","supporting_citations":[],"review_version":1}