REVIEW 3 major objections 3 minor 42 references
On the contribution from the light quarks to $H\to\gamma\gamma , \gamma Z$
T0 review · 3 major / 3 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The light-quark contribution to H→γγ does not vanish: it starts linearly, not quadratically, in the light-quark masses, and shifts the predicted decay rate by about three percent.
desk verdict The qualitative SVV argument is right and refutes the 2025 claim, but the LMD numbers are an estimate and the 'exact' language oversells them. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the scalar-vector-vector three-point function $\Gamma^{abc}_{\mu\nu}(q_1,q_2)=\int d^4x\,d^4y\,e^{iq_1\cdot x}e^{iq_2\cdot y}\langle 0|T\{S^a(0)V^b_\mu(x)V^c_\nu(y)\}|0\rangle$, built from colour-singlet quark bilinears in three-flavour QCD, whose flavour structure matches the singlet and nonsinglet components of the Higgs and vector currents. In the chiral limit this correlator is an order parameter of chiral symmetry breaking, and the short-distance expansions (3.4) and (3.6) show a leading behaviour proportional to the quark condensate, which proves it is non-vanishing. To reach the physical decay point the paper uses the lowest-meson-dominance (LMD) ansatz $F_{\mathrm{LMD}}(q_1^2,q_2^2,q_3^2)=[a+b(q_1^2+q_2^2)+c\,q_3^2]/[(q_1^2-M_V^2)(q_2^2-M_V^2)(q_3^2-M_S^2)]$, a single-pole representation of the form factor. Matching to the operator-product expansion fixes $b=c=-\tfrac12\langle\bar q q\rangle_0$, and a sum rule over $e^+e^-$ widths fixes $a$; the resulting expression (4.5) interpolates from the deep-Euclidean regime to the physical kinematics $q_1^2=q_2^2=0$, $q_3^2=M_H^2$ (and $q_2^2=M_Z^2$ for the $Z$ channel), and it is this interpolation that produces the numerical enhancement.
What would settle it
A direct test would be a lattice QCD computation of the scalar-vector-vector three-point function with physical quark masses, including the flavour-singlet (disconnected) quark contractions, extrapolated to the Higgs kinematics $q_1^2=q_2^2=0$, $q_3^2=M_H^2$ and to $q_2^2=M_Z^2$ for the $Z$ channel: if the form factor came out zero at those points, or if the amplitude scaled quadratically with $m_q$ rather than linearly, the paper's central claim would be false. A complementary experimental check is a percent-level measurement of the $H\to\gamma\gamma$ rate: agreement with the standard perturbative prediction to better than the predicted 3% downward shift would rule out the LMD enhancement.
Extended reading notes
Core claim
The central claim is that the light-quark contribution to the $H\to\gamma\gamma$ and $H\to\gamma Z$ amplitudes does not vanish and is linear in the light-quark masses $m_u$, $m_d$, $m_s$, not quadratic as in the usual perturbative treatment. The object that carries the argument is the $\langle SVV\rangle$ three-point function, the correlator of a colour-singlet scalar quark density with two vector currents, which in the chiral limit (all quark masses set to zero) is an order parameter of spontaneous chiral symmetry breaking: the operator-product expansion gives a leading short-distance behaviour proportional to the quark condensate $\langle\bar q q\rangle_0$, so the correlator cannot be identically zero even though it vanishes at every order of perturbation theory. The decay amplitudes factor a single power of the light-quark mass times this non-vanishing hadronic function, and that is the origin of the linear-mass behaviour. Evaluated in the lowest-meson-dominance approximation to the $N_c\to\infty$ limit of QCD, with parameters fixed by the short-distance constraints, a sum rule, and lattice inputs, the summed light-quark effective couplings for $H\to\gamma\gamma$ come out larger by about three orders of magnitude than the one-loop perturbative values, giving $\Gamma(H\to\gamma\gamma)_{\mathrm{LMD}}=0.97\,\Gamma(H\to\gamma\gamma)_{\mathrm{pQCD}}$, while the $H\to\gamma Z$ contribution is barely modified.
Load-bearing premise
The result assumes that the lowest-meson-dominance form factor, whose parameters are fixed by the short-distance expansion and a sum rule, faithfully represents the true non-perturbative correlator at the physical decay point — two on-shell photons (or a photon and a $Z$) with the Higgs on shell; if the true form factor vanished or behaved very differently at those kinematics, the predicted enhancement, and possibly the linear-mass conclusion itself, would fail.
Editorial extensions
If this is right
- If the claim holds, the light-quark contribution to $H\to\gamma\gamma$ is roughly a thousand times larger than the one-loop perturbative estimate, moving the predicted rate to $\Gamma(H\to\gamma\gamma)_{\mathrm{LMD}}=0.97\,\Gamma(H\to\gamma\gamma)_{\mathrm{pQCD}}$; this shift becomes visible once branching-ratio precision reaches the percent level.
- The $H\to\gamma Z$ amplitude gains no comparable enhancement: the light-quark couplings grow only modestly and the rate prediction is essentially unchanged, so this channel will not discriminate the effect.
- The non-perturbative treatment removes the renormalization-scheme and scale ambiguity of the perturbative light-quark contribution, because the final couplings depend only on ratios of quark masses and on mass-times-condensate products, which are renormalization-group invariant.
- Within the large-$N_c$/LMD description the light-quark contribution is purely real; generating imaginary parts would require subleading $1/N_c$ effects from multi-particle intermediate states, which the paper leaves for future work.
Reading between the lines
- A lattice QCD evaluation of the $\langle SVV\rangle$ form factor at the exact Higgs kinematics, or at small momenta with controlled extrapolation, would settle the dispute directly; the paper uses lattice inputs for masses and condensates but does not propose this computation, and varying the light-quark mass on the lattice would test the predicted linear scaling.
- Because the enhancement is driven by the external scale $M_H$ (the $q_3^2$ term in the LMD numerator), the same linear-mass mechanism should apply to heavier CP-even scalars in extensions of the standard model: their light-quark-induced $\gamma\gamma$ corrections would not decouple like a perturbative $m_q^2$ term.
- The same reasoning applies to any amplitude built from a scalar density and two vector currents in QCD, so other precision observables sensitive to such matrix elements may hide similar non-perturbative linear-quark-mass effects that purely perturbative estimates miss.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the contribution of the light quarks (u, d, s) to the Standard Model amplitudes for H→γγ and H→γZ, using non-perturbative QCD. The authors argue that the relevant object is the ⟨SVV⟩ three-point function, which they show via short-distance OPE (Eqs. (3.4), (3.6)) does not vanish identically in the chiral limit. They then evaluate this correlator in the lowest-meson-dominance (LMD) approximation to the large-Nc limit, obtaining a form factor (Eq. (4.5)) whose parameters are fixed by the OPE, a sum rule, and external inputs. The resulting light-quark contribution is linear in the quark masses, not quadratic, and for H→γγ is larger in magnitude by about three orders than the lowest-order perturbative evaluation, changing the total H→γγ rate by about 3% (Eq. (5.9)). For H→γZ the effect is small. The paper claims that the non-vanishing and linear-mass behavior are exact consequences of QCD short-distance properties.
Significance. If the quantitative results held, the paper would overturn the usual neglect of light-quark contributions in Higgs decays and correct a recent claim (ref. [9]) that these contributions vanish at the non-perturbative level. The identification of the ⟨SVV⟩ correlator as the correct object and the demonstration that it is not identically zero are solid and important. The paper is careful to use external data (e+e- widths, lattice quark condensates, PDG/FLAG masses) rather than fitting to the target amplitudes, so there is no circularity. The predicted few-percent effect on H→γγ is potentially relevant for future precision measurements at HL-LHC and FCC-ee. However, the central quantitative claim relies on the LMD model, whose validity at the physical kinematics is the main weakness; the exactness claimed in the abstract and Section 6 goes beyond what the OPE alone establishes.
major comments (3)
- [Abstract and Section 6] The statements that the light-quark contribution is linear in the quark masses at the physical point are presented as exact ('These statements are exact since they rest only on the identification of the ⟨SVV⟩ correlators ... and on their short-distance properties in QCD'). This overreaches the evidence. The OPE results (3.4) and (3.6) constrain the correlator only in deep-Euclidean or mixed deep-Euclidean regimes; they do not determine the form factor F at q1^2=q2^2=0, q3^2=M_H^2 (or q2^2=M_Z^2), where two of the momenta are on shell and q3 is timelike. A correlator that is nonzero in the Euclidean region could in principle vanish at these specific physical kinematics. The linear-in-mass behavior of the physical amplitudes therefore rests on the LMD ansatz of Eq. (4.2), not on an exact QCD statement. Please revise the abstract and Section 6 to clearly distinguish the exact non-vanishing of the SVV correlator from the model-dependent evaluation at the physical point.
- [Section 4, Eqs. (4.2)-(4.5)] The LMD form factor retains only the lowest vector and scalar poles. In the large-Nc limit, radial excitations in each channel are not parametrically suppressed (each meson contributes at order 1 in 1/Nc). Their contribution to FLMD at q3^2=M_H^2, where the scalar pole is evaluated far from its mass, is not estimated. Since the numerical enhancement in Table 1 and Eq. (5.9) is dominated by the OPE-determined b and c terms evaluated at the physical point, the quantitative results should be framed as an estimate within a specific model, with an estimate of the truncation uncertainty or a justification for why higher resonances are negligible at these kinematics. As written, the paper's quantitative predictions (e.g., the three-order enhancement and Γ_LMD=0.97Γ_pQCD) are not established beyond the LMD ansatz.
- [Section 5, Table 1] The comparison between the LMD and pQCD columns mixes amplitudes with different analytic structure: the LMD results are purely real because of the pole approximation, while the pQCD results are complex because of physical thresholds. For the u and d quarks the imaginary part in pQCD is comparable to or larger than the real part. The statement of an 'increase ... by at least three orders of magnitude' should specify whether it refers to |C| or Re C, and the interpretation of the enhancement should account for the fact that the two evaluations belong to different approximations (large-Nc poles versus perturbative cuts). Please clarify this point and discuss whether a direct comparison of absolute values is meaningful.
minor comments (3)
- [Section 5, Eq. (5.1) and footnote] The quoted value c~ = 4.6(0.8)·10^-3 appears inconsistent with the footnote stating that the central value is '0.52 instead of 0.46'; please clarify the units and whether the table uses the updated value.
- [Throughout] There are several typographical errors: 'runnung' before Eq. (5.3), 'regine' in Section 2, and 'an the one hand' in Section 6. The title also renders 'H→γγ, γZ' with a missing space in 'toH → γγ, γZ' in the header.
- [Section 3, Eq. (3.4)] The OPE formula (3.4) would benefit from an explicit statement of the kinematic conditions (all q_i^2 large and Euclidean) and the order of the neglected terms; the current text says 'deep-Euclidean' but does not specify how large the momenta must be for the leading term to dominate.
Circularity Check
No significant circularity: the central non-vanishing and linear-mass claims rest on OPE and chiral-symmetry arguments, and the LMD parameters are fixed by external inputs rather than by the target H decay amplitudes.
full rationale
The derivation chain is: (i) identify the light-quark contribution with the <SVV> three-point function; (ii) use the OPE to show the correlator is nonzero in the chiral limit; (iii) approximate the form factor with the LMD ansatz whose parameters are fixed by short-distance constraints, a sum rule, and external hadronic inputs; (iv) evaluate numerically. The non-vanishing and linear-in-mass claims do not reduce to a fitted parameter: the quark-mass factor is factored out explicitly in eqs. (2.5) and (2.6), and the form factors are evaluated in the chiral limit, so linearity follows if the form factor does not vanish, which is supported by the OPE results (3.4) and (3.6) with the quark condensate as the leading term. The LMD parameters are not fitted to H->gamma-gamma or H->gamma-Z amplitudes: c_tilde comes from e+e- partial widths (eq. 5.1), m_q<q_bar q>_0 from lattice QCD (eq. 5.2), quark masses from PDG/FLAG (eq. 5.3), and b and c from the OPE constraints (4.3). The self-cited sum rule (4.4) fixes only parameter a, and the paper itself notes that a is numerically subdominant: the numerator of FLMD is dominated by the term proportional to (q^2+M_H^2), larger by three to four orders of magnitude than the term proportional to c_tilde. Thus no load-bearing self-citation chain forces the result. The LMD ansatz is openly an approximation, and Section 6 explicitly qualifies the quantitative conclusions ('should certainly be taken cum grano salis'), so the model dependence is acknowledged rather than disguised as a derivation. The claim in Section 6 that the non-vanishing/linear statements are 'exact' may overstate the support for the on-shell value, but that is a correctness or model-extrapolation concern, not circularity. No equation was found in which a fitted input is renamed a prediction or in which the target quantity is defined into the input.
Assumptions & free parameters
free parameters (3)
- M_V, lowest vector meson mass =
770 MeV
- M_S, lowest scalar meson mass =
1400 MeV
- Tilde c, sum-rule constant =
4.6(0.8) x 10^-3
assumptions (4)
- domain assumption Chiral limit: the three-point functions are evaluated with mu = md = ms = 0 after factoring out the quark masses
- domain assumption Spontaneous chiral symmetry breaking occurs, so the quark condensate <qqbar>_0 is nonzero
- ad hoc to paper The deep-Euclidean OPE results (3.4) and (3.6) control the form factor at physical kinematics
- domain assumption Large-Nc limit with lowest-meson dominance (single-pole ansatz) describes the SVV form factor
Cite this review
Pith. "Pith review of On the contribution from the light quarks to $H\to\gamma\gamma , \gamma Z$." pith.science (2026). https://pith.science/paper/JP4TD2HF
@misc{pith2026250722551,
author = {Pith},
title = {Pith review of: On the contribution from the light quarks to $H\to\gamma\gamma , \gamma Z$},
year = {2026},
howpublished = {\url{https://pith.science/paper/JP4TD2HF}},
note = {Machine review of arXiv:2507.22551}
}
abstract
This Letter addresses the contribution from the light $u , d , s$ quarks to the amplitudes for the decay modes of the Brout-Englert-Higgs scalar boson $H$ into two photons ($H\to\gamma\gamma$) or to a photon and a neutral weak gauge boson ($H\to\gamma Z$), taking into account the non-perturbative aspects of QCD. Contrary to a recent claim, the contribution from the light quarks does not vanish. Rather, it is shown that, in contrast to the perturbative evaluations usually considered in the literature, this contribution to the amplitudes starts with a term that is linear, and not quadratic, in the masses of the light quarks, thus pointing toward a sizeable enhancement of their contribution to $H \to \gamma \gamma$.
Reference graph
Works this paper leans on
-
[9]
A. I. Hern´ andez-Ju´ arez, R. Gait´ an and R. Martinez,Light quark contributions to Higgs decays, [arXiv:2502.01668 [hep- ph]]
-
[1]
J. R. Ellis, M. K. Gaillard and D. V. Nanopoulos, Nucl. Phys. B 106, 292 (1976)
work page 1976
-
[2]
M. A. Shifman, A. I. Vainshtein, M. B. Voloshin and V. I. Zakharov, Sov. J. Nucl. Phys. 30, 711-716 (1979)
work page 1979
-
[3]
M. B. Gavela, G. Girardi, C. Malleville and P. Sorba, Nucl. Phys. B 193, 257 (1981)
work page 1981
-
[4]
R. N. Cahn, M. S. Chanowitz and N. Fleishon, Phys. Lett. B 82, 113 (1979)
work page 1979
-
[5]
L. Bergstrom and G. Hulth, Nucl. Phys. B 259, 137 (1985) [erratum: Nucl. Phys. B 276, 744 (1986)]
work page 1985
-
[6]
K. H. Phan, L. T. Hue and D. T. Tran, PTEP 2021, no.9, 093B05 (2021) [arXiv:2103.14248 [hep-ph]]
work page Pith review arXiv 2021
-
[7]
M. Aiko, J. Braathen and S. Kanemura, Eur. Phys. J. C 85, 489 (2025) [arXiv:2307.14976 [hep-ph]]
arXiv 2025
Show all 42 references
- [8]
-
[10]
Okubo, Phys
S. Okubo, Phys. Lett. 4, 14 (1963); G. Zweig, CERN Report No. 8419/TH412 (1964); J. Iizuka, Prog. Theor. Phys. Suppl. 37, 21 (1966)
1963
-
[11]
Spira, A
M. Spira, A. Djouadi and P. M. Zerwas, Phys. Lett. B 276, 350 (1992)
1992
-
[12]
Gehrmann, S
T. Gehrmann, S. Guns and D. Kara, JHEP 09, 038 (2015) [arXiv:1505.00561 [hep-ph]]
2015 arXiv
-
[13]
Bonciani, V
R. Bonciani, V. Del Duca, H. Frellesvig, J. M. Henn, F. Moriello and V. A. Smirnov, JHEP 08, 108 (2015) [arXiv:1505.00567 [hep-ph]]. 4The ⟨P V V⟩ three-point function would become relevant in cases, like two-Higgs doublet models, where the scalar sector contains also a J P = 0...
2015 arXiv
-
[14]
’t Hooft, Nucl
G. ’t Hooft, Nucl. Phys. B 72, 461 (1974)
1974
-
[15]
Witten, Nucl
E. Witten, Nucl. Phys. B 160, 57 (1979)
1979
-
[16]
Moussallam and J
B. Moussallam and J. Stern, Chiral symmetry aspects of the scalars, arXiv:hep-ph/9404353 [hep-ph]
-
[17]
Kadav´ y, K
T. Kadav´ y, K. Kampf and J. Novotn´ y, JHEP10, 142 (2020) [arXiv:2006.13006 [hep-ph]]
2020 arXiv
-
[18]
B. L. Ioffe and A. V. Smilga, Nucl. Phys. B 232, 109 (1984)
1984
-
[19]
G. S. Bali, F. Bruckmann, M. Constantinou, M. Costa, G. Endrodi, S. D. Katz, H. Panagopoulos and A. Schafer, Phys. Rev. D 86, 094512 (2012) [arXiv:1209.6015 [hep-lat]]
2012 arXiv
-
[20]
P. Ball, V. M. Braun and N. Kivel, Nucl. Phys. B 649, 263 (2003) [arXiv:hep-ph/0207307 [hep-ph]]
2003 arXiv
-
[21]
Peris, M
S. Peris, M. Perrottet and E. de Rafael, JHEP 05 (1998), 011 [arXiv:hep-ph/9805442 [hep-ph]]
1998 arXiv
-
[22]
M. F. L. Golterman and S. Peris, Phys. Rev. D 61 (2000), 034018 [arXiv:hep-ph/9908252 [hep-ph]]
2000 arXiv
-
[23]
Knecht, S
M. Knecht, S. Peris, M. Perrottet and E. de Rafael, Phys. Rev. Lett. 83, 5230 (1999) [arXiv:hep-ph/9908283 [hep-ph]]
1999 arXiv
-
[24]
Knecht, B
M. Knecht, B. Moussallam and J. Stern, Nucl. Phys. B 429, 125 (1994) [arXiv:hep-ph/9402318 [hep-ph]]
1994 arXiv
-
[25]
Knecht and A
M. Knecht and A. Nyffeler, Eur. Phys. J. C 21, 659 (2001) [arXiv:hep-ph/0106034 [hep-ph]]
2001 arXiv
-
[26]
Navas et al
S. Navas et al. [Particle Data Group], Phys. Rev. D 110, 030001 (2024)
2024
-
[27]
J. A. Oller and E. Oset, Phys. Rev. D 60, 074023 (1999) [arXiv:hep-ph/9809337 [hep-ph]]
1999 arXiv
-
[28]
Cirigliano, G
V. Cirigliano, G. Ecker, H. Neufeld and A. Pich, JHEP 06, 012 (2003) [arXiv:hep-ph/0305311 [hep-ph]]
2003 arXiv
-
[29]
Bijnens and G
J. Bijnens and G. Ecker, Ann. Rev. Nucl. Part. Sci. 64, 149 (2014) [arXiv:1405.6488 [hep-ph]]
2014 arXiv
-
[30]
Knecht, S
M. Knecht, S. Narison, A. Rabemananjara and D. Rabetiarivony, Phys. Lett. B 787, 111 (2018) [arXiv:1808.03848 [hep-ph]]
2018 arXiv
-
[31]
Aoki et al
Y. Aoki et al. [Flavour Lattice Averaging Group (FLAG)], Eur. Phys. J. C 82, 869 (2022) [arXiv:2111.09849 [hep-lat]]
2022 arXiv
- [32]
-
[33]
McNeile, C
C. McNeile, C. T. H. Davies, E. Follana, K. Hornbostel and G. P. Lepage, Phys. Rev. D 82, 034512 (2010) [arXiv:1004.4285 [hep-lat]]
2010 arXiv
-
[34]
Durr et al
S. Durr et al. [BMW], JHEP 08, 148 (2011) doi:10.1007/JHEP08(2011)148 [arXiv:1011.2711 [hep-lat]]
2011 arXiv
-
[35]
Durr et al
S. Durr et al. [BMW], Phys. Lett. B 701, 265 (2011) [arXiv:1011.2403 [hep-lat]]
2011 arXiv
-
[36]
Bazavov, C
A. Bazavov, C. Bernard, C. DeTar, X. Du, W. Freeman, S. Gottlieb, U. M. Heller, J. E. Hetrick, J. Laiho and L. Levkova, et al. PoS LA TTICE2010, 083 (2010) [arXiv:1011.1792 [hep-lat]]
2010 arXiv
-
[37]
Blum et al
T. Blum et al. [RBC and UKQCD], Phys. Rev. D 93, 074505 (2016) [arXiv:1411.7017 [hep-lat]]
2016 arXiv
-
[38]
Bruno et al
M. Bruno et al. [ALPHA], Eur. Phys. J. C 80, 169 (2020) [arXiv:1911.08025 [hep-lat]]
2020 arXiv
-
[39]
A. I. Hern´ andez-Ju´ arez, R. Gait´ an and R. Martinez, Phys. Rev. D111, 015001 (2025) [arXiv:2405.03094 [hep-ph]]
2025 arXiv
-
[40]
ATLAS Coll., Updated projections for Higgs boson measurements usingH → γγ , Zγ and µµ decays and for the combined measurement of Higgs boson couplings with the ATLAS detector at the HL-LHC, ATL-PHYS-PUB-2025-014 (2025)
2025
-
[41]
Benedikt et al., CERN physics reports, CERN-FCC-PHYS-2025-0002, Geneva, 2025
Future Circular Collider Feasibility Study Report Volume 1: Physics, Experiments, Detector, preprint edition edited by M. Benedikt et al., CERN physics reports, CERN-FCC-PHYS-2025-0002, Geneva, 2025. Available online: https://cds.cern.ch/record/292819
2025
-
[42]
Atti´ eet al
D. Atti´ eet al. [Linear Collider Vision], A Linear Collider Vision for the Future of Particle Physics, arXiv:2503.19983 [hep-ex]. 10
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.