Riding the Seesaw: what Higgsstrahlung may reveal about massive neutrinos

We investigate if the projected high-precision measurements of the cross section of the Higgsstrahlung process e+e-→Zh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+ e^- \rightarrow Zh$$\end{document} at a future electron–positron collider can be utilised to indirectly probe the fermionic Seesaw models. We consider the two centre-of-mass energies s=240\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=240$$\end{document} GeV and 365 GeV, and compare the collider reaches to constraints from electroweak observables, probes of lepton flavour universality and the existing and prospective bounds from searches for lepton flavour violation. For the analysis we assume the limit of an exactly conserved lepton-number symmetry. We find that while any appreciable correction to the Higgsstrahlung cross section is already strictly constrained in the Type-I Seesaw model, effects of up to O(10%)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {O}}(10\%)$$\end{document} are possible within Type-III Seesaw.


Introduction
Upon the discovery of the Higgs boson in 2012 [1; 2], the particle content as proposed within the Standard Model (SM) was confirmed to exist.Precise determinations of the properties of the Higgs boson, as well as in other sectors of physics, are a major component of the proposed program of next-generation lepton colliders.There are five proposals for so-called electron-positron "Higgs factories" [3]: the Circular Electron Positron Collider (CEPC) [4][5][6], the International Linear Collider (ILC) [7][8][9][10], the Future Circular Collider (FCC-ee) [11][12][13][14], the Compact Linear Collider (CLIC) [15][16][17] and the Cool Copper Collider (C 3 ) [18][19][20].The designation "Higgs factory" mainly refers to the stage of a future e + e − collider being run at a centre-of-mass energy of roughly √ s = 240-250 GeV, where the integrated cross section of the Higgsstrahlung process peaks and so dominates over all Higgs production mechanisms.Other stages involve the operation as a "Z factory" at the Z pole for √ s = m Z , close to the W + W − production threshold at √ s = 2m W , and close to or at the t t threshold for √ s = 350-370 GeV, as well as potential upgrades for runs at even higher centre-of-mass energies [21].
The incompleteness of the SM has been corroborated upon the measurement of neutrino oscillations by the Super-Kamiokande [22] and SNO collaborations [23; 24], which evidence the need to introduce new physics in order to account for neutrino masses.A particularly straightforward way to do so is via postulating the existence of sterile neutrinos, several searches for which at future lepton colliders have been proposed.Most of the sensitivity studies focus on direct production of sterile neutrinos in electroweak processes and Higgs production [25][26][27][28][29][30][31][32][33], see also Refs.[34][35][36][37] for reviews on heavy neutral leptons at colliders.Similarly, electroweak triplet fermions at colliders have been studied in, e.g., Refs.[38; 39].The anticipated high precision attainable at future lepton colliders also serves as a motivation to consider virtual corrections, see e.g.Ref. [40].In Refs.[41; 42], the contribution of sterile neutrinos to the triple-Higgs coupling was studied.
In this work, we investigate the potential of electron-positron colliders to test the fermionic Seesaw scenarios [43][44][45][46][47][48], given the sub-percent precision to which these colliders are projected to measure the cross section of the Higgsstrahlung process e + e − → Zh.Achieving this precision is possible due to the so-called "recoil method" that refers to the selection of Higgsstrahlung events solely through measuring the four-momenta of the decay products of the Z boson which recoils against the Higgs boson.This method is in principle applicable for any Higgs decay mode and thus allows for a model-independent reconstruction of the Higgs-boson mass, see for instance Ref. [3].The Higgsstrahlung process is also well-understood in the SM.Recently, two-loop electroweak corrections to the SM cross section were calculated in Refs.[49; 50]; see also Refs.[51][52][53][54][55][56][57].
We consider the process at two benchmark centre-of-mass energies, √ s = 240 GeV and 365 GeV.As argued for in Ref. [58], the smaller cross section at larger s can be partly compensated for by a higher instantaneous luminosity which scales approximately linearly with s.Moreover, the additional boost of the Z and h bosons allows for a better separation of the respective jets and therefore a more precise measurement of σ Zh ×BR(H → X).Further advantages mentioned are the immediate access to top-pair production as well as the e + e − → Hν e ν e process via W + W − fusion which enables a precise determination of the Higgs-boson width, and the possibility to measure the Higgs-boson mass with a precision similar to the one at √ s = 240 GeV.In the fermionic Seesaw models, the smallness of the absolute mass scale of the light neutrinos generically requires the neutrino Yukawa couplings to be tiny for TeV-scale sterile states, and consequently lepton-number conserving processes are suppressed [59][60][61] and beyond the reach of the proposed lepton colliders.In contrast, symmetry-protected Seesaw models [59; 62-72] offer the possibility to entertain relatively light singlets or triplets around the TeV scale without the need to assume tiny Yukawa couplings only.Essentially, the new interaction states are assigned units of lepton number such that it is (approximately) conserved, which implies that some elements of the Dirac and Majorana mass matrices can be sizeable, whereas the remaining ones must be comparatively suppressed.The smallness of active neutrino masses is then guaranteed via the proportionality to these small entries, and does not only rely on an overall suppression through the large mass scale.In addition, no fine-tuned cancellation between (a priori unrelated) elements of the mass matrices are needed.
The remaining paper is organised as follows.In Section 2 we introduce the relevant theoretical framework.We discuss Higgsstrahlung and other relevant processes which are sensitive to the same parameters in Section 3. Our results are summarised and conclusions are drawn in Section 4. In the Appendices, more technical details and useful formulae are collected.

Theory Framework
Firstly, to set our conventions for notation, signs, and normalisation, we state the electroweak (EW) and leptonic parts of the SM Lagrangian: with the gauge-covariant derivative in terms of the Pauli matrices σ a with a = 1, 2, 3.

SMEFT
In the Standard Model Effective Field Theory (SMEFT) we extend the SM Lagrangian with a set of dimension-five and -six operators O i : In our convention the Wilson coefficients C i are dimensional; for later use we also define the dimensionless variants Ĉi ≡ C i × TeV 2 of the dimension-six coefficients.The conventional basis of operators is the so-called Warsaw basis [73].For the convenience of the reader, we reproduce the most relevant operators for the fermionic Seesaw models (additional operators relevant to Higgsstrahlung are also listed in Table 2): with H = iσ 2 H * , and C the charge conjugation matrix satisfying In addition to directly contributing to physical processes, the SMEFT operators also modify the relations between observed quantities and SM parameters, as well as the relations between different parameters.The result of these effects is that the SM parameters will be numerically different in SMEFT compared to the SM.For example, the measured value of the Fermi constant, ĜF , in muon decay will differ from its standard expression G F = 1/( √ 2v 2 T ), v T being the Higgs VEV in SMEFT (see e.g.Ref. [74]), by HL,11 + C (3) When computing the shift in a quantity such as a cross section from its SM value, one must therefore take care to include the contributions from these parameter shifts.Concretely, if σ SM (g i ) is a cross section computed in the SM written as a function of the parameters g i , then where ∆σ Direct captures direct contributions of new operators, and δg i is the shift in the parameter g i .The values of the parameter shifts critically depend on one's choice of input parameters; in this work we use (α, m Z , G F ) as our electroweak inputs.A thorough review of such shifts may be found for instance in Ref. [75], and we provide additional exposition in Appendix A.

