{"id":"108220af-e52d-4275-9077-a4539da292d8","arxiv_id":"2607.19495","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"First measurement of the ortho-positronium→3γ Dalitz plot over nearly the entire phase space, consistent with QED.","lead":"Physicists used the J-PET plastic-scintillator detector to record, for the first time, the full angular pattern (Dalitz plot) of photons from ortho-positronium decaying into three gamma rays. The measured pattern agrees with standard QED predictions to within a few percent across most of the phase space.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Efficiency and misreconstructed-signal corrections derived from LO QED MC can bias extracted Dalitz plot toward LO, weakening the QED-consistency claim.","rationale":"The reader's weakest assumption matches the main soft spot: reliance on LO-generated MC for efficiency and misreconstructed background makes the QED-consistency conclusion partially circular. This is load-bearing because NLO corrections are at the 1.5–2.5% level, comparable to the quoted statistical/systematic uncertainties, and the systematic studies do not vary the generator model. The first Dalitz-plot measurement is still credible, so the paper merits conditional acceptance pending a model-independent cross-check or a stated caveat; no change to the reader's verdict.","tokens_in":9545,"tokens_out":4709,"duration_ms":48751,"concrete_test":"Reanalyze the same experimental data using an efficiency map and misreconstructed-signal template generated from the NLO O(α) matrix element of Ref. [25] instead of LO. Compare the resulting corrected Dalitz plot bin-by-bin with the LO-corrected result; if any bin shifts by more than the quadrature sum of the quoted statistical and systematic uncertainties, the LO-based correction introduces a significant bias and the consistency claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5 corrects the measured Dalitz distribution using an efficiency map obtained from Monte Carlo events generated with the LO QED matrix element (Eq. 2), and Eq. (4) subtracts a misreconstructed-signal template generated from the same LO model. If the true o-Ps→3γ distribution differs from LO—e.g. by the O(α) corrections shown in Fig. 2 to be 1.5–2.5% in the accessible region—the per-bin acceptance and the misreconstructed fraction can shift. Since the fit uses the same LO template for both signal components, the extracted distribution is pulled toward LO wherever misreconstructed or acceptance corrections are non-negligible. The quoted systematic uncertainties were estimated by varying selection cuts, not by varying the generator model, so this potential bias is not included in the 2–3% systematic error. The headline consistency with both LO and NLO QED is therefore at least partly a consequence of the analysis procedure; the first-measurement claim remains, but the QED test is not independent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports the first measurement of the Dalitz plot for ortho-positronium decay into three photons, performed with the J-PET detector. Using a 22Na source and a porous XAD-4 target, the authors identify 33 million o-Ps→3γ candidates. After event selection based on hit multiplicity, time-over-threshold, decay-plane distance, timing, and angular correlations, a one-dimensional fit (Eq. 4) separates reconstructed signal, misreconstructed signal, and background, with the two signal templates taken from Monte Carlo generated with the leading-order QED matrix element. The final angular Dalitz plot is corrected by an efficiency map obtained from the same Monte Carlo and is stated to be consistent with both leading-order and next-to-leading-order QED predictions, with about 3% statistical and 2-3% systematic uncertainty over most of the phase space.","tokens_in":9809,"tokens_out":6172,"duration_ms":70356,"significance":"If the result holds, this is the first nearly complete experimental map of the differential three-photon decay distribution of ortho-positronium, a quantity that has been poorly explored since the early angular measurements of Mills and Berko. The detector description, the large event sample, and the explicit treatment of backgrounds are strengths. However, the QED-consistency conclusion is weakened by the fact that the acceptance/efficiency correction and the misreconstructed-signal template are both derived from the leading-order theory under test. The measurement itself is valuable, but the paper currently overstates the independence of the QED test.","major_comments":[{"comment":"The acceptance/efficiency correction in §5 and the misreconstructed-signal template subtracted in the fit of Eq. (4) are both computed from Monte Carlo generated with the LO QED matrix element (Eq. 2; Fig. 8a). The final Dalitz plot is therefore corrected with a model-dependent map. If the true o-Ps→3γ distribution differs from LO, for example by the O(α) corrections shown in Fig. 2 to be 1.5–2.5% in the accessible region, the per-bin average efficiency and the misreconstructed fraction can shift. Because the same amplitude A is used for reconstructed and misreconstructed templates in Eq. (4), a deviation from LO can be partially absorbed by the fit normalization. The systematic uncertainties are estimated only by varying selection cuts; no variation of the generator model or closure test is reported. Thus the statement that the result is 'consistent with both the leading-order and next-","section":"§5 and Eq. (4)"},{"comment":"No quantitative comparison with theory is provided. The text says the distribution is in agreement with LO QED (Sec. 6) and consistent with NLO (abstract), but the paper reports no χ², residual distribution, or confidence level for data versus either LO or NLO predictions. Since the bin-by-bin uncertainties are about 3% statistical and 2–3% systematic, a residual plot with uncertainties (or a data table) is necessary to judge whether the agreement is good or merely not excluded. This is central to the consistency claim.","section":"§5/§6"},{"comment":"The statistical model in Eq. (4) is not standard as written. The denominator contains A²(N_sig_i + N_mis_i) + B² N_bcg_i. If N_sig_i and N_bcg_i are template counts scaled by A and B, the expected Poisson variance of the model contribution is A N_sig_i + B N_bcg_i, not A² times the template. If the templates are normalized to unit integrals, the meaning of A,B and this variance should be stated. As written, the χ² value and the fitted uncertainties depend on an arbitrary convention. Please clarify or correct the expression; otherwise the quoted statistical uncertainty of the extracted distribution is not well defined.","section":"Eq. (4)"}],"minor_comments":[{"comment":"The y-axis label says 'NLO−NLO+O(α)' while the caption/text describe the difference between LO and LO+O(α). Please make the notation consistent.","section":"Fig. 2"},{"comment":"The right subpanels of Fig. 8 show systematic uncertainties for 'regions 1 and 2', but these regions are not defined in the text. Please define the boundaries of the regions or refer to the green region explicitly.","section":"§5, Fig. 8"},{"comment":"There is a typo: 'Kracow' should be 'Kraków'.","section":"§2"},{"comment":"For a first measurement of this kind, it would be helpful to provide a data table of the final Dalitz plot values (and uncertainties) as a supplement or in HEPData, so that the result can be used by theorists.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The first-measurement claim is likely sound and the collaboration has a strong record in positronium instrumentation. The main risk is over-interpretation of the QED consistency test because the corrections are generated from the theory being tested. I do not see this as a reason for rejection, but it must be addressed in revision, along with a quantitative comparison and clarification of the fit statistic. If the authors add model-dependence checks and a residual/χ² analysis, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper really does measure the angular Dalitz distribution for o-Ps→3γ for the first time, and the experimental effort is impressive. But the headline claim—consistency with LO and NLO QED—is not a test of QED. The efficiency map and the misreconstructed-signal template are both generated from the LO QED matrix element, and the final corrected distribution is exactly the LO template shape once the fit does its job. So the agreement is circular.\n\nWhat's genuinely new: 33M o-Ps→3γ candidates, a 15°-binned Dalitz plot over almost the full phase space, and the demonstration that J-PET can select and reconstruct this decay with 0.03% overall efficiency and 84% purity. The detector work, the TOT-based photon ID, and the background rejection using angular correlations are credible and well described. The first-measurement claim deserves to stand, but with a caveat: it is a first measurement of the raw event distribution folded with a LO-based acceptance correction, not a model-independent map of the decay.\n\nThe soft spot is central. The fit in Eq. (4) uses N_sig and N_mis from LO MC. After subtracting the best-fit background and misreconstructed components and dividing by the LO efficiency, the resulting 'measured' plot is mathematically pinned to the LO shape. Deviations in the data would show up as a bad chi-square, but no goodness-of-fit is reported, and the systematic errors are from cut variations only—they don't include the generator-model dependence. That means the 2–3% systematic is an underestimate of the model uncertainty, and the claimed agreement with NLO (1.5–2.5% differences) cannot be tested by this analysis.\n\nThe fix is straightforward: repeat the efficiency correction and template fitting using an NLO generator or a phase-space-like flat generator, and quote the difference as a systematic. Without that, the paper overstates what it demonstrates.\n\nThe paper is worth a serious referee because the first-measurement claim is valuable and the detector achievement is real. But the QED consistency claim should be removed from the abstract or heavily qualified, and the model-dependence systematic must be added. I'd send it to review with a request for major revision, not desk reject.","headline":"First Dalitz plot measurement is real, but the QED agreement is baked into the analysis via LO-based corrections.","tokens_in":10479,"tokens_out":4360,"would_cite":false,"duration_ms":45183,"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":"This paper presents the first measurement of the Dalitz plot for ortho-positronium decay into three photons, across nearly all of the allowed phase space.","keywords":["ortho-positronium","Dalitz plot","three-photon annihilation","QED test","J-PET detector","positronium decay","differential decay rate","angular distribution"],"falsifier":"Re-analyze the same 33 million events using an efficiency map and a misreconstructed-signal template generated from the next-to-leading-order QED matrix element instead of the leading-order one; if any bins of the corrected Dalitz plot shift by more than the quoted systematic uncertainties, the extraction is not independent of the theoretical input. Alternatively, a bin-by-bin deviation pattern in the corner regions that persists under independent variations of all selection cuts would show the true distribution does not match the assumed QED shape.","tokens_in":9447,"feed_emoji":"⚛️","tokens_out":3858,"duration_ms":39504,"temperature":0.7,"pith_summary":"The paper claims to have measured, for the first time, the Dalitz plot for the three-photon decay of ortho-positronium over almost the full allowed kinematic region. Using a large-acceptance plastic scintillator detector and 33 million identified decay events, the authors map the differential decay distribution as a function of the angles between the three photons, with statistical uncertainties around 3% and systematic uncertainties of 2–3% across most of the plot. The measured distribution agrees with the leading-order QED prediction and with the next-to-leading-order corrected prediction within these uncertainties. A sympathetic reader would care because the Dalitz plot is the full differential observable for a fundamental three-body QED decay; until now only three narrow angular configurations had been measured. If correct, this opens a new experimental window on testing higher-order QED corrections in a bound-state system.","feed_headline":"First Dalitz plot maps positronium's three-photon decay","feed_subtitle":"A scintillator detector measured the full angular distribution and found agreement with QED at the few-percent level.","key_machinery":"The central object is the Dalitz plot in angular representation — a two-dimensional histogram of decay-event density over two independent angular variables that fully specify the three-photon final state. The mechanism that carries the argument is the event-selection and correction pipeline: three-photon hits are selected by timing, energy-deposition, and angular-correlation cuts; a template fit subtracts misreconstructed signal and background; and a Monte Carlo efficiency map, built by simulating events generated from the leading-order QED amplitude and passing them through the same selection, corrects the measured distribution bin by bin. The comparison of the corrected plot to leading-ord","core_discovery":"On its own terms, the paper reports the first experimental determination of the Dalitz plot for ortho-positronium decay into three photons, covering almost the entire available phase space except for low-energy-photon configurations with relative angles larger than about 170 degrees. After subtracting misreconstructed signal and background using a template fit and correcting with a Monte Carlo efficiency map, the resulting angular distribution is consistent with both the leading-order QED matrix element and the O(alpha)-corrected prediction. The measured distribution contains 33.42 million identified o-Ps to 3-gamma events and has statistical uncertainty of about 3% and systematic uncertaint","pith_inferences":["Because the efficiency map and the misreconstructed-signal template are both generated from Monte Carlo events produced with the leading-order QED matrix element, the correction procedure will pull the measured Dalitz plot toward the leading-order shape if the true distribution differs bin-by-bin; the agreement with QED should therefore be read as a consistency check rather than a fully model-inde","The same template-subtraction analysis could be rerun with a shape model that allows free deviations—for example, a bin-by-bin amplitude factor or an NLO-shape template—turning the measurement into a genuine fit for higher-order or symmetry-violating coefficients.","If the projected order-of-magnitude statistical improvement is achieved, the 1.5–2.5% difference between LO and NLO predictions shown in the paper's comparison becomes resolvable, making the Dalitz plot a sensitive, multidimensional probe of O(alpha) QED corrections in positronium decay.","The angular representation chosen here is well matched to a detector that measures photon directions with roughly 1-degree resolution but has an energy threshold near 30 keV; lowering that threshold in future detectors could extend the plot into the currently missing corner and directly test the expected suppression of the matrix element at the kinematic boundary."],"forward_implications":["The measurement provides the first nearly complete map of the o-Ps to 3-gamma differential decay rate, replacing the three discrete angular points measured in 1967.","The data set a reference for future tests of higher-order QED corrections, since the difference between leading-order and O(alpha) predictions in the region of interest is about 1.5–2.5% and the current statistical uncertainty is about 3%.","The result validates the J-PET detector concept for multiphoton final-state measurements, which is relevant to positronium imaging and discrete-symmetry searches that rely on correlated photon kinematics.","The measurement quantifies the unmeasured corner of the Dalitz plot (low-energy photons with relative angles above 170 degrees), showing where detector thresholds currently limit sensitivity and where future detectors could extend coverage.","Within uncertainties, the agreement with both leading-order and next-to-leading-order QED predictions supports the use of the standard Ore-Powell amplitude as a reliable input for positronium-related simulations and medical-imaging applications."],"fun_headline_variants":["First Dalitz plot of positronium's three-photon decay","J-PET maps positronium decay into three photons","First experimental Dalitz plot for ortho-positronium decay","Positronium's three-photon decay: first Dalitz plot map","Dalitz plot for positronium decay measured for first time"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The measured Dalitz plot is corrected with an efficiency map, and the background is subtracted using templates, both built from Monte Carlo events generated with the leading-order QED matrix element; if the true differential distribution differs from leading order in a way that changes the per-bin acceptance, the correction pulls the measurement toward the predicted shape, so the 'consistent with QED' conclusion is at least partly a consequence of the analysis procedure rathe","fun_headline_variants_meta":{"raw":{"variants":["First Dalitz plot of positronium's three-photon decay","J-PET maps positronium decay into three photons","First experimental Dalitz plot for ortho-positronium decay","Positronium's three-photon decay: first Dalitz plot map","Dalitz plot for positronium decay measured for first time"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1113,"prompt_tokens":619,"completion_tokens":494,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":363,"completion_tokens_details":{"reasoning_tokens":418}},"tokens_in":363,"tokens_out":494,"duration_ms":4889,"temperature":1.0,"reasoning_tokens":418,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T12:33:45.574642+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-analyze the same 33 million events using an efficiency map and a misreconstructed-signal template generated from the next-to-leading-order QED matrix element instead of the leading-order one; if any bins of the corrected Dalitz plot shift by more than the quoted systematic uncertainties, the extraction is not independent of the theoretical input. Alternatively, a bin-by-bin deviation pattern in the corner regions that persists under independent variations of all selection cuts would show the true distribution does not match the assumed QED shape.","supporting_citations":[],"review_version":1}