REVIEW 3 major objections 5 minor 3 cited by
Three-loop induced neutrino mass model in a non-invertible symmetry
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper proposes that a non-invertible symmetry FR(6) forbids the tree-level Yukawa coupling that would give neutrinos mass, generating it at one loop through new fermions L' and a singlet S0, so neutrino masses are a three-loop effect.
desk verdict The central mechanism is invalid: the Eq. (18) interaction terms violate both U(1)_Y and the FR(6) charges listed in Table II, so the loop-induced Yukawa and three-loop neutrino mass are never defined. 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 central object is the non-invertible symmetry $F_R(6)$, defined as a $\mathbb{Z}_2$ gauging of $\mathbb{Z}_6$, with fusion rules $\epsilon\otimes\epsilon=1\oplus\sigma$, $\sigma\otimes\sigma=1\oplus\sigma$, $\rho\otimes\rho=1$, $\epsilon\otimes\sigma=\epsilon\oplus\rho$, $\epsilon\otimes\rho=\sigma$, $\sigma\otimes\rho=\epsilon$. The charge assignments in Table II put the new fields ($L'_L$ and $S_0$) in nontrivial representations so that the tree-level operator that would give $y_\eta$ is not invariant. The one-loop diagram with $L'$ and $S_0$, controlled by the couplings $f$, $g$, and the $\mu$ term, then regenerates $y_\eta$, and this same $y_\eta$ feeds into the standard neutrino mass formula of the base model.
What would settle it
Multiply the $F_R(6)$ charges of the fields in each operator of Eqs. (18) and (19) using the fusion rules in Eq. (15) and check whether the identity representation $1$ appears in the product; if any operator's product contains only non-trivial representations (such as $\epsilon$ or $\rho$), then $y_\eta$ is not forbidden at tree level by this symmetry, and the three-loop mechanism is not realized.
Extended reading notes
Core claim
The central claim is that a non-invertible symmetry can do the work that a discrete parity does in the inert-doublet radiative seesaw model: it forbids the Yukawa coupling $y_\eta$ at tree level, but because the symmetry is non-invertible it is dynamically broken at one loop, so $y_\eta$ is generated radiatively. The new fields $L'$ and $S_0$, together with the scalar trilinear coupling $\mu|\eta|^2 S_0$, produce $y_\eta$ through a one-loop diagram, and since the neutrino mass formula of the base model is already one-loop in $y_\eta$, the physical neutrino mass matrix is a genuine three-loop quantity. The same scalar doublet $\eta$ participates both in generating $y_\eta$ and in the neutrino mass loop, so the model has a single inert doublet doing two jobs.
Load-bearing premise
The mechanism works only if the tree-level interaction terms in Eqs. (18) and (19) are genuinely invariant under the $F_R(6)$ charges assigned in Table II; if any of these operators transforms in a non-trivial representation rather than the identity, the symmetry would not forbid $y_\eta$ at tree level and the loop-generation chain would not start.
Editorial extensions
If this is right
- If the mechanism is correct, neutrino masses are predicted to be three-loop suppressed, which naturally explains their smallness without invoking very heavy right-handed neutrinos.
- The model predicts a positive contribution to the muon anomalous magnetic moment (because the photon attaches to charged fermions rather than charged bosons in the new loop), opposite in sign to the minimal inert-doublet version.
- The fermion $N_{R1}$ cannot be the dominant dark matter because the radiatively generated $y_\eta$ is at most $10^{-4}$, far below the required order-one coupling; dark matter must instead be the scalar $S_0$ or the real inert doublet component $\eta_R$.
- The allowed mixing angles $\theta_{12,23,13}$ show distinctive shapes, and the dark matter mass is bounded from above ($m_S \lesssim 800$ GeV for $S_0$ DM, $m_R \simeq 534 \pm 8.5$ GeV for $\eta_R$ DM), giving testable predictions for future direct-detection and collider searches.
- All viable parameter points satisfy the $\Sigma m_\nu \leq 120$ meV bound from cosmology, and the model gives specific upper limits on the effective neutrinoless double-beta decay mass $m_{ee}$.
Reading between the lines
- A direct check of the $F_R(6)$ invariance of the interaction terms in Eqs. (18) and (19), using the fusion rules of Eq. (15), is a necessary consistency test; if any of these operators carries a nontrivial representation, the tree-level forbidding of $y_\eta$ would not hold and the mechanism would be ill-defined.
- The same fusion-rule technique could be applied to other unwanted tree-level couplings (for instance, higher-dimensional lepton-number-violating operators) to push their appearance to even higher loops, generalizing the idea beyond the neutrino sector.
- Because the non-invertible symmetry is dynamically broken at loop level, the model effectively links the smallness of neutrino mass to the loop order of symmetry breaking rather than to a hierarchy in Yukawa couplings, a pattern that could be probed by detecting the new fields $L'$ and $S_0$ at colliders.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a three-loop radiative neutrino mass model based on the Ma model, adding vector-like lepton doublets L' and a singlet scalar S0, and imposing a non-invertible symmetry FR(6) in the class of a Z2 gauging of Z6. The direct neutrino Yukawa coupling y^eta is forbidden at tree level by FR(6), and the authors claim that y^eta is generated at one loop through the couplings f, g, and the scalar cubic term mu, making the neutrino mass a three-loop effect. The rest of the paper performs numerical scans for neutrino oscillation data, lepton flavor violation, muon g-2, and two scalar dark matter candidates, S0 and eta_R, for normal and inverted hierarchies.
Significance. The general idea of using a non-invertible symmetry to forbid a coupling at tree level while allowing it radiatively is timely and potentially interesting. The paper gives explicit fusion rules, a full field-content table, a one-loop formula for the induced Yukawa coupling, and relic density calculations using micrOMEGAs, which is a strength. However, the central charge assignment is internally inconsistent: the very operators on which the loop mechanism depends are not invariant under the symmetry that defines the model. As written, the tree-level Lagrangian already violates FR(6), so the advertised mechanism cannot be constructed and the numerical results do not constrain a valid model.
major comments (3)
- [Sec. II, Table II and Eq. (18)] The tree-level Lagrangian is not invariant under FR(6) with the charges given in Table II. Using the fusion rules in Eq. (15), the coupling g_ab \bar{L}'_{La} N_{Rb} \tilde{\eta} transforms as \sigma\otimes\rho\otimes\sigma = (\sigma\otimes\rho)\otimes\sigma = \epsilon\otimes\sigma = \epsilon\oplus\rho, which does not contain the identity. Similarly, h_ab \bar{L}'^c_{Ra} N_{Rb} \eta transforms as 1\otimes\rho\otimes\sigma = \epsilon, and the vector-like mass term M_{L'} \bar{L}'_L L'_R transforms as \sigma\otimes 1 = \sigma. None of these operators is invariant. In particular, the g term is non-invariant for any choice of the L'_R and S0 charges, because the charges of L'_L = \sigma, N_R = \rho, and \eta = \sigma are fixed by the table. The one-loop formula in Eq. (20) therefore has no basis in the symmetric Lagrangian, and the model's central mechanism cannot be defined.
- [Sec. II, Table II; Eq. (18)] The hypercharge assignments in Table II are inconsistent with the 4-component notation used in Eq. (1). With Y(L'_L) = -1/2 and Y(L'_R) = +1/2, the Dirac mass term M_{L'} \bar{L}'_L L'_R carries hypercharge (+1/2) + (+1/2) = +1, and the coupling f \bar{L}_L L'_R S0 also carries net hypercharge +1; the analogous f' term is likewise not U(1)_Y invariant. If the authors intended two-component Weyl notation, then the table's hypercharges and the explicit bar and charge-conjugation symbols in Eq. (18) must be rewritten consistently. As it stands, the loop diagram behind Eq. (20) is not gauge-invariant.
- [Sec. IV, Eqs. (22)-(23)] The numerical analysis treats f and g as free input couplings and uses the Casas-Ibarra parametrization to reproduce neutrino oscillation data. But because these couplings are not part of the FR(6)-invariant Lagrangian, the scans in Sec. IV describe an effective theory in which the non-invertible symmetry has been removed by hand. The claim that the observed neutrino parameters follow from a three-loop mechanism generated by the symmetric model is therefore not supported by the presented calculation.
minor comments (5)
- [Eq. (20)] The second logarithm in Eq. (20) is written with M_{L'a} inside the argument, while the prefactor sums over the index b; this should be M_{L'b} for consistency.
- [Eq. (13)] The quoted value \Delta a_\mu \simeq (39 \pm 64) \times 10^{-1} has an implausible normalization; a factor of 10^{-11} or 10^{-10} is presumably intended.
- [After Eq. (8)] The sentence defining C_{21}, C_{31}, and C_{21} repeats C_{21} twice and uses 0.1784 versus 0.173648; one of these should presumably be C_{32}.
- [Table II] The FR(6) row in Table II appears to contain only seven entries for eight fields, leaving the charges of S0 and possibly eta ambiguous. Since these charges are the basis of the whole construction, the table must be corrected to assign a charge to every field explicitly.
- [Sec. II, Eq. (18)] The notation for the L' fields should be clarified: if L'_L and L'_R are 4-component chiral projections, the vector-like mass term requires equal hypercharges; if they are two-component Weyl spinors, the appearance of \bar{L}'_L and L'^c_R in the same equation is a mixed notation that needs to be defined.
Circularity Check
No significant circularity: the neutrino-sector fit is an input-constrained parameter scan, not a prediction, and no load-bearing claim reduces to its own input.
full rationale
The derivation chain is not circular. The FR(6) charges, fusion rules, and the one-loop formula for y_eta in Eq. (20) are assumed or computed ingredients; they are not derived from the numerical outputs. The closest candidate for circularity is Eqs. (21)-(23), where f and g are re-expressed via Casas-Ibarra from the observed neutrino masses and PMNS matrix, but this is an input-constrained parametrization, not a prediction of D_nu or U. The paper explicitly summarizes its numerical work as demonstrating 'allowed space for our input parameters' rather than as a parameter-free prediction. The LFV, muon g-2, and dark-matter observables are computed after imposing the neutrino fit and can exclude points, so no fitted quantity is renamed as a prediction and no result is forced by construction. The self-citations to non-invertible-symmetry papers supply background formalism, but the central loop calculation is self-contained once the charge assignments are accepted. A separate consistency issue exists: the terms in Eq. (18) may violate the U(1)_Y and FR(6) assignments in Table II, which would be an internal-consistency flaw rather than a circularity, and is not counted in this score.
Assumptions & free parameters
free parameters (7)
- m_S =
[0.1-1000] GeV (scanned)
- m_eta =
[m_S-1000] GeV (scanned)
- M_L' =
[m_S-1000] GeV (scanned)
- M_N =
[m_S-1000] GeV (scanned)
- mu =
[0-1000] GeV (scanned)
- O_N angles (3 complex) =
|theta| in [0,pi], arguments in [0,2pi]
- lightest neutrino mass and Majorana phases =
D_nu1(nu3) <= 1 eV, phases in (-pi,pi)
assumptions (3)
- domain assumption The FR(6) fusion rules in Eq. (15) are taken as the correct non-invertible symmetry algebra.
- domain assumption The non-invertible symmetry is dynamically broken at loop level, so operators that are not invariant can be generated radiatively.
- ad hoc to paper The Yukawa and scalar terms in Eq. (18) and Eq. (19) are invariant under the charge assignments in Table II.
invented entities (3)
-
Non-invertible symmetry FR(6)
-
Vector-like lepton doublets L'
-
Singlet scalar S0
Cite this review
Pith. "Pith review of Three-loop induced neutrino mass model in a non-invertible symmetry." pith.science (2026). https://pith.science/paper/X7MOEUJ3
@misc{pith2026250716198,
author = {Pith},
title = {Pith review of: Three-loop induced neutrino mass model in a non-invertible symmetry},
year = {2026},
howpublished = {\url{https://pith.science/paper/X7MOEUJ3}},
note = {Machine review of arXiv:2507.16198}
}
abstract
We propose a new type of radiatively induced neutrino masses at three-loop level based on the Ma model, introducing a non-invertible symmetry in the class under a ${\mathbb Z_2}$ gauging of ${\mathbb Z_6}$ symmetry and adding three isospin doublet vector-like fermions $L'$ and singlet boson $S_0$. Under this symmetry, the Yukawa interactions directly related to the neutrino masses are not allowed at tree-level. However it is allowed at one-loop level due to $L'$ and $S_0$ as well as $\eta$, which is no longer invariant under this symmetry. Therefore, the symmetry is dynamically broken. Intriguingly, $\eta$ plays important roles in contributing to both the radiative matrices $y^\eta$ and $m_\nu$. After constructing our model, we show some numerical analyses to satisfy the lepton flavor violations, muon anomalous magnetic dipole moment, and a boson dark matter candidate $S_0$ or $\eta_R$ for the cases of normal hierarchy and inverted hierarchy. Then, we demonstrate allowed space for our input parameters.
Figures
Figures from the paper (6 more)
Forward citations
Cited by 3 Pith papers
-
Radiative lepton model in a non-invertible fusion rule
A radiative lepton mass model with a Z2-gauged Z5 non-invertible fusion rule can fit neutrino oscillation data and predicts charged-lepton EDMs that indirectly bound 0νββ to ≲28–33.5 meV.
-
Non-Invertible Selection Rules on Heterotic Non-Abelian Orbifolds
Non-Abelian orbifolds in heterotic string theory produce non-invertible coupling selection rules, since twisted sectors are labeled by conjugacy classes whose products contain multiple classes and yield characteristic...
-
A natural realization of inverse seesaw model in a non-invertible selection rule
A Z3 Tambara-Yamagami fusion rule is used to build an inverse seesaw model with radiatively generated Majorana masses and a dark matter candidate, fitted to neutrino oscillation data.
Reference graph
Works this paper leans on
-
[1]
Y. Choi, H. T. Lam, and S.-H. Shao, Phys. Rev. Lett.129, 161601 (2022), 2205.05086
arXiv 2022
-
[2]
C. Cordova, S. Hong, S. Koren, and K. Ohmori, Phys. Rev. X14, 031033 (2024), 2211.07639
arXiv 2024
-
[3]
DM phenomenology For the calculation of the relic density fromS 0 DM, we simply neglect all Yukawa couplings in the calculation, and hence, the following results are common both for the NH and IH cases, as long as the contributions from these Yukawa couplings are enough small. We choose the relevant parameters as mη = [mS −10 3] GeV, m S = [0−800] GeV,{κ ...
- [4]
- [5]
- [6]
- [7]
-
[8]
T. Kobayashi, Y. Nishioka, H. Otsuka, and M. Tanimoto, JHEP05, 177 (2025), 2503.09966
arXiv 2025
Show all 45 references
-
[9]
Suzuki and L.-X
M. Suzuki and L.-X. Xu (2025), 2503.19964. 17
2025
- [10]
- [11]
- [12]
-
[13]
Kobayashi, H
T. Kobayashi, H. Okada, and H. Otsuka (2025), 2505.14878
2025
- [14]
-
[15]
J. Dong, T. Jeric, T. Kobayashi, R. Nishida, and H. Otsuka (2025), 2507.02375
2025 arXiv
- [16]
- [17]
- [18]
-
[19]
Dijkgraaf, E
R. Dijkgraaf, E. P. Verlinde, and H. L. Verlinde, Commun. Math. Phys.115, 649 (1988)
1988
-
[20]
Kobayashi, S
T. Kobayashi, S. Raby, and R.-J. Zhang, Nucl. Phys. B704, 3 (2005), hep-ph/0409098
2005 arXiv
-
[21]
Kobayashi, H
T. Kobayashi, H. P. Nilles, F. Ploger, S. Raby, and M. Ratz, Nucl. Phys. B768, 135 (2007), hep- ph/0611020
2007
-
[22]
F. Beye, T. Kobayashi, and S. Kuwakino, Phys. Lett. B736, 433 (2014), 1406.4660
2014 arXiv
- [23]
-
[24]
J. J. Heckman, J. McNamara, M. Montero, A. Sharon, C. Vafa, and I. Valenzuela (2024), 2402.00118
2024 arXiv
- [25]
-
[26]
J. Dong, T. Kobayashi, R. Nishida, S. Nishimura, and H. Otsuka (2025), 2504.09773
2025 arXiv
- [27]
- [28]
-
[29]
J. A. Casas and A. Ibarra, Nucl. Phys. B618, 171 (2001), hep-ph/0103065
2001 arXiv
-
[30]
Aghanim et al
N. Aghanim et al. (Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], 1807.06209
2020 arXiv
-
[31]
Abe et al
S. Abe et al. (KamLAND-Zen) (2024), 2406.11438
2024
-
[32]
Aker et al
M. Aker et al. (KATRIN), Science388, adq9592 (2025), 2406.13516
2025
-
[33]
D. P. Aguillard et al. (Muon g-2) (2025), 2506.03069
2025
- [34]
-
[35]
J. Kubo, E. Ma, and D. Suematsu, Phys. Lett. B642, 18 (2006), hep-ph/0604114
2006 arXiv
-
[36]
P. A. R. Ade et al. (Planck), Astron. Astrophys.594, A13 (2016), 1502.01589
2016 arXiv
-
[37]
Hambye, F
T. Hambye, F. S. Ling, L. Lopez Honorez, and J. Rocher, JHEP07, 090 (2009), [Erratum: JHEP 05, 066 (2010)], 0903.4010
2009 arXiv
-
[38]
Esteban, M
I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J. a. P. Pinheiro, and T. Schwetz, JHEP12, 216 (2024), 2410.05380
2024 arXiv
-
[39]
R. L. Workman et al. (Particle Data Group), PTEP2022, 083C01 (2022)
2022
-
[40]
B´ elanger, F
G. B´ elanger, F. Boudjema, A. Goudelis, A. Pukhov, and B. Zaldivar, Comput. Phys. Commun.231, 173 (2018), 1801.03509
2018 arXiv
-
[41]
Belanger, A
G. Belanger, A. Mjallal, and A. Pukhov, Eur. Phys. J. C81, 239 (2021), 2003.08621. 18
2021 arXiv
-
[42]
Belanger et al., JHEP02, 042 (2022), 2111.08027
G. Belanger et al., JHEP02, 042 (2022), 2111.08027
2022 arXiv
-
[43]
Alguero, G
G. Alguero, G. Belanger, S. Kraml, and A. Pukhov, SciPost Phys.13, 124 (2022), 2207.10536
2022 arXiv
-
[44]
Alguero, G
G. Alguero, G. Belanger, F. Boudjema, S. Chakraborti, A. Goudelis, S. Kraml, A. Mjallal, and A. Pukhov, Comput. Phys. Commun.299, 109133 (2024), 2312.14894
2024 arXiv
- [45]
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.