LEFT
To describe physics taking place at scales below the electroweak scale we utilise Low-Energy Effective Field Theory (LEFT).The relevant part of the Lagrangian for purely leptonic transitions reads in the Jenkins-Manohar-Stoffer (JMS) basis [76] L We define the covariant derivative in QED as in D µ = ∂ µ + iQeA µ , following Ref.[76].For semi-leptonic neutral-current transitions, the relevant part of the Lagrangian is To obtain the LEFT Wilson coefficients for the Seesaw models introduced below, we utilise the software package DsixTools [77] to (i) compute the renormalisation group (RG) running of the SMEFT coefficients between the Seesaw scale and the electroweak scale, µ = m Z , (ii) match the SMEFT and LEFT coefficients, and (iii) run the LEFT coefficients to the low scale µ = 5 GeV.As there are no sizeable contributions to quark-field operators in the Seesaw models, we assume that this procedure captures the main contributions from RG running in LEFT, and further effects at lower scales do not appreciably change the results.See Appendix B for more details.

Type-I
In the Type-I Seesaw model [43][44][45][46][47] the SM Lagrangian is extended by adding n ν right-handed sterile neutrinos {ν Ri } nν i=1 (so that there are a total of 3 + n ν neutrinos), accompanied by a new Yukawa interaction to generate Dirac neutrino masses, as well as Majorana mass terms for the ν R : Coefficient Type-I Type-III Table 1: SMEFT Wilson coefficients obtained from matching the Type-I and Type-III Seesaw models at the scale µ = M X [78][79][80].In order to properly account for the stringent bounds from the non-observation of lepton flavour violation, the electroweak dipole operators C eB and C eW (last two rows) are matched to one-loop order, while all other operators are matched at tree level.
where the conjugate fields ν c are defined as ν c ≡ γ 0 Cν * [81].After electroweak symmetry breaking we are able to express the combined mass terms as the matrix equation where we refer to √ 2 as the Dirac mass matrix, and M ν is the Majorana mass matrix.Matching this theory onto SMEFT at the scale µ = M ν yields the effective operators collected on the left side of Table 1.

Type-III
In the Type-III Seesaw model [48] the SM Lagrangian is extended by adding n Σ right-handed weak fermion triplets {Σ a Ri } n Σ i=1 with vanishing hypercharge, a new Yukawa interaction to generate Dirac neutrino masses, and Majorana mass terms for the Σ R [69; 82]: where a = 1, 2, 3 is elided from all but the Yukawa term.For a fixed triplet generation i, the eigenstates of electric charge are given by the combinations [69] In a manner completely analogous to the Type-I Seesaw model, after electroweak symmetry breaking we obtain the neutrino mass matrix where we refer to √ 2 as the Dirac mass matrix, and M Σ is the Majorana mass matrix.The states Σ ± Ri instead mix into the charged leptons.Matching this theory onto SMEFT at the scale µ = M Σ yields the effective operators collected on the right side of Table 1.

Conserved Lepton-Number Symmetry
In this work, we study symmetry-protected versions of the fermionic Seesaw models, wherein a lepton-number (LN) symmetry decouples the physics of neutrino masses from the phenomenology associated with the conservation of LN [59; 69-72].Without loss of generality, we fix n ν = n Σ = 2, that is, we focus on the case of two heavy fermion singlet or triplet interaction states which is consistent with at least two massive active neutrinos, as is dictated by neutrino oscillation data.The heavy fermion states are assigned 1 and −1 unit of LN, respectively.After electroweak symmetry breaking, the Dirac mass matrix is given by and the Majorana mass matrix reads where and µ 1,2 are dimensionless parameters.We parametrise the mixing of the SM neutrino ν i with the fermion singlet ν Ri , or the neutral component Σ 0 Ri of the fermion triplet in terms of the dimensionless ratios which are equal to the active-sterile mixing angles in the small-mixing approximation, that is, if O (v/M X ) 3 effects are neglected.For simplicity, we will refer to θ e also as the "electron(-flavour) mixing angle", and to |θ e | as "electron(-flavour) mixing", and equivalently for the other flavours.Light neutrino masses are then proportional to and µ 2 which break LN: The limit µ 2 = 0 and µ 1 = = 0 is referred to as inverse Seesaw [62; 83; 84], and = 0 and µ 1,2 = 0 is commonly known as linear Seesaw [85; 86]. 1e adopt the LN-conserving limit = µ 1,2 = 0 with non-zero M X in this work, which results in massless active neutrinos and a heavy Dirac neutrino of mass M X , and hence assume the textures for both X = ν and Σ.In this way, we neglect the phenomenological implications of LN violation, and instead focus on LN-conserving effects. 2 Note that one may add a further singlet or triplet with vanishing LN such that which supports three massive active neutrinos if one departs from the LN-conserving limit.Still, the additional state trivially decouples from the phenomenology. 3 Phenomenology We choose the benchmark value M X = 1 TeV for the masses of the new interaction states in our analysis.To our knowledge, this is consistent with all performed direct searches for heavy neutral leptons at colliders, see for instance Ref. [37] for a recent overview.In Ref. [88], for sterile neutrinos of a mass M ν ≈ 1 TeV the constraint |θ e | ∼ |θ µ | O(1) was derived via a search for the signature of three charged leptons with any combination of electron and muon flavours.Ref. [89] reports the constraint |θ µ | 2 O(0.1) for TeV-scale Majorana neutrinos, based on a search for same-sign dimuon final states, see also Ref. [61].The bound M Σ 0 ≥ 910 GeV was derived in a recent study [90] which focuses on leptonic final states and takes into account earlier ATLAS results.
In our phenomenological discussion we consider the following observables: • the relative shift ∆σ/σ 0 in the Higgsstrahlung cross section from its SM prediction, • the effective leptonic weak mixing angle sin 2 (θ lept w,eff ) and the W -boson mass m W , • the ratios g X µ/e and g X τ /µ of leptonic gauge couplings as probes of lepton flavour universality (LFU), and the ratios R(K 3 ) and R(V us ), and • the branching ratios of the LFV processes µ → eγ, µ → 3e, τ → eγ, τ → 3e, and the ratios of the µ − e conversion rates over the muon capture rate in different target nuclei.
In Tables 4, 5, 6 and 7, the theoretical expressions for these observables in the fermionic Seesaw models are listed as functions of the mixing angles θ e , θ µ , θ τ , as defined in Eq. ( 16), for a matching scale µ = M X = 1 TeV.While our expressions hold for complex Yukawa couplings, we do not consider CP violation in the analysis and treat the mixings as real numbers which are only constrained from existing bounds on LN-conserving processes.If effects from RG running in SMEFT above the matching scale are neglected, one may naïvely interpret the results also for a larger mass and appropriately rescale the couplings.Nonetheless, we also computed the relevant expressions for µ = M X = 10 TeV, which is commented on in Section 3.5.
The following discussions are supported by plots in the θ e -θ µ plane, as well as plots in the θ e -θ τ plane for the LFU and LFV observables. 4The third mixing angle is fixed to the benchmark values θ τ = 10 −2 and θ µ = 10 −6 ; these choices are most transparently justified (at least for Type-III) by Figure 6, which depicts the most competitive constraints for both models.We have explicitly checked that there are no appreciable changes in the resulting phenomenology if these values are tuned smaller or even zero, apart from the fact that the LFV bounds become weaker and eventually vanish.In each of these plots we produce exclusion regions which either reflect the current bounds at 2σ for the electroweak and LFU observables, or the upper limits on the LFV processes at 90% C.L.
Figure 1: The tree-level cross section for the Higgsstrahlung process in the SM as a function of the centre-of-mass energy √ s.

