Gravitational footprints of massive neutrinos and lepton number breaking
Andrea Addazi,1, 2, ∗ Antonino Marcianò,1, † António P. Morais,3, ‡ Roman Pasechnik,4, § Rahul Srivastava,5, ¶ and José W. F. Valle5, ∗∗
1Department of Physics & Center for Field Theory and Particle Physics, Fudan University, 200433 Shanghai, China 2College of Physics, Sichuan University, Chengdu, 610065, China 3Departamento de Física, Universidade de Aveiro and CIDMA, Campus de Santiago, 3810-183 Aveiro, Portugal, EU 4Department of Astronomy and Theoretical Physics, Lund University, 221 00 Lund, Sweden, EU 5AHEP Group, Institut de Física Corpuscular – C.S.I.C./Universitat de València, Parc Científic de Paterna. C/ Catedrático José Beltrán, 2 E-46980 Paterna (Valencia) - SPAIN
We investigate the production of primordial Gravitational Waves (GWs) arising from First Order Phase Transitions (FOPTs) associated to neutrino mass generation in the context of type-I and inverse seesaw schemes. We examine both “high-scale” as well as “low-scale” variants, with either explicit or spontaneously broken lepton number symmetry U(1)L in the neutrino sector. In the latter case, a pseudo-Goldstone majoron-like boson may provide a candidate for cosmological dark matter. We find that schemes with softly-broken U(1)L and with single Higgs-doublet scalar sector lead to either no FOPTs or too weak FOPTs, precluding the detectability of GWs in present or near future measurements. Nevertheless, we found that, in the majoron-like seesaw scheme with spontaneously broken U(1)L at finite temperatures, one can have strong FOPTs and non-trivial primordial GW spectra which can fall well within the frequency and amplitude sensitivity of upcoming experiments, including LISA, BBO and u-DECIGO. However, GWs observability clashes with invisible Higgs decay constraints from the LHC. A simple and consistent fix is to assume the majoron-like mass to lie above the Higgs-decay kinematical threshold. We also found that the majoron-like variant of the low-scale seesaw mechanism implies a different GW spectrum than the one expected in the high-scale seesaw. This feature will be testable in future experiments. Our analysis shows that GWs can provide a new and complementary portal to test the neutrino mass generation mechanism.
I. INTRODUCTION and neutrinoless double beta decay (0ν2β) searches, neu- trino masses can be probed through Charged Lepton Fla- vor Violation (CLFV) experiments at the high intensity Non-zero neutrino masses constitute one of the most ro- and/or high energy frontier [21–23]. Moreover, neutrino bust evidences for new physics [1–3]. Despite great efforts mass generation can leave signatures at high-energy col- over the last two decades to underpin the origin of neu- liders like the Large Hadron Collider (LHC) [24–27]. trino mass, the basic underlying mechanism remains as elusive as ever. Small neutrino masses can be generated The detection of Gravitational Waves (GWs) by the in many ways, both for Majorana [4,5] and Dirac [6–11] LIGO team has opened an entirely novel method to probe arXiv:1909.09740v2 [hep-ph] 2 May 2020 neutrinos. Here, we focus on the various variants of the the underlying new physics associated to neutrino mass popular type-I seesaw mechanism for Majorana neutri- generation. It was advocated that the spectrum of pri- nos [12–15]. We consider both high- and low-scale [16–20] mordial GWs, potentially measurable at the currently realizations, with explicit or spontaneous lepton number planned GW interferometers, may represent an impor- violation, in which SU(3)c SU(2)L U(1)Y singlet neu- tant cutting-edge probe for new physics. This follows ⊗ ⊗ trinos act as neutrino mass mediators. Besides oscillation from the fact that these interferometers can be sensitive enough to measure the echoes of the possible First Order Phase Transitions (FOPTs), which might have happened in the past cosmological history [28]. ∗ [email protected] † [email protected] In this letter, we focus on possible gravitational foot- ‡ [email protected] § [email protected] prints of the various variants of the popular type-I and ¶ rahulsri@ific.uv.es inverse seesaw mechanisms for Majorana neutrinos. The ∗∗ valle@ific.uv.es relevant part of the minimal type-I seesaw Lagrangian is 2 given by majoron variants of minimal type-I and inverse seesaw Type−I c c c one should replace = Yν LHν¯ + Mν ν + h.c. (1) LYuk T Here, L = (ν, l) are the SM lepton doublets, H is the M Yσ vσ/√2 , µ Yσ vσ/√2 (5) → → SM Higgs doublet, νc are the three SM singlet “right- handed” neutrinos. The 3 3 matrices Yν and M are the in Eq. (17) and (3), respectively. An additional attractive × Yukawa coupling and the νc mass matrix, respectively. feature of majoron models is the existence of a pseudo Due to the Pauli principle the latter is symmetric. No- Nambu-Goldstone boson commonly dubbed as majoron. tice that, for brevity, we omit family indices through- The latter carries odd charge under Z2 thus providing out this letter. Notice also that the mass term explic- a good [39–41], and testable [42, 43] dark matter can- didate. In the standard majoron seesaw schemes, the itly breaks the lepton number symmetry U(1)L to its Z2 subgroup. The electroweak (EW) symmetry is broken by majoron mass is considered to be small, of order keV, for the vacuum expectation value (vev) of the Higgs field, i.e. it to be a suitable warm dark matter candidate. How- H = vh/√2, generating the light neutrino masses ever, the current stringent constraints on invisible Higgs h i 2 boson decay modes [26, 27, 44, 45] put severe limita- Type−I vh T −1 mν = Yν M Yν . (2) tion on how large h σσ coupling could be. On the 2 other hand, the strength→ of the cosmological phase tran- The lightness of the left-handed neutrinos is then as- sition is expected to be strongly correlated with the size cribed to the heaviness of the “right-handed” isosinglet of the Higgs-majoron quartic coupling. It remains an partners e.g. for Y (1),M (1014) GeV, one gets ∼ O ∼ O open question whether it is possible to reconcile the cur- mν (0.1) eV. ∼ O rent LHC bounds on the invisible Higgs decays in the in- Another popular realization of this idea is the “low-scale” verse seesaw scenario featuring a keV-scale majoron dark variant, in which two gauge singlet fermions νc and S are matter with the existence of strong EW FOPTs yielding added sequentially to the SM particle content [16–20]. potentially observable GWs signals. The template of these schemes has exact conservation of lepton number and, as a result, strictly massless neutri- In this work, we consider the case where the lepton num- nos. Yet flavor is violated to a potentially large degree, ber symmetry is broken also explicitly, but softly. This subjected only to constraints from weak interaction pre- way the majoron can pick up a mass. If this is larger cision observables, such as universality tests [29–33]. To than a half of the Standard Model (SM) Higgs boson this template one adds a small seed of lepton number vi- mass, mh/2 62.5 GeV, the invisible Higgs decays will ' olation, leading to nonzero neutrino mass. One example be kinematically forbidden, avoiding the stringent con- is the so-called “inverse seesaw” mechanism, where the straints on the Higgs-majoron quartic coupling. In such smallness of the neutrino mass is linked to the breaking scenario the heavy majoron-like state can provide a vi- able candidate for cold dark matter (CDM). The same of the lepton number symmetry U(1)L to its Z2 sub- group, through the so called µ-term. The relevant part type of scalar CDM scenarios with exact Z2 parity have of Lagrangian in this case is given by been broadly studied in singlet extensions of the SM [46– 56], which share very similar properties with the scalar Inverse c c = Yν LHν¯ + Mν S + µSS + h.c., (3) LYuk sector of the model under consideration. where µ is also a 3 3 symmetric matrix. The light neutrino mass is then× given by 2 Inverse vh T T −1 −1 m = Y M µM Yν . (4) ν 2 ν II. SCALAR SECTOR Note that small neutrino masses are “protected”, since mν 0 as the lepton number symmetry gets restored by → having µ 0 [16–20]. In this case there can be sizable The scalar sector of the majoron inverse seesaw model unitarity violation→ in neutrino propagation [34–36]. has been extensively studied in the literature, including the perturbative unitarity and stability of the scalar po- For both high- and low-scale seesaw, one can have spon- tential, as well as the electroweak precision tests and the taneous breaking of , leading to the so-called U(1)L Z2 bounds on Higgs-majoron couplings [27, 57]. The scalar majoron variants of the seesaw→ [17, 37, 38]. This is ac- potential is written as follows complished by adding the SM singlet scalar σ, which car- ries two units of lepton number charge. Then σ vσ 2 † † 2 2 † † 2 h i ≡ 0(Φ, σ) = µΦΦ Φ + λΦ(Φ Φ) + µσσ σ + λσ(σ σ) spontaneously breaks U(1)L 2, leading to a dynami- V Z † † 1 2 2 → +λΦσΦ Φσ σ + µ σ + h.c. , (6) cal explanation of the small neutrino masses. To get the 2 b 3 with Φ and σ given by model parameters. In the numerical analysis of GW spec- tra below we do not enforce the majoron mass mA to be 0 1 G + iG 1 small (e.g. at a keV scale) and treat it as a free parameter Φ = , σ = (φσ + σR + iσI ) , √2 φh + h + iη √2 instead. (7)
0 where h, η, G, G , σR, σI are real scalars. These latter III. GRAVITATIONAL WAVES FROM FOPTS fields represent quantum fluctuations about the classi- cal mean-fields φα = φh, φσ . In the zero-temperature { } limit, the mean-fields approach the corresponding vevs, In order to characterize the features of the GWs origi- i.e. φh,σ(T = 0) vh,σ, where vh = 246 GeV is nating from FOPTs in seesaw schemes, we calculate the ≡ the SM Higgs vev. Besides, one of the physical CP- strenght of the phase transition α at the bubble nucle- even scalar states, with a small or no mixing with σR, ation temperature Tn typically defined through the trace is identified with the SM-like Higgs boson with mass, anomaly as [58, 59] mh1 mh = 125 GeV. The last (soft) term appearing in ≡ 1 h Tn ∂Vi ∂Vf i Eq. (6) implements the explicit breakdown U(1)L Z2, → α = Vi Vf , (12) and hence provides a pseudo-Goldstone mass to the imag- ργ − − 4 ∂T − ∂T inary part of the field σ known as majoron 1. I with
In what follows, we discuss further implications of the 2 π 4 majoron seesaw scenario in both versions of the scalar ργ = g∗ T (13) 30 n sector, with explicit (vanishing vσ) and with sponta- 2 neous (vσ = 0) lepton-number U(1)L Z2 symmetry being the energy density of the radiation medium at the 6 → breaking, for physics of cosmological EW FOPTs and ex- bubble nucleation epoch found in terms of the number amine the associated GWs spectra. In the first version, of relativistic d.o.f.’s. g∗ 106.75 [60–63]. Above, Vi ' no mixing occurs so that and Vf are the values of the effective potential in the 2 symmetric and broken phases just before and after the 2 2 2 2 2 λσhvh m = 2λhv , m = µ + µ + , (8) h1 h h2 σ b 2 transition takes place, respectively. Another key quantity 2 to calculate the GW spectrum is the inverse time-scale 2 2 2 λσhvh mA = µσ µb + , (9) of the phase transition, which, in units of the Hubble − 2 β parameter H, reads as for the SM Higgs boson, CP-even and CP-odd (majoron) scalars, respectively. In the second version, the physical ! β ∂ Sˆ3 CP-even states acquire masses, = Tn , (14) H ∂T T 2 2 Tn 2 2 2 λσvσ λhvh m = λhv + λσv − , (10) h1,h2 h σ ∓ cos 2θ where Sˆ3 is the Euclidean action. In this work, we do in terms of h-σR mixing angle θ, while the majoron gets not consider the case of runaway bubbles and use the a pseudo-Goldstone mass, formalism outlined in Ref. [63] to calculate the spectrum m2 m2 = 2µ2 , µ2 < 0 . (11) of primordial GWs. A ≡ σI − b b While in the latter case, a very light majoron (compared For the case of non-runaway nucleated bubbles, the inten- to the EW scale vh) would imply setting an equally small sity of the GW radiation grows with the ratio ∆vn/Tn, µb parameter, in the former case, Eq. (9), this limit re- where lies on a very strong fine tuning between µσ, µb and λσh ∆vφ = vf vi , φ = h, σ (15) n | φ − φ| defines the difference between the VEVs of the initial 1 One could also add other other explicit breaking terms such as (metastable) and final (stable) phases at the bubble nu- † σΦ Φ but here we stick to the simplest possibility of mass terms. cleation temperature Tn. The quantity ∆vn/Tn also of- 2 More precisely, here we refer to a mixture of explicit and spon- fers a measure of the strength of the phase transition taneous U(1)L breakings in the scalar sector of the model since alternative to α. However, while the latter is of com- the soft scalar mass term µb 6= 0 is always present in the scalar potential (6) to provide a non-zero mass to the physical majoron mon use in the context of GWs, the former is typically state. This should not be confused with the neutrino sector where referred to in the context of electroweak baryogenesis. In U(1)L is always spontaneously broken by a small or large vσ vev this work, we will consider both quantities on the same in the considered majoron versions of the seesaw schemes. footing. 4
From the discussion in Ref. [59, 64], it follows that bubble sizable FOPTs at the quantum level since the fermion wall collisions do not provide an efficient way of produc- thermal loop contributions generate highly suppressed ing GWs in the models of interest to us here. As a result, terms to the effective potential in the high-temperature GWs originate mainly from two sources: expansion. Indeed, we find relatively weak FOPTs for many points in parameter space in this case and the cor- I. Magnetohydrodynamic (MHD) turbulence; 2 peak responding GW “intensity” parameter h ΩGW lies far be- II. Sound shock waves (SW) of the early Universe low the sensitivity of any conceivable experiment3. plasma, generated by the bubble’s violent expan- sion. In the majoron inverse seesaw scenarios, i.e. when the scalar sector of the SM is extended by incorporating an These contributions arise over transient times in the early additional complex SM-singlet scalar state, the situation Universe and get subsequently “redshifted” by the ex- changes dramatically. Indeed, it is well known that the pansion. To a present observer this appears as a cosmic presence of additional SM scalar singlets significantly en- gravitational stochastic background. Intuitively, one ex- hances the strength and multiplicity of the FOPTs, and pects that from any of these leading order contributions, in some cases leads to potentially detectable GW spectra a high wall velocity is necessary to generate detectable (see e.g. Refs. [68–72]). The presence of an additional GWs. In our numerical analysis, performed with the help classical field φσ = 0 coupled substantially to the Higgs of the CosmoTransitions package [65], we have consid- boson strongly affects6 the shape of the effective potential ered supersonic detonations such that the bubble wall at non-zero temperature T allowing for a richer pattern of velocity maximizes the GW peak amplitude and is above EW FOPTs. Note, in the considered majoron seesaw this the Chapman-Jouguet velocity defined as happens in both variants with explicit and spontaneous q U(1)L breaking at T = 0 discussed above in Sec.II. 1 2 2 vJ = cs + α + α , (16) 1 + α 3 The values of the coupling constants used in our numer- ical analysis satisfy the conservative bounds provided by with √1 being the speed of sound. Besides, in cs = 3 tree-level perturbativity, λi < 4π and Yi < √4π, for our results the SW contribution dominates the peak fre- the quartic and Yukawa couplings,| | respectively| | 4. Since quency and the peak amplitude. Furthermore, the state the tree-level potential also receives both quantum and of the art expressions derived in Ref. [63] do not account finite temperature corrections, we only considered values for MHD-turbulence effects due to large theoretical un- certainties. Therefore, in the remainder of this work, 9 we will not take such effects into consideration. Note, 10− LISA for certain parameter configurations one also expects se- 0.1325 12 quential phase transition patterns leading to multi-peak 10− BBO GWs spectra studied for the first time in Refs. [66, 67]. 0.1300 peak GW
15 DECIGO-corr σh
Ω 10− λ 2
h 0.1275 u-DECIGO A. Seesaw-induced GWs spectra 10 18 − 0.1250
u-DECIGO-corr 21 To investigate the eventual occurrence of phase transi- 10− 6 4 2 0 10− 10− 10− 10 tions, the standard way is to incorporate into the ef- fpeak [Hz] fective potential the tree-level zero temperature compo- nents, the Coleman-Weinberg corrections, the full one- FIG. 1: The GW spectrum as a function of λσh for the case loop finite-temperature corrections, as well as the Daisy of spontaneously broken U(1)L symmetry. No solutions consistent re-summation. It is worth noticing that within the type-I with the LHC bound on invisibly-decaying Higgs were found. Other model parameters are fixed as: m = 1 keV, m = 591 GeV, seesaw mechanism with explicitly broken lepton number, A h2 vσ = 858 GeV, Yσ,1 = 1.20, and Yσ,2 = 1.66. no FOPTs are obtained. The heavy isosinglet neutrinos practically decouple at the EW scale, and do not alter the nature of the EW phase transition. In contrast, in the non-majoron inverse seesaw mechanism (i.e. without 3 We emphasize that variations of the Yukawa coupling Yν in the adding a SM singlet scalar) the singlet neutrinos lie closer range 1 ÷ 10 cannot allow the detectability of the GW signal in to the EW scale, and can have a sizable coupling to the the type-I seesaw scenario. Higgs boson. However, even by adding a large number of 4 For a more involved analysis of the perturbativity constraints in singlet neutrino species it is impossible to generate any a scalar extension of the SM, see e.g. Ref. [73] 5 within the ranges λi < 5 and Yi < 3.5, even more U(1)L breaking scenario at zero temperature, with van- | | | | conservative bounds in a zero-temperature theory. ishing vσ providing a tiny value of µ in the inverse seesaw mechanism according to Eq. (5). We find a large num- Due to the current LHC constraints on invisible Higgs ber of points with strong FOPTs that generate the GW decays [26, 27, 44, 45] one has a bound λ < 0.03 the σh peak-amplitudes well within the projected LISA sensitiv- Higgs-majoron quartic coupling in the case of∼ light keV- ity, with typical values of Y below unity. scale majoron. Under this assumption in our numerical σ scan we did not find any point with an EW FOPT that 3 is strong enough for a potential observability of the re- LISA 10 sulting GW spectra. This is illustrated in Fig1. 10− 2 There is a strong correlation of the peak-amplitude with BBO 14 peak GW λσh value such that requiring the latter to be very small 10− σ Y Ω DECIGO-corr 2 makes the GW signals well below the reach of LISA or h u-DECIGO 18 1 even the planned BBO and DECIGO missions. Note, this 10− is the case for both considered versions of the majoron
u-DECIGO-corr inverse seesaw model, with explicit and spontaneous lep- 22 10− 6 4 2 0 ton number symmetry breaking in the scalar sector at 10− 10− 10− 10 T = 0. Therefore, we conclude that the standard keV- fpeak [Hz] scale (warm) majoron dark matter scenario associated with the inverse seesaw mechanism cannot be probed by FIG. 2: The GW spectrum as a function of the Yukawa Yσ cou- the GW astrophysics in the current simplest formulation. pling in the case of softly-broken U(1)L symmetry (i.e. vσ = 0). Order one variation of Yσ correspond to several order of magnitude For this reason, from now on we only consider the case variations in the GW power spectrum. Other model parameters are of heavy majoron kinematically closing the mA > mh/2 fixed as λσ = 0.37, λσh = 1.07, M = 239.4 GeV, mh2 = 154.6 GeV invisible Higgs decay channel and thus enabling us to and mA =369.9 GeV. consider larger values of λσh that ensure the existence of the strong FOPTs in the model under consideration. As is typical in models with several scalars, due to the presence of two classical fields φh, φσ in the effec- tive potential at finite temperatures,{ besides} a plenty of single-step FOPTs one also finds double-step and, B. Inverse seesaw with majoron: small vσ case in some rare cases, even triple-step FOPTs for a given parameter space point and at well-separated nucleation Let us now consider the case of a genuine inverse seesaw temperatures. One could naturally expect the presence with majoron and very small singlet VEV vσ, effectively of several peaks in the corresponding GW spectrum as- generating an equally small µSS term in the Lagrangian sociated with each FOPTs in such a chain of transitions. (3). As discussed above, this mechanism offers a dynami- Notice also that the analyses of GW spectra for the cal explanation for light neutrino masses, whose scale can multi-step FOPT scenarios involving EW phase transi- be attributed to the (tiny) scale of spontaneous breaking tions require particular care, as discussed for instance in of the lepton-number U(1)L symmetry in the neutrino Ref. [78]. sector (while being softly-broken in the scalar sector). Multi-peak configurations in the GW spectrum occur In order to understand the role of heavy neutrino in the very frequently in the inverse seesaw with majoron. This generation of the GW spectra in the majoron inverse see- fact is further highlighted in Fig.3. Indeed, the double- saw scenario, we study the sensitivity of the GW peak- peak feature of the GW spectrum is a generic prediction 2 peak amplitude h ΩGW originated by the EW FOPTs with of our model, that can arise for many parameter choices, respect to the variation of the Yukawa coupling Yσ. As as shown in Fig. 3a. Configurations with larger peak shown in Fig.2 an order one variation in the Yukawa cou- multiplicities are also possible, as seen in Fig. 3b, where pling reflects into violent variations by several orders of the color denotes the peak number, 1 (blue), 2 (cyan) magnitude in the GWs spectrum, with all the other pa- or 3 (red). Such a rich structure of the GW spectrum rameters fixed. The results are shown together with the is favoured for relatively large quartic couplings involv- projected sensitivities expected in LISA, and the planned ing σ. From Fig. 3b we also see that the GW spectra u-DECIGO and BBO missions [62, 63, 74, 75]. We have with three peaks are rarer than single or double GW- taken the u-DECIGO sensitivity curves from Ref. [76], peak spectra. We also note that a significant fraction whereas the sensitivities of other experiments are taken of the single peak cases are potentially testable at LISA from Ref. [77]. Here, for simplicity we have used the soft and BBO. However, finding a well-resolved double- or 6
double-peak signatures, these scenarios may suffer from 13 10− LISA larger underlined theoretical uncertainties so we decided BBO not to discuss them here. 16 DECIGO-corr 10− u-DECIGO 7 10− peak GW
Ω 19 u-DECIGO-corr 9 2 10− 10− h 11 10− 22 BBO 10− 13 LISA 10− DECIGO-corr
GW 15
25 Ω 10− 10− 3 2 1 0 1 2 10− 10− 10− 10 10 h 17 10− fpeak [Hz] u-DECIGO 19 10− (a) Selected double-peak scenarios within the LISA and BBO (T , ∆vh, ∆vσ) = (122.0, 23.0, 0.0), (92.0, 23.0, 0.0) GeV 21 n n n sensitivity ranges. The two ends of each line represent the h σ u-DECIGO-corr 10− (Tn, ∆vn, ∆vn) = (77.0, 235.0, 114.0), (183.0, 235.0, 114.0) GeV h σ location of the peaks of the double-peak GW spectrum. The (Tn, ∆vn, ∆vn) = (293.0, 0.0, 2941.0) GeV 23 two maxima in each double-peak GW spectra are joined by a 10− 5 4 3 2 1 0 10− 10− 10− 10− 10− 10 straight line, in order to easily identify the peaks associated f [Hz] with each other. FIG. 4: Inverse-seesaw-with-majoron benchmark GW spectra in 9 LISA 10− peak 1 the scenario with soft U(1)L symmetry breaking (i.e. with vanish- ingly small vσ(T = 0) → 0). The green curve represents the case BBO peak 2 with single-step FOPT, while the blue and red curves correspond 16 DECIGO-corr peak 3 10− u-DECIGO to double-step FOPTs whose characteristics are given in TablesI andII. For double-step transitions, the (Tn, vn) pairs in each peak u-DECIGO-corr peak GW h f i 23 are ordered from low to high frequencies. Here, ∆vn = |v − v | Ω h h
2 10− σ f i
h and ∆vn = |vσ − vσ| at a given nucleation temperature T = Tn.
30 10− i i f f Peak Id Tn vh, vσ → vh, vσ α β/H 37 Green 1 293 (0, 0) → (0, 2941) 0.5 4.9 10− 7 4 1 2 5 10− 10− 10− 10 10 Red 1 183 (0, 114) → (235, 0) 7.7 × 10−4 7.2 × 104 fpeak [Hz] Red 2 77 (0, 114) → (235, 0) 0.1 231 Blue 1 122 (193, 0) → (216, 0) 1.4 × 10−2 3.1 × 103 (b) Scatter plot showing the number of peaks for given model −3 4 parameter choices. Notice the appearance of double- and even Blue 2 92 (193, 0) → (216, 0) 9.4 × 10 3.0 × 10 triple-peak features. TABLE I: Phase transition parameters for the three curves in Fig.4. In “peak Id” column, the numbering of multi-step scenarios FIG. 3: The multi-peak feature arising from different phase tran- is ordered from high to low nucleation temperature Tn, given in sitions in the cosmological history of the Universe is very generic i units of GeV. The vevs before the transition (vh,σ) and after the in the inverse seesaw with majoron. f transition (vh,σ) are given in units of GeV. triple-peak feature where more than a single peak could Curve mh2 mA λσh λσ M Yσ be observable in a measurement appears to be very chal- Green 236 708 1.7 5 × 10−3 380 2 lenging. No single point with such a feature has been Red 192 970 2.3 1.5 93 0.1 found by our scans. As seen in Fig. 3a in most cases both Blue 325 169 4 2.7 158 0.1 peak-amplitudes occur below the projected sensitivities of any planned measurements rendering the observability TABLE II: Model parameters for the three curves in Fig.4. Mass of the corresponding scenarios very remote. In a subset parameters are given in units of GeV. of cases, only the largest peak has been found in the sen- sitivity ranges of the planned BBO and u-DECIGO mis- In order to understand the characteristic features of the sions, while the second peak typically lies outside of the GW spectra in the majoron inverse seesaw, in Fig.4 we reach of future detectors. This is a direct consequence of show the GW energy density spectrum obtained for dis- the not-so-large quartic couplings that are restricted to tinct nucleation temperatures Tn. Here, we have depicted be smaller than five. While larger values of the quartic three benchmark scenarios (for the case of vanishing vσ) couplings could generate well-separate and measurable originating from the single-step FOPT (green line) as well 7
7 10− as double-step FOPTs i.e. with two consecutive strong h σ (Tn, ∆vn, ∆vn) = (53.0, 119.0, 786.0), (134.0, 33.0, 3.0) GeV h σ 9 (Tn, ∆vn, ∆vn) = (82.0, 166.0, 693.0), (149.0, 0.0, 353.0) GeV phase transitions (red and blue lines). The correspond- 10− h σ (Tn, ∆vn, ∆vn) = (53.0, 115.0, 786.0), (136.0, 32.0, 2.0) GeV ing values for the model parameters are given in Tab.I 10 11 and Tab.II. − BBO 13 LISA 10− DECIGO-corr The green curve represents a single-peak scenario with a GW 15
Ω 10− 2
single very strong U(1)L phase transition, ∆vσ(Tn)/Tn h ' 17 10− 10, and softly-broken lepton symmetry. It features a very u-DECIGO
19 strong FOPT but not EW one. In fact, this peak is 10− a probe of U(1)L breaking at finite temperatures while 21 10 u-DECIGO-corr the associated EWPT is very weak or even of second − 23 10− 5 4 3 2 1 0 order. This is possible due to relatively large values of a 10− 10− 10− 10− 10− 10 f [Hz] quartic (portal) coupling, λσh 2, and majoron-neutrino ' Yukawa coupling, Yσ 2, which make the m/T -ratio ' (a) The expected GW spectra. sizable. Hence, the cubic (m/T )3 terms in the thermal expansion can produce a potential barrier between both LISA 12 vacua, inducing this type of transitions. 10− BBO
On another hand, the other two benchmark points, which 16 DECIGO-corr 10− u-DECIGO
have the double-peak feature, have strong EWPT for peak GW
Ω u-DECIGO-corr both peaks in the red curve and for the smaller, almost 2 20
h 10− invisible, peak of the blue curve. Besides this, the blue 24 curve has no U(1)L breaking in any of the minima of the 10− effective potential at the corresponding nucleation tem- 28 peratures whereas the red curve exhibits a strong U(1)L 10− 4 2 0 2 10− 10− 10 10 breaking for the higher peak and a weak U(1)L breaking fpeak [Hz] for the hidden peak. Similarly to the single-peak sce- nario, the observable peaks in the red and blue curves (b) Scatter plot showing typical double-peak scenarios. are generated due to large λσh and λσ couplings. Note, the single-peak case (green line) lies well within the LISA FIG. 5: Gravitational footprints of “fake” low-scale seesaw with range [62, 63, 74], while only the largest peaks in the majoron. In both plots we take vσ ∼ O(100) GeV – O(1) TeV. double-step FOPT scenarios are within the planned sen- sitivity range of the BBO (red line) and u-DECIGO-corr C. Type-I seesaw with majoron: large v (blue line) measurements [75], respectively. σ