Lorentz invariance violation in the neutrino sector: a joint analysis from big bang nucleosynthesis and the cosmic microwave background

TL;DR

This paper constrains Lorentz invariance violation in the neutrino sector by combining big bang nucleosynthesis and cosmic microwave background data, revealing strong limits on possible deviations from standard physics.

Contribution

It provides the first joint analysis of BBN and CMB data to constrain Lorentz invariance violation in neutrinos, breaking degeneracies with dark matter density.

Findings

Primordial helium-4 abundance constrains Lorentz violation parameters.

Planck 2015 data combined with helium measurements yields strong limits.

Degeneracy between Lorentz violation and dark matter density is broken.

Abstract

We investigate constraints on Lorentz invariance violation in the neutrino sector from a joint analysis of big bang nucleosynthesis and the cosmic microwave background. The effect of Lorentz invariance violation during the epoch of big bang nucleosynthesis changes the predicted helium-4 abundance, which influences the power spectrum of the cosmic microwave background at the recombination epoch. In combination with the latest measurement of the primordial helium-4 abundance, the Planck 2015 data of the cosmic microwave background anisotropies give a strong constraint on the deformation parameter since adding the primordial helium measurement breaks the degeneracy between the deformation parameter and the physical dark matter density.

Tables

Table 1. Table 1: Mean values and marginalized 68% confidence level for the deformation parameter and other cosmological parameters, derived from the Planck data in combination with baryon acoustic oscillation data, the JLA sample of Type Ia supernovae and the measurement of the helium-4 mass fraction. For comparison with BBN consistency, fixing Y p = 0.24 subscript 𝑌 p 0.24 Y_{\rm p}=0.24 is also considered.

Model	Planck2015+JLA+BAO		Planck2015+JLA+BAO+$Y_p$
	$Y_p=0.24$ (fixed)	BBN consistency	BBN consistency
${\Omega }_b{h}^2$	$0.02127\pm 0.00046$	$0.02212\pm 0.00025$	$0.02227\pm 0.00021$
${\Omega }_c{h}^2$	$0.1124\pm 0.0061$	$0.1163\pm 0.0040$	$0.1197\pm 0.0018$
$100{\theta }_{\mathrm{MC}}$	$1.04181\pm 0.00098$	$1.04173\pm 0.00101$	$1.04090\pm 0.00050$
${\tau }_{\mathrm{re}}$	$0.071\pm 0.018$	$0.076\pm 0.018$	$0.077\pm 0.018$
$n_s$	$0.9583\pm 0.0069$	$0.9679\pm 0.0049$	$0.9659\pm 0.0044$
$\ln \left({10}^{10}{A}_s\right)$	$3.064\pm 0.038$	$3.084\pm 0.035$	$3.087\pm 0.035$
$H_0$(km $s^{-1}$${\mathrm{Mpc}}^{-1}$)	$65.92\pm 1.76$	$67.25\pm 1.00$	$68.01\pm 0.63$
${\sigma }_8$	$0.840\pm 0.014$	$0.846\pm 0.016$	$0.841\pm 0.015$
Age(Gyr)	$14.191\pm 0.35$	$13.937\pm 0.19$	$13.767\pm 0.068$
${\Omega }_{\Lambda }$	$0.691\pm 0.008$	$0.694\pm 0.007$	$0.693\pm 0.007$
${\Omega }_m$	$0.309\pm 0.008$	$0.306\pm 0.007$	$0.307\pm 0.007$
$\xi$	$0.101\pm 0.091$	$0.042\pm 0.049$	$-0.002\pm 0.014$

Equations

$$E^2=m^2+p^2+\xi p^2,$$

$$n_{\nu }$$

$${\rho }_{\nu }$$

$$P_{\nu }$$

$$X_i=\frac{n_i}{n_B},i=n,p,{}^{2}H,\cdots .$$

$$\dot{X_i}=\sum _{j,k,l}{N}_i\left({\Gamma }_{kl\to ij}\frac{X_l^{N_l}{X}_k^{N_k}}{N_l!N_k!}-{\Gamma }_{ij\to kl}\frac{X_i^{N_i}{X}_j^{N_j}}{N_i!N_j!}\right),$$

