UCI–TR–2025–23

Scale-independent relations between neutrino mass parameters

Mu-Chun Chen,  Shaheed Perez  and  Michael Ratz

Department of Physics and Astronomy, University of California, Irvine, CA 92697-4575 USA

Theories of flavor operate at various scales. Recently it has been pointed out that in the context of modular flavor symmetries certain combinations of observables are highly constrained, or even uniquely fixed, by modular invariance and holomorphicity. We find that even in the absence of supersymmetry these combinations are surprisingly immune against quantum corrections.

1 Introduction

Theories of flavor accommodate, or even predict, fermion masses, mixing angles and $\mathcal{CP}$ phases, which constitute a significant fraction of the standard model (SM) parameters. The scale of new physics underlying the corresponding models, which we will denote by $\Lambda_{\mathrm{flavor}}$, generally are different from scales at which experimental measurements are made. This means that quantum corrections to the model predictions have to be taken into account. This raises the question of whether there are predictions that do not depend on the scale $\Lambda_{\mathrm{flavor}}$ at which the model is defined.

In the context of modular flavor symmetries [Feruglio:2017spp] (for reviews see e.g. [Feruglio:2019ybq, Almumin:2022rml, Kobayashi:2023zzc, Ding:2023htn, Ding:2024ozt, Feruglio:2025ztj]) it has recently been pointed out that there are certain combinations of entries of the Weinberg operator are independent of the modulus $\tau$ [Chen:2024otk]. In addition, these combinations are known to be renormalization group (RG) invariant at 1-loop [Chang:2002yr]. This latter statement holds both in SM and minimal supersymmetric standard model (MSSM).

The purpose of this analysis is to discuss the impact of quantum corrections on the above invariance in the absence of supersymmetry (SUSY). This is motivated also by the recent proposal of non-holomorphic modular flavor symmetries [Qu:2024rns, Ding:2024inn, Qu:2025ddz], in non-supersymmetric setups.

2 Neutrino masses described by the Weinberg operator

We consider scenarios in which neutrino masses are described by the Weinberg operator. In the supersymmetric context, the superpotential of the lepton sector is then given by

$$\mathscr{W}_{\mathrm{lepton~mass}}=Y_{e}^{gf}\,L_{g}\cdot H_{d}\,E_{f}+\frac{1}{2}\kappa_{gf}\,L_{g}\cdot H_{u}\,L_{f}\cdot H_{u}\;.$$ (2.1)

Here, the superfields $L^{f}$ and $E^{f}$ denote the three generations of the $\text{SU}(2)_{\mathrm{L}}$ charged lepton doublets and singlets, respectively. The flavor indices $f$ and $g$. The superfields $H_{u/d}$ stand for the MSSM Higgs doublets. In (2.1), “$\cdot$” indicate contractions with the Levi–Civita symbol. $m_{\nu}=v_{u}^{2}\,\kappa$ is the neutrino mass matrix, with $\kappa$ being the effective neutrino mass operator. Finally, $Y_{e}$ denotes the charged lepton Yukawa couplings. In models based on modular flavor symmetries, $\kappa$ and $Y_{e}$ are given in terms of the modular forms.

In the SM amended by the Weinberg operator, the lepton masses are described by

$$\mathscr{L}_{\mathrm{lepton~mass}}=-Y_{e}^{gf}\,\overline{\ell_{\mathrm{L},g}}e_{\mathrm{R}\,g}\cdot\phi-\frac{1}{4}\kappa_{gf}\ell^{g}\cdot\phi\,\ell^{f}\cdot\phi+\text{h.c.}\;.$$ (2.2)

Here, $\ell_{\mathrm{L},f}$ denote the lepton doublets, $e_{\mathrm{R}\,g}$ the right-handed charged leptons, and $\phi$ the SM Higgs.

Apart from the charged lepton masses, $m_{f}=y_{f}\,v_{\mathrm{EW}}$ with $v_{\mathrm{EW}}$ denoting the vacuum expectation value (VEV) of the electroweak Higgs $\phi$, the lepton sector has 9 flavor parameters,

$$\{\xi_{i}\}=\{m_{1},m_{2},m_{3},\theta_{12},\theta_{13},\theta_{23},\delta,\varphi_{1},\varphi_{2}\}\;.$$ (2.3)

Out of these parameters, two mass squared differences and the mixing angles $\theta_{ij}$ have been measured with relatively good precision, see e.g. [Esteban:2024eli]. On the other hand, the absolute neutrino mass scale and the Dirac phase $\delta$ are subject to constraints but not determined precisely. We currently do not know whether neutrinos are Majorana fermions, and thus have no knowledge of the values of the Majorana phases $\varphi_{i}$.

