The $N_{\rm eff}$ at CMB challenges $U(1)_{X}$ light gauge boson scenarios

The Cosmic Microwave Background (CMB) radiations play a crucial role in elucidating the dynamics of the early universe as well as shaping our present understanding of the large-scale structure of the universe. Observations of the CMB, including temperature anisotropies, polarization patterns, matter-energy distribution, reionization phenomena in the early universe, are consistent with both the standard $\Lambda$CDM cosmology and the Standard Model of particle physics (SM) [1]. One of the interesting cosmological parameters $N_{\rm eff}$, which represents the number of relativistic degrees of freedom at the CMB scale ($T\sim$ eV), serves an important role in understanding the dynamics of the thermal history during the epoch $T\sim$ MeV to eV. At high temperature of the universe, photons and neutrinos shared the same temperature. As the universe cooled down ($T\lesssim 2$ MeV) and the interaction rate fell below the Hubble expansion rate ($H$), neutrinos decoupled from photons, resulting in the formation of two separate baths i.e. neutrino bath and photon bath (combined with electrons). At temperature $T\lesssim m_{e}$, the two baths evolve with different temperature [1]. The parameter $N_{\rm eff}$ is then parameterized by the ratio of energy densities in the neutrino and photon baths. In the particle content of the SM, one can expect the neutrino decoupling event to take place around $T\sim 2$ MeV, resulting in $N_{\rm eff}^{\rm SM}=3.046$ at the CMB [2, 3]. Note the slight deviation in $N_{\rm eff}$ from the expected three light neutrino degrees of freedom in SM ($3~{}\nu_{L}$) is due to non-instantaneous neutrino decoupling, finite temperature QED corrections, and neutrino flavor oscillations [2, 3].

Observations of the CMB by the Planck satellite can also measure the parameter $N_{\text{eff}}$ ¹¹1Measurements of baryon acoustic oscillations (BAO), large-scale structure (LSS) datasets, and other cosmological probes contribute to determining the value of $N_{\rm eff}$ [4, 5].. The latest Planck 2018 data provides a precise measurement of $N_{\rm eff}$ at the time of the CMB to be $2.99^{+0.34}_{-0.33}$ with $95\%$ confidence level (CL), under the assumption of standard $\Lambda$CDM cosmology [5]. The upper limit of the Planck data suggests an additional contribution (apart from SM) to $N_{\rm eff}$, extending up to $0.28$ ($\Delta N_{\rm eff}=N_{\rm eff}^{\rm obs}-N_{\rm eff}^{\rm SM}$), providing the hint for BSM physics.

The presence of beyond standard model (BSM) light mediators, interacting with either the photon or neutrino bath at temperatures relevant to neutrino decoupling, introduces additional contributions to the radiation energy density, leading to an increase in $N_{\rm eff}$ [6, 7, 8] ²²2Note additional contributions to $N_{\rm eff}$ may also come from various other sources like extra radiation [9, 10, 11, 12], models involving early dark energy [10, 13], and relativistic decaying dark matter [14].. The observations of $N_{\rm eff}$ impose constraints on the BSM light mediators, as presented by the red shaded region in the cartoon Fig.1. In certain mass regions of these mediators, $N_{\rm eff}$ turns out to be a severe constraint compared to other existing bounds. $U(1)_{X}$ models are among such widely studied BSM paradigms with light mediators and their phenomenological consequences have been explored in several contexts³³3See for example Ref. [15, 16, 17, 18, 19, 20, 21, 22] and the references therein.. In this work, we consider light gauge boson ($Z^{\prime}$) originating from different types of $U(1)_{X}$ gauge extended scenarios, which naturally involve interactions with neutrinos and electron-positron pairs (either at the tree level or loop-induced). The BSM interactions in the presence of light $Z^{\prime}$ can potentially modify the late-time dynamics between the photon and neutrino baths, prolonging the process of neutrino decoupling and increasing the $N_{\rm eff}$ value. As a result, the upper limit on $N_{\rm eff}$ ($3.33$ with $95\%$ CL) derived from the CMB can be utilized to constrain both the mass of the light gauge boson ($M_{Z^{\prime}}$) and its corresponding gauge coupling ($g_{X}$).