$$n+{\nu }_e$$

$$n+e^{+}$$

$$n$$

$$\Gamma =\left[1-\frac{3}{8}\xi -\frac{3\left(C_V^2-C_A^2\right)}{4\left(C_V^2+3C_A^2\right)}\xi \right]{\left(1+\xi \right)}^{-3/2}{\Gamma }^{\left(0\right)},$$

$$Y_{\mathrm{p}}^{\mathrm{FIT}}=0.2334+0.52{\Omega }_b{h}^2+0.258\xi ,$$

$$N_{\mathrm{eff}}=3.046(1+\xi {)}^{-3/2}.$$

$${\rho }_r=N_{\mathrm{eff}}\frac{7}{8}{\left(\frac{4}{11}\right)}^{4/3}{\rho }_{\gamma }.$$

$$\delta {\rho }_{\nu }$$

$$\delta P_{\nu }$$

$$({\rho }_{\nu }+P_{\nu }){\theta }_{\nu }$$

$$({\rho }_{\nu }+P_{\nu }){\sigma }_{\nu }$$

$$\Psi (\vec{k},\hat{n},q,\tau )=\sum _{\ell =0}^{\infty }{\left(-i\right)}^{\ell }\left(2\ell +1\right){\Psi }_{\ell }(\vec{k},q,\tau )P_{\ell }(\hat{k}\cdot \hat{n}),$$

$${\dot{\Psi }}_{\ell }+\left(1+\xi \right)\frac{q}{ϵ}\frac{k}{2\ell +1}\left[\left(\ell +1\right){\Psi }_{\ell +1}-\ell {\Psi }_{\ell -1}\right]$$

$$+\left({\delta }_{2\ell }\frac{1}{15}\dot{h}+{\delta }_{2\ell }\frac{2}{5}\dot{\eta }-{\delta }_{0\ell }\frac{1}{6}\dot{h}\right)\frac{d\ln f_0}{d\ln q}=0,$$

$${\Psi }_{{\ell }_{\max }+1}\approx \frac{\left(2{\ell }_{\max }+1\right)ϵ}{\left(1+\xi \right)qk\tau }{\Psi }_{{\ell }_{\max }}-{\Psi }_{{\ell }_{\max }-1}.$$

$$F(\vec{k},\hat{n},\tau )$$

$$=\sum _{\ell =0}^{\infty }(-i{)}^{\ell }(2\ell +1)F_{\ell }(\vec{k},\tau )P_{\ell }(\hat{k}\cdot \hat{n}).$$

$${\dot{\delta }}_{\nu }$$

$${\dot{\theta }}_{\nu }$$

$$2{\dot{\sigma }}_{\nu }$$

$${\dot{F}}_{\ell }$$

$$F_{{\ell }_{\max }+1}\approx \frac{\left(2{\ell }_{\max }+1\right)}{\sqrt{1+\xi }k\tau }{F}_{{\ell }_{\max }}-F_{{\ell }_{\max }-1}.$$

$$h$$

$$\eta$$

$${\delta }_{\nu }$$

$${\theta }_{\nu }$$

$${\sigma }_{\nu }$$

$$F_3$$

$${\Psi }_0$$

$${\Psi }_1$$

$${\Psi }_2$$

$${\Psi }_{\ell }$$

$$\{{\Omega }_b{h}^2,{\Omega }_c{h}^2,{\theta }_{\mathrm{MC}},{\tau }_{\mathrm{re}},n_s,A_s,\xi \},$$

Full text

I Introduction

Neutrino oscillations have shown that there are small but nonzero mass squared differences between three neutrino mass eigenstates (see Ref. GonzalezGarcia:2007ib and the reference therein), which imply the existence of physics beyond the standard model of particle physics. However, neutrino oscillation experiments cannot tell us the overall mass scale of neutrinos. Fortunately, cosmology provides a promising way to determine or constrain the total mass of neutrinos by the gravitational effects of massive neutrinos which can significantly change the CMB power spectrum due to lensing Kaplinghat:2003bh ; Lesgourgues:2005yv ; Lesgourgues:2012uu ; Ade:2015xua and the formation of large-scale structures Hu:1997mj ; Ringwald:2004np , therefore, alter the cosmic microwave background (CMB) anisotropies and the large-scale structure distribution of matter (see Wong:2011ip for a review). Recent Planck 2015 data combined with some low-redshift data give an upper limit on the total mass of neutrinos at $95\%$ confidence level, $\Sigma m_{\nu}<0.19$ eV Ade:2015xua .