3 Lepton flavor parameters and quantum corrections

Equations˜2.1 and 2.2 contain the Weinberg operator,

$$\mathscr{L}_{\kappa}=-\frac{1}{4}\kappa_{gf}\ell^{g}\cdot\phi\,\ell^{f}\cdot\phi+\text{h.c.}\;.$$ (3.1)

$\kappa$ is a symmetric matrix of mass dimension $-1$.

Throughout this study, we will work in a basis in which $Y_{e}$ is diagonal and positive,

$$Y_{e}=\operatorname{diag}(y_{e},y_{\mu},y_{\tau})\quad\text{with }y_{f}>0\text{ for }f\in\{e,\mu,\tau\}\;.$$ (3.2)

In this basis, all the renormalizable interactions in the lepton sector are diagonal in flavor space.

3.1 Invariants

In the basis chosen as given in (3.2), we define the invariants

$$I_{fg}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\frac{(m_{\nu})_{ff}\,(m_{\nu})_{gg}}{\bigl((m_{\nu})_{fg}\bigr)^{2}}=\frac{\kappa_{ff}\,\kappa_{gg}}{\bigl(\kappa_{fg}\bigr)^{2}}\;,$$ (3.3)

where no summation over the flavor indices $f$ and $g$ is implied. We are interested in quantum corrections to these combinations. In order to obtain the second equality in (3.3), we have to assume that the normalizations of the three lepton doublets coincide at a given scale. This can be achieved in bottom-up [Chen:2021prl] and top-down models [Li:2025bsr]. The focus of this study is on the RG stability of the $I_{fg}$ (3.3).

A key feature of these expressions is that they can be entirely expressed in terms of observable flavor parameters. Explicitly,

$$I_{12} =\frac{a_{0}\left[\widetilde{m}_{1}(c_{23}s_{12}+\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{12}s_{13}s_{23})^{2}+\widetilde{m}_{2}(c_{12}c_{23}-\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}s_{12}s_{13}s_{23})^{2}+m_{3}c_{13}^{2}s_{23}^{2}\right]}{\left[\widetilde{m}_{1}c_{12}(c_{23}s_{12}+\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{12}s_{13}s_{23})-\widetilde{m}_{2}s_{12}(c_{12}c_{23}-\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}s_{12}s_{13}s_{23})-\mathrm{e}^{\thinspace\mathrm{i}\thinspace\delta}m_{3}s_{13}s_{23}\right]^{2}}$$ (3.4a)

$$I_{13} =\frac{a_{0}\left[\widetilde{m}_{1}(\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{12}c_{23}s_{13}-s_{12}s_{23})^{2}+\widetilde{m}_{2}(\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{23}s_{12}s_{13}+c_{12}s_{23})^{2}+m_{3}c_{13}^{2}c_{23}^{2}\right]}{\left[\widetilde{m}_{1}c_{12}(\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{12}c_{23}s_{13}-s_{12}s_{23})+\widetilde{m}_{2}s_{12}(\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{23}s_{12}s_{13}+c_{12}s_{23})-\mathrm{e}^{\thinspace\mathrm{i}\thinspace\delta}m_{3}c_{13}c_{23}\right]^{2}}$$ (3.4b)

$$I_{23} =\left[m_{3}c^{2}_{13}s^{2}_{23}+\widetilde{m}_{1}\left(c_{23}s_{12}+\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{12}s_{13}s_{23}\right)^{2}+\widetilde{m}_{2}\left(c_{12}c_{23}-\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}s_{12}s_{13}s_{23}\right)^{2}\right]$$

$$\qquad{}\times\frac{4\left[m_{3}c^{2}_{13}c^{2}_{23}+\widetilde{m}_{2}\left(\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{23}s_{12}s_{13}+c_{12}s_{23}\right)^{2}+\widetilde{m}_{1}\left(\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}c_{12}c_{23}s_{13}-s_{12}s_{23}\right)^{2}\right]}{\left[(\widetilde{m}_{1}a_{1}+\widetilde{m}_{2}a_{2})-m_{3}\sin(2\theta_{23})c^{2}_{13}\right]^{2}}\;,$$ (3.4c)

$s_{ij}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\sin\theta_{ij}$, $c_{ij}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\cos\theta_{ij}$, $t_{ij}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\tan\theta_{ij}$, and