In the context of $N_{\rm eff}$, various types of gauged $U(1)_{X}$ scenarios have been explored in the existing literature. These scenarios aim to explain the excess of $N_{\rm eff}$ observed in the CMB, involving $U(1)_{L_{\mu}-L_{\tau}}$ [7], $U(1)_{B-L}$ [8] or, generic $U(1)_{X}$ model [23]. A comprehensive discussion on $N_{\rm eff}$ within the generic $U(1)_{X}$ scenario is presented in Ref. [23], categorizing $U(1)_{X}$ models into two classes: $(i)$ those with tree-level coupling of $Z^{\prime}$ with electron and $(ii)$ those without tree-level coupling of $Z^{\prime}$ with electron. The study also highlights that the varying $U(1)_{X}$ charges for $e,\nu_{e,\mu,\tau}$, corresponding to various gauged extended scenarios, notably impact the value of $N_{\rm eff}$. This work is an extension of Ref. [23], where we consider a variety of $U(1)_{X}$ models, such as $U(1)_{B-L}$, $U(1)_{L_{i}-L_{j}}(i\neq j)$, $U(1)_{B-3L_{i}}$, and $U(1)_{B_{i}-L_{j}}$ etc., and study their impact on $N_{\rm eff}$, assuming the scenario where $Z^{\prime}$ was initially ($T>M_{Z^{\prime}}$) in thermal bath (which requires $g_{X}\gtrsim\mathcal{O}(10^{-9})$). The primary objective of this work is to identify the exclusion regions in the $M_{Z^{\prime}}-g_{X}$ plane for each gauged extended model, guided by the upper bound on $N_{\rm eff}$ derived from CMB observations [5].

The gauged extended models are also motivated by their potential detectability across a broad mass range, from eV to TeV, as sketched in Fig.1. For gauge boson $Z^{\prime}$ with TeV or sub TeV masses, the most severe constraints are derived from collider experiments like LHC [24, 25, 26, 27] and LEP [28]. Whereas, in the mass range spanning from a few MeV to a few GeV, experiments carried out by LHCb [29], Belle[30], BaBar [31], and fixed target experiments [16, 32] are notably more effective in providing constraints. In the sub GeV mass region, with $g_{X}\gtrsim\mathcal{O}(10^{-6})$, the low energy neutrino scattering experiments (neutrino electron scattering [18], neutrino-nucleus scattering [33] etc.), provide most stringent constraints. However, in the mass range where $M_{Z^{\prime}}\lesssim\mathcal{O}$(MeV) and $g_{X}\lesssim\mathcal{O}(10^{-6})$, direct searches become less sensitive. Furthermore, as we will show, CMB observations in terms of $N_{\rm eff}$ play a significant role in constraining the parameter space for $M_{Z^{\prime}}\lesssim\mathcal{O}({\rm MeV})$. In this work, we identify the CMB exclusion region with $95\%$ CL for each $U(1)_{X}$ scenario and compare it to the other aforementioned existing constraints in the $M_{Z^{\prime}}-g_{X}$ plane. We highlight that in certain regions of $M_{Z^{\prime}}$ and coupling $g_{X}$, observations of $N_{\rm eff}$ in the CMB provide more stringent constraints compared to other observations. We present a schematic diagram in Fig.1, where various existing constraints are summarised by different color patches based on the mass and coupling. Our key focus in this study is the red-shaded region, which represents the exclusion region derived from $N_{\rm eff}$ observations in the CMB.

The paper is organized in the following way. In section 2, we provide a brief overview of the evolution of $N_{\rm{eff}}$ in the presence of a light $Z^{\prime}$, with a focus on its dependence on the gauge extension. Then we present the exclusion regions in the $M_{Z^{\prime}}-g_{X}$ plane for various $U(1)_{X}$ models, determined by the upper bound on $N_{\rm{eff}}$ derived from CMB observations in section 3. Finally, in section 4 we summarize our results and conclusions.

Before going to the $N_{\rm eff}$ analysis of light $Z^{\prime}$ emerging from various $U(1)_{X}$ symmetries we present a brief and generic overview of the effective $Z^{\prime}$ models in this section. The discussion in this section is based on our previous work [23] where we have presented a comprehensive analysis of $N_{\rm eff}$ for a generic $U(1)_{X}$ model. For a generic $U(1)_{X}$ model, we assume that SM quark doublets ($Q_{i}$) and singlets ($u_{i},d_{i}$), as well as lepton doublets ($L_{i}$) and singlets ($\ell_{i}$) are charged under this $U(1)_{X}$ symmetry⁴⁴4Note that for now we have left the $U(1)_{X}$ charges of SM particles as free parameters, whose values depend on the details of the $U(1)_{X}$ symmetry, see Table 1.. We also need a SM singlet scalar $\sigma$ carrying $U(1)_{X}$ charge to generate the mass of $Z^{\prime}$ ($M_{Z^{\prime}}$) by acquiring a vacuum expectation value (VEV) to break the $U(1)_{X}$ symmetry. As pointed out in Ref. [23] that to affect the $\nu_{L}$ decoupling $Z^{\prime}$ has to be in the MeV scale. To keep $M_{Z^{\prime}}$ very light compared to the SM $Z$ gauge boson ($M^{2}_{Z^{\prime}}\ll M^{2}_{Z}$), we consider $U(1)_{X}$ charge of the SM Higgs doublet $\mathbb{X}_{\Phi}=0$ and hence $M^{2}_{Z^{\prime}}$ does not get any contribution from the SM vev. The charge assignments are shown in Table 1.