Another possible signal of new physics is Lorentz invariance violation in the neutrino sector. The observed neutrino oscillations may originate from a combination of effects involving neutrino masses and Lorentz invariance violation Coleman:1997xq ; Coleman:1998ti ; Kostelecky:2003xn ; Kostelecky:2003cr ; Kostelecky:2011gq ; Motie:2012qj . Lorentz invariance is a fundamental symmetry in the standard model of particle physics. Although present experiments have confirmed Lorentz invariance to a good precision Diaz:2013wia , it may be broken in the early Universe when energies approach the Planck scale. The standard model itself is believed to be a low-energy effective theory of an underlying unified theory. Lorentz invariance violation has been explored in quantum gravity AmelinoCamelia:2002dx , loop quantum gravity Alfaro:1999wd , non-commutative field theory Carroll:2001ws , and doubly special relativity theory Magueijo:2001cr . Various searches for Lorentz invariance violation have been performed with a wide range of systems Colladay:1998fq ; Kostelecky:2008ts .

The cosmological consequences of neutrinos with the Coleman-Glashow type Coleman:1998ti dispersion relation have been studied in Ref. Guo:2012mv . As can be expected, the Lorentz-violating term affects not only the evolution of the cosmological background but also the behavior of the neutrino perturbations. The former changes the expansion rate prior to and during the epoch of photon-baryon decoupling, which alters heights of the first and second peaks of the CMB temperature power spectrum, while the latter alters the shape of the CMB power spectrum by changing neutrino propagation. Since these two effects can be distinguished from a change in the total mass of neutrinos or in the effective number of extra relativistic species, CMB data have been proposed as a probe of Lorentz invariance violation in the neutrino sector.

Moreover, since the Lorentz-violating term influences the abundances of the light elements by altering the energy density of the Universe and weak reaction rates prior to and during the Big bang nucleosynthesis (BBN) epoch, BBN provides a promising probe of Lorentz invariance violation in the neutrino sector in the early Universe Guo:2013zwa . In particular the BBN-predicted abundance of helium-4 is very sensitive to the Lorentz-violating term. It is well known that the primordial helium-4 abundance plays a non-negligible role at the recombination epoch because it has an influence on the number density of free electron. Therefore, Lorentz invariance violation affects the CMB power spectrum by changing the BBN-predicted abundance of helium-4. A reasonable way to test Lorentz invariance violation in the neutrino sector with CMB data is to take account of the BBN-predicted helium-4 abundance as a prior rather than fixing the helium-4 abundance.

In this paper we investigate constraints on Lorentz invariance violation in the neutrino sector from a joint analysis of BBN and CMB. In combination with the latest measurement of the primordial helium-4 abundance, the Planck 2015 data of the CMB anisotropies give a strong constraint on the deformation parameter. Adding the primordial helium measurement can break effectively the degeneracy between the deformation parameter and the physical dark matter density.

This paper is organized as follows. In Sect. II we parameterize Lorentz invariance violation in the neutrino sector. In Sect. III, we derive the weak reaction rates in the Lorentz-violating extension of the standard model and calculate the BBN prediction of the helium-4 abundance. In Sect. IV, we derive the Boltzmann equation for neutrinos in the synchronous gauge and calculate the CMB power spectrum. In Sect. V, we place constraints on the deformation parameter using the Planck 2015 data in combination with the latest measurement of the primordial helium-4 abundance. Sect. VI is devoted to conclusions.

II Deformed dispersion relation

At a phenomenological level, the deformed dispersion relation can be characterized by a power series of momentum. In this paper, we consider a simple instance of the dispersion relation for neutrinos constructed in Coleman:1998ti , which can be parameterized as follows:

[TABLE]

