{"id":"a76aab88-1b00-4b35-be25-5d08a0324596","arxiv_id":"2506.16077","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A simulation-based analysis estimates that quantum entanglement in LHC Z-pair production can be detected at 2-3.75 sigma with current and future luminosities.","lead":"This paper uses quantum state tomography on simulated LHC events to argue that the spin states of Z boson pairs produced in proton collisions are entangled. It projects that this entanglement could be seen at about 2 sigma with existing LHC data and up to 3.75 sigma with the High-Luminosity upgrade.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 2σ/3.75σ sensitivities are not supported: the significance uses the biased mean of the quadratic concurrence estimator rather than a null-hypothesis separation, and the pseudo-experiment statistics are inconsistent with the stated event counts.","rationale":"The reader correctly identified the statistical treatment as the weakest point. My analysis sharpens this: even if the linear coefficient estimators were unbiased, the quadratic concurrence estimator in Eq. 9 is not, and the paper's significance definition (mean over width) does not test against the separable null. This is the single most load-bearing concern because the abstract's quantitative claim is exactly those significances. The pseudo-experiment count inconsistency compounds the problem by leaving the quoted uncertainties unverifiable. I do not object to the physics idea; the qutrit tomography formalism is standard, and the leading-order study is a reasonable starting point. But in its current form the sensitivity projection is not credible. A corrected statistical treatment with an explicit null hypothesis, a closure test, and correctly sized pseudo-experiments could make the claim testable. Therefore the reader's REJECT verdict stands unchanged.","tokens_in":12454,"tokens_out":14895,"duration_ms":164456,"concrete_test":"Run a dedicated closure and null test: (i) from the same MadGraph sample, generate 1000 pseudo-experiments at L = 300 fb⁻¹ and L = 3000 fb⁻¹ with the correct expected event counts (Tab. I), resampling with replacement from a parent sample large enough (or generating 12.7M/126.9M events); (ii) build a separable reference state with the same single-particle polarizations but with c²_MB ≤ 0 (e.g., by taking a product of single-particle density matrices) and pass it through the identical reconstruction; (iii) compare the signal mean, the null mean, and the null 95% quantile. If the null mean is positive and the separation (signal mean - null mean)/std falls below 2σ/3.75σ, or if the closure test on the signal shows a bias exceeding the quoted uncertainty, the abstract's claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim (abstract; Figs. 5–6) depends entirely on the statistical interpretation of the reconstructed lower bound c²_MB from Eq. 9. Two related defects invalidate the quoted significances. First, the 'significance' plotted in Fig. 6 is the mean of the pseudo-experiment distribution divided by its standard deviation, with no reference to the null hypothesis of a separable state. Because c²_MB is a quadratic function of the estimated coefficients A_i, B_i, C_ij, its expectation value is shifted from the true value by a variance-dependent bias (Jensen-type). This is visible in the paper's own numbers: the low-luminosity mean in Fig. 6 is high while the converged truth-level value is about 0.375. A separable state (true c²_MB ≤ 0) will also yield a positive reconstructed mean at low statistics, so the reported 2σ at L ≈ 300 fb⁻¹ measures bias, not entanglement, and does not establish a 2σ separation from the null. Second, the pseudo-experiment description is internally inconsistent: 1000 pseudo-experiments with the event counts of Tab. I would require 5.8M (137 fb⁻¹), 12.7M (300 fb⁻¹), and 126.9M (3000 fb⁻¹) events, not the stated one million total. Unless the resampling procedure is specified, the quoted standard deviations and significances are unvalidated. A closure test and a null test are required before the sensitivity claim can be accepted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies quantum entanglement in pp -> ZZ -> 4l at 13 TeV. The authors parametrize the two-Z spin density matrix with Gell-Mann matrices, reconstruct the polarization and correlation coefficients from simulated lepton angular distributions via quantum state tomography, and compute the lower bound of concurrence c^2_MB (Eq. 9). Using MadGraph5 aMC@NLO at LO and MadSpin, they generate simulated events, split them into 1000 pseudo-experiments, and report that the reconstructed c^2_MB converges to about 0.375 with increasing luminosity. They claim a 2 sigma sensitivity with LHC Run 2+3 luminosity (300/fb) and 3.75 sigma with HL-LHC luminosity (3/ab).","tokens_in":12731,"tokens_out":10145,"duration_ms":116781,"significance":"If the numerical sensitivity claim were reliable, the paper would provide a valuable projection of an entanglement measurement in massive gauge-boson pair production at the LHC, complementing the existing top-quark entanglement measurements. The theoretical framework is standard and the use of quantum state tomography for spin-1 systems is well motivated. However, the central statistical argument has serious flaws: the quoted significance is not a test against a separable null hypothesis, the pseudo-experiment construction is internally inconsistent with the stated event counts, and the background and detector effects are not quantitatively treated. These issues directly affect the abstract's main claim, so the significance of the result cannot be assessed from the present manuscript.","major_comments":[{"comment":"The quoted significance is computed as the mean of the pseudo-experiment distribution of c^2_MB divided by its standard deviation, with no reference to the null hypothesis of a separable state. Because c^2_MB is a quadratic function of the reconstructed coefficients A_i, B_i, C_ij (Eq. 9), its expectation value is positively biased at finite statistics by a Jensen-type variance contribution. The paper's own Fig. 4 shows low-luminosity reconstructed values larger than one while the converged high-luminosity value is about 0.375, which is clear evidence of this bias. A separable state with zero true correlations will also give a positive reconstructed mean at low statistics, so the reported 2 sigma at 300/fb may be measuring bias rather than entanglement. A closure test with an injected separable density matrix and a null-hypothesis significance test are required before the sensitivity claim can be accepted.","section":"Section III, Eq. (9) and Figs. 5-6"},{"comment":"The pseudo-experiment description is internally inconsistent. The paper states that one million events are divided into 1000 pseudo-experiments, each containing events matching the total reachable events at a given luminosity. However, Table I gives 5793 events at 137/fb, 12687 events at 300/fb, and 126870 events at 3000/fb. One thousand pseudo-experiments would require 5.8 million, 12.7 million, and 126.9 million events, respectively, not one million total. Unless a resampling procedure (e.g., bootstrap with replacement) is explicitly described, the standard deviations and significances in Figs. 5 and 6 are not reproducible or validated.","section":"Section III, Table I and the text describing one million events"},{"comment":"The background rejection argument is not quantitative. The text asserts that background events do not matter in the signal region based on the invariant-mass distribution, but no selection criteria, signal-region definition, or residual background fraction is given. Moreover, only triboson backgrounds (ZZZ, WWZ, WZZ) are considered; Z+jets and ttbar backgrounds, which can also produce four leptons, are not addressed. Since the sensitivity claim assumes a background-free sample, the impact of background contamination on the quoted significance is not established.","section":"Section III, Fig. 2 and the background discussion"},{"comment":"All results are obtained from parton-level leptons without any detector simulation. There is no treatment of lepton reconstruction efficiency, acceptance, momentum resolution, or isolation requirements, which are essential for estimating the event yield and angular reconstruction quality at the LHC. For a paper whose central claim is an LHC sensitivity projection, this is a significant limitation; either a detector-level study, or at minimum a smearing/acceptance model, is needed to support the abstract's quantitative claims about Run 2+3 and HL-LHC data.","section":"Section III (all numerical results)"}],"minor_comments":[{"comment":"The expression 2Tr[ρ]^2 in Eq. (8) appears to be a typo; consistency with Eq. (9) requires 2Tr[ρ^2]. For a normalized state, 2(Tr ρ)^2 = 2, which would not lead to Eq. (9). The authors should correct the formula and re-derive the numerical expressions.","section":"Eq. (8)"},{"comment":"The axis label in Fig. 4 is truncated (\"|cos( )|\") and the text states that values of C^2 larger than one are \"unphysical.\" For a two-qutrit pure state, C^2 can reach 4/3, so values moderately above one are not necessarily unphysical; this statement should be revised.","section":"Fig. 4 and surrounding text"},{"comment":"The sentence \"The functions given in Eq. A2 and together with the matrix Aj i A3 in Appendix can be used to extract...\" is garbled. The notation A_j^i in Eq. (13) is not defined in the main text, and the sentence should be rewritten for clarity.","section":"Section II C, paragraph after Eq. (12)"},{"comment":"The paper should specify exactly how the one million generated events are converted into pseudo-experiments: whether events are resampled with replacement, whether the same events are reused for different luminosity points, and how the quoted standard deviations are computed. This is needed for reproducibility.","section":"Section III, pseudo-experiment description"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely and interesting topic, and the underlying physical framework (tomography of ZZ spin states) is sound. However, the statistical analysis as presented does not support the abstract's sensitivity claims. I recommend major revision rather than rejection because the flaws are, in principle, addressable: a null-hypothesis test, a closure test, an explicit and consistent pseudo-experiment procedure, and a detector-level or smeared treatment could materially improve the paper. If the authors cannot provide these, the quantitative claims should be removed and the paper reframed as a purely qualitative proof-of-concept."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is a Monte Carlo sensitivity study for detecting quantum entanglement in ZZ pairs at the LHC via qutrit quantum state tomography. The headline claim—2σ at 300/fb and 3.75σ at 3/ab—is not supported by the analysis as written. The significance is computed as the mean of the reconstructed concurrence lower bound c²_MB divided by its standard deviation, with no reference to the null hypothesis of a separable state. Since c²_MB is a quadratic function of the estimated spin coefficients, its expectation value carries a positive variance-induced bias. The paper's own Fig. 4 shows this: the 300/fb band exceeds 1, which is unphysical for a concurrence bounded by 1. That bias, not entanglement, is what the 2σ number measures. The pseudo-experiment description is also internally inconsistent: 1000 pseudo-experiments with the event counts in Table I would require 5.8M, 12.7M, and 126.9M events for 137/fb, 300/fb, and 3/ab, not the stated one million. Unless a resampling scheme is spelled out, the quoted uncertainties are not reproducible.\n\nWhat the paper does well is set up a concrete LO pipeline—MadGraph plus MadSpin—for reconstructing the ZZ spin density matrix from lepton angular distributions, including the symmetrization needed for identical Z bosons. That is a legitimate extension of the methodology, and the reconstruction coefficients in Fig. 3 look plausible. The qualitative result that SM ZZ production yields entangled pairs is almost certainly correct and consistent with the earlier work in [50], which weakens the novelty claim. The authors should state explicitly what their analysis adds beyond [50].\n\nSmaller issues: there is no detector simulation or NLO corrections; that's fine for an idealized projection, but it should be labeled as such rather than presented as a realistic experimental sensitivity. The background rejection argument, based solely on m4l shapes, is too quick to be convincing, though it's not load-bearing.\n\nBottom line: this is a serious effort at a timely topic, but the central numerical claims are not credible until the authors run a closure test and a null-hypothesis test, and clean up the pseudo-experiment statistics. I would not publish as is. I would, however, send it to peer review with a request for major revision—the methodology is worth refereeing rather than desk rejecting.","headline":"ZZ entanglement sensitivity study with a solid LO pipeline but unsupported significance claims due to a biased quadratic estimator and inconsistent pseudo-experiment counts.","tokens_in":13265,"tokens_out":4591,"would_cite":false,"duration_ms":50105,"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":"ZZ pairs at the LHC are entangled: simulated data show a positive concurrence lower bound, reaching 3.75σ at HL-LHC luminosity.","keywords":["quantum entanglement","concurrence","ZZ production","quantum state tomography","LHC","spin density matrix","Gell-Mann matrices","pseudo-experiments"],"falsifier":"Run a closure test on the same 1000 pseudo-experiments: generate events from a known density matrix, reconstruct its C², and compare the mean to the known value. If the reconstructed mean exceeds the maximum possible value of one at 137 fb⁻¹, the estimator is biased and the quoted 2σ and 3.75σ significances do not hold.","tokens_in":12215,"feed_emoji":"🔗","tokens_out":5504,"duration_ms":63056,"temperature":0.7,"pith_summary":"This paper argues that the ZZ pairs produced in proton-proton collisions at the LHC carry a quantum-entangled spin state, and that the entanglement can be seen in data. Using quantum state tomography on simulated pp→ZZ→4ℓ events at 13 TeV, the authors reconstruct the two-qutrit spin density matrix from lepton angular distributions and compute the lower bound of concurrence C²_MB. They find C²_MB converges to about 0.375 as luminosity grows, with significance 2σ for combined Run 2+3 data and 3.75σ for the HL-LHC sample of 3 ab⁻¹. If correct, this would make di-boson production a new experimental platform for quantum information tests in high-energy collisions.","feed_headline":"ZZ entanglement reaches 3.75σ in simulated LHC data","feed_subtitle":"Tomography of pp→ZZ→4ℓ finds a positive concurrence bound, opening a new collider test of quantum mechanics.","key_machinery":"The central object is the two-qutrit spin density matrix of the ZZ pair, parametrized with Gell-Mann matrices and reconstructed through quantum state tomography from the angular distribution of the decay leptons. The entanglement observable is the analytic lower bound of concurrence for mixed states, C² = 2(Tr[ρ²] - Tr[ρ_A²] - Tr[ρ_B²]), expressed in terms of the extracted polarization and correlation coefficients. The extraction uses Wigner P-symbol angular functions, and statistical sensitivity is assessed by dividing one million simulated events into 1000 pseudo-experiments at three luminosity benchmarks.","core_discovery":"At leading order in the Standard Model, the spin state of the on-shell ZZ system produced in pp collisions is not a classical mixture: the reconstructed lower bound of concurrence is positive, meaning the two Z bosons are entangled. The authors obtain this by writing the density matrix with Gell-Mann matrices, extracting the polarization and correlation coefficients A_i, B_i, C_ij from the angular distribution of the final leptons, and applying the analytic lower bound C² = -4/9 - 6∑A_i² - 6∑B_i² + 8∑C_ij². Averaged over 1000 pseudo-experiments, C²_MB approaches 0.375, and the statistical significance grows from about 2σ at 137-300 fb⁻¹ to 3.75σ at 3 ab⁻¹. The paper also notes that at lower luminosity the reconstructed value can exceed the physical maximum of one, which it attributes to statistical fluctuation.","pith_inferences":["A closure test on the reconstruction, injecting events from a known density matrix, would show whether the physically impossible values above one seen at low luminosity are a bias or a fluctuation; if biased, the quoted significances would need downward revision.","Because the extraction uses leading-order angular distributions, higher-order QCD and electroweak corrections could shift C²_MB, so the HL-LHC sensitivity estimate likely overstates the reach until those corrections are folded in.","The same lower-bound observable could be used to search for anomalous ZZ couplings, since such couplings would alter the correlation coefficients and hence the measured concurrence."],"forward_implications":["A positive measurement of C²_MB in real HL-LHC data would establish entanglement between massive gauge bosons without relying on Bell-inequality assumptions.","The same tomography pipeline can be applied to other boson-pair channels, such as WW and ZH, to compare their quantum-correlation content.","The convergence of C²_MB with luminosity provides a direct, background-suppressed measurement of Standard Model ZZ spin correlations at 13 TeV.","If the 3.75σ significance holds with detector-level effects included, the ZZ channel becomes a second high-energy platform for entanglement, complementing top-quark pairs."],"supporting_citations":[{"why":"Provides the normalized differential cross-section formula (Eq. 10) connecting the decay-lepton angles to the joint spin density matrix.","marker":"[66]"},{"why":"Supplies the analytic lower bound on concurrence for mixed bipartite states used as the entanglement observable.","marker":"[65]"},{"why":"Provides the Gell-Mann basis used to parametrize the two-qutrit spin density matrix.","marker":"[62]"},{"why":"Supplies the Monte-Carlo event generator used to produce the simulated pp→ZZ signal events.","marker":"[69]"},{"why":"Provides the spin-correlation treatment of Z decays in the simulation.","marker":"[70]"},{"why":"Defines the Wigner P angular functions used to extract the polarization and correlation coefficients from lepton directions.","marker":"[73]"},{"why":"Also defines the angular functions and extraction formalism used in the tomography.","marker":"[74]"}],"fun_headline_variants":["ZZ pairs show quantum entanglement in simulated LHC data","Entanglement in ZZ production: up to 3.75σ at HL-LHC","Z bosons entangled: 3.75σ significance from LHC sim","Quantum entanglement in ZZ collisions: 3.75σ in sim"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claimed significances assume that the reconstructed spin-correlation coefficients are unbiased estimates of the true values; because the same procedure produces physically impossible values above one at low luminosity, that assumption has not been demonstrated.","fun_headline_variants_meta":{"raw":{"variants":["ZZ pairs show quantum entanglement in simulated LHC data","Entanglement in ZZ production: up to 3.75σ at HL-LHC","Z bosons entangled: 3.75σ significance from LHC sim","Quantum entanglement in ZZ collisions: 3.75σ in sim"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001159,"raw_usage":{"total_tokens":4796,"prompt_tokens":938,"completion_tokens":3858,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":3780}},"tokens_in":554,"tokens_out":3858,"duration_ms":33316,"temperature":1.0,"reasoning_tokens":3780,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:44:42.478399+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a closure test on the same 1000 pseudo-experiments: generate events from a known density matrix, reconstruct its C², and compare the mean to the known value. If the reconstructed mean exceeds the maximum possible value of one at 137 fb⁻¹, the estimator is biased and the quoted 2σ and 3.75σ significances do not hold.","supporting_citations":[{"cited_title":"Breaking down the entire spectrum of spin correlations of a pair of particles involving fermions and gauge bosons,","cited_arxiv_id":null,"evidence_quote":"Provides the normalized differential cross-section formula (Eq. 10) connecting the decay-lepton angles to the joint spin density matrix."},{"cited_title":"Observable entanglement measure for mixed quantum states,","cited_arxiv_id":null,"evidence_quote":"Supplies the analytic lower bound on concurrence for mixed bipartite states used as the entanglement observable."},{"cited_title":"Symmetries of baryons and mesons,","cited_arxiv_id":null,"evidence_quote":"Provides the Gell-Mann basis used to parametrize the two-qutrit spin density matrix."},{"cited_title":"The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations,","cited_arxiv_id":null,"evidence_quote":"Supplies the Monte-Carlo event generator used to produce the simulated pp→ZZ signal events."},{"cited_title":"Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations,","cited_arxiv_id":null,"evidence_quote":"Provides the spin-correlation treatment of Z decays in the simulation."},{"cited_title":"Bell inequalities and quantum entanglement in weak gauge boson production at the lhc and future colliders,","cited_arxiv_id":null,"evidence_quote":"Defines the Wigner P angular functions used to extract the polarization and correlation coefficients from lepton directions."},{"cited_title":"Quantum state tomography, entanglement detection and bell violation prospects in weak decays of massive particles,","cited_arxiv_id":null,"evidence_quote":"Also defines the angular functions and extraction formalism used in the tomography."}],"review_version":1}