SM Tree-Level Contribution
The tree-level differential cross section in the SM is well-known, and is given by [91; 92] dσ where unpolarised beams are assumed, with the spin-averaged matrix element.Here θ is the angle between the incoming electron and outgoing Z boson, and is the relevant Källèn function.The corresponding integrated cross section is The dependence of σ 0 on √ s is depicted in Figure 1.It peaks around a centre-of-mass energy of √ s ≈ 245 GeV.

Corrections in SMEFT
Including corrections from new physics, we write where δ|M t | 2 denotes the effect of parameter shifts in SMEFT to the tree-level cross section, as discussed in Section 2.1.1,and 2 Re M * t M c is the interference term of the tree-level amplitude with corrections from new operators.The explicit result reads5 where the parameter shifts and coefficients d i are presented below, and the form factors F i may be found in Ref. [92].
The full resulting fractional shift for the cross section is Here the parameter shifts in the (α, m Z , G F ) input scheme are HL,11 + C HL,11 , and (28b) where we adopt the shorthands s w ≡ sin θ w and c w ≡ cos θ w , and δ ĜF is defined in Eq. ( 5).The integrated form factors f i are and their corresponding coefficients d i are HL,11 + C HL,11 ), and (30c) The diagrams giving rise to the d i 's are depicted in Figure 2. Note that d 1 and f 1 are absent here, as we have elected instead to absorb their contribution into δg ZZh .The dimension-six SMEFT operators which constitute these corrections are listed in Table 2.We comment that the d i 's are zero at tree level in the Seesaw models, apart from d 4 in Type-III.
Table 2: Dimension-6 SMEFT operators in the Warsaw basis [73] which enter the correction to the Higgsstrahlung cross section.Here

Discussion
An essential part of the program of most of the proposed next-generation lepton colliders is to be run as "Higgs factories" at a centre-of-mass energy of √ s = 240 GeV (CEPC, FCC-ee) or 250 GeV (ILC), and to perform a scan of the t t production threshold in the range √ s = 350-380 GeV.Measurements of the Higgsstrahlung cross section are foreseen at both stages for all collider proposals considered herein, apart from CLIC which is envisioned to directly run at √ s = 380 GeV in its initial stage.Therefore, we evaluate the relative shift ∆σ/σ 0 of the Higgsstrahlung cross section, computed in Eq. ( 27), at √ s = 240 GeV and 365 GeV.Regarding the precision of the measurement, we assume the two benchmark values of 0.5% and 1.0% for √ s = 240 GeV, and 1.0% for √ s = 365 GeV.This is representative of the results of several conducted studies of the attainable precision, which are collected in Table 3.
The shifts in the cross section as functions of the mixing angles are listed for both Seesaw models in Table 4.To aid in the following discussion, we moreover find for and for where the Wilson coefficients are evaluated at µ = √ s, respectively, and as a reminder, we define the dimensionless Wilson coefficients by Ĉ = C × TeV 2 .These approximate expressions deviate from the exact results, presented in Table 4, by maximally 5% in either model.
As the couplings of the Z boson to charged leptons are not directly altered at tree level in the Type-I Seesaw model, σ(e + e − → Zh) is predominantly modified via the shift in the Fermi constant, Eq. ( 5), which enters through the shifts δg ZZh , δg L , and δg R in Eq. (28).While the contributions of electron and muon mixing are of fairly similar magnitudes, ∆σ/σ 0 turns out to be slightly more sensitive to the latter.This is due to a partial cancellation of the coefficient d 4 in Eq. (30c) (which acquires a nonzero value from RG running) against the dominant contribution from the Fermi constant for electron mixing.As the corresponding form factor f 4 scales with s, FCC-ee 5 + − 240 0.5% [14] 1.5 + − , q q, ν ν 365 0.9% [14] ILC 1.35 Table 4: Shifts of the Higgsstrahlung cross section at different centre-of-mass energies in terms of the mixing angles, when M X = 1 TeV.For the sake of comparison, we have also computed them for the larger centre-of-mass energy √ s = 500 GeV.
the resulting sensitivity to electron mixing shrinks even further at higher energies.Consequently, if sterile neutrinos are to be searched for via precision Higgs measurements, we do not expect running a next-generation lepton collider at higher centre-of-mass energies to reveal much for Type-I.Contrariwise, in Type-III the couplings of the Z boson to left-handed charged leptons are modified at tree level which results in a sizeable contribution to d 4 .That is, in an EFT language, the model induces the effective four-point interaction eγ µ P L eZ µ h depicted on the right in Figure 2 which interferes with the tree-level contribution to Higgsstrahlung in the SM.As the latter is suppressed by s due to the Z-boson propagator, this results in a very pronounced sensitivity of the ratio ∆σ/σ 0 to electron mixing which approximately scales with s.Consequently, if enough luminosity can be attained to compensate for smaller statistics, fermion triplets may well be searched for in Higgsstrahlung measurements at larger centre-of-mass energies.Note that the contributions from electron mixing could in principle be (partly) cancelled by large muon mixing, however, we will find that this scenario is tightly constrained by existing phenomenological bounds.From Table 4   The W -boson mass listed in the PDG [98] corresponds to the mass parameter in a Breit-Wigner distribution with a mass-dependent width.The SM prediction for the effective leptonic weak mixing angle is taken from Table II in Ref. [96].The model predictions are obtained from matching onto SMEFT at the Seesaw scale µ = M X = 1 TeV and running to the electroweak scale µ = m Z .