where $E$ is the neutrino energy, $m$ the neutrino mass, $p=\left(p^{i}p_{i}\right)^{1/2}$ the magnitude of the 3-momentum, and $\xi$ the deformation parameter characterizing the Lorentz symmetry violation.

In the spatially flat Friedmann-Lemaintre-Robertson-Walker metric with the scale factor $a$, the number density $n_{\nu}$, energy density $\rho_{\nu}$ and pressure $P_{\nu}$ for neutrinos with (1) are given by

[TABLE]

The phase space distribution for neutrinos is $f_{0}(q)=g_{s}\left[1+\exp\left(\epsilon/T_{0}\right)\right]^{-1}$, where $g_{s}=2$ is the number of spin degrees of freedom, $T_{0}$ the neutrino temperature today, $q$ the magnitude of the comoving 3-momentum, and $\epsilon=\sqrt{m^{2}a^{2}+\left(1+\xi\right)q^{2}}$ the comoving energy.

III BBN prediction

The abundances of the light elements produced during the BBN epoch depend on the competition between the nuclear and weak reaction rates and the expansion rate of the Universe. In the standard cosmological scenario and in the framework of the electroweak standard model, the dynamics of this phase is controlled by only one free parameter, the baryon to photon number density.

We consider $N$ species of nuclides whose abundances $X_{i}$ are the number densities $n_{i}$ normalized with respect to the baryon number density $n_{B}$,

[TABLE]

Their evolutions are ruled by the following Boltzmann equations,

[TABLE]

where $\Gamma$ denotes the reaction rate and $N_{i}$ is the number of nuclide $i$ involved in the reaction.

In the case of Lorentz invariance violation in the neutrino sector, the changes of the reaction rates of the following weak reactions have to be taken into account:

[TABLE]

which determine the neutron-to-proton ratio when the baryons become uncoupled from the leptons. The abundances of the light elements depend on the neutron-to-proton ratio at the onset of BBN. As an example, let us compute the reaction rate of the process (7) in the Lorentz-violating extension of the standard model Colladay:1998fq ; Lambiase:2005kb . The reaction rate is Guo:2013zwa

[TABLE]

where $C_{V}$ and $C_{A}$ are the vector and axial coupling of the nucleon, $\Gamma^{\left(0\right)}$ is the standard reaction rate per incident nucleon derived in Ref. Lopez:1998vk . The first prefactor on the right-hand side of Eq. (10) arises from the neutrino propagator and the $e\nu W$ coupling in the Lorentz-violating extension of standard model. The second prefactor comes from the statistical distribution for neutrinos. As for the other processes in (7)-(9), their reaction rates can be simply derived by properly changing the statistical factors and the delta function determined by the energy conservation for each reaction. Therefore, the corrections to the conversion rate of neutron into proton and its inverse rate are the same as in (10).

In order to calculate the abundance of the light elements produced during the BBN epoch, we modified the publicly available PArthENoPE code Pisanti:2007hk to appropriately incorporate the Lorentz-violating term in the neutrino sector. As shown in Guo:2013zwa , the BBN-predicted abundance of helium-4 is sensitive to the deformed parameter $\xi$ due to the fact that it is mainly determined by the neutron-to-proton ratio, which is related to the weak reaction rate, neutrino number density and expansion rate of the Universe. In our analysis we shall focus on the helium-4 abundance that influences the CMB power spectrum.

Although in our analysis we use an interpolation method to obtain the value of $Y_{\mathrm{p}}$ on a grid of points in $(\Omega_{b}h^{2},\xi)$ space (see Sect. V), it is possible to describe the dependencies of the BBN-predicted helium-4 abundances upon the two parameters by simple, linear fits which, over their ranges of applicability, are accurate to a few percent. Following Kneller:2004jz , we get the following linear fit for $Y_{\mathrm{p}}$ versus $\Omega_{b}h^{2}$ and $\xi$:

[TABLE]

over our adopted range in $\Omega_{b}h^{2}$ and $\xi$, as shown in Fig. 1. We note that this fit is a very good approximation over the adopted parameter ranges, which is useful in facilitating studies of the viability of various options for non-standard physics and cosmology, prior to undertaking detailed BBN calculations

