{"id":"6ba54549-0db8-4392-bf17-dfdfa08e2eb5","arxiv_id":"2607.26662","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In the 2D quantum Ising ferromagnet, a domain-wall interface relaxes inside a lightcone with a universal linear magnetization ramp near the center, and the order parameter at any fixed point decays as 1/t in the rough phase.","lead":"This paper derives the late-time magnetization profile of a domain-wall interface in the two-dimensional quantum Ising model, finding a universal linear ramp inside the lightcone and a 1/t decay at any fixed point in the rough phase. It also explains why the rough interface phase breaks down as the system moves away from criticality.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (33) rests on the imported form-factor residue (16) and the unverified vanishing of regular terms; no independent check is provided.","rationale":"The paper's central claim is the lightcone profile (33) and the resulting 1/t decay of the order parameter at fixed x. The derivation is coherent and the low-density argument is physically motivated, but the single most load-bearing input is the singular part of the order-parameter form factor, Eqs. (15)–(16), imported from ref. [7] without independent derivation. The reader identified exactly this as the weakest assumption. My assessment agrees: the residue formula is crucial but unverified in this manuscript, and the treatment of regular terms, while plausible, also depends on an oddness argument that is not checked. The paper itself notes that numerical verification in d=2 is difficult and relies on a recent work for the roughening transition, but does not test Eq. (33). Therefore the appropriate verdict remains CONDITIONAL, unchanged from the reader's verdict: the analytic derivation is internally consistent, but the conclusion should be accepted conditionally on an independent verification of the form-factor input or a direct numerical check of the predicted profile. No internal inconsistency or fatal flaw was found that would justify rejection.","tokens_in":10596,"tokens_out":25625,"duration_ms":268895,"concrete_test":"Perform a tensor-network (e.g., TDVP) real-time simulation of the 2D transverse-field Ising model on a cylinder of circumference L_y=8, with a Gaussian domain-wall initial state of width W0=2, at h/h_c≈0.85 (rough phase). Measure ⟨σ^x(x,t)⟩ for x=0,1,2,3 and times t up to ~10. Check that t⟨σ^x(x,t)⟩/⟨σ^x⟩+ approaches a constant A_f x, with A_f matching Eq. (24) computed from the chosen f, and that the profile collapses onto a straight line for |x|/t≪1. Failure of this scaling collapse or a slope deviating from A_f indicates that the imported residue (16) or the cancellation of regular terms is invalid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result, Eq. (33), is obtained from the form-factor expansion (15) with the pole residue c^{(-1)} = −2iML/N ⟨σ⟩+ quoted from ref. [7], and from the assertion that regular terms (k≥0) do not contribute to the |x|≪t profile. This input is not rederived in the present setting. The residue contains the macroscopic factor L/N, which is surprising for a connected form factor of a local operator and is not explained. The k=0 term is O(1/t) at fixed x and is therefore of the same order as the linear term A_f x/t; it is discarded only by an oddness argument that forces c^{(0)}=0, but the necessary cancellation is not independently verified. Neither Eq. (16) nor c^{(0)}=0 is checked against an independent derivation or a numerical simulation. If the residue is incorrect, the step heights outside the lightcone would not match the initial condition, and if c^{(0)}≠0, the 1/t coefficient at fixed x would differ from the predicted A_f x. Thus the dissolution of the interface and the universal linear profile are conditional on an unverified, imported form-factor property.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the unitary time evolution of a domain-wall initial condition in the ferromagnetic phase of the 2D transverse-field Ising model. The initial state interpolates between the two broken-symmetry ground states, with a generic left-right symmetric profile and translation invariance in y. Using the asymptotic quasiparticle basis and a form-factor decomposition of the order parameter, the authors derive the large-time limit of ⟨σ^x(x,t)⟩. The central result is Eq. (33): outside the lightcone the magnetization retains the bulk values ±⟨σ⟩_+, while for |x|≪t it grows linearly as \\tilde A_f x/t, with amplitude set by the initial condition. Consequently, at every fixed x the order parameter decays as 1/t, so the interface dissolves. The same formalism is used to argue for a rough phase at small quasiparticle density and for a breakdown of that phase away from criticality.","tokens_in":10921,"tokens_out":10032,"duration_ms":114967,"significance":"If the derivation is complete, this is a valuable extension of the d=1 quench results to two dimensions: it makes a sharp, falsifiable prediction for the lightcone profile of the 2D Ising ferromagnet and identifies which features are universal. The strength of the paper lies in the clean stationary-phase/contour argument in §2.2, the generic treatment of the initial condition (which survives only through one amplitude), and the probabilistic interpretation in terms of quasiparticle trajectories. However, the central result relies on the form-factor residue (16) imported from ref. [7] and on the vanishing of the regular term c^(0)=0; neither is derived here, and no direct numerical verification of Eq. (33) is supplied. These are not cosmetic issues, because the k=0 term is of the same order as the predicted 1/t decay. The result is therefore conditional on external inputs.","major_comments":[{"comment":"The low-momentum expansion and the pole residue c^(−1)_σ = −2iML/N ⟨σ⟩_+ are imported from ref. [7] and not rederived in the present setting. This is the load-bearing input for the entire late-time profile: through Eqs. (21)–(25), the step heights and the linear central region are direct consequences of this residue. The macroscopic factor L/N in a connected matrix element of a local operator is surprising and is not explained; the limiting order L,N→∞ with fixed ratio is only specified later, in Eq. (27). I request a self-contained derivation of Eqs. (15)–(16) in the 2D context, or an independent check (e.g., a cylinder numerical simulation initialized in the same class, or a direct form-factor computation).","section":"§2.2, Eqs. (15)–(16)"},{"comment":"The discard of the regular part of the form factor reduces to the statement c^(0)_σ=0. This is load-bearing. For fixed x, the k=0 contribution is O(1/t), i.e. of the same order as the predicted A_f x/t term; it is not a subleading correction. The paper's only justification is the footnote that the exact d=1 form factor has c^(0)=0, which is not a statement about the 2D model. The additional appeal to oddness of ⟨σ(x,t)⟩ is not sufficient as written: oddness constrains the full expectation value, not individual contributions to an asymptotic expansion, unless all other O(1/t) even terms are proved absent. Please provide a derivation of c^(0)_σ=0, or a calculation of the full O(1/t) coefficient, for the 2D Ising model.","section":"§2.2, after Eq. (23)"},{"comment":"The treatment of the M2≈1.8M1 species is qualitative. The statement that the off-diagonal (species 1↔2) form factor has no pole is plausible because annihilation requires equal masses, but the conclusion that its regular part does not contribute for |x|≪t again relies on the unproved c^(0)=0 property for mixed-species matrix elements. Since a regular low-momentum limit would produce an O(1/t) contribution at fixed x, the final Eq. (33) is not fully established without this check. A short low-momentum analysis of the mixed form factor, or a numerical estimate of its low-momentum limit, is needed.","section":"§3, second-quasiparticle-species paragraph"},{"comment":"The manuscript explicitly notes that numerical verification in d=2 is difficult and cites only the recent Ref. [25] for the roughening transition. Given that Eqs. (15)–(16) and c^(0)=0 are imported from earlier work, the lack of any direct check of Eq. (33) is a serious gap. The Monte Carlo check in ref. [7] concerns the static interface width, not the time-dependent one-point function. I am not requesting numerics as a substitute for a proof, but the authors should either provide the missing analytic derivations or present a numerical test of the x/t scaling and the 1/t decay.","section":"§3, last paragraphs and Conclusion"}],"minor_comments":[{"comment":"The text says 'This allows us to rewrite (21)' but the displayed equation being rewritten is (17). Please correct the reference.","section":"§2.2, Eq. (19)"},{"comment":"Typos: 'm odel' in the running title, and 'appearence' in the Figure 1 caption.","section":"Title, Fig. 1"},{"comment":"The notation σ^x(x,t) uses x both as the operator label (longitudinal Pauli matrix) and as the spatial coordinate. Consider writing σ^L(x,t) or using bold x for the spatial vector, to avoid confusion.","section":"§3, Eqs. (32)–(33)"},{"comment":"The limits N→∞ and L→∞ are introduced informally. Since Eq. (27) later uses a fixed ratio N/L=κ/ξ, the order of limits should be stated explicitly.","section":"§2.1–§2.2"},{"comment":"The relation N/L=κ/ξ is stated without derivation in the main text; the reference to ref. [7] appears later. An explicit citation at Eq. (27) and a definition of κ would improve readability.","section":"§2.3, Eq. (27)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of cond-mat.stat-mech and the central prediction is interesting. The main editorial concern is the amount of weight placed on a previous result from the same group (ref. [7]) without a self-contained derivation or independent check. I would ask the authors to make that input transparent and either derive it or support it numerically before publication. I do not regard this as a reason for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a genuine step beyond the 1D analysis: it handles the 2D interface with an infinite number of quasiparticles along y, and the connectedness-structure argument plus the stationary-phase evaluation are clean and convincing. The main prediction, Eq. (33), says the order parameter decays as 1/t at any fixed point inside the lightcone, with a universal linear ramp in x/t whose slope is the only memory of the initial condition. That is a new analytical result and is physically significant.\n\nThe soft spots are real but not disqualifying. The singular residue c^{(-1)} = -2iML/N ⟨σ⟩+ is quoted from [7] rather than derived here. The L/N factor looks odd at first, but it is just the inverse density N/L = κ/ξ, so the residue is actually finite and controlled; the paper would be easier to trust if it said that explicitly. The k=0 term, which would also contribute at O(1/t) at fixed x, is set to zero by an oddness argument. That is plausible because the initial state is odd, but it is not independently verified. If c^(0) were nonzero, the 1/t coefficient would shift. The multi-species discussion is qualitative, though the conclusion that cross-species terms cannot produce the pole is sound. Finally, there is no numerical check of Eq. (33). That is not a fatal objection for a paper of this type, but it keeps the result conditional.\n\nOverall, the central mechanism holds up. The paper deserves a serious referee. I would send it to review, asking for a derivation or a more careful justification of the residue and the vanishing of c^(0), and ideally a numerical check of the linear ramp.","headline":"First 2D analytical interface profile, but the load-bearing form-factor residue is imported and the vanishing regular term is only argued by symmetry.","tokens_in":11357,"tokens_out":3477,"would_cite":true,"duration_ms":38723,"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":"At large times in the rough phase of the two-dimensional quantum Ising model, an interface between the two ferromagnetic ground states dissolves: the magnetization at any fixed point decays as 1/t, and the late-time profile inside the light","keywords":["two-dimensional quantum Ising model","interface dynamics","order parameter","form factors","quasiparticles","roughening transition","lightcone","spontaneous symmetry breaking"],"falsifier":"Simulate the real-time evolution of the 2D transverse-field Ising model from a domain-wall initial state and measure ⟨σ^x(x,t)⟩ at fixed x for h well inside the claimed rough phase: the prediction is a 1/t decay with an x/t-linear profile at intermediate x. Alternatively, compute the pole residue c^{(-1)}_σ directly from the lattice form factor; if it differs from −2iML/N⟨σ⟩₊, or if a regular term contributes at order x/t, Eq. (33) fails.","tokens_in":10543,"feed_emoji":"🧲","tokens_out":5784,"duration_ms":61198,"temperature":0.7,"pith_summary":"This paper analyzes how a sharp interface between the two ferromagnetic ground states of the two-dimensional transverse-field Ising model evolves in time. It claims that at late times, inside the lightcone, the magnetization profile becomes a universal linear ramp in x/t whose only memory of the initial condition is one amplitude, while outside the lightcone the two bulk values remain. At any fixed spatial point the order parameter therefore decays to zero as 1/t, meaning the interface dissolves and its width grows without bound—the quantum analogue of interfacial roughening. The same quasiparticle argument predicts where this rough phase breaks down as the system moves away from criticality.","feed_headline":"Magnetization decays as 1/t as 2D quantum Ising interface dissolves","feed_subtitle":"Inside the lightcone the late-time profile is universal up to one amplitude; the interface roughens and its width grows without bound.","key_machinery":"The central mechanism is the low-momentum singular part of the one-particle form factor of the order parameter: F^{σ,c}_1(p|q) ≈ c^{(-1)}_σ /(p_x − q_x) with c^{(-1)}_σ = −2iML/N ⟨σ⟩₊ for p_y = q_y and p, q → 0. This kinematical pole, combined with stationary-phase dominance of small momenta, turns the many-quasiparticle integral into a single residue integral over rapidities; the result is the step function outside the lightcone and the linear ramp inside. Regular terms in the form-factor expansion are shown not to contribute in |x| ≪ t.","core_discovery":"The paper claims that for the 2D transverse-field Ising ferromagnet, starting from any left-right symmetric, y-translation-invariant initial state interpolating between the two ferromagnetic ground states, the large-time longitudinal magnetization in the rough phase is given by Eq. (33): −⟨σ_x⟩₊ for x < −t, A_f ⟨σ_x⟩₊ x/t for |x| ≪ t, and ⟨σ_x⟩₊ for x > t. Consequently, lim_{t→∞} ⟨σ_x(x,t)⟩ = 0 at every fixed x. The profile is universal in the inner region except for the amplitude A_f, which encodes the initial state; the second quasiparticle species only renormalizes this amplitude. The same analysis yields a mechanism for breakdown of the rough phase away from criticality.","pith_inferences":["Because the argument uses only the kinematical pole of the one-particle form factor, the same linear-ramp profile should appear for other local operators whose form factors carry the same pole, with only the amplitude changing—a prediction the paper does not spell out.","The predicted 1/t decay at fixed x is a concrete signature for quantum simulators with two-dimensional spin arrays: a domain-wall initial state should show local magnetization decaying to zero with a 1/t tail above the roughening field, and a different behavior below it.","The factorization of the profile into x/t times a state-dependent amplitude suggests a scaling form ⟨σ(x,t)⟩ ≈ t^{-1} F(x/t) in the rough phase; checking this data collapse in numerics would test universality more sharply than a single time slice.","If the smooth phase exists, the transition should be visible as a change in the decay exponent or the appearance of persistent oscillations in ⟨σ(x,t)⟩, since the pole-driven argument would no longer apply."],"forward_implications":["At any fixed spatial point, the magnetization tends to zero as 1/t, so the sharp interface dissolves and the interfacial width grows linearly with time.","For |x| much smaller than t, the magnetization profile is a universal linear ramp, ⟨σ^x⟩/⟨σ^x⟩₊ ≈ A_f x/t, and the details of the initial interpolation survive only in the single amplitude A_f.","Outside the lightcone x = ±t the magnetization keeps its initial bulk values, so information about the initial state propagates at the quasiparticle speed of light.","A second stable quasiparticle species (mass about 1.8 times the lightest) does not alter the profile except through a redefinition of A_f, because the pole term is absent in off-diagonal channels.","The same quasiparticle-density argument identifies when the rough phase breaks down: as the correlation length shrinks away from criticality, the interface becomes smooth below some h_r > 0."],"fun_headline_variants":["Interface meltdown in 2D quantum Ising model","2D Ising interface fades, magnetization dies as 1/t","Quantum Ising interface dissolves with universal decay","Universal 1/t decay in 2D Ising interface dynamics","Rough phase collapse: 2D Ising interface dissolves"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the magnetization's two-quasiparticle matrix element has exactly the quoted 1/(p_x − q_x) pole with the quoted residue in the 2D Ising ferromagnet; the paper takes this from earlier work and relies on all regular terms vanishing in the inner region. If that pole strength is wrong or the regular terms contribute, the predicted profile collapses.","fun_headline_variants_meta":{"raw":{"variants":["Interface meltdown in 2D quantum Ising model","2D Ising interface fades, magnetization dies as 1/t","Quantum Ising interface dissolves with universal decay","Universal 1/t decay in 2D Ising interface dynamics","Rough phase collapse: 2D Ising interface dissolves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000105,"raw_usage":{"total_tokens":829,"prompt_tokens":657,"completion_tokens":172,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":401,"completion_tokens_details":{"reasoning_tokens":88}},"tokens_in":401,"tokens_out":172,"duration_ms":2683,"temperature":1.0,"reasoning_tokens":88,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T11:17:51.444547+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the real-time evolution of the 2D transverse-field Ising model from a domain-wall initial state and measure ⟨σ^x(x,t)⟩ at fixed x for h well inside the claimed rough phase: the prediction is a 1/t decay with an x/t-linear profile at intermediate x. Alternatively, compute the pole residue c^{(-1)}_σ directly from the lattice form factor; if it differs from −2iML/N⟨σ⟩₊, or if a regular term contributes at order x/t, Eq. (33) fails.","supporting_citations":[],"review_version":1}