Electroweak Sector
In this section, we introduce the shifts in the weak mixing angle and the mass of the W boson.
The expressions obtained for the fermionic Seesaw models in terms of the mixing angles are listed in Table 5.

Weak Mixing Angle
In SMEFT the weak mixing angle is modified in accordance with (see e.g.Ref. [74]) where the shift in the Fermi constant δ ĜF is defined in Eq. ( 5).There are numerous ways to extract the weak mixing angle from data; the most precise determination is that of the effective leptonic weak mixing angle s 2 w,eff ≡ sin 2 (θ lept w,eff ) at LEP [97], achieved via measurements of the left-right asymmetry factor Apart from the general shift in Eq. ( 33), we must also take into account the fact that a modification of the Z couplings to charged leptons will directly affect the extraction of s 2 w,eff from A f .Incorporating the "direct" shifts to these couplings, we find ≈ s 2 w,SM + 0.020 Ĉ(3) HL,11 + Ĉ(3) HL,22 − 0.005 where the right-hand side is evaluated at the scale µ = m Z .Contributions from C He are not sourced at tree level in either fermionic Seesaw model and thus we neglect them in the approximate expression above.

W -Boson Mass
The shift incurred in SMEFT is [74] which approximately evaluates to at the scale µ = m Z .Evidently, the predicted W -boson mass is mainly sensitive to O HL as the major correction is induced via the modified extraction of the Fermi constant in the fermionic Seesaw models.

Discussion
The constraints arising from the electroweak observables s 2 w,eff and m W are illustrated in Figure 3.For Type-I Seesaw one immediately notices that observing a deviation in the Higgsstrahlung cross section is already in direct conflict with the determination of s 2 w,eff , with a 0.5% shift at √ s = 240 GeV suffering from a 2.9σ tension, and a 1% shift excluded at ∼ 6σ.While this tension can be reduced by turning up θ τ (see the expression in Table 5), LFU constraints discussed in Section 3.3 below preclude this from occurring.Thus, the Type-I Seesaw model is unlikely to be a viable minimal SM extension that can be probed in precision Higgs measurements, unless a significant reduction of the statistical uncertainty of these measurements can be attained.
In the Type-III Seesaw model, the effects of the triplets contributing to the extraction of the Fermi constant and from direct contributions to the leptonic gauge couplings largely cancel out, which equally applies to electron and muon flavour, see Eq. ( 36).Thus, s 2 w,eff acts as a rather weak constraint on the Type-III Seesaw model, and is in fact most relevant for tau-flavour mixing, implying |θ τ | 0.06 at 2σ.
The existing tension between the SM prediction for the W -boson mass m W and the larger experimental world average is however exacerbated in Type-III, leading to a much stronger constraint.In contrast, as C HL is induced with equal magnitude, but opposite signs via matching at the Seesaw scale in the two models under consideration, the tension can be alleviated in Type-I; if uncertainties are ignored, the current world average is reproduced for |θ e | 2 + |θ µ | 2 ≈ 0.051, and the CDF measurement [99] for |θ e | 2 + |θ µ | 2 ≈ 0.097.A detectable shift in σ(e + e − → Zh) could in principle only be induced in the latter case.
Testing the Type-III Seesaw model via Higgsstrahlung measurements is generally compatible with the current constraint arising from m W .Still, it disfavours detectable ∆σ/σ 0 induced via muon mixing, and clearly prefers contributions from electron mixing.The bound from m W can be expected to become more competitive when the CDF measurement will be included in the experimental average in the future.

Lepton Flavour Universality
In SMEFT the W ν coupling is directly altered due to O (3) This modifies (semi-)leptonic decays mediated by W bosons, which, among other effects, drives the predictions for LFU ratios away from 1.In the following we consider ratios constituted by the leptonic decays π → ν, K → ν and τ → ν ν.Furthermore, we utilise the extraction of the Cabibbo-Kobayashi-Maskawa (CKM)-matrix element V us from the semi-leptonic decay K 3 : K → π ν and from nuclear beta decays, together with the assumption of unitarity of the CKM matrix.

Ratios of Leptonic Gauge Couplings
We consider the ratios of leptonic gauge couplings collected in Table 6 which are probes of LFU.
For the first three listed therein, [102; 103] where (g µ /g e ) X is the ratio of leptonic gauge couplings extracted from the ratio Γ(X → µν)/Γ(X → eν) for X = π, K, and Γ(τ → µν ν)/Γ(τ → eν ν) for X = . 6Leptonic W -boson decays can also be used to derive constraints on the ratios of leptonic gauge couplings as in g W µ/e ≡ (g µ /g e ) W and g W τ /µ ≡ (g τ /g µ ) W ; still, we do not list them in the table, as the corresponding bounds are merely weaker versions of other constraints.We similarly have Table 6: Current constraints on and model predictions for the LFU ratios taken into account in this work.For R(V us ), the given experimental value refers to the case of N f = 2 + 1 + 1 dynamical quark flavours in the lattice simulations from which the relevant decay constants are extracted.The case of N f = 2 + 1 quark flavours gives rise to a less competitive bound [98].The model predictions are obtained from matching onto SMEFT at the Seesaw scale µ = M X = 1 TeV and running to the electroweak scale µ = m Z .In the right column, the upper sign refers to Type-I Seesaw, whereas the lower sign stands for Type-III Seesaw.
In our Seesaw models the predicted deviation of the considered LFU ratios from 1 is at leading order proportional to ±(|θ µ | 2 − |θ j | 2 ) where j = e or τ , and the sign depends on j as well as which Seesaw model is considered, see Table 6.The derived constraints therefore give rise to hyperbolic contours in the figures presented in this section.Clearly, if the data favours a ratio to be, say, smaller than 1, predicting it to be larger than 1 will lead to a tighter constraint.This then translates into one of the mixing angles being slightly more stringently bounded than another one, with the roles reversed if the other Seesaw model is considered.Consequently, if the contribution from a specific mixing angle accommodates the data well in one model, the other model will necessarily increase the tension with the SM.