Since the neutrino species share the temperature with the photons prior to neutrino decoupling, the prediction for the effective number of neutrino species is

[TABLE]

The total relativistic energy density $\rho_{r}$ in neutrinos is given in terms of the photon density $\rho_{\gamma}$ by

[TABLE]

IV CMB power spectrum

The Lorentz-violating term affects not only the evolution of the cosmological background but also the behavior of the neutrino perturbations. Following Ref. Ma:1995ey we shall derive the Boltzmann equation of the distribution function for neutrinos with (1) to linear order in the perturbations. In the synchronous gauge, the perturbed energy density $\delta\rho_{\nu}$, pressure $\delta P_{\nu}$, energy flux $\theta_{\nu}$, and shear stress $\sigma_{\nu}$ are, respectively, given by

[TABLE]

where $\Psi_{l}$ are the expansion coefficients of the perturbed neutrino distribution function expanded in a series of Legendre polynomials

[TABLE]

where $\tau$ is the conformal time. These perturbations $\Psi_{l}$ satisfy the following Boltzmann equations:

[TABLE]

in the synchronous gauge, where dots denote derivatives with respect to the conformal time, $h$ and $\eta$ are the two scalar modes of the metric perturbations in the Fourier space $k$. To avoid the reflections from high-$\ell$ equations by simply set $\Psi_{\ell}=0$ for $\ell>\ell_{\rm max}$, such a Boltzmann hierarchy is effectively truncated by adopting the following scheme Ma:1995ey :

[TABLE]

Making a transformation $q\rightarrow\sqrt{1+\xi}q$, we note that the number density, energy density and pressure of neutrinos are proportional to $(1+\xi)^{-3/2}$, which means that increasing $\xi$ decreases the number density, energy density and pressure of neutrinos. Moreover, there is a factor $\sqrt{1+\xi}$ in the second term of the Boltzmann equations, which comes from neutrino propagation.

For the massless neutrinos, the Boltzmann equations are simplified by setting $\epsilon=\sqrt{1+\xi}q$. The number of variables can be reduced by integrating out the $q$ dependence in the neutrino distribution function so that

[TABLE]

Using the orthonormality and recursion relation of the Legendre polynomials, from the Boltzmann equations we derive the evolution equations of the perturbed energy density $\delta_{\nu}=\delta\rho_{\nu}/\rho_{\nu}$, energy flux $\theta_{\nu}$, shear stress $\sigma_{\nu}$ independent of the momentum $q$,

[TABLE]

The truncation scheme for massless neutrinos is

[TABLE]

In the synchronous gauge the adiabatic initial conditions for the metric perturbations and massless neutrinos in the cold dark matter frame are

[TABLE]

where $C$ is a dimensionless constant determined by the amplitude of the fluctuations from inflation, $R_{\nu}\equiv\rho_{\nu}/(\rho_{\gamma}+\rho_{\nu})$, and $\omega\equiv\Omega_{m}H_{0}/\sqrt{\Omega_{\gamma}+\Omega_{\nu}}$, which corresponds to the matter contribution to the total energy density of the Universe. In the radiation-dominated era the massive neutrinos are relativistic. The initial condition of $\Psi_{0}$ for massive neutrinos is related to $\delta_{\nu}$. Then, using Eqs. (19) and ignoring the mass terms in the differential equations of high-order moment, we get the initial conditions of $\Psi_{\ell}$ for massive neutrinos:

[TABLE]

In order to compute the theoretical CMB power spectrum, we modified the Boltzmann CAMB code Lewis:1999bs to appropriately incorporate the Lorentz-violating term in the neutrino sector. As pointed out in Guo:2012mv , the effects of the Lorentz-violating term on the CMB power spectrum are distinguished from a change in the total mass of neutrinos or in the effective number of extra relativistic species Hu:1995fqa ; Bashinsky:2003tk ; Hamann:2007pi ; Hou:2011ec ; Guo:2012jn . Therefore, the measurements of the CMB anisotropies provide a cosmological probe of Lorentz invariance violation in the neutrino sector. Moreover, we calculate the matter power spectrum today for $\xi=-0.1,0,0.1$. From Fig. 2 we find that increasing $\xi$ can enhance the matter power spectrum at small scales. It is well known that variances in the helium-4 abundance modify the density of free electrons between helium and hydrogen recombination and therefore influence the CMB power spectrum. However, the effects of the Lorentz-violating term on the BBN-predicted helium-4 abundance were not considered in Guo:2012mv . In the next section, we shall put constraints on the deformation parameter from a joint analysis of BBN and CMB.

