Cosmological constraints on mirror matter parameters

Up-to-date estimates of the cosmological parameters are presented as a result of numerical simulations of cosmic microwave background and large scale structure, considering a flat Universe in which the dark matter is made entirely or partly of mirror matter, and the primordial perturbations are scalar adiabatic and in linear regime. A statistical analysis using the Markov Chain Monte Carlo method allows to obtain constraints of the cosmological parameters. As a result, we show that a Universe with pure mirror dark matter is statistically equivalent to the case of an admixture with cold dark matter. The upper limits for the ratio of the temperatures of ordinary and mirror sectors are around 0.3 for both the cosmological models, that show the presence of a dominant fraction of mirror matter, $0.06 \lesssim \Omega_{\rm mirror} h^2 \lesssim 0.12$.


I. INTRODUCTION
Since a missing mass in the Universe was pointed out in the 1930s, astrophysical evidences for dark matter have been accumulating, increasingly confirming its presence at all cosmological scales. But even if the processes of structure formation can draw a picture of its main features, the elementary nature of dark matter remains unknown.
Modern cosmology provides powerful tools for testing dark matter: Big Bang Nucleosynthesis (BBN), the Cosmic Microwave Background (CMB) and the Large Scale Structure (LSS) power spectra can be reproduced in numerical simulations to better discriminate between the different classes of candidates, one of them giving good agreements with all these observations being known as Cold Dark Matter (CDM). Similarly, mirror dark matter can also account for the cosmological observations, and it is of primary interest to determine whether or not it gives a better agreement than CDM, or if it is equivalent.
Mirror matter was originally proposed by Lee and Yang [1], and further considered by several authors [2][3][4][5], in order to restore the parity symmetry of the Lagrangian of the Standard Model. The most natural way to do so is to add to the existing Lagrangian its parity-symmetric counterpart, so that the whole Lagrangian is invariant under the parity transformation, each part transforming into the other. This corresponds to reintroduce all the known fields with the same coupling constants, but with opposite parities. We therefore end up with a new sector of particles, called mirror sector, which is an exact duplicate of the ordinary sector, but where ordinary par- * paolo.ciarcelluti@gmail.com † quentin.wallemacq@ulg.ac.be ticles have left-handed interactions, mirror particles have right-handed interactions [5]. As a consequence, the three gauge interactions act separately in each sector, the only link between them being gravity. Because mirror baryons, just like their ordinary counterparts, are stable and can be felt only through their gravitational effects, the mirror matter scenario provides an ideal interpretation of dark matter. Its particularity is that it is a self-interacting candidate 1 , but without any new parameter at the level of particle physics.
Since it was introduced, mirror matter has been widely studied and its compatibility with experimental and observational constraints has been verified [6][7][8][9][10][11]. Some potential indications of its existence came from the observation of neutron stars [12,13], but also from the directsearch experiments, for which mirror matter gives one of the few possible explanations [14,15]. Previous analytical and numerical studies on CMB and LSS power spectra [8,[16][17][18][19][20] have given encouraging results, but have only limited the parameter space of mirror dark matter. Here, we explore it in detail using the fast numerical code CAMB, in order to quantify the compatibility of mirror matter with cosmological observations and to obtain constraints on its parameters.
In Section II, we recall the important cosmological quantities and epochs in presence of mirror matter. We describe briefly how we modified the numerical codes CAMB and cosmoMC in order to incorporate the mirror components into the evolution of the Universe in Section III. We then present, in Section IV, the 1-σ constraints on the cosmological parameters coming from the CMB and LSS data for different dark matter compositions, i.e. pure CDM, mixed CDM and mirror matter, and pure mirror matter, and we compare the corresponding bestfit models. We show also the temporal evolution of the density perturbations for the different components of the Universe.

II. BASIC COSMOLOGY WITH MIRROR MATTER
Even if the ordinary and mirror sectors are characterized by identical Lagrangians, and obey therefore the same physical laws, their macroscopic realizations are not necessarily the same. The differences between the evolutions of the two sectors are parametrized by two cosmological free parameters: the ratio x of the temperatures of the ordinary and mirror cosmic background radiations; the relative amount β of mirror baryons compared to the ordinary ones.
where T (T ′ ), Ω b (Ω ′ b ), and S (S ′ ) are respectively the ordinary (mirror) photon temperature, cosmological baryon density (normalized to the critical density of the Universe), and entropy per comoving volume [6].
There are several components that contribute to the total present energy density Ω tot : the energy density of relativistic species (radiation) Ω r , the energy density of non-relativistic species (matter) Ω m , and the energy density of the vacuum (cosmological constant or dark energy) Ω Λ . According to the observations of the CMB anisotropies, Ω tot = Ω m + Ω r + Ω Λ ≈ 1, meaning that the Universe today is almost flat. In presence of mirror matter, mirror components are present in both radiation and matter energy densities, and the latter is expressed by where the term Ω DM includes the contributions of any other possible dark matter particles but mirror baryons. It can be shown that during BBN, the mirror species γ ′ , e ± ′ and ν ′ e,µ,τ , respectively for mirror photon, electrons, positrons and neutrinos, bring a contribution to the relativistic degrees of freedom equivalent to an effective number of extra neutrino families ∆N ν ≃ 6.14 x 4 . In view of the current bounds on ∆N ν [21], this corresponds to an upper limit x 0.7, which means that the temperature of the mirror sector was smaller than that of the ordinary one at the epoch of nucleosynthesis, T ′ < T .
Due to the separate conservation of entropies in the two sectors, this initial temperature difference holds throughout the expansion of the Universe, so that the cosmological key epochs take place at different redshifts in the two sectors, happening earlier in the mirror sector than in the ordinary one [6,17]. The relevant epochs for structure formation are the matter-radiation equality (MRE), the matter-radiation decouplings (MRD) in the ordinary and mirror sectors, and the photon-baryon equipartitions in each sector, occurring respectively at the redshifts z eq , z dec , z ′ dec , z bγ and z ′ bγ . The MRE is common to both sectors and, in presence of mirror matter, happens at the redshift while the MRDs and photon-baryon equipartitions in each sector are respectively related by and since T ′ dec ≃ T dec up to small corrections to Eq. (4). Because x > 0, the value of z eq obtained in the mirror scenario is always smaller than in the standard case, while the upper bound x 0.7 ensures that z ′ dec > z dec and z ′ bγ > z bγ , showing that the MRD and the photon-baryon equipartition occur earlier in the mirror sector.
By identifying z eq and z ′ dec from Eqs. (3) and (4), one obtains a reference value x eq [6,17] under which the mirror photons decouple during the radiation-dominated era, with the consequence that the primordial perturbations evolve, in the linear regime, in a way that is very similar to the standard CDM case. Also, analytical and numerical studies on CMB and LSS power spectra [8,[16][17][18][19][20]22] have already shown, by comparing qualitatively the results with the observations, that for relatively cold mirror sectors (x 0.3), the dark matter of the Universe can be fully realized by mirror baryons, while for higher x (x 0.3) mirror baryons and CDM would form an admixture.

III. THE MODIFIED NUMERICAL TOOLS CAMB AND COSMOMC
We modified the publicly available code CAMB [23] for the simulation of the anisotropies of the microwave background and the large scale structure of the Universe, together with its Markov Chain Monte Carlo sampling tool cosmoMC [24], in order to include the mirror components. CAMB is used in its most standard mode, i.e. with adiabatic initial conditions for the perturbations, in linear regime, in a flat Universe (Ω tot = 1), with an equation of state of vacuum p Λ = −ρ Λ (w Λ = −1), without any massive neutrinos and with the standard number of neutrino families N eff = 3.046.
We defined the necessary mirror variables and included them into the calculations of all the relevant quantities related to the evolution of the background. In particular, we considered both ordinary and mirror matter when describing the gravitational interactions. As the mirror particles obey the same physical laws as ordinary ones, we doubled the equations related to the evolution of the perturbations of the background, adding the corresponding variables for the perturbations of mirror energy densities and peculiar velocities.
The recombinations are calculated separately in each sector, using two times per model the same code REC-FAST [25] present in CAMB. This considerably increases the computational time of a model, especially for small values of x. Indeed, the recombination of mirror hydrogen (as well as of mirror helium) scales as x −1 , as stated by Eq. (4), so that the integration of the equations giving the ionized fractions of mirror hydrogen and helium has to be performed over a wider range of redshifts. However, the program is still fast enough to calculate the huge number of models required for a Monte Carlo analysis of the parameter space.
Compared with previous numerical studies [8,[16][17][18][19], we have used here an updated estimate of the primordial composition of mirror particles from Refs. [6,[26][27][28][29], checking that the models obtained using CAMB with this more accurate treatment of mirror BBN are consistent with the previous ones.
We performed Markov Chain Monte Carlo (MCMC) analyses of the parameter space constituted by the standard cosmological parameters plus the two mirror ones using cosmoMC. To this aim, we added the mirror parameters to the parameter list of cosmoMC and linked them to the equivalent ones in CAMB. We end up with an eight-dimensional set of cosmological parameters for which we adopt flat priors and broad distributions, as summarized in Table I. Ω b h 2 and Ω cdm h 2 are respectively the baryon and cold dark matter densities, x and β are the mirror photon temperature and mirror baryon density relative to the corresponding ordinary quantities, θ s is the ratio of the sound horizon to the angular diameter distance at decoupling, τ is the reionization optical depth, n s is the scalar spectral index and A s is the scalar fluctuation amplitude. The upper limit on x is set by the BBN limit mentioned in Section II. CAMB also calculates derived parameters such as the matter and dark energy densities Ω m and Ω Λ , the redshift of reionization z re , the Hubble parameter h, the age of the Universe in Gyr and the density fluctuation

IV. RESULTS
This modified version of CAMB can be used in three configurations of dark matter: standard CDM, pure mirror matter, and mixed mirror-CDM matter. We considered all these three cases and performed for each of them two analyses, one using the CMB data only and the other using the CMB data combined with the LSS ones. The standard CDM case, with the same assumptions and priors, serves as a reference model.
The CMB datasets come from the WMAP7 team [30], together with the APT [31] and SPT [32] observations. The former provides the acoustic oscillations of the Cosmic Microwave Background on degree scales with limited cosmic-variance precision while the other two give accurate power spectra at higher l's. For the LSS, we included the power spectrum from the SDSS-DR7 luminous red galaxy sample [33] limited to k 0.2h Mpc −1 , corresponding to sufficiently large length scales, where non-linear clustering and scale-dependent galaxy biasing don't have to be taken into account.
In Table II, we show the 1-σ constraints on the parameters obtained from the different dark matter compositions and cosmological tests. The density probability of the parameter x was found to be almost flat in the low-x region and sharply decreasing at higher x. For that reason, we chose to give the upper limits on that parameter at the 95% confidence level. The best-fit models obtained by using the CMB data only and both CMB and LSS data From Table II, we directly note that the primary cosmological parameters, except the CDM energy density, do not vary significantly from one kind of model to the other and for both analyses (CMB and CMB+LSS). On the other hand, some derived parameters are more perturbed, as the total matter density, that increases at the expenses of the dark energy density in models using CMB data only. This is coupled to the decrease of the Hubble parameter, which still falls within the current constraints. The non-baryonic matter density is in all cases 5 or 6 times larger than the baryonic density, as usually derived from standard analyses. Turning to the mirror parameters, the 95% c.l. upper limits on x are found to be x < 0.456 for a pure mirror Universe and x < 0.479 for the mixed mirror-CDM scenario, in case of the CMB only analysis. Adding the LSS constraints in the computations significantly lowers these upper lim-its, since we obtain x < 0.297 and x < 0.315 at 95% c.l. respectively for the pure mirror and mixed mirror-CDM cases. This confirms the higher sensitivity on x of the formation of the large scale structure of the Universe already highlighted in previous studies [6,16]. Note that all these allowed intervals for x contain the values that make possible the interpretation of the direct-darkmatter-search experiments by mirror matter [14,15]. Finally, the parameter β, quantifying the presence of mirror matter in the Universe, has similar values with CMB only and CMB combined with LSS. In the pure mirror case, β lies between 5 and 5.5, indicating that mirror models require consistent amounts of mirror matter in order to reproduce the observables. In the mixed mirror-CDM case, we find mirror matter densities that are between 2 and 4 times larger than those of CDM. This suggests that, in a Universe with a dark matter composed of several components, the mirror one would be the dominant part.
The likelihoods of the best-fit models for the three compositions of dark matter have, for each analysis (CMB and CMB+LSS), very close values. Respectively for pure CDM, pure mirror and mixed mirror-CDM, they are − ln(L) = 3772, − ln(L) = 3771 and − ln(L) = 3771 in case of the CMB only analysis, and 3795, 3794 and 3795 in the other case. In view of the difference of one or two degrees of freedom between the models, these values do not show any statistical preference.
In Fig. 3 we show the temporal evolution of density perturbations for CDM, ordinary and mirror baryons, ordinary and mirror photons, computed at two different scales and for our estimated values of the parameters (for x we considered the upper limits) in the two configurations of pure mirror (upper figures) and mixed mirror-CDM (lower figures). In every plot it is visible that the decoupling of mirror baryons and photons happens before the one of the ordinary species, and this is more evident in the right plots. In particular, comparing the top and bottom right plots, one can see the equivalence between CDM and mirror matter for this scale of perturbations and this region of x. The left panels are for the perturbations at a smaller scale. The mentioned CDM-mirror equivalence is at first sight less evident, but considering the final temporal effect (that gives the distribution of the cosmological matter structures that we see today), again the mirror baryons have a role comparable to the CDM, driving the evolutions of the structures, and providing the gravitational seeds where ordinary baryons can fall and accrete.
The fact that we can access only an upper limit on x, together with the equivalence of the different best fits, could suggest a mirror matter with a CDM-like behavior. Indeed, as announced by Eq. (4) and further confirmed by Refs. [6,16,17], the smaller x, the earlier the decoupling of mirror baryons, and the more they behave like cold dark matter at linear scales, the case x → 0 being equivalent to CDM. If such a trend is verified, CDM could find a possible interpretation through mirror mat- ter. Future data on LSS, especially in non-linear regimes where mirror matter shows more marked differences with CDM at non-zero x, should help to discriminate between mirror and CDM models, or confirm their equivalence.

V. CONCLUSION
In summary, we presented up-to-date estimates of the cosmological mirror parameters x and β coming from the observations of the anisotropies of the Cosmic Microwave Background and of the Large Scale Structure of the Universe. The most stringent constraints were obtained by using both the CMB and LSS data, for which we determined that x < 0.297 (95% c.l.) and β = 5.21 ± 0.21 (1σ) in the pure mirror case, and x < 0.315 (95% c.l.) and β = 3.64 ± 0.46 (1σ) in the mixed mirror-CDM case. These parameter ranges contain the values favored by the direct searches for dark matter. On the other hand, we have seen that cosmological models with dark sectors constituted by pure CDM, pure mirror matter and both CDM and mirror matter are equivalent concerning the CMB and LSS power spectra. The upper limits on x together with the equivalence of the different compositions of dark matter may indicate that, if present, mirror matter could behave like CDM. Future data on LSS at nonlinear length scales should help to discriminate between CDM and mirror matter or confirm their equivalence, in which case the latter would be an interpretation of the former.