$$a_{0} \mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\widetilde{m}_{1}c_{12}^{2}+\widetilde{m}_{2}s_{12}^{2}+\mathrm{e}^{2\thinspace\mathrm{i}\thinspace\delta}m_{3}t_{13}^{2}$$ (3.5a)

$$a_{1} \mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\left[\left(s^{2}_{12}-\mathrm{e}^{-2\thinspace\mathrm{i}\thinspace\delta}c^{2}_{12}s^{2}_{13}\right)\sin(2\theta_{23})-\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}\cos(2\theta_{23})\sin(2\theta_{12})s_{13}\right]\;,$$ (3.5b)

$$a_{2} \mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\left[\mathrm{e}^{-\thinspace\mathrm{i}\thinspace\delta}\cos(2\theta_{23})\sin(2\theta_{12})s_{13}+\left(c^{2}_{12}-\mathrm{e}^{-2\thinspace\mathrm{i}\thinspace\delta}s^{2}_{12}s^{2}_{13}\right)\sin(2\theta_{23})\right]\;.$$ (3.5c)

The invariants $I_{fg}$ depend on $m_{1}$, $m_{2}$, $\varphi_{1}$ and $\varphi_{2}$ only via the combinations $\widetilde{m}_{1}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}m_{1}\,\mathrm{e}^{\thinspace\mathrm{i}\thinspace\varphi_{1}}$ and $\widetilde{m}_{2}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}m_{2}\,\mathrm{e}^{\thinspace\mathrm{i}\thinspace\varphi_{2}}$. As $I_{fg}$ are complex, each of them contains two real flavor parameters. This means that, unless there are degeneracies, 6 independent linear combinations out of the 9 flavor parameters $\xi_{i}$ in (2.3) are described by the $I_{fg}$.

Why are we interested in these invariants, $I_{fg}$? There are two main reasons. First of all, they are RG invariant at the 1-loop level [Chang:2002yr], as we shall discuss in more detail in Section˜3.2. Additionally, they turn out to have remarkable properties in the framework of modular flavor symmetries. For instance, in the Feruglio model [Feruglio:2017spp], $I_{12}=-2$ and $I_{13}\,I_{23}=-32$, independently of the value of the modulus [Chen:2024otk]. That is, these invariants carry a large amount of the information on the modular symmetries. As we shall see next, they are not only independent of the modulus but also, for all practical purposes, insensitive to the definition flavor scale $\Lambda_{\mathrm{flavor}}$ of the model.

3.2 Renormalization group equations

In [Chang:2002yr] it has been found that $I_{fg}$ defined in (3.3) are independent of the renormalization scale at one-loop. In the supersymmetric context, one may view this as being a simple consequence of the non-renormalization theorem. Only the wave-function renormalization constants depend on the scale, and the latter cancel in the $I_{fg}$ expressions [Haba:1999ca]. The RG-invariance of $I_{fg}$ thus holds at all loop levels in supersymmetric models. So in the following, we will focus on non-supersymmetric case.

Since $\kappa$ is a symmetric matrix, its renormalization group equation (RGE) has the form

$$\frac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}t}\kappa=\sum_{k}\kappa^{(k)}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\sum_{k}\left(\frac{1}{16\piup^{2}}\right)^{k}\,\left[\alpha^{(k)}\kappa+P^{(k)}\,\kappa+\kappa\bigl(P^{(k)}\bigr)^{{\mathchoice{\raisebox{0.0pt}{$\displaystyle\intercal$}}{\raisebox{0.0pt}{$\textstyle\intercal$}}{\raisebox{0.0pt}{$\scriptstyle\intercal$}}{\raisebox{0.0pt}{$\scriptscriptstyle\intercal$}}}}+Q^{(k)}\,\kappa\,\bigl(Q^{(k)}\bigr)^{{\mathchoice{\raisebox{0.0pt}{$\displaystyle\intercal$}}{\raisebox{0.0pt}{$\textstyle\intercal$}}{\raisebox{0.0pt}{$\scriptstyle\intercal$}}{\raisebox{0.0pt}{$\scriptscriptstyle\intercal$}}}}\right]\;.$$ (3.6)

The $t$-derivative is the logarithmic derivative with respect to the renormalization scale $\mu$,

$$\frac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}t}\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\mu\frac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}\mu}\;,$$ (3.7)