Light Quark Mixing
In Table 6 we furthermore consider ratios of the CKM-matrix element V us extracted from the semi-leptonic kaon decays K µ3 and K e3 , the leptonic kaon decay K µ2 , and nuclear beta decaythese are The dependence of R(K 3 ) on new physics is identical to that of g π µ/e and g K µ/e , of which both can similarly be viewed as CKM ratios, with the latter being commonly denoted by R(K 2 ).In the following we illustrate why R(V us ) is a special case, and refer the reader to Ref. [102], wherein the ratio was originally proposed, for further information.
The extracted value of V us from the leptonic kaon decay K µ2 is where V us is a Lagrangian parameter, the C HL,22 is due to the new direct contribution, Eq. ( 39), and δ ĜF /G F parameterises the shift in the Fermi constant.Similarly, beta decay results in which, as per Ref. [102], is translated to V β us using CKM unitarity of the Lagrangian parameters: From this we derive the ratio HL,22 + C HL,11 .
As is explained in Ref. [102], a crucial feature of R(V us ) is the enhanced sensitivity to new physics due to (V ud /V us ) 2 ≈ 20.Experimental data favours the ratio to be smaller than 1 with a significance between 1σ and 2σ, depending on the number of quark flavours assumed for the calculation of the relevant decay constant, see also Table 6.

Discussion
In Figure 4, the constraints arising from LFU ratios are illustrated.Among those, R(V us ) is very sensitive to muon mixing and highly relevant for the fermionic Seesaw models.For Type-I Seesaw, the ratio is enhanced and thus driven further away from the data.In fact, R(V us ) is in conflict with any visible effect in Higgsstrahlung induced via electron or muon mixing, even more so than s w .The next-to-most competitive bounds on electron and muon mixing arise from g K µ/e and g π µ/e , respectively.
In the case of Type-III Seesaw, the most important constraints stem from g π µ/e and R(V us ).Since requiring a discernible deviation in the Higgsstrahlung cross section mainly translates into a lower bound on the electron mixing angle in this model and muon mixing may thus be tuned arbitrarily small, the obtained upper limits on |θ e | are more relevant in this context.g π µ/e demands |θ e | 0.04 at 2σ, unless one allows for larger muon mixing and thus a cancellation, whereupon the maximal electron mixing angle allowed by LFU bounds can be rendered up to 50% larger.R(V us ) implies |θ µ | 0.05 at 2σ, which holds largely independently of the value chosen for |θ e |.Note that in a vein similar to m W , the bounds from g K µ/e and R(V us ) both constrain muon mixing efficiently enough so that no appreciable cancellations of the contributions to Higgsstrahlung from electron mixing can occur for Type-III Seesaw.
As can be seen in Figure 4b, LFU data constrains tau mixing to |θ τ | 0.06 for either Seesaw model, which arises from g τ /µ in Type-I, and g π τ /µ for Type-III.These constraints can in principle be weakened if |θ µ | is more sizeable; still, large changes are only observed for O(0.1) muon mixing, a scenario which is nonetheless excluded by other observables.While the contributions to R(V us ) (a) LFU constraints in the θ e -θ µ plane.
(b) LFU constraints in the θ e -θ τ plane.To avoid clutter, subdominant constraints like the one from R(K 3 ) are not drawn, but their locations may be inferred from the above plots.from electron and tau mixing may cancel, this does not open up parameter space in the Type-I Seesaw model, as tau mixing itself is too constrained.Further LFU ratios not contained in Table 6 deviate from the SM prediction by (close to) 2σ and thus present moderate anomalies in themselves, see Ref. [100].Explicitly, g K τ /µ is measured to be smaller than 1, while g τ /e and g W τ /e exceed the SM expectation.If the models under consideration ought to accommodate the data on g K τ /µ , one would require large tau mixing in comparison with muon mixing for Type-I Seesaw, with the flavours swapped for Type-III Seesaw.The latter is   unlikely to be realised for scenarios which are testable via Higgsstrahlung measurements, as large ∆σ/σ 0 are likely induced via electron mixing in this model, which then demands muon mixing to be very small due to the bounds arising from LFV, see the following section.Similarly, for Type-III Seesaw, the other two ratios necessitate tau mixing to substantially exceed electron mixing in magnitude, which is not a promising scenario either to be tested in the given context, in particular in light of the bound on BR(τ → 3e).

Lepton Flavour Violation
As is generically the case for models of neutrino mass generation, the Seesaw models predict sizeable rates for flavour-violating decays of charged leptons.These processes have not been observed to date and thus impart stringent bounds on the parameter space, which will likely get refined in the near future due to several ongoing or upcoming experiments, see Table 7.The scales relevant for these LFV decays are within the realm of LEFT and thus it is instrumental to discuss them in terms of the contributions to LEFT operators.Since we focus on the comparison with the sensitivities to the Higgsstrahlung process at colliders, we restrict ourselves to observables involving electronflavoured transitions.A comprehensive investigation of LFV effects in the symmetry-protected Type-I Seesaw model can for instance be found in Ref. [118].We will relegate the explicit matching conditions used in this section to Appendix B.

Radiative Charged-Lepton Decays
The branching ratios for radiative flavour-violating charged-lepton decays read [119] BR where the SMEFT Wilson coefficients on the right are evaluated at the electroweak scale µ = m Z .
In both Seesaw models the one-loop matching contributions to the electromagnetic dipole operator O eγ from the electroweak dipole operators O eB and O eW are of the same order of magnitude as the contributions from O (1) HL and O HL which originate from RG running, see also Eq. ( 75) in Appendix B.

Trilepton Decays
The branching ratio for trilepton decays with identical flavours in the final state is given by [120] 7 In the Type-III Seesaw model, these decays are dominated by the vector operators O V LL ee and O V LR ee , with the flavour change occurring in the left-handed lepton bilinear.These operators receive large contributions from tree-level matching of the Type-III Seesaw model onto SMEFT, and then onto LEFT (see Appendix B).By neglecting all Wilson coefficients apart from C V LX ee,jijj with X = L, R, we thus find where the SMEFT Wilson coefficients on the right are evaluated at the electroweak scale µ = m Z .
In the Type-I Seesaw model, all Wilson coefficients entering the branching ratios for trilepton decays receive contributions from matching onto SMEFT only at loop level.In this case, the branching ratios are relatively more sensitive to the contributions from the electromagnetic dipole operator O eγ , and the above approximations are only accurate to about 20%.