V Joint analysis from BBN and CMB

In our analysis we use a modified version of the publicly available CosmoMC package to explore the parameter space by means of Monte Carlo Markov chains technique Lewis:2002ah . Since the deformation parameter $\xi$ is nearly uncorrelated with the total mass of neutrinos when CMB data are used to detect the signature of Lorentz invariance violence Guo:2012mv , we only consider massless neutrinos with vanishing chemical potentials in our analysis. Our cosmological model is the spatially flat $\Lambda$CDM plus three types of massless neutrino with the deformed dispersion relation (1), which can be described by the following seven parameters:

[TABLE]

where $h$ is the dimensionless Hubble parameter defined by $H_{0}=100h$ km $\rm s^{-1}$ $\rm Mpc^{-1}$, $\Omega_{b}h^{2}$ and $\Omega_{c}h^{2}$ are the physical baryon and dark matter densities relative to the critical density, $\theta_{\mathrm{MC}}$ is an approximation to the ratio of the sound horizon to the angular diameter distance at the photon decoupling, $\tau_{\rm re}$ is the reionization optical depth, $n_{s}$ and $A_{s}$ are the spectral index and amplitude of the primordial curvature perturbations at the pivot scale $k_{0}=0.05$ Mpc*-1*.

We firstly run the PArthENoPE BBN code with two free input parameters, $\Omega_{b}h^{2}$ and $\xi$, to precalculate the primordial helium-4 abundance on a grid of points in $(\Omega_{b}h^{2},\xi)$ space. For each couple of parameter values the helium-4 abundances are then obtained from the grid by two-dimensional cubic spline interpolation, which is adopted to calculate the CMB power spectrum. This method is called “BBN consistency”.

We use the recently released Planck 2015 likelihood code and data, including the Planck low-$\ell$ likelihood at multipoles $2\leq\ell\leq 29$ and Planck high-$\ell$ likelihood at multipoles $\ell\geq 30$ based on pseudo-$C_{\ell}$ estimators Ade:2015xua . The former uses foreground-cleaned LFI $70$ GHz polarization maps together with the temperature map obtained from the Planck $30$ to $353$ GHz channels by the Commander component separation algorithm over 94% of the sky. The latter uses $100$, $143$, and $217$ GHz half-mission cross-power spectra, avoiding the galactic plane as well as the brightest point sources and the regions where the CO emission is the strongest. “Planck 2015” denotes the combination of the low-$\ell$ temperature-polarization likelihood and the high-$\ell$ temperature likelihood.

We use the “Joint Light-curve Analysis” (JLA) sample of type Ia supernovae, which is constructed from the SNLS and SDSS supernova data, together with several samples of low-redshift supernovae Betoule:2012an . Baryon acoustic oscillation (BAO) measurements are another important astrophysical data set, which are powerful to break parameter degeneracies from CMB measurements. We use BAO measurements of $D_{V}/r_{\rm drag}$ from the 6dFGS at $z_{\rm eff}=0.106$, SDSS Main Galaxy Sample at $z_{\rm eff}=0.15$, BOSS LOWZ at $z_{\rm eff}=0.32$, and BOSS CMASS at $z_{\rm eff}=0.57$ Ade:2015xua . Here $D_{V}$ is the effective distance measure for angular diameter distance, $r_{\rm drag}$ is the comoving sound horizon at the end of the baryon drag epoch and $z_{\rm eff}$ is the effective redshift. In our analysis we also use the helium-4 mass fraction measurement of $Y_{\rm p}=0.2449\pm 0.0040$, derived from helium and hydrogen emission lines from metal-poor extragalactic HII regions and from a regression to zero metallicity Aver:2015iza .