i.e. $t=\ln(\mu/\mu_{0})$ with some reference scale $\mu_{0}$. In (3.6), $k$ indicates the loop level, and the $\alpha^{(k)}$ are flavor-independent coefficients. The matrices $P^{(k)}$, $Q^{(k)}$ are composed of the renormalizable couplings of the theory and diagonal,

$$P^{(k)} =\operatorname{diag}\bigl(P^{(k)}_{1},P^{(k)}_{2},P^{(k)}_{3}\bigr)\;,$$ (3.8a)

$$Q^{(k)} =\operatorname{diag}\bigl(Q^{(k)}_{1},Q^{(k)}_{2},Q^{(k)}_{3}\bigr)\;.$$ (3.8b)

At one-loop, $P^{(1)}=C_{e}\,Y_{e}Y_{e}^{\dagger}=C_{e}\,\operatorname{diag}(y_{e}^{2},y_{\mu}^{2},y_{\tau}^{2})$ with $Y_{e}$ being the charged lepton Yukawa matrix (3.2). $C_{e}=-3/2$ in the SM [Antusch:2001ck] and two-Higgs models [Antusch:2001vn], and $C_{e}=1$ in the MSSM [Chankowski:1993tx, Babu:1993qv]. At the 1-loop level in (3.6), there is only one matrix in flavor space, $P^{(1)}$, and we can choose $Q^{(1)}=\mathds{1}$.

The two-loop contribution in (3.6) in the SM is given by [Ibarra:2024tpt]

$$\kappa^{(2)} =\left(\frac{1}{16\piup^{2}}\right)^{2}\left[\alpha^{(2)}\,\kappa+P^{(2)}\,\kappa+\kappa\,\bigl(P^{(2)}\bigr)^{{\mathchoice{\raisebox{0.0pt}{$\displaystyle\intercal$}}{\raisebox{0.0pt}{$\textstyle\intercal$}}{\raisebox{0.0pt}{$\scriptstyle\intercal$}}{\raisebox{0.0pt}{$\scriptscriptstyle\intercal$}}}}+Q^{(2)}\,\kappa\,\bigl(Q^{(2)}\bigr)^{{\mathchoice{\raisebox{0.0pt}{$\displaystyle\intercal$}}{\raisebox{0.0pt}{$\textstyle\intercal$}}{\raisebox{0.0pt}{$\scriptstyle\intercal$}}{\raisebox{0.0pt}{$\scriptscriptstyle\intercal$}}}}\right]\;,$$ (3.9)

where now there is a nontrivial $Q$-matrix,

$$P^{(2)} =\left(-\frac{57}{16}g_{1}^{2}+\frac{33}{16}g_{2}^{2}+\frac{5}{4}T\right)\,Y_{e}Y_{e}^{\dagger}+\frac{19}{4}Y_{e}Y_{e}^{\dagger}\,Y_{e}Y_{e}^{\dagger}\;,$$ (3.10a)

$$Q^{(2)} =\sqrt{2}\,Y_{e}Y_{e}^{\dagger}\;.$$ (3.10b)

Here, $g_{1}$ and $g_{2}$ are the running gauge coupling constants and $T\mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\operatorname{Tr}[Y_{e}Y_{e}^{\dagger}+3Y_{u}Y_{u}^{\dagger}+3Y_{d}Y_{d}^{\dagger}]$, with $Y_{u}$ and $Y_{d}$ being the Yukawa coupling matrices for the up-type quarks and the down-type quarks, respectively.

Analogously to (3.6), we can write the loop expansion of the $I_{fg}$ as

$$I_{fg} \mathchoice{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}=}{\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}=}\frac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}t}I_{fg}=\sum_{k}I_{fg}^{(k)}\;,$$ (3.11)

where

$$I_{fg}^{(k)} =\frac{\kappa_{ff}^{(k)}\,\kappa_{gg}}{\bigl(\kappa_{fg}\bigr)^{2}}+\frac{\kappa_{ff}\,\kappa_{gg}^{(k)}}{\bigl(\kappa_{fg}\bigr)^{2}}-2\frac{\kappa_{ff}\,\kappa_{gg}}{\bigl(\kappa_{fg}\bigr)^{3}}\,\kappa_{fg}^{(k)}\;.$$ (3.12)

Truncating (3.6) at two-loop level, i.e. $k=2$, and inserting this truncation into (3.11) we obtain