µ − e Conversion in Nuclei
As the scalar and gluon operators are suppressed in the fermionic Seesaw models, the µ − e conversion rate takes the simple form [122; 123] where the overlap integrals D, V (p) and V (n) , and muon capture rates ω capt can be found in Refs.[122; 124], and the effective coupling constants are see Appendix B for approximate matching expressions.We are interested in the conversion ratio CR(µ → e), defined as the ratio of the µ−e conversion rate ω conv over the muon capture rate ω capt .For the Type-I and Type-III Seesaw models it approximates to where the Wilson coefficients on the right are evaluated at the low scale µ = 5 GeV, and where the upper, middle and lower entry in the brackets refers to a titanium (Ti), gold (Au) and aluminium (Al) target, respectively.We also include the predictions for the conversion ratios for lead (Pb) and sulfur (S) in Table 7, from which one can infer that the respective current bounds do not impose relevant constraints.As is reflected by the above approximation, µ − e conversion is dominated by contributions from left-handed vector operators in both Seesaw models.This is evident in Type-III where these contributions are sourced at tree level, but also holds in Type-I.The electroweak dipole operator O eγ,12 plays a subdominant role in both models.

Discussion
We find that for the Type-I Seesaw model, as per the absence of tree-level contributions to trilepton decays and µ−e conversion, the most competitive bounds currently arise from the non-observation of µ → eγ and τ → eγ, see Figure 5. Therefore, a detectable shift in σ(e + e − → Zh) would enforce either |θ e | 0.1 and |θ µ | 10 −4 , or vice versa.Still, in the µ − e sector, we expect the limits from BR(µ → 3e) and CR(µ−e) to become more stringent in the future, and further improve this bound by up to two orders of magnitude.The relatively loose current bound from CR(µ − e; Au) is due to cancellations between the effective vector couplings to protons and neutrons.These generically occur for all target materials in both fermionic Seesaw models, but are insignificant for Type-III.GeV and 365 GeV, respectively.We do not depict the bounds from tau-flavoured processes in the θ e -θ µ plots, and vice versa.If included, they would appear as vertical lines with positioning highly dependent on the choice of the third mixing angle, which limits their relevance.
In the case of Type-I, matching at a larger scale reduces the cancellations for CR(µ − e; Au), cf.Section 3.5.See also for instance Ref. [125] for a pertinent discussion.Currently, the nonobservation of τ → eγ enforces |θ τ | 0.  can be expected to become slightly more stringent in light of the future reaches of BR(τ → eγ) and BR(τ → 3e), still, the improvements are expected to be less than an order of magnitude.
In Type-III Seesaw, due to tree-level contributions to the respective pertinent operators, the bounds from BR(µ → 3e) and CR(µ − e) are stricter than that of BR(µ → eγ), and τ → 3e is also more competitive than τ → eγ.Note that the projected sensitivity to BR(µ → eγ) at MEG II cannot even be expected to supersede the current bound on BR(µ → 3e).Moreover, the existing bound on µ − e conversion in gold effectively enforces the muon mixing angle to be smaller than |θ µ | 10 −5 for |θ e | > 10 −2 which is required for an observable deviation of the Higgsstrahlung cross section at a future lepton collider.Thus, the measurements of the cross section are only sensitive to a rather pronounced hierarchy |θ µ /θ e | 10 −3 .A hierarchy |θ e /θ µ | 10 −5 , corresponding to the region visible in the top-left of Figure 5a, is in principle also compatible with the LFV bounds, but still disfavoured by the constraints from electroweak and LFU observables, see the relevant sections above.
Similarly, the non-observation of τ → 3e presently constrains the tau mixing angle to |θ τ | 0.03 if a non-SM signature in σ(e + e − → Zh) is to be attainable, that is, tau mixing should not substantially exceed electron mixing in magnitude.This limit on |θ τ | would be strengthened by an order of magnitude if no decay τ → 3e is observed at Belle II, which will then also necessitate a hierarchy |θ τ /θ e | 0.1 if a chance of detecting ∆σ/σ 0 is to be retained.In a similar vein, if no signals in µ → 3e or µ − e conversion in aluminium are observed in the future, the ratio of the relevant mixing angles will be constrained to be even as small as |θ µ /θ e | 10 −4 or 10 −5 , respectively.

Larger Seesaw Scale
Lastly, we will comment on the scenario with a larger matching scale.In Table 8 we have collected our results for the shifts in the observables considered in our work for the Seesaw scale µ = M X = 10 TeV.
In the case of the Type-III Seesaw model, we find that the results discussed in this section seem to be fairly robust with respect to raising the triplet mass to O(10 TeV) at least.That is, the numerical coefficients entering the expressions for the considered observables typically change by less than 20%.The only notable differences lie in the sensitivity of δm W and ∆σ/σ 0 to tau-flavour mixing at larger centre-of-mass energies, where the respective coefficients grow by a factor of 2-5.Still, it is a subleading effect, as these observables remain much more sensitive to electron-and muon-flavour mixing.Therefore, the results for the Type-III Seesaw model discussed in Section 3 so far, for which M Σ = 1 TeV is assumed, will also approximately hold for (moderately) larger masses.
In contrast, the observable phenomenology of Type-I Seesaw morphs somewhat nontrivially upon raising the Seesaw scale to M ν = O(10 TeV).Most profoundly, the Higgsstrahlung shift ∆σ/σ 0 now experiences a crossing near √ s = 500 GeV, whereupon the dependence on the squared electron mixing angle |θ e | 2 reverses from positive to negative.Additionally, the trilepton decay rates receive a relative numerical boost of 200%, and the µ − e conversion rates are, in general, significantly altered.Specifically, µ − e conversion in gold and lead increase substantially due to the fact that the effective left-handed vector couplings to neutrons increase by a factor larger than 2, while the proton couplings and the dipole operator remain largely unchanged and thus the cancellations are much less efficient.On the contrary, µ − e conversion in aluminium and titanium experience a suppression.This implies in particular that the current bound arising from µ − e conversion in gold is clearly stronger than the one from titanium, unlike the scenario with M ν = 1 TeV.