In Table 1 we list our results. With the combined data of Planck 2015+JLA+BAO, the deformation parameter is estimated to be $\xi=0.042\pm 0.049$ if the helium-4 abundance is set from the BBN prediction. For comparison, if the helium-4 abundance is not varied independently of other parameters, with the same data set we find $\xi=0.101\pm 0.091$, which indicates that the constraint of $\xi$ is not improved compared to the result derived in Guo:2012mv from the 7-year WMAP Komatsu:2010fb in combination with lower-redshift measurements of the expansion rate. The reason is that the estimated values of cosmological parameters are biased by fixing $Y_{\rm p}=0.24$. The gray and red regions in Fig. 3 show 68% and 95% contours in the $\xi–\Omega_{c}h^{2}$ plane setting $Y_{\rm p}=0.24$, derived from Planck 2015+JLA+BAO and from WMAP7+$H_{0}$+BAO Guo:2012mv , respectively. We see that the Planck 2015 data marginally favor a positive $\xi$ compared to the 7-year WMAP data. The addition of the latest measurement of the helium-4 abundance leads to a strong constraint on the deformation parameter $\xi=-0.002\pm 0.014$. As shown in Figure 3, adding the $Y_{\rm p}$ data breaks the degeneracy between the deformation parameter and the physical dark matter density. This result shows that no signal of Lorentz invariant violation is detected by the joint analysis of CMB and BBN. Figure 4 shows the constraints in the $\xi-\Omega_{b}h^{2}$ plane, $\xi-H_{0}$ plane, $\xi-n_{s}$ plane and $\xi-A_{s}$ plane from Planck 2015+JLA+BAO.

VI Conclusions and discussions

We have studied the cosmological constraints on Lorentz invariance violation in the neutrino sector by a joint analysis of CMB and BBN. Instead of fixing the value of the helium-4 abundance, we have applied the BBN prediction of the helium-4 abundance determined by $\xi$ and $\Omega_{b}h^{2}$ to calculate the CMB power spectra. Using the Planck 2015 data in combination with the JLA sample of type Ia supernovae and the BAO feature, we have put a constraint on the deformation parameter. Adding the measurement of the helium-4 abundance breaks the degeneracy between the deformation parameter and the physical dark matter density and so improves the constraint by a factor 3.

As found in Guo:2012mv , the deformation parameter is nearly uncorrelated with the total mass of neutrinos $\Sigma m_{\nu}$ when WMAP data are used to constrain $\xi$. Planck 2015 data is much more sensitive to the total mass of neutrinos than WMAP data. Therefore, we consider the effects of the neutrino masses on our results. Figure 5 shows 68 and 95% contours in the $\xi-\Sigma m_{\nu}$ plane for Planck 2015 data combined with lower-redshift measurements. We see that $\xi$ is nearly uncorrelated with $\Sigma m_{\nu}$ from Fig. 5.

Compared to cosmological bounds on the Lorentz-violating coefficient, observations of high-energy astrophysical neutrinos give stronger constraints. As listed in Table XIII of Ref. Kostelecky:2011gq , the coefficient is constrained down to ${\cal O}(10^{-9})$ from the time-of-flight measurements, under the assumption that neutrino oscillations are negligible. Cohen and Glashow have argued that the observation of neutrinos with energies in excess of 100 TeV and a baseline of at least 500 km allows us to deduce that the Lorentz-violating parameter is less than ${\cal O}(10^{-11})$ Cohen:2011hx . The present work provides a new way to probe the signal of Lorentz invariance violation in the early Universe, which can in principle be used to constrain $\xi$ in the sterile neutrino sector Hamann:2010bk ; Hamann:2011ge . Although present cosmological data are too weak to yield competitive constraints, future measurements of CMB and BBN offer prospects for placing stringent constraints.

Acknowledgements.

Our numerical analysis was performed on the “Era” of Supercomputing Center, Computer Network Information Center of Chinese Academy of Sciences. This work is supported in part by the National Natural Science Foundation of China Grants No.11575272, No.11690021, No.11690022, No.11335012, No.11375247, No.11435006 and No.91436107, in part by the Strategic Priority Research Program of the Chinese Academy of Sciences, Grant No.XDB23030100 and No.XDB21010100, and by Key Research Program of Frontier Sciences, CAS.