$$\frac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}t}I_{fg} =\frac{(\kappa_{ff}^{(1)}+\kappa_{ff}^{(2)})\,\kappa_{gg}}{\kappa_{fg}^{2}}+\frac{\kappa_{ff}\,(\kappa_{gg}^{(1)}+\kappa_{gg}^{(2)})}{\kappa_{fg}^{2}}-2\frac{\kappa_{ff}\,\kappa_{gg}}{\kappa_{fg}^{3}}\,(\kappa_{fg}^{(1)}+\kappa_{fg}^{(2)})$$

$$=\left[\frac{\kappa_{ff}^{(1)}\,\kappa_{gg}}{\kappa_{fg}^{2}}+\frac{\kappa_{ff}\,\kappa_{gg}^{(1)}}{\kappa_{fg}^{2}}-2\frac{\kappa_{ff}\,\kappa_{gg}}{\kappa_{fg}^{3}}\,\kappa_{fg}^{(1)}\right]+\left[\frac{\kappa_{ff}^{(2)}\,\kappa_{gg}}{\kappa_{fg}^{2}}+\frac{\kappa_{ff}\,\kappa_{gg}^{(2)}}{\kappa_{fg}^{2}}-2\frac{\kappa_{ff}\,\kappa_{gg}}{\kappa_{fg}^{3}}\,\kappa_{fg}^{(2)}\right]$$

$$\mathchoice{=\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}}{=\mathrel{\hbox to0.0pt{\raisebox{1.03334pt}{$\cdot$}\hss}\raisebox{-1.03334pt}{$\cdot$}}}{=\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\small$\cdot$}\hss}\raisebox{-1.20554pt}{\small$\cdot$}}}{=\mathrel{\hbox to0.0pt{\raisebox{0.34444pt}{\tiny$\cdot$}\hss}\raisebox{-1.20554pt}{\tiny$\cdot$}}}I_{fg}^{(1)}+I_{fg}^{(2)}\;.$$ (3.13)

Then, we can calculate $I_{fg}^{(k)}$ as

$$I_{fg}^{(k)} =\frac{\kappa_{ff}\,\kappa_{gg}}{(16\piup^{2})^{k}\kappa_{fg}^{2}}\left[\alpha^{(k)}+2P_{ff}^{(k)}+\bigl(Q_{ff}^{(k)}\bigr)^{2}+\alpha^{(k)}+2P_{gg}^{(k)}+\bigl(Q_{gg}^{(k)}\bigr)^{2}\right.$$

$$\hphantom{{}=\frac{\kappa_{ff}\,\kappa_{gg}}{(16\piup^{2})^{k}\kappa_{fg}^{2}}}\quad{}\left.\vphantom{\bigl(Q_{gg}^{(k)}\bigr)^{2}}-2\bigl(\alpha^{(k)}+P_{ff}^{(k)}+P_{gg}^{(k)}+Q_{ff}^{(k)}Q_{gg}^{(k)}\bigr)\right]$$

$$=\frac{\kappa_{ff}\,\kappa_{gg}}{(16\piup^{2})^{k}\kappa_{fg}^{2}}\left(Q_{ff}^{(k)}-Q_{gg}^{(k)}\right)^{2}\;.$$ (3.14)

At 1-loop, $Q_{ff}^{(1)}=0$, so there is no correction to $I_{fg}$ at this order. At 2-loop, $Q^{(2)}$ is a diagonal matrix with diagonal elements given by

$$Q_{ff}^{(2)}=\sqrt{2}\,(Y_{e}Y_{e}^{\dagger})_{ff}=\sqrt{2}\,y_{f}^{2}\;,$$ (3.15)

where we use $y_{f}>0$. Therefore, $\dot{I}_{fg}$ up to 2-loop using equation (3.14) is explicitly given by

$$\frac{\mathop{}\!\mathrm{d}I_{fg}}{\mathop{}\!\mathrm{d}t}=\frac{2\bigl(y_{f}^{2}-y_{g}^{2}\bigr)^{2}}{(16\piup^{2})^{2}}I_{fg}\;.$$ (3.16)

Specifically,

$$\frac{\mathop{}\!\mathrm{d}I_{12}}{\mathop{}\!\mathrm{d}t} =\frac{2\bigl(y_{e}^{2}-y_{\mu}^{2}\bigr)^{2}}{(16\piup^{2})^{2}}I_{12}\;,$$ (3.17a)

$$\frac{\mathop{}\!\mathrm{d}I_{13}}{\mathop{}\!\mathrm{d}t} =\frac{2\bigl(y_{e}^{2}-y_{\tau}^{2}\bigr)^{2}}{(16\piup^{2})^{2}}I_{13}\;,$$ (3.17b)