Note that to generate the SM quark and charged lepton masses, we take the following $U(1)_{X}$ charge assignments: $\mathbb{X}_{Q_{i}}=\mathbb{X}_{u_{i}}=\mathbb{X}_{d_{i}}$ and $\mathbb{X}_{L_{i}}=\mathbb{X}_{\ell_{i}}$. It is also worth highlighting that although the $U(1)_{X}$ charges of quark doublet should be the same as that of the up and down quark singlets, the charges may differ across different generations i.e. $\mathbb{X}_{Q_{1}}\neq\mathbb{X}_{Q_{2}}\neq\mathbb{X}_{Q_{3}}$ as shown in Table 1.

Besides the particle mentioned in Table 1, we also need right-handed neutrino (RHN) $\nu_{R}$ for anomaly cancellation and light neutrino mass generation [34, 23, 35, 15]. In this work, we consider only Majorana type mass models where, $\nu_{R}$ are also too heavy to affect $\nu_{L}$ decoupling [36, 37] and we ignore their contribution⁵⁵5In Dirac-type mass models, $\nu_{R}$ are relativistic at MeV temperature and can significantly alter $N_{\rm eff}$ [9].. Apart from the minimal particle content of Table 1, most of the popular $U(1)_{X}$ models studied in literature may also contain additional BSM particles which are typically much heavier than MeV scale and will not take part in our analysis [15].

The $\nu_{L}$ decoupling is governed by only the weak interaction processes in SM scenario and usually takes place at around $T\sim$ MeV temperature [6]. The only particles relevant at that temperature are $e^{\pm},~{}\gamma$ and $\nu_{i}~{}(i\equiv e,\mu,\tau)$. The scenario changes drastically in the presence of any light (mass $\sim\mathcal{O}(1)$ MeV) BSM particles interacting with either of these three sectors [7, 23]. Thus one can easily anticipate that light $U(1)_{X}$ gauge boson $Z^{\prime}$ (with mass, $M_{Z^{\prime}}\sim\mathcal{O}(1)$ MeV) interacting with $\nu_{i}$ and $e^{\pm}$ will significantly alter the $N_{\rm eff}$ at the time of CMB. In this work we will assume all $3$ $\nu_{i}$ acquire the same temperature ($T_{\nu_{L}}$) and consider the case where $Z^{\prime}$ was initially in a thermal bath ($g_{X}\gtrsim 10^{-9}$) [23]. $N_{\rm eff}$ is parametrised in terms of the temperature ratios of $\nu_{i}$ sector and $\gamma$ bath [1],

$N_{\rm eff}$ at the time of CMB is determined by computing the temperature ratio $\frac{T_{\nu_{L}}}{T_{\gamma}}$ at $T_{\gamma}=T_{\rm CMB}$. Note that, throughout this paper, whenever we say $N_{\rm eff}$, we refer $N_{\rm eff}$ at the time of CMB. To evaluate $T_{\nu_{L}}/T_{\gamma}$ we solve the following coupled equations ⁶⁶6Further details can be found in Appendix A-C of Ref. [23].,

where, $\rho_{r},P_{r}$ and $T_{r}$ stand for the energy density, pressure density, and temperature of species $r$. On the other hand, $\frac{\delta\rho_{a\to b}}{\delta t}$ signify the energy transfer rate from bath $a$ to $b$, deduced by integrating the collision terms [7, 23]. As mentioned earlier, assuming same temperature between 3 $\nu_{i}$ sectors, we denote $\rho_{\nu_{L}}=\sum_{i=e,\mu,\tau}\rho_{\nu_{i}}$, $\frac{\delta\rho_{\nu_{L}\to e}}{\delta t}=\sum_{i=e,\mu,\tau}\frac{\delta\rho_{\nu_{i}\to e}}{\delta t}$ and $\frac{\delta\rho_{Z^{\prime}\to\nu_{L}}}{\delta t}=\sum_{i=e,\mu,\tau}\frac{\delta\rho_{Z^{\prime}\to\nu_{i}}}{\delta t}$. Throughout the paper, by $\nu_{L}$ we refer all $3$ generation of $\nu_{i}~{}(i\equiv e,\mu,\tau)$ as a whole. We do not discuss the energy transfer rates here as they are already elaborated in detail in Ref. [23].

