{"id":"7b6854da-d9bc-4c2d-951d-ef274456dbf8","arxiv_id":"2512.13809","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Born-averaged measurement-induced entanglement in Tomonaga-Luttinger liquids is equivalent to averaging over conformal boundary conditions, yielding universal cumulants and a bimodal heavy-tailed distribution.","lead":"Using conformal field theory and a replica trick, this paper derives formulas for the full statistics of measurement-induced entanglement in one-dimensional quantum critical systems, going beyond the previously known average. The authors find universal scaling laws for all cumulants and show the distribution is bimodal with heavy tails, verified against exact numerics on the XX chain.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (33) reduces every measurement outcome to a constant Dirichlet shift; non-conformal outcomes are asserted negligible, but heavy tails make them capable of controlling the ζ→0 cumulant scaling.","rationale":"The paper's strongest claim — the equivalence in Eq. (39) and the universal scaling Eq. (40) — is exactly as the reader describes. The most load-bearing condition is that the microscopic Born average over all lattice measurement outcomes collapses, at low energy, to an average over conformal boundary conditions. The paper itself flags this in Sec. 3.1: the ferromagnetic outcome is known not to flow to a conformal boundary condition, and the authors argue only that such contributions should be small. My stress test sharpens this: the derivation of the winding function in Sec. 3.2 effectively uses only the zero mode of the outcome profile, so the non-conformal question is not merely a boundary RG subtlety but an assumption built into the form of Eq. (25). Moreover, the predicted distribution is heavy-tailed, so the 'bad' outcomes cannot be dismissed on rarity grounds without a quantitative estimate of their Born weight and entanglement. The replica-limit analytic continuation is a real but secondary concern; it is standard in this literature and is not where the argument is most fragile. The proposed concrete test avoids the replica limit and the inaccessible ζ→0 window: for integer n,k1,k2, Eq. (33) is a concrete prediction for a lattice quantity that can be computed exactly or with Monte Carlo on the free-fermion chain. If that check passes, the conditional status is substantially strengthened; if it fails, the central equivalence is broken. Given the reader already assigned CONDITIONAL for essentially this reason, my read does not change the verdict.","tokens_in":28199,"tokens_out":10145,"duration_ms":93552,"concrete_test":"For the XX chain, use the free-fermion correlation-matrix update rules of Appendix C to compute the lattice generalized replica partition function Z_A(k1,k2) = E_m[(tr ρ_{m,A}^n)^{k1}(tr ρ_m)^{k2}] for integer values (n,k1,k2) = (2,2,1) and (3,1,1). Evaluate this by exact enumeration of all outcomes for L = 24–48 (with |B| ≤ ~20) and by Monte Carlo sampling for L = 64–96, comparing with the RHS of Eq. (33) over ζ ∈ [0.01, 0.5]. A deviation beyond a few percent (after finite-size extrapolation) would falsify the Born-average-to-conformal-BC replacement before any replica limit is taken; agreement would independently support Eq. (39), including the role of the winding contribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central equivalence, Eq. (39), depends on replacing the Born average over microscopic outcomes by an integral over conformal boundary conditions. The critical step is in Sec. 3.2: after the rotation, the outcome-dependent factor is dropped using ∑_m e^{-S[...]} ≡ 1, leaving only the winding function W beyond the survival of an integral over the zero-mode shift δφ. However, the classical action used in Eq. (25), S_{C(1)}[φ_cl,m] = (g/2)(δm)^2/h, is the action for a constant boundary-value difference; it contains no dependence on spatial variations m(θ). Thus the derivation implicitly assumes every measurement outcome is a constant Dirichlet condition, not a general profile. The paper acknowledges in Sec. 3.1 that outcomes such as the ferromagnetic state do not flow to conformal boundary conditions, but it only asserts that their contributions are 'small in practice' — no bound or RG estimate is given. This is load-bearing because the predicted distribution P(S_m) is heavy-tailed; rare outcomes with S_m near log 2 carry a ζ^{g/2}/√log(1/ζ) tail in the ζ→0 limit, the same order as the leading cumulant scaling in Eq. (40). If a small fraction of non-conformal outcomes produces near-maximal entanglement with non-negligible Born weight, κ_l for l≥2 could acquire leading-order corrections. The numerical evidence — good agreement for κ_2, κ_3 and P(S_m) at ζ=0.02, L≤600 — does not probe the ζ→0 regime; the authors state that the scaling window of Eq. (40) is beyond accessible numerics. Therefore the central universal-statistics claim is conditional on an unquantified and untested assumption about the lattice–continuum outcome map.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the statistics of measurement-induced entanglement (MIE) after projective charge measurements on disjoint intervals of a one-dimensional Tomonaga-Luttinger liquid. Using a generalized replica partition function Z_A(k1,k2) = <(tr ρ_{m,A}^n)^{k1} (tr ρ_m)^{k2}>, the authors map the problem to a compact free boson on a cylinder, perform a rotation that isolates the measurement dependence in a single replica, and sum over winding sectors. The central result, Eq. (39), states that the cumulant-generating function of the Born-averaged MIE equals that of the 'forced' MIE averaged over conformal boundary conditions with weight p(δφ) ∝ Z_{C(1),δφ}. From this they derive the ζ→0 scaling Eq. (40), with regimes separated at nl=1/2, and the full distribution P(S_m), which is bimodal with a log-normal tail near S_m=0 and a 1/sqrt(log 2 - S_m) tail near S_m=log 2. They also define disorder-induced entanglement (DIE) and find 1/log(1/ζ) scaling. Numerical checks on the XX chain for κ2, κ3, P(S_m), and DIE show good agreement at moderate sizes and at ζ=0.02.","tokens_in":28673,"tokens_out":14067,"duration_ms":125825,"significance":"If the central equivalence Eq. (39) is correct, this is a substantial advance: it provides parameter-free, conformal-field-theoretic predictions for the full outcome-averaged entanglement statistics in a generic critical 1D system, including heavy tails and a rare Bell-pair-generating outcome. The absence of fitted parameters and the quantitative agreement of several independent numerical observables (Figs. 3a,b, 4, 5) are genuine strengths. The universality claim, however, rests on two non-rigorous steps: the assertion that non-conformal microscopic measurement outcomes are negligible in the Born average, and the analytic continuation of the winding function before the replica limit. The paper explicitly acknowledges the former (Sec. 3.1) but provides no bound, and the numerical verification of the asymptotic scaling Eq. (40) is admittedly beyond the accessible regime. These issues are load-bearing but, in my view, addressable within the manuscript's scope.","major_comments":[{"comment":"The derivation reduces every lattice measurement outcome to a constant Dirichlet shift. The classical action used in Eq. (25), S_C(1)[φ_cl,m] = (g/2)(δm)^2/h, depends only on a constant boundary-value difference, not on the spatial profile m(θ). The sum over m in Eq. (24) is therefore effectively a sum over constant boundary conditions only. The paper acknowledges in Sec. 3.1 that outcomes such as |↑↑↑…⟩ do not flow to conformal boundary conditions, but only states that their contributions are 'small in practice'; no RG estimate, operator-content argument, or numerical bound is supplied. This is load-bearing because the predicted P(S_m) and the ζ→0 scaling Eq. (40) are controlled by rare near-maximal outcomes; if a small Born-weight fraction of non-conformal outcomes yields S_m near log 2, the leading-order cumulant scaling could be modified. The numerical checks at ζ=0.02 and L≤600 do n","section":"Secs. 3.1–3.2, Eq. (25)"},{"comment":"The analytic continuation of the winding function W from Eq. (31) to Eq. (32) is not justified. The derivation in Appendix B uses Poisson resummation, completing the square, and a Dirac-delta Fourier representation for positive integers k1 and k2; the final expression involves fractional powers of sums (the k1 and k2 powers of theta functions), and the replica limit k→0 is then taken after this continuation. This is a standard replica-trick heuristic, but it is not a derivation unless a uniqueness or continuity argument is supplied. Because the replica limit is the entire content of the paper, this step needs either a clearer justification (e.g., via the theta-function representation being entire in the exponent) or an explicit check: expand both sides of Eq. (31)/(32) to order k^2 along k1=k, k2=-nk and verify against the stated Eq. (37). At minimum, the paper should flag this as a work","section":"Sec. 3.2 and Appendix B"},{"comment":"The statement that 'taking the replica limit in (17) with the derived generalized replica partition function (30) reproduces precisely (37)' is not demonstrated. Eq. (37) is the first nontrivial test of the central equivalence and the basis for the general cumulant formula Eq. (39), so the intermediate algebra should be shown (or sketched in an appendix). This is not purely cosmetic: the prefactors and normalization in Eqs. (30) and (32) can contribute to derivatives along k1=k, k2=-nk, and the reader cannot verify that the κ2 formula is free of such contributions.","section":"Sec. 4.1, Eq. (37)"},{"comment":"The derivation of Eq. (40) computes the scaling of moments and then asserts the same scaling for cumulants. For a heavy-tailed distribution this is plausible, but it should be justified: cumulants are linear combinations of products of lower moments, and one must check that no cancellation removes the leading term. A small-ζ numerical check of the moment/cumulant ratio, even at the smallest reachable ζ, would strengthen the claim. The paper's own statement (Sec. 4.2) that the scaling window is 'beyond the numerical window accessible in this work' is an important caveat that should be stated more prominently; Fig. 3(c) is a theoretical extension, not a numerical verification.","section":"Sec. 4.2 and Appendix A"}],"minor_comments":[{"comment":"Typo: 'indciates' should be 'indicates'. Also 'à priori' should be italicized and without a space in French ('a priori').","section":"Sec. 1"},{"comment":"The symbol δm is not defined before use. Please define it as the difference of the (constant) boundary values on the two cylinder boundaries, and explicitly note that the spatial-profile dependence has been dropped at this stage; this will make the approximation in the main text transparent.","section":"Eq. (25)"},{"comment":"The prefactor sqrt((nk1+k2+1)g/(2πh)) in Eq. (32) is not derived in the main text. The naive completion-of-square calculation in Appendix B gives different powers of sqrt(2πg/h) and sqrt(2πgn/h); please clarify how these are absorbed and whether they affect the replica limit.","section":"Eq. (32)"},{"comment":"There appears to be a typo in the expression involving 'ε g π logζ' around Eq. (59); the sign and placement of π should be checked against Eq. (58).","section":"Sec. 4.4.2, Eq. (59)"},{"comment":"The caption says 'we see the scaling Eq. (40)' for the theoretical curves in panel (c). Rephrase to make clear that this is a prediction extended into the asymptotic regime, not a numerical observation.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is clearly within scope and the results are potentially important. My main concern is not about the internal consistency of the CFT calculation but about the unproven suppression of non-conformal measurement outcomes, which is acknowledged but not quantified. Given that the predicted distribution is heavy-tailed, this is not a minor caveat. I recommend asking for a concrete diagnostic or bound, and for the analytic-continuation step to be either justified or explicitly flagged as a replica-trick assumption. I do not recommend rejection: the numerical agreement at moderate ζ is good, there are no fitted parameters, and the requested items are likely fixable within the manuscript's scope. The heavy reliance on the companion paper [69] is acceptable but should be made more self-contained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is the full statistics, not just the mean: closed forms for all cumulants of MIE and the full distribution, plus universal ζ→0 scaling forms (ζ^{g/2}/√log with a transition at nl = 1/2) and a distinct prediction for disorder-averaged entanglement (DIE). The two-parameter replica trick and the representation of the Born average as a weighted average over conformal boundary conditions are real advances. There are no fitted parameters, and the numerical checks on the XX chain — second and third cumulants and P(S_m) at ζ = 0.02, L up to 600 — agree well with the theory. This is a serious piece of work, not a routine extension of the authors' earlier mean-MIE paper.\n\nThe main soft spot is exactly the one the authors acknowledge in Sec. 3.1 and then partially wave away. The derivation of the classical winding contribution, Eq. (25), uses an action for a constant boundary-value difference; spatial variations of the measurement outcome m(θ) never enter. So the sum over all microscopic outcomes is effectively replaced by a sum over constant Dirichlet shifts. The paper asserts that non-conformal outcomes such as the ferromagnetic state have small weight, but no bound or RG argument is given. That matters more here than in the mean calculation, because the predicted distribution is heavy-tailed: rare outcomes with S_m near log 2 carry a ζ^{g/2}/√log weight, the same order as the leading cumulant scaling in Eq. (40). If a small fraction of non-conformal outcomes produces near-maximal entanglement with non-negligible Born weight, the higher cumulants could receive leading-order corrections. The numerical agreement at ζ = 0.02 is reassuring but does not probe the ζ→0 regime, which the authors say is beyond accessible system sizes.\n\nThe replica-limit analytic continuation in Appendix B follows a known method and is plausible, though not fully rigorous. I do not see that as a deal-breaker; the physics is clear and the structure of the result is sensible. The scaling forms are checked analytically within the paper but not against numerics in the scaling window — a limitation the authors disclose.\n\nBottom line: this deserves a serious referee. The central claim is conditional on a stated but unquantified assumption, and the paper would be stronger with any quantitative handle on the non-conformal contributions — an RG estimate, a bound, or a cleverer numerical probe. But the derivation is coherent, the results are concrete and falsifiable, and there is no circularity. Send it to peer review, ask the referees to press on the non-conformal outcomes, and expect the authors to have something to say. I would cite this in my own work.","headline":"First analytic result for the full distribution of measurement-induced entanglement in Tomonaga-Luttinger liquids, with a real but unquantified assumption about non-conformal outcomes that should be pressed in review.","tokens_in":29071,"tokens_out":3231,"would_cite":true,"duration_ms":33154,"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 full statistics of measurement-induced entanglement in one-dimensional Tomonaga-Luttinger liquids are universal: Born-averaging over measurement outcomes is equivalent, at low energy, to averaging over conformal boundary conditions weig","keywords":["measurement-induced entanglement","Tomonaga-Luttinger liquid","conformal field theory","Born average","boundary conditions","replica trick","entanglement statistics","Rényi entropy"],"falsifier":"Compute in an interacting gapless chain (e.g., XXZ at Δ≠0) the second and third cumulants of MIE at small cross-ratio and check the predicted scaling ζ^{g/2}/√log(1/ζ) and the crossover at nl=1/2; alternatively, directly compute the contribution of the ferromagnetic outcome to the Born average at small ζ and test whether it vanishes compared with conformal outcomes.","tokens_in":1578,"feed_emoji":"⚛️","tokens_out":3033,"duration_ms":62298,"temperature":0.7,"pith_summary":"What this paper tries to establish: after projective measurements on part of a one-dimensional quantum critical system, the Born-averaged post-measurement entanglement—and all its higher cumulants—are universal and conformally invariant. The key mechanism is a replica trick that turns the average over exponentially many measurement outcomes into an average over conformal boundary conditions of the compact free boson, with each boundary condition weighted by its partition function. As a result, all cumulants share a universal scaling in the maximal-separation limit, with a crossover at nl = 1/2. The full distribution is bimodal with heavy tails. The predictions are confirmed numerically on the XX chain.","feed_headline":"Born-averaged entanglement statistics reduce to one conformal rule","feed_subtitle":"In Tomonaga-Luttinger liquids, outcome weights become partition functions, making every cumulant universal.","key_machinery":"The central object is the generalized replica partition function Z_A(k1,k2) and the conformal map from the n-sheeted cylinder to a finite cylinder with length h(ζ). A reflection rotation isolates the measurement dependence in a single replica, producing a 'winding function' that sums over compact-boson windings. Poisson resummation and analytic continuation convert this into the integral representation that becomes the boundary-condition average.","core_discovery":"The central claim is that the generalized replica partition function Z_A(k1,k2) = (tr ρ^n_{m,A})^{k1} tr ρ_m^{k2} equals, at low energies, an integral over conformal boundary conditions δφ: Z_A ~ ∫ dδφ Z_{C(1),δφ} (Z_{C(n),δφ})^{k1} (Z_{C(1),δφ})^{k2}. Consequently, the l-th cumulant of the post-measurement entanglement equals the l-th cumulant of the forced MIE under p(δφ) ∝ Z_{C(1),δφ}. In the maximal-separation limit ζ→0, cumulants scale as ζ^{g/2}/√log(1/ζ) for n > 1/(2l), as ζ^{g/2} at n = 1/(2l), and as ζ^{2gnl(1−nl)} for n < 1/(2l). The full distribution is bimodal with heavy tails.","pith_inferences":["One could test the recipe in other critical models (e.g., Ising, Potts) where measurements may or may not flow to conformal boundary conditions; the present derivation provides a natural template.","The log-normal tail suggests that experimental single-shot measurements will often find entanglement far from the mean; sample-to-sample fluctuations may be large.","The crossover at nl = 1/2 could be measured by tuning the Rényi index in cold-atom or ion-trap implementations of monitored critical chains."],"forward_implications":["The full statistics of MIE in TLLs are universal: the mean, variance, skewness, and all higher cumulants are fixed by the Luttinger parameter and the cross-ratio.","For large enough Rényi index, every cumulant decays as ζ^{g/2}/√log(1/ζ) when the unmeasured regions are far apart, and the crossover at nl=1/2 is a sharp signature of Born-averaged physics.","The tail near S_m = log 2 implies a finite (if small) probability of generating a Bell pair across the measured region, so a critical state can act as a quantum wire.","Uniform (quenched-disorder) averaging produces different scaling, 1/log(1/ζ), independent of the Rényi index and Luttinger parameter—so Born weighting matters."],"fun_headline_variants":["Conformal weights govern all measurement-induced entanglement moments","Measurement-induced entanglement: universal stats via CFT","Born averaging maps measurement outcomes to conformal boundary conditions","Luttinger liquids: measurement entanglement becomes universal and bimodal","Universal cumulants of measurement-induced entanglement in Luttinger liquids"],"cache_read_input_tokens":30336,"weakest_assumption_plain":"Every microscopic measurement outcome can be replaced, at low energies, by a conformally invariant Dirichlet boundary condition, with non-conformal outcomes such as the ferromagnetic state contributing negligibly to the Born average.","fun_headline_variants_meta":{"raw":{"variants":["Conformal weights govern all measurement-induced entanglement moments","Measurement-induced entanglement: universal stats via CFT","Born averaging maps measurement outcomes to conformal boundary conditions","Luttinger liquids: measurement entanglement becomes universal and bimodal","Universal cumulants of measurement-induced entanglement in Luttinger liquids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000279,"raw_usage":{"total_tokens":1496,"prompt_tokens":747,"completion_tokens":749,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":491,"completion_tokens_details":{"reasoning_tokens":685}},"tokens_in":491,"tokens_out":749,"duration_ms":7141,"temperature":1.0,"reasoning_tokens":685,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T16:20:34.778489+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute in an interacting gapless chain (e.g., XXZ at Δ≠0) the second and third cumulants of MIE at small cross-ratio and check the predicted scaling ζ^{g/2}/√log(1/ζ) and the crossover at nl=1/2; alternatively, directly compute the contribution of the ferromagnetic outcome to the Born average at small ζ and test whether it vanishes compared with conformal outcomes.","supporting_citations":[],"review_version":1}