$$\frac{\mathop{}\!\mathrm{d}I_{23}}{\mathop{}\!\mathrm{d}t} =\frac{2\bigl(y_{\mu}^{2}-y_{\tau}^{2}\bigr)^{2}}{(16\piup^{2})^{2}}I_{23}\;.$$ (3.17c)

Interestingly, these results show that if $I_{fg}$ vanishes at some scale, it will stay zero at all scales. Since the coefficients on the right-hand sides of (3.17) are real, this statement applies separately to the real and imaginary parts of the $I_{fg}$. That is, if $\operatorname{Re}I_{fg}$ or $\operatorname{Im}I_{fg}$ vanishes at some scale, it will remain zero at all scales.

Given the hierarchy $y_{\tau}\gg y_{\mu}\gg y_{e}$, we see that the corrections of $I_{12}$ are even more suppressed than the RG effects on $I_{13}$ and $I_{23}$. Even the latter are basically RG stable. Since $y_{\tau}\sim 10^{-2}$, the coefficient is of the order $10^{-10}$. Multiplying this by $\ln(\Lambda_{\mathrm{flavor}}/v_{\mathrm{EW}})$ still leads to RG effects at most of the order $10^{-8}$. This means that, for all practical purposes, the $I_{fg}$ are invariant under the renormalization group in the SM, and thus not sensitive to the flavor scale $\Lambda_{\mathrm{flavor}}$.

Using equations (3.17), we can estimate benchmark values to the quantum corrections. Choosing $y_{e}\approx 2\times 10^{-6}$, $y_{\mu}\approx 5\times 10^{-4}$, and $y_{\tau}\approx 7\times 10^{-3}$, the correction to $I_{fg}$ using Equation˜3.16 is given by

$$\Delta I_{fg}\approx\frac{2\bigl(y_{f}^{2}-y_{g}^{2}\bigr)^{2}}{(16\piup^{2})^{2}}I_{fg}\Delta t\;.$$ (3.18)

We have used a modified version of REAP [Antusch:2005gp] to verify that, when running the invariants at two loop over a few orders of magnitude, they remain practically unchanged. A more detailed numerical study will be presented elsewhere.

Let us comment on two-Higgs doublet models. Usually one imposes symmetries to make sure that the charged leptons only couple to one of the Higgs doublets in order to avoid flavor changing neutral currents [Weinberg:1976hu, Glashow:1976nt, Paschos:1977ay], cf. the discussion in [Antusch:2001vn]. In these models $y_{\tau}$ may be of the order unity. Even in this case, the corrections (3.17) remain well below the percent level.

3.3 Limitations

In our analysis, we have focused on the case in which the model gives rise to SM, MSSM or a 2HDM below its definition scale. If there are additional renormalizable couplings that are sensitive to specific lepton flavors, our analysis may no longer apply. Studying such scenarios is beyond the scope of this work.

4 Summary

Motivated by the analytic properties of certain combinations of the neutrino mass matrix $I_{fg}$ in the context of modular flavor symmetries, we have studied the stability of these expressions under the renormalization group. While the $I_{fg}$ receive corrections at the two-loop level, for all practical purposes they remain RG invariant in the SM and 2HDM, i.e. in the absence of SUSY. This leads to predictions that are insensitive to the scale $\Lambda_{\mathrm{flavor}}$ at which the model is defined. The conclusions drawn from the analytical properties of the $I_{fg}$ can therefore be confronted to data without the need of a detailed renormalization group analysis. In other words, experimental measurements can directly probe ultraviolet (UV) physics.

Acknowledgments

We are indebted to Michael Schmidt for help with the REAP package. S.P. acknowledges support from the APS Bridge program.

2HDM

two-Higgs doublet model

BSM

beyond the standard model

BU

bottom-up

CG

Clebsch–Gordan

EFT

effective field theory

FCNC

flavor changing neutral current

FI

Fayet–Iliopoulos [Fayet:1974jb]

GS

Green–Schwarz [Green:1984sg]

GUT

Grand Unified Theory

IO

inverted ordering

IR

infrared

LEET

low-energy effective theory

LHC

Large Hadron Collider

MIHO

modular invariant holomorphic observables

MSSM

minimal supersymmetric standard model

NO

normal ordering

OPE

operator product expansion

QFT

quantum field theory

RG

renormalization group

RGE

renormalization group equation

SM

standard model

SUSY

supersymmetry

TD

top-down

UV

ultraviolet

VEV

vacuum expectation value

VVMF

vector-valued modular form