REVIEW 3 major objections 3 minor 31 references
At the FCC-ee, invisibly decaying light scalars produced with a Z boson could be discovered up to about 80 GeV, with expected cross-section limits as low as 10^-2 fb.
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-02 11:12 UTC pith:JC7DXA4G
load-bearing objection Solid FCC-ee sensitivity study with a public model and clean analysis chain, but the benchmark parameters contradict the paper's own mixing equations, so the discovery-reach claim is unsupported as written. the 3 major comments →
Search for Invisibly Decaying Light Scalars at the FCC-ee
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the e+e- -> Z h1 (h1 -> invisible) process at sqrt(s) = 240 GeV and L = 10.8 ab^-1 can be measured with expected 95% CL upper limits on sigma(e+e- -> q qbar h1) x B(h1 -> invisible) of roughly 10^-2–10^-1 fb for scalar masses below the Z mass and 0.1–1 fb in the 80–120 GeV range. The analysis uses a singlet-scalar toy model in which a new scalar h1 mixes with the 125 GeV Higgs boson and decays invisibly to a 5 GeV dark-matter scalar. Event generation includes interference between the new scalar and the Standard Model Higgs, which matters when their masses are close. With the benchmark mixing parameter sin alpha = 0.985, the expected significance exceeds 5-sigma for
What carries the argument
The load-bearing observable is the Z-boson recoil mass, M_recoil = sqrt(s + m_jj^2 - 2 E_jj sqrt(s)), built from the hadronic Z decay; a true invisible scalar of mass M creates a peak in this variable. The analysis combines preselection cuts on the Z mass, jet topology, and missing momentum with a boosted-decision-tree classifier trained on jet kinematics and the recoil mass, then converts the classifier output into expected upper limits using an asymptotic profile-likelihood treatment with a 10% background uncertainty. The underlying model is a scalar singlet mixing with the Standard Model Higgs (mixing angle alpha), with the new scalar decaying invisibly to a dark-matter candidate, and the
Load-bearing premise
The reach numbers assume the benchmark scalar model is physically viable—in particular, that a 125 GeV state with a roughly 6.9 MeV width and about 7.5% invisible-plus-cascade branching is the observed Higgs boson, and that the low-mass, strongly mixed scalar states used for the discovery claim are not already ruled out by earlier collider searches.
What would settle it
A decisive check is to take the LEP-era decay-mode-independent bounds on e+e- -> Z S (S -> invisible) and evaluate them at the benchmark masses and couplings: if those bounds exclude the benchmark production rates of roughly 13–20 fb after scaling to LEP center-of-mass energies, the claimed discovery reach is already excluded. Independently, compare the assumed Higgs total width of about 6.9 MeV with the measured total width of the 125 GeV Higgs, which is far smaller; that comparison alone would settle whether the benchmark model is viable.
If this is right
- At 240 GeV and 10.8 ab^-1, the FCC-ee could exclude invisible-scalar production cross-sections down to about 10^-2–10^-1 fb for scalar masses below the Z mass, and 0.1–1 fb in the 80–120 GeV range.
- For scalar–Higgs mixing angles near 0.985, the expected significance exceeds 5σ for masses up to about 80 GeV, making discovery plausible in a single run.
- The reported upper limits are nearly independent of the mixing angle over 0.973–0.995, so they can serve as approximate model-independent constraints on this production topology.
- The analysis includes interference between the new scalar and the Standard Model Higgs, so its sensitivity statement remains meaningful where the two resonances overlap near 125 GeV.
- For low masses (15–60 GeV), the scalar peak in the recoil-mass distribution is well separated from backgrounds, which is why the strongest sensitivity appears at the low end of the mass range.
Where Pith is reading between the lines
- Editorially: the limit curves, expressed as sigma x B, can be rescaled by any model's production cross-section and branching fraction, so the sensitivity numbers carry over to dark-photon, axion-like, or hidden-sector scenarios with the same Z-associated topology—even though the paper's benchmark is a singlet-plus-dark-matter toy model.
- Editorially: a natural extension is to search for the scalar's visible decays using the same recoil-mass spectrum; the present invisible analysis's strong low-mass reach suggests an analogous visible search would extend the programme without needing additional luminosity.
- Editorially: FCC-ee's proposed high-luminosity run at the Z pole would multiply Z-associated production rates, so recasting this analysis there could push the cross-section reach below 10^-2 fb for low-mass scalars.
- Editorially: the sensitivity dips near the Z and Higgs masses arise from ZZ and ZH backgrounds; a future analysis with a finer background decomposition or alternative jet definitions might recover some of the lost sensitivity, but that is a suggestion beyond the paper's claims.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a simulation-based projection for the process e+e− → Z(→jj) h1, with h1 decaying invisibly, at the FCC-ee (√s = 240 GeV, L = 10.8 ab−1). The authors introduce a simplified model of a real scalar singlet mixed with the SM Higgs plus a scalar dark-matter candidate, generate signal and background events with MadGraph/Pythia/Delphes using the IDEA detector parametrization, and analyze them with a cut-based strategy and an XGBoost BDT. Expected significances are computed with the Asimov formula and expected 95% CL upper limits on σ(e+e−→Zh1) × B(h1→invisible) are derived with CMS Combine. The main claims are: sensitivities of ~10−2–10−1 fb below the Z mass, 0.1–1 fb in the 80–120 GeV range, and discovery reach for scalar masses up to 80 GeV depending on the mixing angle.
Significance. If the results are correct, the paper would provide a useful, technically detailed projection for a key FCC-ee search channel. The analysis chain is standard and largely reproducible: the UFO model is publicly available, the event generation uses well-established tools, the BDT is carefully described, and the limit-setting uses the CMS Combine machinery. The treatment of h1–h2 interference is a positive feature. There is no fitted-to-target circularity. However, the benchmark model on which the signal normalization and the discovery-reach claim rest is internally inconsistent with the stated mixing convention. This undermines the central quantitative claims as currently presented, and the lack of any LEP/OPAL constraint overlay leaves the statement that 'existing experimental constraints do not yet exclude these states' unsupported.
major comments (3)
- [II, Eq. (3), Table I] The benchmark is internally inconsistent with the mixing convention. Equation (3) gives h2 = −sinα φS + cosα φ, so with sinα = 0.985, cosα ≈ 0.173 and h2 has a doublet component of only cos^2α ≈ 0.03. Its SM-like width would then be ≈ 0.03 × 4.1 MeV ≈ 0.12 MeV. Adding the BSM contribution implied by Table I (B_cascade ≈ 7.5% of Γ_h2 = 6.85 MeV, i.e. ≈ 0.51 MeV) gives a total width ≈ 0.6 MeV, not 6.85 MeV. Conversely, to obtain Γ_h2 ≈ 6.85 MeV with only 7.5% BSM width requires SM-like partial widths summing to ≈ 6.3 MeV, which is impossible for cos^2α = 0.03. Similarly, with sinα = 0.985, h1 is mostly doublet (component ≈ 0.985), so σ(e+e−→Zh1) at √s = 240 GeV should be of order 0.97 × σ_SM(ZH) × BR(Z→jj), i.e. roughly 100 fb, whereas Table I reports 13–20 fb. The quoted cross-sections are of the size expected for a mostly-singlet h1 (doublet component ≈ 0.17). Thus the parameter input an
- [II, §VI, abstract] The paper does not apply LEP/OPAL constraints, despite citing OPAL [5] and stating in §II that 'we do not impose all constraints'. With the stated sinα = 0.985, h1 is mostly doublet, so LEP e+e−→Zh1 searches, including invisible decays, would be directly relevant and likely exclude most of the 15–80 GeV mass range where the paper claims discovery reach. The abstract statement that 'existing experimental constraints do not yet exclude these states' is therefore not established. The authors should either overlay the relevant LEP limits on the benchmark points or demonstrate explicitly that the surviving parameter space still yields the quoted cross-sections and significances.
- [IV, Fig. 6] The significance curves in Fig. 6 are computed using the cross-sections from Table I, which are the same cross-sections that are inconsistent with the mixing convention in Eq. (3). Even if the kinematic shapes used for the BDT are approximately model-independent, the absolute signal normalization is not. After correcting the mixing convention, the signal yields, and therefore the significances and the 'up to 80 GeV discovery' conclusion, will change. The model-independent upper-limit part of the paper may survive, but the discovery-reach claim needs to be re-derived with a physically consistent benchmark.
minor comments (3)
- [Table II] The units in Table II appear inconsistent for the rows labeled 'ννH(→ ...)': for example, 'e+e−→ννH(→bb)' is listed as 26.7 pb, which is three orders of magnitude larger than the expected ZH production with Z→νν and H→bb (of order 20–30 fb). The numbers appear to be in fb rather than pb. Please correct the table header or the values, and verify that the analysis uses the correct normalizations.
- [Fig. 2] The caption for Fig. 2 refers to 'the simulation for e+e−→Z(jj)h1(χχ) with sinα = 0.985' in black dotted lines, while the text says this is obtained using h1 mediation only. Please clarify whether this curve is generated with the full interference or with only the h1 propagator, since the two are compared.
- [IV, Eq. (12)] For M_S = 25 GeV the pmiss threshold in Eq. (12) is 10 GeV, which is the same as the preselection cut in Table III. This makes the additional cut redundant for that mass point. Check whether the threshold is intentionally 10 GeV or should be a slightly larger value.
Circularity Check
No significant circularity: the expected limits are extracted from simulated signal/background samples via a profile-likelihood fit; benchmark inputs are not fitted to the claimed limit.
full rationale
The central result is an expected-upper-limit projection, not a parameter extraction. Signal events are generated from a toy model (Eq. 1) with fixed inputs (sinα=0.985, Mχ=5 GeV, M2=125 GeV, v=246 GeV), with benchmark points in Table I selected to satisfy external constraints (B(h2→invisible)<10%, Γ_h2∈[4,8] MeV, perturbativity |λi|≤4π). The 95% CL limits are then obtained with the CMS Combine AsymptoticLimits machinery from MVA score distributions, with signal normalized to a reference value, so the quoted σ×B upper limits are statistically extracted rather than equal to an input. The 'discovery reach' statements are model-dependent projections, which is standard and not circular. Self-citations (e.g., [17], [26]) motivate the toy model and one benchmark value, but the analysis is not forced by them: the mixing angle is varied (sinα=0.973, 0.985, 0.995) and the limits are shown to be largely independent of it. The paper explicitly notes in Section II that 'we do not impose all constraints that would in principle be viable for a UV complete theory'; this admitted limitation, and any internal inconsistency in the benchmark, are correctness/external-constraint concerns rather than circularity—no equation is defined in terms of the target and no fitted parameter is renamed as a prediction.
Axiom & Free-Parameter Ledger
free parameters (5)
- Mixing angle sin α =
0.985 (varied 0.973, 0.995)
- Dark matter mass Mχ =
5 GeV
- Portal coupling λχ =
per-benchmark values, e.g. -5.0e-4 for 15-55 GeV and -2.5e-3 for 65-120 GeV (Table I)
- Singlet VEV vS =
per-benchmark values, e.g. 2936-8514 GeV (Table I)
- pmiss cut thresholds =
10-36 GeV depending on MS (Eq. 12)
axioms (5)
- domain assumption Eq. (1) scalar potential with a softly broken Z2 and the specified mixing matrix is a physical realization of a light scalar plus dark matter model.
- ad hoc to paper h2 in Eq. (3) can be identified with the observed 125 GeV Higgs while sin α = 0.985.
- domain assumption Existing experimental constraints (OPAL, LEP, LHC) do not exclude the scanned benchmark points.
- domain assumption The 10% background systematic and 1% luminosity uncertainties are representative.
- domain assumption Narrow-width approximation with Γ_h1/M1 < 15%.
invented entities (2)
-
Light scalar singlet h1
no independent evidence
-
Scalar dark matter candidate χ
no independent evidence
read the original abstract
We investigate the production of invisibly decaying light scalars in association with hadronically decaying $Z$ bosons at the Future Circular Collider-ee at a centre-of-mass energy $\sqrt{s}=240$ GeV. Several new physics models predict the existence of these low-mass scalar states, while the existing experimental constraints do not yet exclude these states. We study the low-mass scalar based on a simplified extension of the Standard Model, introducing an additional scalar singlet and a scalar dark matter candidate. The analysis is performed for a set of new scalars with mass in the range $(15, 120)$ GeV, by employing a selection-based strategy complemented with Multivariate Analysis techniques to discriminate the signal from background. The expected upper limits on the production cross-section times the branching fraction of the new scalars decaying invisibly are evaluated as a function of the scalar mass. We find that sensitivities of $\sim 10^{-2}$--$10^{-1}$~fb are achievable for scalar masses below the $Z$ boson mass, while sensitivities of $0.1$--$1$~fb are obtained in the mass range 80--120 \GeV. Depending on the mixing angle, novel scalars with masses up to 80 \GeV are within discovery reach.
Figures
Reference graph
Works this paper leans on
-
[4]
T. Robens and R. Santos,BSM: Extended Scalar Sectors, 7, 2025.2507.21910. [5]OPALcollaboration, G. Abbiendi et al.,Decay-Mode Independent Searches for New Scalar Bosons with the OPAL Detector at LEP,Eur. Phys. J. C27(2003) 311–329, [hep-ex/0206022]
arXiv 2025
-
[6]
G. Ripellino, M. Vande Voorde, A. Gallén and R. Gonzalez Suarez,Searching for Long-Lived Dark Scalars at the FCC-ee,JHEP06(2025) 143, [2412.10141]
Pith/arXiv arXiv 2025
-
[7]
A. Bal, E. Curtis, A.-M. Magnan, B. Maier, T. Robens and N. Wardle,Search for additional scalar bosons within the Inert Doublet Model in a final state with two leptons at the FCC-ee, Eur. Phys. J. C85(2025) 891, [2504.12178]
arXiv 2025
-
[8]
J. Kalinowski, W. Kotlarski, T. Robens, D. Sokolowska and A. F. Zarnecki,Exploring Inert Scalars at CLIC,JHEP07(2019) 053, [1811.06952]. 17 [9]CLICdpcollaboration, K. Mekala, A. F. Zarnecki, B. Grzadkowski and M. Iglicki,Sensitivity to Invisible Scalar Decays at CLIC,Eur. Phys. J. Plus136(2021) 160, [2002.06034]. [10]CLICdpcollaboration, J. Klamka and A. ...
Pith/arXiv arXiv 2019
-
[11]
Y. Wang, M. Berggren and J. List,ILD Benchmark: Search for Extra Scalars Produced in Association with aZBoson at √s= 500GeV,2005.06265. [12]ILD Concept Groupcollaboration, M. Berggren, M. T. N. P. de Vera, B. Brudnowski, K. Zembaczyński and A. F. Żarnecki,Search for Exotic Scalars at the International Linear Collider,PoSICHEP2024(2025) 080
Pith/arXiv arXiv 2005
-
[13]
B. Brudnowski, K. Zembaczyński and A. F. Żarnecki,Prospects for light exotic scalar measurements at the e+e−Higgs factory,EPJ Web Conf.315(2024) 01026, [2409.19761]
Pith/arXiv arXiv 2024
-
[14]
B. Brudnowski, K. Zembaczyński and A. F. Żarnecki,Searches for light exotic scalar decays at the e+e− Higgs factory,Eur. Phys. J. Plus141(2026) 392, [2604.22572]
Pith/arXiv arXiv 2026
-
[15]
T. Robens,Low mass scalars ate+e− colliders, inInternational Workshop on Future Linear Colliders, 4, 2026.2604.27541
arXiv 2026
-
[16]
de Blas et al.,Focus topics for the ECFA study on Higgs / Top / EW factories,2401.07564
J. de Blas et al.,Focus topics for the ECFA study on Higgs / Top / EW factories,2401.07564
-
[17]
Altmann et al.,ECFA Higgs, electroweak, and top Factory Study, vol
J. Altmann et al.,ECFA Higgs, electroweak, and top Factory Study, vol. 5/2025 ofCERN Yellow Reports: Monographs. 6, 2025, 10.23731/CYRM-2025-005. [18]FCCcollaboration, M. Benedikt et al.,Future Circular Collider Feasibility Study Report: Volume 1, Physics, Experiments, Detectors,Eur. Phys. J. C85(2025) 1468, [2505.00272]. [19]FCCcollaboration, M. Benedikt...
Pith/arXiv arXiv 2025
-
[21]
Robens,Interference effects in new physics searches,2602.00256
T. Robens,Interference effects in new physics searches,2602.00256
-
[22]
IDEA—FCC-ee Winter 2023 Monte-Carlo Production
FCC Collaboration, “IDEA—FCC-ee Winter 2023 Monte-Carlo Production.” https://fcc-physics-events.web.cern.ch/fcc-ee/rec/winter2023/IDEA, 2023
2023
-
[23]
T. Robens and T. Stefaniak,Status of the Higgs Singlet Extension of the Standard Model after LHC Run 1,Eur. Phys. J. C75(2015) 104, [1501.02234]
Pith/arXiv arXiv 2015
-
[24]
T. Robens and T. Stefaniak,LHC Benchmark Scenarios for the Real Higgs Singlet Extension of the Standard Model,Eur. Phys. J. C76(2016) 268, [1601.07880]
Pith/arXiv arXiv 2016
-
[25]
A. Ilnicka, T. Robens and T. Stefaniak,Constraining Extended Scalar Sectors at the LHC and beyond,Mod. Phys. Lett. A33(2018) 1830007, [1803.03594]
Pith/arXiv arXiv 2018
-
[26]
T. Robens,Overview on Low Mass Scalars ate+e− Facilities – Theory, inInternational Workshop on Future Linear Colliders, 2024.2409.19657. [27]ATLAScollaboration, G. Aad et al.,Combination of searches for invisible decays of the Higgs boson using 139 fb−1 of proton-proton collision data at s=13 TeV collected with the ATLAS experiment,Phys. Lett. B842(2023) ...
arXiv 2024
-
[29]
N. D. Christensen and C. Duhr,FeynRules - Feynman rules made easy,Comput. Phys. Commun.180(2009) 1614–1641, [0806.4194]
Pith/arXiv arXiv 2009
-
[30]
Darmé et al.,UFO 2.0: the ‘Universal Feynman Output’ format,Eur
L. Darmé et al.,UFO 2.0: the ‘Universal Feynman Output’ format,Eur. Phys. J. C83(2023) 631, [2304.09883]
Pith/arXiv arXiv 2023
-
[31]
J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer et al.,The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations,JHEP07(2014) 079, [1405.0301]. 18
Pith/arXiv arXiv 2014
-
[32]
Bierlich et al.,A Comprehensive Guide to the Physics and Usage of PYTHIA 8.3,SciPost Phys
C. Bierlich et al.,A Comprehensive Guide to the Physics and Usage of PYTHIA 8.3,SciPost Phys. Codeb.2022(2022) 8, [2203.11601]. [33]IDEA Study Groupcollaboration, M. Abbrescia et al.,The IDEA Detector Concept for FCC-ee,2502.21223. [34]DELPHES 3collaboration, J. de Favereau, C. Delaere, P. Demin, A. Giammanco, V. Lemaître, A. Mertens et al.,DELPHES 3: A M...
Pith/arXiv arXiv 2022
-
[35]
W. Kilian, T. Ohl and J. Reuter,WHIZARD: Simulating Multi-Particle Processes at LHC and ILC,Eur. Phys. J. C71(2011) 1742, [0708.4233]
Pith/arXiv arXiv 2011
-
[36]
T. Sjostrand, S. Mrenna and P. Z. Skands,PYTHIA 6.4 Physics and Manual,JHEP05(2006) 026, [hep-ph/0603175]. [37]Key4hepcollaboration, A. Sailer et al.,The Key4hep software stack: Beyond Future Higgs factories,J. Phys. Conf. Ser.3206(2026) 012043, [2312.08151]
Pith/arXiv arXiv 2006
-
[38]
S. D. Ellis and D. E. Soper,Successive Combination Jet Algorithm for Hadron Collisions,Phys. Rev. D48(1993) 3160–3166, [hep-ph/9305266]
Pith/arXiv arXiv 1993
-
[39]
Catani, Y
S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock and B. R. Webber,New Clustering Algorithm for Multi-Jet Cross-Sections ine+e− Annihilation,Phys. Lett. B269(1991) 432–438
1991
-
[40]
Catani, Y
S. Catani, Y. L. Dokshitzer, M. H. Seymour and B. R. Webber,Longitudinally Invariantkt Clustering Algorithms for Hadron–Hadron Collisions,Nucl. Phys. B406(1993) 187–224
1993
-
[41]
M. Cacciari, G. P. Salam and G. Soyez,FastJet User Manual,Eur. Phys. J. C72(2012) 1896, [1111.6097]
Pith/arXiv arXiv 2012
-
[42]
C. Helsens, E. Perez, M. Selvaggi, V. Volkl, L. Forthomme and J. Munch Torndal, HEP-FCC/FCCAnalyses: v0.11.0, may, 2025. 10.5281/zenodo.15528870
-
[43]
T. Chen and C. Guestrin,XGBoost: A Scalable Tree Boosting System, inProceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, pp. 785–794, ACM, 2016.1603.02754. DOI
arXiv 2016
-
[44]
H. Voss, A. Hocker, J. Stelzer and F. Tegenfeldt,TMVA, the Toolkit for Multivariate Data Analysis with ROOT,PoSACA T(2007) 040
2007
-
[45]
Brun and F
R. Brun and F. Rademakers,ROOT: An Object-Oriented Data Analysis Framework,Nucl. Instrum. Meth. A389(1997) 81–86
1997
-
[46]
G. Cowan, K. Cranmer, E. Gross and O. Vitells,Asymptotic formulae for likelihood-based tests of new physics,Eur. Phys. J. C71(2011) 1554, [1007.1727]. [47]CMScollaboration, A. Hayrapetyan et al.,The CMS Statistical Analysis and Combination Tool: Combine,Comput. Softw. Big Sci.8(2024) 19, [2404.06614]
Pith/arXiv arXiv 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.