REVIEW 2 major objections 6 minor 1 cited by
This paper projects that a Tera-Z electron-positron collider can constrain the rare decay Z→ννγ to branching ratios near 10^-9, several orders of magnitude below the current LEP bound.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-04 07:59 UTC pith:NI5WXP2K
load-bearing objection Useful Tera-Z projection for Z→ννγ, but the missing SM-EFT interference term and some typos make the headline BR limits less secure than they look. the 2 major comments →
Search for new physics effects in νbar{ν}γ production at a Tera-Z factory
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
On its own terms, the paper establishes a sensitivity projection: with √s=91.2 GeV and L_int=150 ab^-1, the e+e−→ννγ channel, after cuts requiring a photon with Eγ>4 GeV, missing transverse energy /ET>4 GeV with /ET≤Eγ, and missing-energy significance S_ET>16, yields an expected 3σ upper limit on BR(Z→ννγ) of 1.73×10^-9 for the dimension-6 coupling κ/Λ^2 (and 1.77×10^-9 for the dimension-8 coupling α8/Λ^4) when systematic uncertainty is neglected. The same analysis reaches 5σ evidence for κ/Λ^2≲1 TeV^-2 and α8/Λ^4 on the order of 10 TeV^-4. With 5% systematics the limits degrade to about 10^-7, still well below the LEP bound. The paper treats the EFT signal simply as an additive contribution
What carries the argument
The central object is an effective Lagrangian for the Zννγ vertex: a dimension-6 term with CP-odd and CP-even couplings κ1/Λ^2 and κ2/Λ^2, and a dimension-8 term with coupling α8/Λ^4, added to the Standard Model. The authors implement this vertex in a Monte Carlo event generator with detector simulation, then use three kinematic handles—the photon energy Eγ, the missing transverse energy /ET, and its significance S_ET=/ET^2/E_γ^T—to separate the monophoton signal from the dominant e+e−→ννγ and e+e−→ℓ+ℓ−γ backgrounds. The cut sequence and the significance formula in Eq. (9) carry the whole projection: they convert the difference between signal-plus-background and background alone into the rep
Load-bearing premise
The load-bearing premise is that the anomalous Zννγ signal and the Standard Model ννγ background can be added independently, with no interference between the EFT coupling and the SM one-loop amplitude; if that interference is significant, the projected branching-ratio limits change.
What would settle it
Compute, for a benchmark such as κ/Λ^2=5 TeV^-2 or α8/Λ^4=15 TeV^-4, the complete e+e−→ννγ matrix element containing both the SM and EFT amplitudes, apply the same selection (Eγ>4 GeV, /ET>4 GeV, /ET≤Eγ, S_ET>16), and compare the event count with the paper's independent S+B sum. A difference comparable to the statistical uncertainty would mean the quoted limits need revision; a difference much smaller would support them.
If this is right
- A Tera-Z run can test the one-loop Standard Model prediction for Z→ννγ (7.16×10^-10) for the first time, since the projected 3σ reach at δ_sys=0 is about 1.7×10^-9, within a factor of roughly 2.4 of the SM value.
- New-physics couplings κ/Λ^2 and α8/Λ^4 can be bounded at the TeV^-2 / TeV^-4 scale, probing new physics scales far beyond direct reach.
- A simple cut-based monophoton selection is sufficient to suppress the huge e+e−→ℓ+ℓ−γ background by roughly four orders of magnitude while keeping about 5% of the signal.
- Even at 5% systematic uncertainty, the projected limits (≈3.6×10^-7) improve on the LEP bound of ≈10^-6, so the Tera-Z program has discovery potential even in a conservative scenario.
Where Pith is reading between the lines
- Not addressed in the paper: the EFT signal and SM ννγ background are generated separately and summed (S+B) without computing the interference between the new-physics vertex and the SM amplitudes, even though they share the same final state. If that interference is non-negligible, the quoted κ/Λ^2 and α8/Λ^4 limits—and the BR values in Table II—could shift appreciably.
- A testable extension is to run the same analysis at the off-peak energies shown in the paper (87.7 and 93.9 GeV): the ratio of event rates would separate the resonant Z-pole signal from non-resonant backgrounds and directly check the EFT signal hypothesis.
- A further extension: because the dimension-6 Lagrangian contains both a CP-odd (κ1) and a CP-even (κ2) coupling, a photon angular or polarization analysis could distinguish these two couplings rather than only constraining the sum κ1^2+κ2^2, which the current integrated rate cannot do.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a search for anomalous Zννγ couplings through e+e−→ννγ production at the FCC-ee/CEPC Tera-Z stage (√s=91.2 GeV, L_int=150 ab^-1). The authors implement dimension-6 (κ1,κ2/Λ^2) and dimension-8 (α8/Λ^4) operators in a UFO model, generate signal and SM background events with MadGraph/Pythia/Delphes, and apply a cut-based selection using Eγ, /ET, and S_ET. From significance curves they derive projected upper limits on the couplings and, via the analytic width relations of Eqs. (7)–(8), on BR(Z→ννγ), obtaining values of order 10^-9 at δ_sys=0 and a few times 10^-7 at δ_sys=5% (Table II). The paper claims this is an improvement of several orders of magnitude over the LEP bound and demonstrates the potential of Tera-Z to test rare Z decays.
Significance. If the projected sensitivities are correct, the analysis would be a useful addition to the Tera-Z physics case: the Z→ννγ channel is not strongly constrained by LEP, and a future Z-pole run could test loop-level SM predictions. The use of a full simulation chain (MadGraph+Pythia+Delphes) with a detailed cut-flow table is a strength, and the analytic decay-width expressions provide a transparent mapping between couplings and branching ratios. However, the central quantitative claims are not yet supported because two key ingredients are missing: (i) interference between the EFT amplitude and the SM e+e−→ννγ amplitude, which is not modeled, and (ii) an explicit charged-lepton veto, despite a large surviving l+l−γ background in Table I. Both are standard, calculable additions and could change the limits in Fig. 9 and Table II substantially. The paper's broad conclusion is plausible, but the specific numbers should be treated as preliminary until these omissions are addressed.
major comments (2)
- [Sec. II, Eq. (9)] The signal (Fig. 4a) and the SM e+e−→ννγ background (Fig. 4b–d) are generated separately and combined as S+B in Eq. (9), but both amplitudes have identical external states. The physical rate is |M_SM + M_EFT|^2, and the interference term 2Re(M_SM M_EFT*) is absent from the analysis. No helicity, CP, or kinematic argument is given for its vanishing; it is in principle computable in the same MadGraph setup. Since the coupling limits in Fig. 9 and the BR limits in Table II are extracted from these S values, this omission is load-bearing. The authors should generate the complete e+e−→ννγ process including both SM and EFT contributions, or present an explicit amplitude-level argument showing the interference vanishes after integration over the selected phase space.
- [Table I, Sec. II] The cut-flow table reports 5.37×10^6 surviving l+l−γ events after the final S_ET>16 requirement, comparable to the 8.17×10^6 SM ννγ background. The text never describes a charged-lepton veto or isolation requirement; the baseline selection is only Nγ>0. Since l+l−γ events contain two visible charged leptons, they cannot mimic a single-photon + missing-energy signature unless both leptons are lost. A standard lepton veto would suppress most of this background. As written, the background composition, significance curves (Fig. 9), and resulting Table II limits are not representative of a genuine ννγ selection. A documented lepton veto must be applied, or the survival mechanism for these events must be demonstrated.
minor comments (6)
- [Table II] Several entries appear to have typographical issues: the α8/Λ4 3% δsys, 5σ entry reads '3.72×10^7' and should presumably be 10^-7; also the α8/Λ4 5% δsys 3σ and 5σ limits (6.18×10^-7 and 6.19×10^-7) are almost identical, which is suspicious since a 5σ limit should be weaker than a 3σ limit. Please check and correct.
- [Sec. I, Eqs. (7)–(8)] The integration limits for the decay width are described as 'minimum value for the photon energy greater than half of the beam energy' and maximum equal to M_Z. The relation between these limits and the analysis cuts (Eγ>4 GeV) is not transparent; clarify whether the width used for Figure 3 and Table II matches the phase space accepted by the experimental selection.
- [Fig. 7] The left-panel caption states κ/Λ^2=5.0 TeV^-4, while the text and other figures use TeV^-2; this is presumably a typo. Please correct.
- [Sec. II, Table I] The text names e+e−→γγγ as a subleading background, but Table I does not include it in the cut-flow. Either include it in the table or explain why it is negligible after the selection.
- [Abstract and conclusion] The claim of 'several orders of magnitude' improvement over LEP is only valid for δ_sys=0. At δ_sys=5%, the limits in Table II are a few times 10^-7, which is only a factor of about 3 below the LEP bound of 10^-6. Please qualify this statement in the abstract and conclusion.
- [Sec. III] The limits on BR(Z→ννγ) in Table II are model-dependent translations of coupling limits via Eqs. (7)–(8), not direct, model-independent measurements. The text should state this more explicitly to avoid the impression that the analysis constrains the SM branching ratio itself.
Circularity Check
No significant circularity: the BR limits are a deterministic translation of coupling limits within the same EFT model, and the SM prediction is an external input.
full rationale
The derivation chain is: (i) adopt EFT operators from the literature (Eqs. 1-6); (ii) generate signal and SM background events with MadGraph/Pythia/Delphes; (iii) apply kinematic cuts and compute statistical significance via Eq. (9); (iv) read coupling limits from the significance curves in Fig. 9; (v) convert those coupling limits to BR(Z->nu nu gamma) limits using the analytic width relations in Eqs. (7)-(8). Step (v) is a deterministic one-to-one map within the same EFT model; it does not fit a branching ratio to data and rename it as a prediction. The SM branching ratio 7.16e-10 is taken from the external Ref. [14] and is not used to construct the exclusion. There are no load-bearing self-citations: the EFT operators and SM calculation are cited from other authors, and no uniqueness theorem or ansatz is imported from the present authors' prior work. The possible concern that the signal and SM background are added without interference in Eq. (9) is a physics/completeness issue affecting sensitivity estimates, not a circularity of the derivation chain. Therefore the circularity score is 0.
Axiom & Free-Parameter Ledger
free parameters (4)
- κ1/Λ^2 and κ2/Λ^2 (dimension-6 Zννγ couplings) =
Scanned 1–10 TeV^-2; 3σ/5σ limits implied by Table II but not quoted in TeV^-2
- α8/Λ^4 (dimension-8 coupling) =
Scanned 5–40 TeV^-4; 5σ sensitivity around 10 TeV^-4 for δsys=0
- Systematic uncertainty δsys =
0%, 1%, 3%, 5%
- Analysis cut thresholds =
Eγ>4 GeV, /ET>4 GeV, /ET≤Eγ, S/ET>16
axioms (4)
- domain assumption The effective operator basis in Eqs. (1)–(6), taken from Refs. [15,16], is the complete set of Zννγ couplings at dimension 6 and 8.
- ad hoc to paper Interference between the anomalous Zννγ amplitude and the SM e+e−→ννγ background is negligible.
- domain assumption The simulated SM backgrounds, including l+l−γ after the missing-energy cuts, are modeled correctly by MadGraph+Pythia+Delphes with the IDEA card.
- domain assumption The SM one-loop Z→ννγ amplitude is negligible for the new-physics limits considered.
read the original abstract
Rare decays of the Z boson provide a sensitive probe for physics beyond the Standard Model (SM). This study investigates the $e^{+}e^{-} \to Z \to \nu\bar{\nu}\gamma$ process within the context of the Tera-Z programmes at future colliders such as the FCC-ee and CEPC. The SM predicts a one-loop branching ratio of $7.16 \times 10^{-10}$ for $Z \to \nu\bar{\nu}\gamma$, a value four times smaller than the current experimental limit from the LEP. To explore this window for new physics, we parameterize anomalous $Z\nu\bar{\nu}\gamma$ interactions using an Effective Field Theory framework, considering both dimension-6 and dimension-8 operators. A detailed simulation is performed by generating signal and background events with MadGraph, modeling particle showers with Pythia, and simulating detector effects with Delphes. The analysis employs key kinematic variables-including the photon energy ($E_\gamma$), missing transverse energy ($\not{E}_T$), and the missing transverse energy significance ($S_{\not{E}_T}$) to isolate the signal. The results yield upper limits on the anomalous couplings, from which we infer branching ratios for $Z \to \nu\bar{\nu}\gamma$ on the order of $10^{-9}$. This represents a significant improvement of several orders of magnitude over the LEP sensitivity. Consequently, this study demonstrates the unique potential of the Tera-Z runs not only to test the SM loop-level predictions with unprecedented precision but also to tightly constrain or reveal new anomalous interactions.
Figures
Forward citations
Cited by 1 Pith paper
-
A framework for missing-energy searches with anomalous light vectors
For light U(1)X vectors coupled to anomaly-sensitive currents, Wess-Zumino coefficients are fixed by anomaly matching, which lets all missing-energy searches (K, B, D, Z) be mapped onto one (m_X,g_X) plane with minima...
Reference graph
Works this paper leans on
-
[1]
M. A. Perez, G. Tavares-Velasco and J. J. Toscano, Int. J. Mod. Phys. A19(2004), 159-178 [arXiv:hep- ph/0305227 [hep-ph]]
arXiv 2004
-
[2]
R. Mohanta, Phys. Rev. D71(2005), 114013 doi:10.1103/PhysRevD.71.114013 [arXiv:hep-ph/0503225 [hep-ph]]
Pith/arXiv arXiv 2005
-
[3]
D. N. Gao and X. Gong, Chin. Phys. C47(2023) no.4, 043106 [arXiv:2208.12652 [hep-ph]]
Pith/arXiv arXiv 2023
-
[4]
A. Blondel, J. Gluza, S. Jadach, P. Janot, T. Riemann, A. Akhundov, A. Arbuzov, R. Boels, S. Bon- darenko and S. Borowka,et al.doi:10.23731/CYRM-2019-003 [arXiv:1809.01830 [hep-ph]]
Pith/arXiv arXiv 2019
-
[5]
V. Maura, B. A. Stefanek and T. You, JHEP10(2025), 022 [arXiv:2412.14241 [hep-ph]]
Pith/arXiv arXiv 2025
-
[6]
L. Allwicher, M. McCullough and S. Renner, JHEP02(2025), 164 [arXiv:2408.03992 [hep-ph]]
Pith/arXiv arXiv 2025
-
[7]
J. B. Guimarães da Costaet al.[CEPC Study Group], [arXiv:1811.10545 [hep-ex]]
-
[8]
Abadaet al.[FCC], Eur
A. Abadaet al.[FCC], Eur. Phys. J. ST228(2019) no.2, 261-623
2019
-
[9]
Benediktet al.[FCC], [arXiv:2505.00272 [hep-ex]]
M. Benediktet al.[FCC], [arXiv:2505.00272 [hep-ex]]
-
[10]
Chenget al.[CEPC Physics Study Group], [arXiv:2205.08553 [hep-ph]]
H. Chenget al.[CEPC Physics Study Group], [arXiv:2205.08553 [hep-ph]]
-
[11]
The proposed Z-pole operation plan of FCC-ee comprises data-taking at three center-of-mass energies: 87.7 GeV, 91.2 GeV, and 93.9 GeV
indicate that the CEPC could now deliver over3×1012 Z bosons within a two-year period. The proposed Z-pole operation plan of FCC-ee comprises data-taking at three center-of-mass energies: 87.7 GeV, 91.2 GeV, and 93.9 GeV. These energies are chosen to correspond to half- integer spin tunes, thereby enabling precise energy calibration through resonant depol...
-
[12]
Abdallahet al.[CEPC Study Group], Radiat
W. Abdallahet al.[CEPC Study Group], Radiat. Detect. Technol. Methods8(2024) no.1, 1-1105 [erratum: Radiat. Detect. Technol. Methods9(2025) no.1, 184-192] [arXiv:2312.14363 [physics.acc- ph]]
arXiv 2024
-
[13]
Acciarriet al.[L3], Phys
M. Acciarriet al.[L3], Phys. Lett. B412(1997), 201-209
1997
-
[14]
Abreuet al.[DELPHI], Z
P. Abreuet al.[DELPHI], Z. Phys. C74(1997), 577-586
1997
-
[15]
J. M. Hernandez, M. A. Perez, G. Tavares-Velasco and J. J. Toscano, Phys. Rev. D60(1999), 013004 [arXiv:hep-ph/9903391 [hep-ph]]
Pith/arXiv arXiv 1999
-
[16]
D. W. Dusedau and J. Wudka, Phys. Lett. B180(1986), 290-294
1986
-
[17]
M. Maya, M. A. Perez, G. Tavares-Velasco and B. Vega, Phys. Lett. B434(1998), 354-357 [arXiv:hep- ph/9804366 [hep-ph]]
arXiv 1998
-
[18]
C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer and T. Reiter, Comput. Phys. Commun. 183(2012), 1201-1214 [arXiv:1108.2040 [hep-ph]]. 16
Pith/arXiv arXiv 2012
-
[19]
A. Alloul, N. D. Christensen, C. Degrande, C. Duhr and B. Fuks, Comput. Phys. Commun.185(2014), 2250-2300 [arXiv:1310.1921 [hep-ph]]
Pith/arXiv arXiv 2014
-
[20]
J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli and M. Zaro, JHEP07(2014), 079 [arXiv:1405.0301 [hep-ph]]
Pith/arXiv arXiv 2014
-
[21]
T. Sjöstrand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, P. Ilten, S. Mrenna, S. Prestel, C. O. Ras- mussen and P. Z. Skands, Comput. Phys. Commun.191(2015), 159-177 [arXiv:1410.3012 [hep-ph]]
Pith/arXiv arXiv 2015
-
[22]
de Favereauet al.[DELPHES 3], JHEP02(2014), 057 [arXiv:1307.6346 [hep-ex]]
J. de Favereauet al.[DELPHES 3], JHEP02(2014), 057 [arXiv:1307.6346 [hep-ex]]. 17
Pith/arXiv arXiv 2014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.