It is worth mentioning that the $\nu_{i}\bar{\nu_{i}}Z^{\prime}$ coupling and the $e^{+}e^{-}Z^{\prime}$ coupling are the most crucial parameters in the BSM processes that affect the $N_{\rm eff}$. However, interactions involving heavier charged leptons and quarks are insignificant here, as their energy densities are suppressed at temperatures nearing $\sim 1$ MeV. When $Z^{\prime}$ has coupling to both $\nu_{L}$ and $e^{\pm}$ ($\mathbb{X}_{1}\neq 0$) , both the following BSM processes affect $\nu_{L}$ decoupling: (i) $Z^{\prime}$ decaying to both $e^{+}e^{-}$ and $\nu_{i}\bar{\nu_{i}}$ and (ii) scattering process $\nu_{i}\bar{\nu_{i}}\to e^{+}e^{-}$ mediated by $Z^{\prime}$. Any increase in the effective coupling ($\mathbb{X}_{1}g_{X}$) dictating these BSM processes enhances the BSM contribution, thus leading to an increment in $N_{\rm eff}$. On the other hand in the absence of tree level $Z^{\prime}e^{+}e^{-}$ vertex $(\mathbb{X}_{1}=0)$, all the number density of $Z^{\prime}$ finally gets diluted to $\nu_{L}$ bath (with $100\%$ branching ratio) irrespective of the coupling strength ($\mathbb{X}_{2/3}g_{X}$). Hence for such models ($\mathbb{X}_{1}=0$), the final $N_{\rm eff}$ becomes independent of effective couplings $\mathbb{X}_{2,3}g_{X}$. Thus following our argument in Ref. [23] the light $Z^{\prime}$ models can be broadly classified into two cases: $\mathbb{X}_{1}\neq 0$ and $\mathbb{X}_{1}=0$, depending on whether $Z^{\prime}$ has tree level coupling with electron or not (see Fig.2).

It is worth highlighting that even in the absence of tree level $Z^{\prime}e^{+}e^{-}$ interaction ($\mathbb{X}_{1}=0$), induced coupling with $e^{\pm}$ may exist which varies depending on the specific model [38]. We do not delve into the details of induced coupling which may arise from kinetic mixing or can be loop induced. We take the following effective $Z^{\prime}e^{+}e^{-}$ interaction Lagrangian:

with a benchmark value $\epsilon=-\frac{g_{X}}{70}$, which is very commonly used for $L_{\mu}-L_{\tau}$ models (with $\mathbb{X}_{1}=0$)[39, 7]. Needless to say that with the increase in $M_{Z^{\prime}}$, the energy density of $Z^{\prime}$ gets suppressed and hence $N_{\rm eff}$ decreases for both of the above mentioned two cases.

The excess in $N_{\rm eff}$ at CMB in the presence of $Z^{\prime}$ can be used to constrain the parameter space, as the deviation from SM predicted value of $N_{\rm eff}$ is also very precisely constrained from Planck 2018 [5]. We translate this constrain from $N_{\rm eff}$ in the parameter space of $M_{Z^{\prime}}-g_{X}$ plane for previously mentioned two classes of generic $U(1)_{X}$ models in our previous work [23]. We extend this exercise for all popular $U(1)_{X}$ models by performing exhaustive numerical scans on the parameter space (for $N_{\rm eff}$), in the next section.

In this section, we will discuss some phenomenological consequences of various well motivated $U(1)_{X}$ symmetries. As mentioned earlier, while doing this we consider that the corresponding gauge boson $Z^{\prime}$ was initially in thermal equilibrium and the RHNs are heavy enough that they do not take part in $\nu_{L}$ decoupling. In order to estimate $N_{\rm eff}$, we find the temperature ratio $T_{\nu_{L}}/T_{\gamma}$ at $T_{\gamma}=T_{\rm CMB}$ by solving eq.(2)-(4) and then plug into eq.(1). As described earlier the value of $N_{\rm eff}$ is strongly dependent on parameters: $g_{X}$ and $M_{Z^{\prime}}$, and hence the upper limit on $N_{\rm eff}$ by Planck 2018 [5] leads to strong constraint on the same parameter space. Before going into the details of the results we summarize the existing laboratory-based constraints on light $Z^{\prime}$ in Table-2. The detailed discussion on existing constraints can be found in Ref. [16, 17, 18, 19, 21, 20]. In addition to the constraints mentioned in the table, the light $Z^{\prime}$ may also encounter astrophysical constraints like SN1987A [40, 41, 42]. However, this supernova constraint differs depending on the specific gauge model and requires a detailed analysis which is beyond the scope of this work.