Summary
We have computed the correction to the tree-level cross section of the Higgsstrahlung process e + e − → Zh in the LN-conserving limit of the Type-I and Type-III Seesaw models, and compared several benchmark sensitivities of next-generation lepton colliders to existing and prospective constraints from electroweak precision measurements, and LFU and LFV probes.Summary plots in the θ e -θ µ and θ e -θ τ planes are presented in Figure 6.
As a major result, we found that existing data on the effective leptonic weak mixing angle and LFU observables preclude substantial corrections to the Higgsstrahlung cross section for Type-I Seesaw.The most likely signature of this model at a future lepton collider is therefore the absence of a detectable deviation from the SM prediction, at least if no further new physics modifying the electroweak and LFU sectors is introduced.For Type-III Seesaw, the current constraints (at 2σ) leave genuinely viable parameter space that can be probed at an e + e − Higgs factory.Figure 7 provides a magnified view of this region.Concretely, for a centre-of-mass energy √ s = 240 GeV the largest permitted shift in the Higgsstrahlung cross section is ∼5%; at √ s = 365 GeV it is ∼12%.The viable region in Type-III is isolated by three main considerations.Firstly, the nonobservation of LFV tightly constrains any scenario with sizeable mixing of heavy fermion singlets or triplets with two lepton flavours.These constraints are particularly strong for Type-III Seesaw, which induces tree-level contributions to trilepton decays and µ−e conversion.Indeed, a detectable deviation in σ(e + e − → Zh) already necessitates a sizeable hierarchy between θ e and θ µ which will become more pronounced if signals of LFV remain elusive in the future.(The situation is similar for Type-I Seesaw, where in the absence of contributions to LFV at tree level, the radiative decays µ → eγ and τ → eγ are more important.)Secondly (and thirdly), both the W-boson mass m W and LFU data currently disfavour detectable corrections induced via muon mixing at the level of  Focusing on the viable region in Figure 7, the constraints arising from the LFU ratio g π µ/e and m W are similarly competitive and provide the most stringent constraint on electron mixing in the Type-III Seesaw model, with |θ e | 0.04 at 2σ.Note also that in the region of parameter space where ∆σ/σ 0 is detectable, tau-flavour mixing is more strictly constrained by the current bound on BR(τ → 3e) than by measurements of the weak mixing angle or pion decays.As is expected from Since fermion triplets induce a tree-level contribution to e + e − → Zh which is not mediated via the s channel as in the SM, deviations from the cross section induced via electron-flavour mixing grow approximately with s.In that sense, if the drop in statistics can be compensated by higher luminosity, the Type-III Seesaw model motivates precision measurements of the Higgsstrahlung process at higher centre-of-mass energies as well, whereas this is not indicated for Type-I Seesaw.
Overall, we have corroborated the expectation that a rich interplay of neutrino, Higgs, electroweak and flavour physics is to be expected for Seesaw models at low energies, and demonstrated the benefit of measuring the Higgsstrahlung cross section at multiple centre-of-mass energies for the Type-III Seesaw model.One may extend the research conducted in this work along two major avenues.Firstly, although the list of processes which we consider in the analysis captures a wide range of phenomenology of the fermionic Seesaw models, it is not exhaustive.In particular, taking into account observables sensitive to angular distributions for Higgs physics [126] as well as a comprehensive global fit in the electroweak sector should help to further differentiate between the low-energy signatures of Seesaw models.Secondly, since we relied on the assumption of an exactly conserved LN symmetry on the Lagrangian level, existing data on lepton mixing and the mass hierarchies in the neutrino sector was per definition not incorporated.While we expect the implications of explicit breakings of lepton number for the induced low-energy phenomenology to be small in general, any viable model of neutrino mass generation eventually needs to be tested against them.Lastly, we leave similar studies for different models of neutrino mass generation for future work.

A A Brief Overview of Parameter Shifts in SMEFT
Here we supplement this work with a somewhat pedagogical overview aimed at explaining how to arrive at the parameter shifts presented in the main text.This overview is not entirely selfcontained, and we recommend the reader refer to Section 5 of Ref. [75], which provides an excellent guided tour of much upon which this exposition relies.We will match their notation with the exception of the gauge fields and their couplings, as their redefinition due to O HW and O HB is just a matter of bookkeeping.(Practically, this means we will continue to write e.g.g 2 and Z µ instead of ḡ2 and Z µ respectively.)Herein we work to order (v/Λ) 2 , where Λ is the cutoff scale of the effective theory, equal to M X in the Seesaw models.
Part 1: Effective Parameters We denote with a bar, e.g.ḡ, an effective parameter that appears in place of its unbarred form in the SMEFT Lagrangian after all effective operators have been expanded out and the dust settles.For example, the Z-boson part of the gauge-covariant derivative, which in the SM reads becomes in SMEFT where [75] ḡZ To denote the shifts of these parameters from their SM expressions we write Note that δḡ Z /g Z = δḡ Z /ḡ Z at order (v/Λ) 2 (and similarly so for every other parameter), so we will consistently opt to write the former, which uses less ink.
To reiterate the notation and have a referenceable equation, the shift for a general parameter g is written ḡ = g + δḡ.
There is one exception to this notational rule in v T , which should be written v for consistency, but isn't by convention.In this case we have v T = v + δv.Part 2: Input Parameters We denote with a hat, e.g.α, a parameter directly measured, or derived from measured values using tree-level SM relations.For a parameter directly measured we have Here there is an additional contribution, labelled δ α, which comes from other diagrams contributing to the measurement process, and which cannot be absorbed into a redefinition of the parameter as above.Such diagrams typically arise at the loop order or at O (v/Λ) 4 , which we neglect, therefore leaving δ α = 0.The sole exception to this is the Fermi constant, which acquires treelevel contributions to δ ĜF from the diagrams in Figure 8.For a parameter computed from the input parameters using the tree-level SM relations, ĝ ≡ g(α i ) (one such example would be computing the elementary charge from the measured value of the fine structure constant via ê = √ 4π α), we instead have ĝ = g + i ∂g ∂α i (δ ᾱi + δ αi ) ≡ g + δĝ.
As will soon become apparent, we will be most interested in computing this total shift, Part 3: Shifts in Observables We now come to the climax of this exposition.Suppose now that SMEFT represents the "true" theory, and that we wish to compute the correction to a quantity such as a cross section which has been calculated under the mistaken assumption that the SM was the correct theory: ∆σ = σ SMEFT − σ SM .In addition to direct contributions from new operators, there will be an indirect contribution from the shifts of the SM parameters.
To arrive at this conclusion carefully, convince yourself that the correct cross section should be written in terms of the barred parameters, σ(ḡ i ), as-after all-these parameters are the ones extant in the "true" Lagrangian.In the SM, however, it's evidently the case that both δḡ i = 0 and δĝ i = 0,8 so one sees no issue with setting g i = ĝi and computing the cross section as σ(ĝ i ).Explicitly, σ SMEFT = σ(ḡ i ) + ∆σ Direct and σ SM = σ(ĝ i ), (65) and so ∆σ = ∆σ Direct + i ∂σ ∂g i g i =ĝ i (ḡ i − ĝi ) ( = ∆σ Direct + i ∂σ ∂g i g i =ĝ i δg i . It is worth emphasising that Eq. ( 66) applies only when computing a correction to an SM prediction, which we do so in Sections 3.1 and 3.2 of the main text.However, as in Sections 3.3 and 3.4, it does not apply when one is only interested in the SMEFT prediction, σ SMEFT .

