REVIEW 1 major objections 5 minor 29 references
Structure of leptonic Yukawa couplings in the Zee model
T0 review · 1 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read The Zee model's neutrino mass matrix satisfies a single identity that fixes five of the nine lepton Yukawa couplings.
desk verdict Solid but incremental Zee-model parametrization; the central identity is known, and the claimed five/four split needs a genericity condition on f_ij. 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 skew-symmetric $3\times 3$ Yukawa matrix $F$ of the Zee model, whose three independent entries define a pseudovector $u_i = \epsilon_{ijk} f_{jk}/2$. Because $F u = 0$, the master formula $m^\nu = \kappa(F m_\ell Y^\ell + Y^{\ell T} m_\ell F^T)$ implies $u^T m^\nu u = 0$, an identity independent of $Y^\ell$. This identity, together with the five independent entries of $m^\nu$, is what fixes five entries of $Y^\ell$ and leaves four free; it is the engine of the parameter-counting argument.
What would settle it
Compute a Zee-model variant with an extra scalar loop or tree-level neutrino-mass term and show that its physical $m^\nu$ has nonzero $u^T m^\nu u$; or take a neutrino mass matrix from a global fit to oscillation data and check whether Eq. (8) admits any real solution for $f_{12}, f_{13}, f_{23}$ — if no solution exists, the claimed identity cannot hold for the physical $m^\nu$.
Extended reading notes
Core claim
On the paper's own terms: in the Zee model, the one-loop master formula is $m^\nu = \kappa(F m_\ell Y^\ell + Y^{\ell T} m_\ell F^T)$, and the skew-symmetric $F$ defines a pseudovector $u$ with $F u = 0$. Sandwiching the master formula with $u$ gives $u^T m^\nu u = 0$, equivalent to the explicit entry-wise constraint $f_{23}^2 m^\nu_{11}+f_{13}^2 m^\nu_{22}+f_{12}^2 m^\nu_{33}-2f_{13}f_{23}m^\nu_{12}-2f_{12}f_{32}m^\nu_{13}-2f_{21}f_{31}m^\nu_{23}=0$. This relation holds for any $Y^\ell$ and fixes one element of $m^\nu$ once the other five are specified. Solving the master formula for $Y^\ell$ then determines five of its nine entries in terms of $m^\nu$, the charged-lepton masses, and $F$, leaving four undetermined; the paper identifies $Y^\ell_{11}$, $Y^\ell_{12}$, $Y^\ell_{13}$, $Y^\ell_{21}$ as natural free entries and shows they can be set to zero to suppress muonium-antimuonium oscillation and tree-level lepton-flavor-violating decays. The paper then applies this counting to the two-zero texture $B2$, defined by vanishing $m^\nu_{12}$ and $m^\nu_{33}$, and finds a benchmark consistent with oscillation data.
Load-bearing premise
The argument collapses if the one-loop master formula is not the complete and exact source of neutrino masses: the identity $u^T m^\nu u = 0$ follows only when a single loop factor $\kappa$ and a single skew-symmetric $F$ generate the whole mass matrix, so any additional tree-level or radiative contribution would invalidate it.
Editorial extensions
If this is right
- Every Zee-model neutrino mass matrix must obey Eq. (8) regardless of the form of $Y^\ell$, making it a hard consistency condition for model building.
- Once $m^\nu$ and $F$ are chosen, five entries of $Y^\ell$ are fixed, so numerical scans of the model reduce to four free Yukawa entries.
- Setting $Y^\ell_{11}=Y^\ell_{12}=Y^\ell_{13}=Y^\ell_{21}=0$ forbids tree-level muonium-antimuonium oscillation, $\mu\to 3e$, and $\tau\to (3e,\mu e^- e^+)$ while leaving the neutrino sector unchanged.
- With the remaining couplings, the two-zero texture $B2$ gives a muon $g-2$ near the upper edge of the current measured value and $\tau\to e\mu^-\mu^+$ and $\tau\to 3\mu$ rates within reach of projected flavor experiments.
Reading between the lines
- An implicit consequence of the paper's argument is that the same sandwiching trick applies to any radiative model in which the charged-scalar Yukawa matrix is skew-symmetric, so the identity is a general classification tool for one-loop neutrino masses beyond the Zee model itself.
- A natural extension the paper does not pursue is to run all surviving two-zero textures through the same five/four split, testing which combinations of vanishing $Y^\ell$ entries remain consistent with cLFV bounds and the muon $g-2$.
- The parameter counting suggests a fitting strategy: impose Eq. (8) and the five fixed $Y^\ell$ entries as priors, then vary only the four free entries when scanning charged-lepton flavor observables; this would sharpen predictions for upcoming $\mu\to e\gamma$ and tau-decay searches.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies the Zee model and starts from the standard one-loop Majorana mass formula m^nu = kappa (F m_ell Y^ell + Y^{ell T} m_ell F^T), where F is skew-symmetric. The authors define the pseudovector u satisfying F u = 0, derive u^T m^nu u = 0, and write out the resulting component identity, Eq. (8), which relates F and m^nu without involving Y^ell. They then argue that five entries of Y^ell can be determined from m^nu and F, identify Y^ell_{11}, Y^ell_{12}, Y^ell_{13}, Y^ell_{21} as the four undetermined entries, and give simplified inversion formulas, Eq. (10), under the assumptions m_e approx 0 and Y^ell_{21} approx 0. The framework is applied to the two-zero texture B2 (m^nu_{12}=m^nu_{33}=0), with a benchmark point for normal and inverted neutrino mass ordering that satisfies the 3 sigma ranges of NuFit 6.0 and predicts Delta a_mu approx 1 x 10^-9, together with tau decay branching ratios accessible at Belle II.
Significance. If the algebraic identities are correct, the paper provides a compact and potentially useful parametrization of the Zee-model parameter space: the component identity Eq. (8) is a genuine constraint on any Zee-model neutrino mass matrix, and the explicit inversion formulas in Eqs. (9) and (10) can simplify scans over the leptonic Yukawa couplings. The relation itself is not entirely new, however; analogous statements already appear in Refs. [9] and [10], and the Z-Q parametrization of Ref. [9] captures a similar counting of determined versus undetermined parameters. The paper's added value is the explicit component equations and the B2/g-2 application. The numerical benchmark is internally consistent as an illustration, but it is a hand-picked point rather than a statistical fit, and the paper does not supply reproducible code or machine-checked derivations. The central identity is correct, but the universal five/four split of Y^ell requires a genericity condition that the manuscript does not state.
major comments (1)
- [§III.A-B, Eqs. (8)-(10)] The claimed universal five/four split of Y^ell, and in particular the statement that m^nu_{12} is "implicitly encoded in Eq. (8)", are only valid when the relevant coefficients do not vanish. The inversion in Eq. (10) divides by f13 and f23, so it requires f13 f23 != 0. If f13 = 0, Eq. (8) reduces to f23^2 m^nu_{11} + f12^2 m^nu_{33} + 2 f12 f23 m^nu_{13} = 0 and contains no m^nu_{12} at all; in that case Eq. (10a) and Eq. (10d) cannot be used, the set of Y^ell entries that are solved from the displayed equations changes, and m^nu_{12} becomes an independent input rather than an output of Eq. (8). The parametric count of four free entries survives for generic f, but the specific five determined entries are not parameter-independent. The paper should state the genericity condition, specify which components of m^nu are used to solve for which Y^ell entries, and discuss the degenerate cases separately.
minor comments (5)
- [§II, Eq. (3)] The second mixing matrix in Eq. (3) uses H and h for the charged scalar mass eigenstates, which conflicts with the neutral CP-even states H and h introduced in the first matrix. Please use distinct labels for the charged scalars, such as H_1^+ and H_2^+.
- [References [4] and [19]] References [4] and [19] are the same paper (Conlin and Petrov); they should be consolidated into a single citation.
- [§IV, Table I and §V] The summary states that the B2 texture is "fitted to neutrino oscillation data", but the paper shows only one benchmark point constrained by the 3 sigma ranges, without reporting the resulting neutrino mixing parameters such as sin^2 theta_12, sin^2 theta_13, sin^2 theta_23, delta_CP, Delta m^2_21, and Delta m^2_3l, and without any likelihood or pull information. I recommend either performing a real fit or rephrasing this as an illustrative benchmark point.
- [§III.B, Eq. (10)] The abstract says that five entries of Y^ell can be determined directly from m^nu and F, but Eq. (10) shows that only the combinations kappa Y^ell_{ij} are fixed unless the loop factor kappa is treated as known; in Sec. IV kappa is set to 10^-5 by hand. Please qualify the statement to say "given kappa and the charged-lepton masses".
- [§III.B, Eq. (9)] The derivation of the compact relation in Eq. (9d) is not shown and is not transparent from the text. A brief indication of which linear combinations of the component equations produce Eqs. (9a)-(9e), or a supplementary check, would improve readability.
Circularity Check
No circularity: the central neutrino-mass identity and the five/four Yukawa split follow algebraically from the master formula, and the B2 numbers are benchmark postdictions rather than fitted inputs.
full rationale
The load-bearing derivation is self-contained. Equation (7), u^T m_nu u = 0, is obtained directly from the master formula m_nu = kappa(F m_l Y_l + Y_l^T m_l F^T) plus the skew-symmetry of F, which makes both F u = 0 and u^T F = 0; no fitted parameter or external result is needed. Equation (8) is just the component expansion of that scalar identity. The five relations in Eq. (9) are componentwise rearrangements of the same master formula, and Eq. (10) follows after the explicitly stated simplifying choices m_e = 0 and Y_21 = 0. The claim that five entries of Y_l are determined while four remain free is a linear-algebra counting statement following from the rank of these equations, not a quantity that was fitted to data and then renamed as a prediction. In the B2 application, benchmark values of f_ij and Y_l^ij are chosen to satisfy the 3-sigma neutrino oscillation fit, and then Delta a_mu and the tau branching ratios are computed as outputs; the proximity of Delta a_mu to the WP25 1-sigma upper value is a postdiction of the chosen benchmark, not a parameter fitted to that observable. There are no load-bearing self-citations: the cited prior works on the identity and the Z-Q parametrization are external, and the identity is re-derived here. The skeptic's concern about needing nonzero f_13, f_23 and a chosen set of five mass-matrix components is a conditionality/genericity caveat about the stated five/four split, not a circular reduction of the derivation to its own inputs. Hence no circularity is present.
Assumptions & free parameters
free parameters (4)
- κ (loop factor) =
10^-5
- f12, f13, f23 =
See Table I (NO, IO)
- Y22, Y31, Y32, Y33 =
See Table I (NO, IO)
- Y11, Y12, Y13, Y21 =
0
assumptions (5)
- domain assumption The Zee model with a second Higgs doublet and a singly charged singlet scalar, with Lagrangian in Eq. (2), is the correct framework.
- domain assumption F is skew-symmetric with three independent entries.
- domain assumption The neutrino mass matrix is exactly given by mν = κ(F mℓ Yℓ + Yℓ^T mℓ F^T), Eq. (5), with a single loop factor κ and diagonal mℓ.
- domain assumption The two-zero texture B2 (mν12 = mν33 = 0) is compatible with oscillation, cosmology, and 0νββ data.
- ad hoc to paper The four undetermined Yℓ entries (11,12,13,21) can be set to zero without conflicting with other experimental constraints.
Cite this review
Pith. "Pith review of Structure of leptonic Yukawa couplings in the Zee model." pith.science (2026). https://pith.science/paper/PTCVMJAB
@misc{pith2026250818757,
author = {Pith},
title = {Pith review of: Structure of leptonic Yukawa couplings in the Zee model},
year = {2026},
howpublished = {\url{https://pith.science/paper/PTCVMJAB}},
note = {Machine review of arXiv:2508.18757}
}
abstract
The radiative neutrino mass matrix $m^\nu$ in the Zee model depends on leptonic Yukawa couplings $F$ to a singlet scalar and $Y^\ell$ to a new Higgs doublet. Leveraging the skew-symmetric structure of $F$, we derive a unique identity linking $F$ and $m^\nu$ that is explicitly independent of $Y^\ell$. This relation implies that five entries of $Y^\ell$ can, in principle, be determined directly from $m^\nu$ and $F$, while the remaining four can be selected based on phenomenological assumptions. As an illustration, we apply this framework to the two-zero texture $B2$, highlighting its enhancement of the muon $g-2$.
Reference graph
Works this paper leans on
-
[9]
A. C. B. Machado, J. Monta˜ no, P. Pasquini and V. Pleitez, [arXiv:1707.06977 [hep-ph]]
-
[10]
The Singly-Charged Scalar Singlet as the Origin of Neutrino Masses
T. Felkl, J. Herrero-Garcia and M. A. Schmidt, JHEP 05, 122 (2021) [erratum: JHEP 05, 073 (2022)] [arXiv:2102.09898 [hep-ph]]
work page Pith review arXiv 2021
-
[1]
I. Cordero-Carri´ on, M. Hirsch and A. Vicente, Phys. Rev. D 99, no.7, 075019 (2019) [arXiv:1812.03896 [hep-ph]]
arXiv 2019
-
[2]
and HT 2 = (H+, (H0 2 +iA0)/ √ 2). Here, G+ and G0 denote the Goldstone bosons, H0 1,2 are CP-even 3 neutral scalars, and A0 is a CP-odd scalar. If CP-even and CP-odd scalars do not mix, A0 corresponds directly to a mass eigenstate. The CP-even scalars H0 1 and H0 2, as well as the charged scalarsH+ andχ+, can mix. Their relations with the physical mass e...
- [3]
-
[5]
T. A. Chowdhury, J. Heeck, A. Thapa and S. Saad, Phys. Rev. D 106, no.3, 035004 (2022) [arXiv:2204.08390 [hep-ph]]
work page Pith review arXiv 2022
-
[6]
Zee-model predictions for lepton flavor violation
J. Heeck and A. Thapa, Phys. Lett. B 841, 137910 (2023) [arXiv:2303.13383 [hep-ph]]. 10
work page Pith review arXiv 2023
-
[7]
I. Cordero-Carri´ on, M. Hirsch and A. Vicente, Phys. Rev. D 101, no.7, 075032 (2020) [arXiv:1912.08858 [hep-ph]]
arXiv 2020
Show all 29 references
-
[8]
Primulando, J
R. Primulando, J. Julio and P. Uttayarat, Phys. Rev. D 107, no.5, 055034 (2023) [arXiv:2211.16021 [hep-ph]]
2023 arXiv
-
[11]
=κ−1(f12mν 33 + 2f23mν 13) + 2 f 2 23mµYℓ 21 +f12f13meYℓ 13 , (9d) f23(mτYℓ 33−mµYℓ
-
[12]
Herrero-Garc´ ıa, T
J. Herrero-Garc´ ıa, T. Ohlsson, S. Riad and J. Wir´ en, JHEP04, 130 (2017) [arXiv:1701.05345 [hep-ph]]
2017 arXiv
-
[13]
Di Valentino, A
E. Di Valentino, A. Melchiorri and J. Silk, JCAP 01, 013 (2020) [arXiv:1908.01391 [astro- ph.CO]]
2020 arXiv
-
[14]
Aker et al
M. Aker et al. [KATRIN], Science 388, no.6743, adq9592 (2025) [arXiv:2406.13516 [nucl-ex]]
2025
-
[15]
Abe et al
S. Abe et al. [KamLAND-Zen], Phys. Rev. Lett. 130, no.5, 051801 (2023) [arXiv:2203.02139 [hep-ex]]
2023
-
[16]
Davidson and H
S. Davidson and H. E. Haber, Phys. Rev. D 72, 035004 (2005) [erratum: Phys. Rev. D 72, 099902 (2005)] [arXiv:hep-ph/0504050 [hep-ph]]
2005 arXiv
-
[17]
Aghanim et al
N. Aghanim et al. [Planck], Astron. Astrophys. 641, A6 (2020) [erratum: Astron. Astrophys. 652, C4 (2021)] [arXiv:1807.06209 [astro-ph.CO]]
2020 arXiv
-
[18]
Z. z. Xing, Phys. Lett. B 530, 159-166 (2002) [arXiv:hep-ph/0201151 [hep-ph]]
2002 arXiv
-
[19]
Conlin and A
R. Conlin and A. A. Petrov, Phys. Rev. D 102, no.9, 095001 (2020) [arXiv:2005.10276 [hep- ph]]
2020 arXiv
- [20]
-
[21]
These couplings can, in principle, be constrained by the cLFV processes. On the other hand, if the goal is to minimize the number of free parameters to explore the phenomenological implications of the Zee model, one may consider the scenario with Yℓ 11 = Yℓ 12 = Yℓ 13 = Yℓ 21 ...
2025
-
[22]
=κ−1mν 23 +me(f13Yℓ 12 +f12Yℓ 13). (9e) It is evident that Yℓ 11, Yℓ 12, and Yℓ 13, which appear alongside me, could be insensitive to neutrino oscillation parameters, unless the couplings associated with mµ and mτ vanish or are sufficiently small. As shown in Eq. (4), Yℓ cont...
-
[23]
Equivalently, each mν ij is determined by specific Yℓ ij–fij combinations
(10e) Thus, under minimal assumptions, the five entriesYℓ 22,Yℓ 23,Yℓ 31,Yℓ 32, andYℓ 33 associated with fij are fixed by five elements of mν. Equivalently, each mν ij is determined by specific Yℓ ij–fij combinations. It is worth noting that mν 12 is absent from Eq. (10); howe...
-
[24]
Treesukrat, N
W. Treesukrat, N. Supanam and P. Uttayarat, Phys. Rev. D 111, no.7, 075003 (2025) [arXiv:2501.11572 [hep-ph]]
2025 arXiv
-
[25]
P. H. Frampton, S. L. Glashow and D. Marfatia, Phys. Lett. B 536, 79-82 (2002) [arXiv:hep- ph/0201008 [hep-ph]]
2002
-
[26]
Afanaciev et al
K. Afanaciev et al. [MEG II], Eur. Phys. J. C 84, no.3, 216 (2024) [erratum: Eur. Phys. J. C 84, no.10, 1042 (2024)] [arXiv:2310.12614 [hep-ex]]
2024 arXiv
-
[27]
Navas et al
S. Navas et al. [Particle Data Group], Phys. Rev. D 110, no.3, 030001 (2024)
2024
-
[28]
Arora and B
Priya, S. Arora and B. C. Chauhan, [arXiv:2501.00776 [hep-ph]]
-
[29]
Esteban, M
I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J. P. Pinheiro and T. Schwetz, JHEP 12, 216 (2024) [arXiv:2410.05380 [hep-ph]]
2024 arXiv
-
[30]
Aliberti, T
R. Aliberti, T. Aoyama, E. Balzani, A. Bashir, G. Benton, J. Bijnens, V. Biloshytskyi, 11 T. Blum, D. Boito and M. Bruno, et al. [arXiv:2505.21476 [hep-ph]]. 12
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.