An explicit example
We illustrate the above procedure by working out the shift in g ZZh , which appears in the SM Lagrangian as In SMEFT this becomes-with additional complications arising from the fact that the Higgs field acquires a new normalisation, see Ref. [75], and that O HD contributes a Z µ Z µ h term- which is a tree-level relation valid in the SM.In accordance with Eq. ( 62), we then have Plugging in the known shifts 9 [75] 9 One arrives at the shift δ ḠF by writing δ mZ is more involved, as m2 Z additionally receives a direct contribution from OHD and an indirect contribution from OHW B due to the rediagonalisation of Zµ and Aµ.Lastly, δ mZ is zero as there are no tree-level diagrams at order (v/Λ) 2 which contribute to the Z-boson self energy.
and leaving δ ĜF /G F symbolic as it has a cumbersome expression (see Eq. ( 5)), this evaluates to Lastly, we bring it all together as per Eq. ( 64) to obtain which one sees matches Eq. ( 28).We note that the dependence on δv completely cancels.This is a general trend for all shifts we consider -thus, for our purposes, the replacement v → v T is functionally unphysical.

B Approximate Matching Conditions
From DsixTools we numerically get the matching condition for the electromagnetic dipole operator, relevant to LFV decays: HL,ij − 0.27353 C with the dashed entries being irrelevant for our purposes.We are conscious that the scales associated with the decays of taus and muons are smaller than µ = 5 GeV.The largest contributions from RG running at lower scales can be expected to originate from QCD.The only SMEFT operator which one may a priori expect to yield a sizeable contribution to LEFT operators which involve quark fields and mix into O eγ is O (3) LeQu = (Lσ µν e R ) (Qσ µν u R ).Still, this operator is not induced at 1-loop level in the Seesaw models under consideration [78][79][80], and it only yields a tiny sub-percent contribution to O eγ .Therefore, we do not expect to have missed any sizeable effects from RG running below µ = 5 GeV. 10 The same holds for the (semi-)leptonic vector operators mediating trilepton decays and µ − e conversion; we have explicitly checked that the contribution from quark operators is at most 1% in Eqs. ( 77) and (78) below.
For all the matching conditions listed in the following, the respective LEFT Wilson coefficient on the left-hand side is given at µ = 5 GeV, while the SMEFT Wilson coefficients entering on the right-hand side are evaluated at the scale µ = m Z .For the operators relevant to trilepton decays 10 Paraphrasing this, even if the numerical factor multiplying C LeQu (mZ ) itself is so small for the fermionic Seesaw models that we do not expect the resulting (relative) contribution to Ceγ(mµ) to become sizeable in any way.

Figure 3 :
Figure 3: Current constraints from electroweak observables at 2σ, in comparison with projected sensitivities of precision Higgsstrahlung measurements.The red-ruled regions indicate |∆σ/σ| < 0.5% at √ s = 240 GeV.The dot-dashed and dotted red lines are the corresponding 1% contours at √ s = 240 GeV and 365 GeV, respectively.The ivory regions for large θ e,µ in the upper-right of the plots are excluded by measurements of both s w and m W .For Type-I Seesaw, the orange dashed line marks where the current experimental world average for m W is exactly accommodated, and in the orange-ruled region the CDF measurement [99] is explained at 2σ.

Figure 4 :
Figure 4: Current constraints arising from LFU ratios at 2σ in comparison with projected sensitivities of precision Higgsstrahlung measurements.The red-ruled regions indicate |∆σ/σ| < 0.5% at √ s = 240 GeV.The dot-dashed and dotted red lines are the corresponding 1% contours at √ s = 240 GeV and 365 GeV, respectively.
−15 is expected after Phase I.The upper (lower) value listed for the future reach of CR(µ − e; Al) refers to COMET (Mu2e).The model predictions are obtained from matching onto SMEFT at the Seesaw scale µ = M X = 1 TeV, running to the electroweak scale µ = m Z , matching onto LEFT and running to the low scale µ = 5 GeV.
(a) LFV constraints in the θ e -θ µ plane.Though we checked the current bounds arising from µ−e conversion for different target materials, here we display only the most competitive ones for each model, respectively.The PRISM/PRIME[112] proposal may even further extend the new-physics reach of µ − e conversion.(b)LFV constraints in the θ e -θ τ plane.

Figure 5 :
Figure 5: Current constraints from LFV observables depicted at 90% C.L., together with their prospective future reaches, in comparison with projected sensitivities of precision Higgsstrahlung measurements.The red-ruled regions indicate |∆σ/σ| < 0.5% at √ s = 240 GeV.The dot-dashed and dotted red lines are the corresponding 1% contours at √ s = 240GeV and 365 GeV, respectively.We do not depict the bounds from tau-flavoured processes in the θ e -θ µ plots, and vice versa.If included, they would appear as vertical lines with positioning highly dependent on the choice of the third mixing angle, which limits their relevance.

Figure 6 :
Figure 6: Summary plots featuring the most competitive current constraints as well as future reaches considered in this work in comparison with projected sensitivities of precision Higgsstrahlung measurements.The red-ruled regions indicate |∆σ/σ| < 0.5% at √ s = 240 GeV.The dotdashed and dotted red lines are the corresponding 1% contours at √ s = 240 GeV and 365 GeV, respectively.

Figure 7 :
Figure 7: Plots zoomed in on the viable parameter regions in the Type-III Seesaw model.Only the most constraining observables are depicted.The red-ruled regions indicate |∆σ/σ| < 0.5% at √ s = 240 GeV.The dot-dashed and dotted red lines are the corresponding 1% contours at √ s = 240 GeV and 365 GeV, respectively.

( 1 )
HL,ij , where the Wilson coefficients on the right-hand side are evaluated at the scale µ = m Z , repeated indices do not indicate summation, and A ij = 10 −3

Table 3 :
[94]cast (statistical) precision of measurements of the Higgsstrahlung cross section at different proposed next-generation colliders.In the third column, the Z-decay final states taken into account in the respective analysis are given; + − always implies both Z → e + e − and Z → µ + µ − .For the results from Ref.[93]a polarisation (P e − , P e + ) = (−80%, +30%) is assumed, still, this changes the expected Higgsstrahlung event rate by maximally 50% compared to unpolarised beams[58; 94].No polarisation is assumed for the ILC precision for √ s = 350 GeV and hadronic Z-boson decays in Ref.[94], wherein an attainable precision of 1.76% for polarised beams and an integrated luminosity of 0.35 ab −1 is reported as well.

Table 5 :
SM predictions for and current measurements of the electroweak observables considered in this work, together with approximate expressions for their shifts in terms of the mixing angles.

Table 7 :
Current and projected constraints on the LFV observables taken into account in this work, together with the pertinent model predictions.The current bounds hold at 90% C.
,The δĝ ZZh part of the full shift depends on one's choice of input parameters; in the (α, m Z , G F ) ZĜF ,