Introducing the DiracMilne universe
^{1} UPMCCNRS, UMR 7095, Institut d’Astrophysique de Paris, 75014 Paris, France
email: benoitl@iap.fr
^{2} CSNSM, Université ParisSud, CNRSIN2P3, 91405 Orsay, France
email: chardin@csnsm.in2p3.fr
Received: 9 November 2010
Accepted: 7 October 2011
The ΛCDM standard model, although an excellent parametrization of the present cosmological data, contains two as yet unobserved components, dark matter and dark energy, that constitute more than 95% of the Universe. Faced with this unsatisfactory situation, we study an unconventional cosmology, the DiracMilne universe, a matterantimatter symmetric cosmology, in which antimatter is supposed to present a negative active gravitational mass. The main feature of this cosmology is the linear evolution of the scale factor with time, which directly solves the age and horizon problems of a matterdominated universe. We study the concordance of this model to the cosmological test of type Ia supernovae distance measurements and calculate the theoretical primordial abundances of light elements for this cosmology. We also show that the acoustic scale of the cosmic microwave background naturally emerges at the degree scale despite an open geometry.
Key words: cosmology: theory / dark energy / primordial nucleosynthesis / dark matter
© ESO, 2012
1. Introduction
The inflationbased ΛCDM cosmological model represents a great achievement in modern cosmology. Its predictions are in good agreement with a wide varieties of observational cosmological tests, ranging from bigbang nucleosynthesis and cosmic microwave background anisotropies to distances measurements using type Ia supernovae. However, this success comes at an expensive cost as two poorlyunderstood components are required in the theory to provide the concordance of the cosmological model to the observations. These components are dark matter, whose density is estimated to be around ~20% of the critical density, and dark energy, which accounts for as much as ~75% of the energy content of the Universe.
The physical nature of dark energy remains to be elucidated but all the cosmological observations tend to indicate that it behaves in a similar way as a vacuum energy, characterized by an equation of state p = −ρ, where pressure and density are equal in absolute value but opposite. Dark energy therefore appears to be responsible for the acceleration of the expansion of the Universe, which tends to be confirmed notably by type Ia supernovae distance measurements.
Faced with this uncomfortable situation where we understand less than 5% of our Universe, we study the unconventional cosmology of a symmetric universe i.e. containing equal quantities of matter and antimatter, where antimatter is supposed to present a negative active gravitational mass. The main consequence of this provoking hypothesis is that on large scales, above the characteristic length of the matterantimatter distribution, the expansion factor evolves linearly with time, which is reminiscent of the Milne cosmology (Milne 1933).
The main motivation for considering the possibility that antimatter is endowed with a negative active gravitational mass is provided by the study of Carter (1968) of the KerrNewman geometry representing charged spinning black holes. As noted initially by Carter, the KerrNewman geometry with the mass, charge, and spin of an elementary particle such as an electron bears several of the features expected of the real corresponding particle. In particular, the KerrNewman “electron” has no horizon, has automatically the g = 2 gyromagnetic ratio, and has a ring structure of radius equal to half the Compton radius of the electron. In addition, this geometry has charge and mass/energy reversal symmetries that strongly evoke the CP and T matterantimatter symmetry (Chardin 1997; Chardin & Rax 1992): when the nonsingular interior of the ring is crossed, a second R^{4} space is found where charge and mass change sign (Carter 1968). Therefore, starting from an electron of negative charge and positive mass as measured in the first R^{4} space, we find in the second space a “positron” of positive charge and negative mass. The relation of the KerrNewman geometry to Dirac’s equation and therefore to antimatter has been noted by some authors (Arcos & Pereira 2004; Burinskii 2008).
The introduction of negative masses is not a popular idea amongst physicists and here we do not attempt to demonstrate that antimatter has a negative active gravitational mass. We consider this question instead from a cosmological point of view, by studying the properties of this alternative cosmology and deriving necessary conditions that this model must fulfill to comply with observational tests. Nevertheless, it is important to recall that an analogous system exists in condensed matter, where electrons and holes in a semiconductor appear as quasiparticles with positive inertial mass, but where holes have an energy density lower than the “vacuum” constituted by the semiconductor in its ground state and antigravitate in a gravitational field (Tsidil’kovskiĭ 1975).
The gravitational behavior of antimatter has yet to be investigated in laboratory experiments. Gravity experiments on antihydrogen are under way (Drobychev et al. 2007; Pérez et al. 2008) and are expected to provide important insight into this behavior, hopefully within the next decade.
The paper is organized as follows. In Sects. 2 and 3, we present some general properties of the DiracMilne universe. Section 4 is devoted to the study of primordial nucleosynthesis within the framework of the DiracMilne cosmology. In Sect. 5, we study the compliance of the model to type Ia supernovae distance measurements. In Sect. 6, we investigate the geometric implications on CMB anisotropies and baryonic acoustic oscillations.
2. The DiracMilne universe
We consider a model in which spacetime is endowed with a FriedmannRobertsonWalker (FRW) metric, given by the formula (1)where sink = sinχ,χ,sinhχ according to the value k = 1,0, − 1 of the parameter k describing the spatial curvature, and a(t) is the scale factor. For such a FRW metric, the timetime and spacespace components of the Einstein tensor G_{μν} take the forms (2)On scales larger than the characteristic size of the emulsion, the presence of equal quantities of matter with positive mass and antimatter with negative mass nullifies the stressenergy tensor T_{μν}.
Using Einstein equation G_{μν} = 8πGT_{μν}, we therefore have the equivalence (3)The metric of the DiracMilne universe then reads (4)As in the standard ΛCDM cosmology, the DiracMilne universe has a distinctive geometry: while the standard ΛCDM model has a curved spacetime and flat spatial sections, the DiracMilne universe has a flat spacetime and negatively curved spatial sections.
The Friedmann equation for the DiracMilne universe reads (5)where a_{0} is the present value of the scale factor. Integrating this equation enables one to compute the present age of the Universe t_{U}(6)where H_{0} is the Hubble constant. It should be emphasized that this is a strict equality, whereas in the standard ΛCDM model, the age of the Universe is only approximatively equal to 1/H_{0}. It should also be noted that the linear evolution of the scale factor solves the age problem of the Universe (Chaboyer et al. 1998), which was a prime concern before the introduction of dark energy, and does not affect the DiracMilne cosmology.
Using the metric in Eq. (1), the particle horizon, i.e., the distance a photon can travel since the origin, is given by the limit (7)which diverges logarithmically in the case of the DiracMilne universe. This simple relation has profound implications as it signifies that any two given places in space were causally connected in the past. The DiracMilne universe therefore does not have an horizon problem, which removes the main motivation for the introduction of inflation theories.
Therefore, the linear evolution of the scale factor and the property that the spacetime is flat and spatial sections are open naturally solves two major problems in standard cosmology without requiring additional ingredients such as dark energy or inflation, and is an important motivation for studying in detail such a cosmology.
Before presenting the properties of the DiracMilne universe and for the sake of clarity, we list the underlying hypothesis and necessary constraints that have been assumed so far:

1.
The existence of an efficient mechanism for matter andantimatter separation. We also assume that the DiracMilneuniverse is composed of separated domains of matter andantimatter. These universes have been abundantly studied in thecontext of inhomogeneous nucleosynthesis. Since this universeis matterantimatter symmetric, this system is a percolationemulsion, meaning that the probability that a given domainextends to infinity is close to unity. It remains unclear, however,whether there is an efficient mechanism to produce such anemulsion. A fundamental parameter is the typical size of thisemulsion. We define it as the ratio of the volume V to the surface S separating the two phases: L = V/S. At this stage of work, this characteristic size is a free parameter, which is assumed here to be a constant, that is constrained using primordial nucleosynthesis (cf. Sect. 4).

2.
Antimatter has a negative active gravitational mass and there is a gravitational repulsion between matter and antimatter. This hypothesis is a necessity to avoid contact between matter and antimatter after cosmological recombination. Indeed, it has been demonstrated (Ramani & Puget 1976; Cohen et al. 1998) that in the case of a symmetric universe, where antimatter is assumed to have a positive gravitational mass, the global size of domains of matter and antimatter must be approximatively that of the observable universe. Otherwise, annihilations at the frontiers of matter and antimatter domains would generate a diffuse gamma ray emission that would be in contradiction with observational data.

3.
One of the main interests of the DiracMilne universe consists in providing a solution to the horizon problem without resorting to inflation. However, this in turn relies on the hypothesis that the contribution of radiation to the stressenergy tensor is negligible. If we attribute to radiation a positive contribution to the energy, the evolution of the universe will not differ significantly from the standard evolution for all epochs where the universe is radiation dominated, i.e. for redshifts above a few thousand. However, in the presence of equal quantities of particles with positive and negative masses, this might no longer be the case. A first step toward a full demonstration that the stressenergy tensor of radiation averages to zero can be deduced from the papers of Hoyle & Narlikar (1964), in the case studied by Hawking (1965), where the particle masses are of both signs and in equal numbers. In this case, it can be shown that the stressenergy tensor of radiation is of order 1/N compared to its value when all masses are of the same sign, and where N is the number of mass carriers (particles) in a Hubble horizon.
3. Some aspects of the thermal history of the DiracMilne universe
Some simple properties of the DiracMilne universe are common to those of a purely linear cosmology studied in Lohiya et al. (1998), Kaplinghat et al. (1999, 2000), and Sethi et al. (2005). We briefly recall some of these properties.
3.1. Timetemperature relation
Using Eq. (5), it is straightforward to obtain the relation between the age of the Universe and the redshift, hence the temperature (8)where T_{0} is the present temperature of the Universe, as measured by CMB experiments, T_{0} = 2.725 ± 0.001K (Fixsen & Mather 2002).
Fig. 1 Age of the universe as a function of temperature for the DiracMilne universe (full line, red) and for the fiducial ΛCDM cosmology (dashed line, blue). 

Open with DEXTER 
This relation between time and temperature is valid throughout the whole history of the universe and implies that the thermal history of the DiracMilne universe is drastically modified from the evolution in the standard ΛCDM cosmology. Figure 1 represents the age of the Universe as a function of the temperature for the DiracMilne and the ΛCDM models. It can be seen that, at high temperatures, the DiracMilne universe is much older than the corresponding ΛCDM cosmology. For instance, the traditional 1 MeV ~ 1 s approximation for the standard model becomes 1 MeV ~ 3.3 years in the DiracMilne cosmology. As noted in Lohiya et al. (1998), this difference has profound implications for bigbang nucleosynthesis calculations, and is discussed in Sect. 4.
Another temperature of interest is the temperature of the quarkgluonplasma (QGP) transition. Omnès (1972) proposed that matterantimatter separation occurred around that temperature, owing to a putative repulsive interaction between nucleons and antinucleons. The maximum size of a domain of (anti)matter was controlled by the diffusion of neutrons. Aly (1974) found a maximum size of 7 × 10^{4}cm at a temperature of T ~ 330 MeV, which was at this epoch the estimated temperature of the QGP transition. This size was later found to differ from the minimum size a domain should have in order to ensure a production of primordial helium compatible with observations (Combes et al. 1975; Aly 1978a).
The situation is rather different in the DiracMilne universe as the timescale of the QGP transition is much longer. Since the temperature of the transition is estimated today to be around T ~ 170 MeV (Schwarz 2003), it corresponds to an age of 6 × 10^{5}s in the DiracMilne universe, which is a factor ~10^{10} older than in the standard case. This implies that the maximum size to which a domain could possibly grow (assuming the existence of an efficient separation mechanism) is five orders of magnitude higher. For this reason, the DiracMilne universe is far more weakly constrained by observations than the Omnès cosmology.
3.2. Weak interaction decoupling
A fundamental example of the modifications induced by a linear scale factor can be seen in the epoch of decoupling of the weak interactions. This example was analyzed extensively in Lohiya et al. (1998), so we only provide here the main results. In the standard model, weak decoupling occurs at a temperature T ~ 1MeV. In the DiracMilne universe, this decoupling happens at a lower temperature of around T ~ 80keV, because of the slower variation and the lower value of the expansion rate. This effect is illustrated in Fig. 2. Weak interactions control the n ↔ p equilibrium. At low temperatures, this reaction is limited to the free neutron disintegration (green shortdashed line), but at temperatures higher than 80 keV, the equilibrium between proton and neutrons remains possible. The longdashed line represent the proton conversion rate as a function of the temperature for the DiracMilne universe (red) and the standard cosmology (blue). The analytical expressions for these reaction rates come from Wagoner (1969) and Dicus et al. (1982). Weak interactions decouple when the expansion rate becomes higher than the p ↔ n rate. The small difference in the p ↔ n rate between the two cosmologies is caused by a difference in the neutrino background temperature. As the weak interactions decouple at a temperature of T ~ 80keV, neutrinos indeed also decouple from the photon background also at this temperature, but only after the annihilation of most of the electronpositron pairs. This implies, as noted in Lohiya et al. (1998), that the cosmic neutrino background should have a temperature equal to that of the CMB. This constitutes a distinctive feature of the DiracMilne universe.
4. Primordial nucleosynthesis
Fig. 2 Comparison between expansion rates and weak reaction rates. The shortdashed green line represents the neutrontoproton conversion rate. The two solid lines represent the expansion rate of the DiracMilne universe (red) and the standard cosmology (blue). The long dashed lines represent the protontoneutron conversion rates. Weak interactions decouple when these rates become lower than the expansion rate. 

Open with DEXTER 
The primordial nucleosynthesis is a key success of the standard model of cosmology. Theoretical predictions and observations are in good agreement for ^{4}He, ^{3}He, and D abundances, although some discrepancies exist for ^{7}Li (see Steigman 2007; Cyburt et al. 2008 for reviews of standard bigbang nucleosynthesis (SBBN)). Given the tremendous change in the timescale (of a factor ~10^{8} at T = 1MeV), important modifications to the traditional BBN mechanism arise (Lohiya et al. 1998; Kaplinghat et al. 1999, 2000). Primordial nucleosynthesis in the DiracMilne universe can be described as a twostep process: first the thermal and homogeneous production of ^{4}He and ^{7}Li, and second, the production of D and ^{3}He, the latter being one of the main novelties introduced by the matterantimatter symmetry.
4.1. Thermal and homogeneous BBN
The age of the universe at a temperature of T ~ 80keV, below which deuterium is able to survive its photodisintegration, allowing the formation of ^{4}He nuclei, is about 40 years in the DiracMilne universe. This extremely long timescale, relative to the mean lifetime of the neutron, was initially considered as an impossible obstacle to any kind of primordial helium production in a linear cosmology (Kaplinghat et al. 1999). However, it has been pointed out (Lohiya et al. 1998; Kaplinghat et al. 2000) that ^{4}He production is possible but relies on a somewhat different mechanism than the standard one.
As discussed in the previous sections, weak interactions decouple in the DiracMilne universe at the low temperature of T ~ 80keV. This implies that the neutron and proton populations remain in thermal equilibrium until that temperature. The ratio of their densities is therefore regulated by the Boltzmann factor (9)where Q = 1.29MeV is the difference between the neutron and proton mass.
Fig. 3 ^{4}He production rate for the DiracMilne (dashed red line). The rate drops abruptly when weak interactions decouple around T ~ 80keV. As a comparison, ^{4}He production rate for the standard BBN is also presented (dotted blue line). 

Open with DEXTER 
For temperatures of the order of T ~ 80keV, this number is a very small quantity (~10^{7}) but a few neutrons nevertheless combine with ambient protons to form deuterium, which is then incorporated into helium nuclei. As the weak interactions are still efficient, some protons inversebetadecay to neutrons, restoring the equilibrium value. Given the timescale of 40 years, this very slow process ends up in an effective production of ^{4}He nuclei. Figure 3 present the production rate of ^{4}He in both the DiracMilne universe (dashed line) and the standard cosmology (dotted line). We can discern the sharp decline in the production rate when the weak interactions decouple, leading to the nearly total disappearance of neutrons below T ~ 80keV. As a simple argument illustrating how this mechanism leads to the right amount of ^{4}He, we note that the product of the ^{4}He production rate with the Hubble time at production is roughly the same in the two models. The production rates of ^{4}He were obtained by calculating the time derivative of the ^{4}He abundances.
To compute the primordial abundances of the different light elements, we solved the nonlinear systems of firstorder differential equations describing the network of the nuclear reactions (10)where Y_{i} is the abundance of nuclide i, N_{isot} the number of nuclides, N_{reac} the numbers of nuclear reactions included in the network, the rate of reaction i + j → k + l, and the rate of reaction k + l → i + j.
To ensure accurate results, it is necessary to include a large number (over a hundred) of reactions. As the time interval during which BBN takes place is very long, some slow reactions that are ineffective during SBBN must be integrated into the nuclear reaction network for the DiracMilne universe (Lohiya et al. 1998). Compared to previous studies of light element production in linear cosmology (Kaplinghat et al. 2000), we used the nuclear reaction rates provided by Angulo et al. (1999), Coc et al. (2000), Tang et al. (2003), Descouvemont et al. (2004), as well as the usual rates of Wagoner (1969), Dicus et al. (1982), Caughlan & Fowler (1988), Fukugita & Kajino (1990), Rauscher et al. (1994), and Chen & Savage (1999).
Fig. 4 ^{4}He (top) and ^{7}Li (bottom) theoretical abundances for the DiracMilne universe as a function of the baryon density. Here η_{10} = 10^{10}η. Shaded areas correspond to the observational range. 

Open with DEXTER 
The precise theoretical predictions for the primordial abundances strongly depend – as in the standard model – on the baryonic density. The final primordial abundances of ^{4}He and ^{7}Li as a function of the baryonic density, characterized by the ratio η of the number of baryons to the number of photons (η = n_{b}/n_{γ}), are presented in Fig. 4. The shaded areas correspond to observational constraints (Olive & Skillman 2004; Ryan et al. 2000) on the corresponding primordial abundances. The values of η such that ensure that the theoretical value of Y_{p} is compatible with the observations lie in the range (11)As can be seen in Fig. 4 there is, as in the standard model, no value of η that permits the compatibility of both ^{4}He and ^{7}Li with the model predictions. The lowest possible value of η yields a lithium abundance of (12)Although this value is somewhat larger than the one inferred by observations, we note that it is smaller than the predicted value in standard cosmology calculations, namely (Cyburt et al. 2008). The DiracMilne universe clearly does not solve the ^{7}Li problem, but nevertheless alleviates it.
The value of the baryonic density inferred by the ^{4}He observational constraint is almost 15 times higher than the value usually admitted in the framework of standard cosmology. This high baryonic density is an important feature of the DiracMilne universe as it suppresses the need for nonbaryonic dark matter. Estimates of the dynamical mass in various structures all indicate a matter density in the Universe that is much higher than the baryonic density deduced by the BBN in the standard cosmology, which constitutes a strong motivation to postulate the existence of nonbaryonic massive particles. In the DiracMilne universe, however, the BBN predicts a baryonic density comparable to the matter density estimated by different techniques. There is therefore no compulsory need for nonbaryonic dark matter in the DiracMilne universe.
We note, however, that this high value of the baryonic density aggravates the socalled missing baryon problem (Fukugita 2004). In the standard cosmology, this problem is that only approximately half the baryons predicted by SBBN are observationally detected. The other half is currently believed to reside in the warmhot intergalactic medium (WHIM) (Bregman 2007). In the context of the DiracMilne universe, this WHIM would have to be the major source of baryons, although this possibility remains to be investigated. We note that attempts to explain the dynamical behavior of galaxies consisting primarily of cold molecular (baryonic) gas instead of nonbaryonic dark matter have been performed in the past (Pfenniger et al. 1994) and these scenarios would need to be revisited in the context of the DiracMilne model.
Figure 5 presents the evolution of the primordial abundances for the light elements as a function of temperature for a baryonic density of η = 8.8 × 10^{9}. It can be seen that D and ^{3}He are almost totally destroyed during the stage of thermal production of ^{4}He. This result was found to have severe consequences for linear cosmologies (Kaplinghat et al. 2000). Although this is indeed the case for a regular linear cosmology without antimatter domains, the presence of distinct domains of matter and antimatter in the DiracMilne cosmology provides a natural scenario for the production of D and ^{3}He during a second stage of nucleosynthesis.
Fig. 5 Abundances of light elements obtained in the Milne universe with a baryontophoton ratio η = 8.8 × 10^{9}. Both ^{4}He and ^{7}Li are produced at observationally compatible levels (see Fig. 4), but D and ^{3}He are almost totally destroyed by the very slow thermal nucleosynthesis. 

Open with DEXTER 
4.2. Secondary production of D and ^{3}He
Studies of inhomogeneous bigbang nucleosynthesis as well as primordial nucleosynthesis in the presence of both matter and antimatter domains have been conducted since the late 70’s (Combes et al. 1975; Aly 1978a; Witten 1984; Alcock et al. 1987; Applegate et al. 1987; KurkiSuonio & Sihvola 2000; Rehm & Jedamzik 2001). It has been shown that matterantimatter annihilations at the frontiers of the domains can lead to the production of D, ^{3}He, and T (later decaying as ^{3}He) mainly through two channels of production: nucleodisruption ( and reactions) and photodisintegration of ^{4}He nuclei. These studies provide us with the necessary material to compute the amount of deuterium produced by these various mechanisms.
Our purpose in what follows is to demonstrate the possibility of deuterium secondary production. We therefore consider that the emulsion has a static behavior, in the sense that its comoving size is assumed to remain constant. More precise studies investigating the dynamical behavior of the emulsion are beyond the scope of this first paper and will be treated in upcoming studies.
Annihilations at the frontiers of a domain are driven by the diffusion of nuclei towards the frontiers. The photodisintegration of ^{4}He nuclei by energetic photons resulting from electromagnetic cascades induced by the annihilation photons and nucleodisruption are two possible processes that could produce deuterium.
The main quantity to consider is the diffusion length. It represents the average distance over which a (anti)nucleus can diffuse toward the frontier of the domain on a Hubble time. This length gives an absolute lower bound to the size of the domains, as any concentration of (anti)matter smaller than this diffusion length would be annihilated during a Hubble time. The diffusion length is given by Applegate et al. (1987)(13)where D is the diffusion coefficient and t_{H}(T) is the Hubble time at temperature T. Using the diffusion coefficients given in Jedamzik & Rehm (2001) and Sihvola (2001), we computed the comoving diffusion length represented in Fig. 6.
Fig. 6 Comoving diffusion length (solid red line). The blue dashed line is the comoving thermalizing length of deuterium. Comoving Hubble distance is also shown in green. 

Open with DEXTER 
We can distinguish three regimes of diffusion. The first regime for T ≥ 1MeV is regulated by neutron diffusion. Since neutrons are neutral particles, their electromagnetic interactions with other charged particles are very weak. As the temperature decreases, neutrons disintegrate and diffusion is then maintained by protons, which being charged have a much lower diffusion coefficient, which causes a dip in the diffusion length around 100 keV. As the density decreases because of the expansion, the diffusion length gradually increases until the temperature reaches T ~ 50keV. At this temperature, the density is low enough for the mean distance between protons to be larger than the Debye length of electrons. Protons therefore do not behave as free particles but drag electrons along with them, ensuring charge neutrality. These electrons are themselves subject to Thomson drag (Peebles 1993) and thus limit the diffusion of protons (Jedamzik & Rehm 2001). The effective diffusion coefficient of protons is then (14)where σ_{T} is the Thomson crosssection and ρ_{γ}, the photon energy density.
4.2.1. Annihilation rate
In the most general case, the computation of the annihilation rate is a difficult task (Aly 1978b). However, with our simplified approach based on diffusion, the estimation of the annihilation rate is rather straightforward. The annihilation rate is the number of annihilations per unit of time and surface. The quantity of matter (and antimatter) annihilated over a Hubble time per unit of surface is n_{b}L_{diff}, so that the annihilation rate is simply this quantity divided by the Hubble time (15)This simple expression is found to be in good agreement with the one derived in Cohen et al. (1998). For simplicity, and following previous studies of antimatter BBN (KurkiSuonio & Sihvola 2000; Rehm & Jedamzik 2001), we made the hypothesis that hydrodynamic turbulence, which could be produced by energy release near the domain boundary, can be neglected.
4.2.2. Production by photodisintegration of ^{4}He nuclei
The secondary production of light elements by ^{4}He photodisintegration can occur in many scenarios and has been extensively discussed in the literature. In the framework of standard cosmology, this mechanism is known to produce D and ^{3}He nuclei (Ellis et al. 1992; Protheroe et al. 1995).
Protonantiproton annihilations result in the production of neutral pions, which themselves decay to highenergy photons. Depending on the temperature of the background and their energy, these photons can create e^{ + }e^{ − } pairs on CMBphotons. These newly created pairs can also interact with CMB photons and therefore lead to the creation of electromagnetic cascades (Ellis et al. 1992; Protheroe et al. 1995). These cascades stop when the energy of a photon becomes lower than the pair creation threshold , where E_{γ} is the energy of a thermal photon. Owing to the high number of thermal photons, the threshold energy for pair creation is (Ellis et al. 1992). Photons with energy lower than E_{max} but higher than also undergo elastic scattering on background photons (Svensson & Zdziarski 1990). The resulting spectrum of cascaded photons can be parametrized by (Ellis et al. 1992; KurkiSuonio & Sihvola 2000) (16)where is a normalization constant, and E_{0} is the total energy injected.
Photodisintegration reactions (^{4}He(γ,p)^{3}H, ^{4}He(γ,n)^{3}He, and ^{4}He(γ,np)D) require photons with energies higher than the respective threshold energies Q_{4He(γ,np)D} = 26.07MeV, Q_{4He(γ,p)3H} = 19.81MeV, and Q_{4He(γ,n)3He} = 20.58MeV (Cyburt et al. 2003). The existence of these threshold energies implies that photodisintegration is a late process, as it becomes efficient only when E_{max} becomes higher than one of the previous threshold values. This happens at a temperature T_{ph} ~ 0.5keV.
Fig. 7 Quantity of D (solid green) and ^{3}He (dashed red) nuclei produced by ^{4}He nuclei photodisintegration for one annihilation. 

Open with DEXTER 
Highenergy photons responsible for ^{4}He photodisintegration mainly interact with ambient electrons by means of Compton scattering or with ambient nuclei by means of the BetheHeitler process (Jedamzik 2006). Taking this into account, the number of D nuclei produced per annihilation is given by (17)where f(E_{γ}) = n_{p}σ_{BH}(E_{γ},1) + n_{α}σ_{BH}(E_{γ},2) + k(E_{γ})n_{e}σ_{KN}(E_{γ}), σ_{BH} is the BetheHeitler process crosssection, σ_{KN} is the Compton scattering crosssection in the KleinNishina regime (Rybicki & Lightman 1979), and k(E_{γ}) ≈ 1 − 4/3 [ln(2E_{γ}/m_{e}) + 1/2 ] ^{1} is the mean fractional energy loss by Compton scattering (Protheroe et al. 1995). A similar formula exists for ^{3}He. Figure 7 presents the results of this calculation.
When photodisintegration is most effective (around T ~ 100eV), the mean free path of high energy photons relative to these processes is longer than the diffusion length, implying that D and ^{3}He nuclei produced by photodisintegration will be able to survive and add to the overall light element production. More precisely, one can estimate the temperature (or equivalently the redshift) at which newly produced nuclei will not be annihilated later on. It is then necessary to assume that there exists a redshift z_{end} where matter and antimatter cease to annihilate. The hypothesis of gravitational repulsion between matter and antimatter leads to this gravitational decoupling but the exact mechanism, which would be analogous to the separation of electrons and holes in a gravitational field, still needs to be established. We emphasize that this decoupling is a necessary condition for the viability of the DiracMilne model. A theoretical determination of z_{end} appears to be possible by analogy with the electronhole system (Tsidil’kovskiĭ 1975), but we assume here that it is a free parameter that will be later constrained by observations. We define the comoving value, calculated at a redshift z, of the diffusion length at redshift z_{end} to be (18)Any nuclei produced at a distance from the domain boundary smaller than this length will be annihilated by z_{end}, whereas nuclei produced farther away will survive once annihilation stops. We can then determine the redshift z_{ ∗ } below which the mean free path of high energy photons is longer than , where z_{ ∗ } is then the redshift below which D and ^{3}He nuclei produced by photodisintegration will add up to the final abundances.
Multiplying Eq. (17) by the annihilation rate and integrating between the times corresponding to z_{ ∗ } and z_{end}, we get , the number of D nuclei produced by photodisintegration of ^{4}He nuclei (19)We can get rid of the surface term S(z) by dividing Eq. (19) by the total number of baryons in the volume V, n_{b}V, which yields the final D abundance (20)This expression depends only on two variables, z_{end}, the redshift of gravitational decoupling and L_{1keV}, the comoving emulsion size at 1 keV. The values of the final D and ^{3}He abundances as a function of these two parameters are presented on Fig. 8.
Fig. 8 D and ^{3}He final abundance as a function of the size of the emulsion and the redshift of gravitational decoupling. Hatched region is forbidden because of CMB distortions constraints. 

Open with DEXTER 
Electromagnetic cascades initiated by annihilation photons cause the injection of nonthermal energy into the background radiation. In the standard cosmology, these energy injections are tightly constrained, most notably by the FIRAS measurements (Fixsen et al. 1996). However, in the DiracMilne universe, owing to the low value of the expansion rate, the radiative processes that can thermalize virtually any energy injection (Danese & de Zotti 1977; Hu & Silk 1993) decouple at a lower redshift than in the standard case. The detailed calculation of the thermalization of CMB distortions will be presented elsewhere, but the exclusion contour in the (z_{end} − L_{1keV}) plane is shown in Fig. 8.
To ensure that the final D abundance agrees with observations (D/H ~ 3 × 10^{5}, Pettini et al. 2008) and imposing the gravitational decoupling as late as possible, we find that the comoving size of the domains at 1 keV has to be around (21)which corresponds to a comoving size of 7 kpc today. This estimate of the size of the emulsion may seem small compared, for instance, to the typical size of a galaxy, implying that there should be numerous matter and antimatter domains inside a single galaxy. However, we recall that this size is only constrained when gravitational decoupling occurs and that domains, which are in a nonlinear mode of evolution immediately after photon matterradiation decoupling, will have a dynamical behavior after recombination that effectively increases this size (Dubinski et al. 1993; Piran 1997).
As suggested by Fig. 7, the secondary production of light elements by photodisintegration of ^{4}He leads to an overproduction of ^{3}He relative to D. It is indeed usually assumed that D can only be destroyed by stellar processes. On the other hand, ^{3}He can be either produced or destroyed, but the ratio D/^{3}He can only decrease (Sigl et al. 1995). Since this ratio is observed in our Galaxy to be D/^{3}He ~ 1 but predicted to be D/^{3}He ~ 0.1 in the DiracMilne universe, the overproduction of ^{3}He is a priori a strong constraint on the DiracMilne model. However, we note (see Jedamzik 2002 for a review) that D enhancement could occur at high redshift diminishing the efficacy of the ^{3}He constraint.
4.2.3. Production by nucleodisruption
Nucleodisruption has been found to be a significant producer of D and ^{3}He nuclei (Sihvola 2001) within the standard evolution of the scale factor. However, with the slow evolution of the expansion rate in the DiracMilne universe and the hypothesis we made about the spatial repartition of the matter and antimatter domains, this situation changes.
Nuclei produced by nucleodisruption possess a kinetic energy ranging from a few MeV for nuclei to a few tens of MeV for nucleons (Balestra et al. 1988). These newly produced nuclei thermalize by Coulomb scattering on ambient protons and electrons. The thermalization length for D nuclei produced with an energy E_{0} = 10MeV is presented in Fig. 6 (blue dotted line). This distance is always much smaller than the diffusion length, implying that any D nucleus produced by nucleodisruption will finally return towards the annihilation zone and be destroyed there. A possible way to produce a higher fraction of deuterium by nucleodisruption would be to consider small domains of (anti)matter within a larger domain of antimatter (matter). This situation occurs continually in an emulsion, which suffers a redistribution of “domains” when bridges in the emulsion disappear by annihilation. If the dimension of the larger domain is larger than the diffusion length, then an important fraction of the D and ^{3}He produced by nucleodisruption could survive. However, precise calculations of this production require the knowledge of the statistical properties of the spatial distribution of domains, which strongly depends of the separation mechanism. This point should be investigated in future studies of nucleosynthesis in the DiracMilne universe.
5. Type Ia supernovae
In 1998 (Riess et al. 1998; Perlmutter et al. 1999), distance measurements for type Ia supernovae (SNe Ia) revealed that these objects are dimmer than expected if our Universe was correctly described by a decelerating Einsteinde Sitter model. The introduction of a cosmological constant Λ in the field equations of general relativity, which is apt to produce an accelerating expansion, provided an impressive fit to the observational data. Today, SNe Ia are one of the most important cosmological tests and are considered as prime evidence of an acceleration of the expansion. We recall however that the strong evidence of a recent transition between a decelerating phase and an accelerating phase of expansion heavily relies on the prior hypothesis of spatial flatness. Without this hypothesis, the evidence is less clear (Seikel & Schwarz 2008).
The DiracMilne universe has neither acceleration nor deceleration and is therefore equivalent to an open empty universe. In terms of the usual cosmological parameters, this universe corresponds to the combination (22)In this context, the luminosity distance in the DiracMilne universe follows the simple expression (23)It is usually claimed that the empty universe, hence the DiracMilne universe, is strongly inconsistent with the SNe Ia observations. We wish here to elaborate on this statement using the data of the first release of the SNLS collaboration (Astier et al. 2006). The SNLS data consist of two distinct datasets. The highredshift sample, from the SNLS, comprises 71 SNe Ia with redshifts between 0.2 ≤ z ≤ 1.01. The second sample is a lowredshift set, consisting of 44 SNe Ia taken from the literature with redshifts z ≤ 0.15. These data come from different experiments and are therefore possibly subject to different sources of systematic errors.
Following the definition given in Astier et al. (2006), the distance modulus is (24)where M is the absolute magnitude of SNe Ia, α and β are global parameters that link the stretch s and the color c to the distance modulus, and is the apparent magnitude of the supernova. It should be emphasized that, in contrast to the ΛCDM cosmology, there is no cosmological parameter dependence in the DiracMilne luminosity distance. The only degrees of freedom are the nuisance parameters, M, α, and β. Following the procedure described in Astier et al. (2006), we minimize the expression (25)Here, σ(μ_{B}) takes into account measurement errors in the apparent magnitude , stretch, and color parameters derived by the analysis of light curves (Guy et al. 2005), σ_{int} is the socalled “intrinsic” dispersion, which is a parameter introduced to account for SNe Ia being astrophysical objects that naturally have some intrinsic dispersion in their absolute magnitude. However, the value of this parameter is unknown, and in the fitting procedure, σ_{int} is adjusted to ensure that the reduced chisquared is unity.
5.1. Analysis with only the highz sample
We first performed our analysis on the highz sample without including any lowz SN Ia. Without this lowz anchoring, the analysis does not permit us to discriminate between the ΛCDM and the DiracMilne universes. In this respect, we note that the threeyear analysis of SNLS using their data alone (Guy et al. 2010) is consistent to a better than 68% CL with the DiracMilne universe, while the Einsteinde Sitter (EdS) model is clearly excluded. The evidence of an expansion acceleration therefore relies on a comparison between lowz and highz SNe Ia. We also present the results for the EdS model.
Values of the different parameters for the fit using only the 71 highz SNe Ia of the SNLS.
Fig. 9 Residuals of the Hubble diagram for the DiracMilne (a) and b)), flat ΛCDM (c) and d)), and Einstein deSitter (e) and f)) models. The left column represents the residuals obtained by the minimization of the χ^{2} defined by Eq. (25), in which the values of the intrinsic dispersion is adjusted so that χ^{2}/d.o.f. = 1. In the right column, this intrinsic dispersion is fixed to 0. 

Open with DEXTER 
Values of the different parameters for the fit using only the 71 highz SNe Ia of the SNLS.
The results of the analysis of the 71 SNe Ia are presented in Tables 1 and 2. In Table 1, the intrinsic dispersion is fixed to a null value. The analysis is therefore performed using only the measurement errors, thereby giving a stronger weight to SNe Ia with redshifts 0.2 ≤ z ≤ 0.4, which have smaller errors. The total and reduced χ^{2} of the EdS model are much larger than those of the ΛCDM and DiracMilne models, which are in turn similar. In Table 2, the intrinsic dispersion is determined by assuming that the reduced chisquared is unity. The value of the intrinsic dispersion required to renormalize the reduced χ^{2} of the EdS model to unity is here again much larger than those required for the other two models. The conclusion of this first analysis is twofold. First, we confirm that the decelerating Einsteinde Sitter model is extremely unlikely, which should come as no surprise. Second, the DiracMilne and the flat ΛCDM models are almost identical, DiracMilne being in even closer agreement with the data than the flat ΛCDM model. Our analysis stresses that all SNe Ia analysis depend strongly on the use of lowz data to anchor the Hubble diagram.
5.2. Analysis on the full data set
Similarly, we proceeded with our analysis of the full data sample used by SNLS in its oneyear analysis, i.e. using an heterogeneous sample of lowz SNe Ia. Our results are given in Tables 3 and 4. In this case, as expected, the flat ΛCDM provides a closer fit to the data than the DiracMilne universe.
Values of the different parameters for the fit using the complete data set of 115 SNe Ia.
Values of the different parameters for the fit using the complete data set of 115 SNe Ia.
In Fig. 9, we present the residuals of the Hubble diagram for the DiracMilne, the flat ΛCDM, and Einsteinde Sitter models. The lefthand column represents the residuals when the value of the intrinsic dispersion was adjusted to normalize the χ^{2} to 1 per degree of freedom. In the righthand column, we present the residuals obtained when the intrinsic dispersion parameter is fixed to the value σ_{int} = 0. Setting σ_{int} to 0 enables one to consider the “real” measurement errors. It appears that for SNe with redshifts in the range 0.2 ≤ z ≤ 0.8 these errors are smaller than the intrinsic dispersion. The use of such an ad hoc parameter in the analysis may degrade the quality of the data in the redshift interval 0.2 − 0.4, where the quality of the observations is highest. We note that such an analysis with statistical errors only was previously performed by the SNLS threeyear analysis (Guy et al. 2010).
Panels e) and f) present the residuals for the Einsteinde Sitter model. Even though the reduced χ^{2} has been constrained to unity in panel e), the characteristic slope in the residuals illustrates the nonconformity of the Einsteinde Sitter model to the data. The difference between the DiracMilne and the ΛCDM models is however almost negligible, and can hardly be seen by simply examining the residuals.
We note that it is only the use of nearby SNe Ia that enables us to distinguish between the DiracMilne and the ΛCDM models, the latter then being the most likely model at a confidence level of more than 3σ as previously announced (Astier et al. 2006; Kowalski et al. 2008).
It appears that sources of previously unaccounted for systematic errors (Kelly et al. 2010) are present in the nearby sample data set. To investigate this, we determined the constant offset to the apparent magnitude of nearby SNe Ia required for the chisquared for the DiracMilne and the flat ΛCDM to become equal. We found that an offset of is sufficient to ensure that the two models are equally probable. This value should be compared to the budget of systematics errors estimated in recent studies: ΔM = 0.04mag (Kowalski et al. 2008). Therefore, a relatively mild systematic error of 1.5σ for nearby SNe Ia would lead us to favor the DiracMilne universe over the conventional ΛCDM cosmology in the SNe Ia analysis.
6. Other tests
A major result of CMB experiments has been the precise measurement of the position of the first acoustic peak on the degree scale, which seems to imply that the spatial curvature is nearly zero (Komatsu et al. 2011). In the open spatial geometry of the DiracMilne universe, this position would naively be expected at a much smaller angle. The ratio of the angular distances in the two models taken at redshift z ~ 1100, which corresponds to the surface of last scattering surface, is (26)The value of this ratio implies that an astrophysical object at redshift z = 1100 is seen under an angle 169 times smaller in the DiracMilne universe than in the ΛCDM cosmology.
The angular position of the first peak is defined by the angle under which the sound horizon is seen at recombination (27)where χ_{s}(z_{ ∗ }) is the sound horizon, d_{A}(z_{ ∗ }) the angular distance, and z_{ ∗ } the redshift of the last scattering surface. It is of interest to consider the equivalent multipole ℓ_{a} ~ π/θ. The sound horizon is defined as the distance that acoustic waves can travel in the primordial plasma. Taking into account the universe expansion, this distance reads (28)where the speed of sound , R being a corrective factor caused by the presence of baryons (Hu & Sugiyama 1995). Its value can be related to the baryon to photon ratio η by R ≈ 1.1 × 10^{12}η/(1 + z).
The definition of the lower bounds of the integral requires some care as we have shown previously that this integral diverges near the initial singularity. However, the mechanisms of sound generation in the DiracMilne and ΛCDM universes differ radically. In contrast to the generation of sound waves in the standard model where inhomogeneities are produced at the epoch of inflation, it seems reasonable to consider that sound waves in the DiracMilne universe are probably produced by annihilation at the matterantimatter frontiers. It is therefore natural to consider for this value the epoch of QGP transition around T ~ 170MeV. In the absence of inflation, the universe is indeed very homogeneous before this event, and QGP transition seems naturally to be the scale of interest in the DiracMilne cosmology. There is at present a rather general consensus that the QCD transition is of neither first nor second order at null chemical potential, but is probably an analytic crossover (see e.g. Aoki et al. 2006; Bazavov et al. 2009; Endrődi et al. 2011, and references therein). Under this hypothesis, it would be extremely difficult, if not impossible, to understand how significant amounts of matter and antimatter could have survived annihilation in the very slow evolution of the DiracMilne universe. However, this present belief is based on a QCD calculation, whereas the situation of the primordial universe is more complicated, including additional light leptons (electrons, positrons, neutrinos, and antineutrinos). In addition, present calculations use light quark masses for the u and d quarks that are significantly higher than the actual masses of these quarks. Therefore, the observation of large matterantimatter domains would be an extremely useful indication that, in contrast to present expectations, there is a sharp transition allowing the survival of significant regions of antimatter at the QCD transition. It should also be noted that some authors are clearly considering the possibility that baryogenesis occurs at the QCD transition or at a O(100 MeV) temperature (see for example Dolgov 1992 for a review).
Acoustic waves then propagate in the plasma as long as matter and antimatter are in contact, i.e. until the gravitational decoupling, estimated in the previous section at z_{end} ≈ 3 × 10^{4}. With these values, the comoving sound horizon is found to be (29)The expression of the angular position of the first acoustic peak then follows (30)It can be shown that the redshift of the last scattering surface in the DiracMilne universe is a few percent lower than in the standard cosmology, here again due to the late decoupling of the radiative processes leading to recombination. We found that z_{∗} ~ 1040.
Calculating the multipole of the acoustic scale using expression (30), we obtain ℓ_{a} ~ 160. The standard value of this quantity is ℓ_{a} ~ 300 (Spergel et al. 2003). Instead of a discrepancy of a factor ≈ 169 , there is an almost exact compensation between the larger geometrical term, induced by the open geometry of the DiracMilne universe, and the larger sound horizon, caused by the slow evolution of the expansion rate before recombination. Taking into account the numerous approximations in the model, this remarkable coincidence is quite unexpected and represents a fascinating motivation to study the DiracMilne universe in more detail.
However, the sound horizon scale is also imprinted in the largescale structure power spectrum under the form of small oscillations (Eisenstein & Hu 1998) called baryonic acoustic oscillations (BAOs). These oscillations are expected, at least in a first approximation, on the same scale as the sound horizon. As discussed above, the sound horizon is much larger in the DiracMilne universe than in the standard cosmology. Admitting that there are BAOs in the DiracMilne universe, they should be expected on a scale much larger than that of the standard cosmology. The claimed detection of these BAOs on the expected scale (within the standard cosmology) of ~150Mpc (Eisenstein et al. 2005), presently detected at the ~3σ level, if confirmed by ongoing experiments, would therefore provide a strong constraint on the DiracMilne universe.
7. Conclusion
Since the standard ΛCDM model is in good agreement with observations but rather poorly theoretically motivated, we have studied here an alternative cosmological model, the DiracMilne universe. Inspired by the work of Dirac, Kerr, and Carter, this model restores the symmetry between matter and antimatter. Relying on the symmetries of the KerrNewman solutions in general relativity, it makes the hypothesis that particles and antiparticles behave similarly to quasiparticles such as electrons and holes in a semiconductor, and that antimatter has a negative active gravitational mass. A fundamental characteristic of this Universe is the linear evolution of its scale factor, which solves in an elegant way both the problems of the horizon and the age of the universe.
For primordial nucleosynthesis, we have found that the DiracMilne universe is able to produce ^{4}He at an adequate level, while producing ^{7}Li nuclei in proportions admittedly a factor three higher than the observed values but with a smaller disagreement between observations and predictions than the standard cosmology. We have also shown that surface annihilations at the frontiers of matter and antimatter naturally lead to the production of D nuclei, in amounts that are proportional to the inverse of the characteristic size of the emulsion. The main focus of this study has been the production of D, although the assumption that the emulsion has a fixed characteristic size is clearly an approximation. Relaxing this assumption will allow the possibility of the total annihilation of small patches of antimatter inside larger regions of matter. This might lead to a net production of D by nucleodisruption, which has a higher D/^{3}He production ratio (Balestra et al. 1988), hence alleviate the constraints made by the overproduction of ^{3}He nuclei.
The DiracMilne universe does not undergo an accelerated expansion and therefore seems to conflict with the usual interpretation of type Ia distance measurements. We have shown that the DiracMilne universe can nevertheless be reconciled with type Ia supernovae observations if we take into account a reasonable size of systematic errors in the lowz supernovae subset. Finally, and perhaps the most surprising result of this unusual cosmology is that the acoustic scale naturally emerges at the degree scale, despite the open geometry.
These first results are encouraging but some issues have not been addressed in this study. Amongst these is how structure formation occurs in the presence of separate domains of positive and negative mass. The use of numerical simulations will most probably be a necessity as the usual linear approximation does not hold. Immediately after matterradiation decoupling, density contrast is indeed on the order of unity for any distribution of matter with positive mass and antimatter with negative mass.
Acknowledgments
It is a pleasure to acknowledge fruitful discussions with J. Andrea, E. Armengaud, B. Carter, N. Fourmanoit, J. Fric, K. Jedamzik, T. Jolicœur, E. Keihanen, R. Pain, J. Rich, and the members of the SNLS collaboration. We express our special thanks to A. Coc for allowing us to adapt his nucleosynthesis code. We also thank an anonymous referee for his/her questions and comments. Needless to say, these people are not responsible for the errors present in this admittedly provocative but hopefully interesting paper.
References
 Alcock, C., Fuller, G. M., & Mathews, G. J. 1987, ApJ, 320, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Aly, J. J. 1974, A&A, 35, 311 [NASA ADS] [Google Scholar]
 Aly, J. J. 1978a, A&A, 64, 273 [NASA ADS] [Google Scholar]
 Aly, J. J. 1978b, A&A, 67, 199 [NASA ADS] [Google Scholar]
 Angulo, C., Arnould, M., Rayet, M., et al. 1999, Nuclear Phys. A, 656, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Aoki, Y., Endrődi, G., Fodor, Z., Katz, S. D., & Szabó, K. K. 2006, Nature, 443, 675 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Applegate, J. H., Hogan, C. J., & Scherrer, R. J. 1987, Phys. Rev. D, 35, 1151 [NASA ADS] [CrossRef] [Google Scholar]
 Arcos, H. I., & Pereira, J. G. 2004, General Relativity and Gravitation, 36, 2441 [NASA ADS] [CrossRef] [Google Scholar]
 Astier, P., Guy, J., Regnault, N., et al. 2006, A&A, 447, 31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Balestra, F., Bossolasco, S., Bussa, M. P., et al. 1988, Nuovo Cimento A Serie, 100, 323 [Google Scholar]
 Bazavov, A., Bhattacharya, T., Cheng, M., et al. 2009, Phys. Rev. D, 80, 014504 [NASA ADS] [CrossRef] [Google Scholar]
 Bregman, J. N. 2007, ARA&A, 45, 221 [NASA ADS] [CrossRef] [Google Scholar]
 Burinskii, A. 2008, Gravit. Cosmol., 14, 109 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Carter, B. 1968, Phys. Rev., 174, 1559 [NASA ADS] [CrossRef] [Google Scholar]
 Caughlan, G. R., & Fowler, W. A. 1988, Atom. Data Nucl. Data Tables, 40, 283 [NASA ADS] [CrossRef] [Google Scholar]
 Chaboyer, B., Demarque, P., Kernan, P. J., & Krauss, L. M. 1998, ApJ, 494, 96 [NASA ADS] [CrossRef] [Google Scholar]
 Chardin, G. 1997, Hyperfine Interactions, 109, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Chardin, G., & Rax, J.M. 1992, Phys. Lett. B, 282, 256 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, J.W., & Savage, M. J. 1999, Phys. Rev. C, 60, 065205 [NASA ADS] [CrossRef] [Google Scholar]
 Coc, A., Hernanz, M., José, J., & Thibaud, J.P. 2000, A&A, 357, 561 [NASA ADS] [Google Scholar]
 Cohen, A. G., de Rujula, A., & Glashow, S. L. 1998, ApJ, 495, 539 [NASA ADS] [CrossRef] [Google Scholar]
 Combes, F., FassiFehri, O., & Leroy, B. 1975, Ap&SS, 37, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Cyburt, R. H., Ellis, J., Fields, B. D., & Olive, K. A. 2003, Phys. Rev. D, 67, 103521 [NASA ADS] [CrossRef] [Google Scholar]
 Cyburt, R. H., Fields, B. D., & Olive, K. A. 2008, J. Cosmol. AstroPart. Phys., 11, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Danese, L., & de Zotti, G. 1977, Nuovo Cimento Rivista Serie, 7, 277 [Google Scholar]
 Descouvemont, P., Adahchour, A., Angulo, C., Coc, A., & VangioniFlam, E. 2004, Atom. Data Nucl. Data Tables, 88, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Dicus, D. A., Kolb, E. W., Gleeson, A. M., et al. 1982, Phys. Rev. D, 26, 2694 [NASA ADS] [CrossRef] [Google Scholar]
 Dolgov, A. D. 1992, Phys. Rep., 222, 309 [NASA ADS] [CrossRef] [Google Scholar]
 Drobychev, G. Y., Nédélec, P., Sillou, D., et al. 2007, Proposal for the AEGIS experiment at the CERN antiproton decelerator (Antimatter Experiment: Gravity, Interferometry, Spectroscopy), Tech. Rep. SPSCP334. CERNSPSC2007017, CERN, Geneva [Google Scholar]
 Dubinski, J., da Costa, L. N., Goldwirth, D. S., Lecar, M., & Piran, T. 1993, ApJ, 410, 458 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, 560 [NASA ADS] [CrossRef] [Google Scholar]
 Ellis, J., Gelmini, G. B., Lopez, J. L., Nanopoulos, D. V., & Sarkar, S. 1992, Nucl. Phys. B, 373, 399 [NASA ADS] [CrossRef] [Google Scholar]
 Endrődi, G., Fodor, Z., Katz, S. D., & Szabó, K. K. 2011, J. High Energy Phys., 4, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Fixsen, D. J., & Mather, J. C. 2002, ApJ, 581, 817 [NASA ADS] [CrossRef] [Google Scholar]
 Fixsen, D. J., Cheng, E. S., Gales, J. M., et al. 1996, ApJ, 473, 576 [NASA ADS] [CrossRef] [Google Scholar]
 Fukugita, M. 2004, in Dark Matter in Galaxies, ed. S. Ryder, D. Pisano, M. Walker, & K. Freeman, IAU Symp., 220, 227 [Google Scholar]
 Fukugita, M., & Kajino, T. 1990, Phys. Rev. D, 42, 4251 [NASA ADS] [CrossRef] [Google Scholar]
 Guy, J., Astier, P., Nobili, S., Regnault, N., & Pain, R. 2005, A&A, 443, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guy, J., Sullivan, M., Conley, A., et al. 2010, A&A, 523, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hawking, S. W. 1965, Roy. Soc. London Proc. Ser. A, 286, 313 [NASA ADS] [CrossRef] [Google Scholar]
 Hoyle, F., & Narlikar, J. V. 1964, Roy. Soc. London Proc. Ser. A, 282, 191 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., & Silk, J. 1993, Phys. Rev. D, 48, 485 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., & Sugiyama, N. 1995, ApJ, 444, 489 [NASA ADS] [CrossRef] [Google Scholar]
 Jedamzik, K. 2002, Planet. Space Sci., 50, 1239 [NASA ADS] [CrossRef] [Google Scholar]
 Jedamzik, K. 2006, Phys. Rev. D, 74, 103509 [NASA ADS] [CrossRef] [Google Scholar]
 Jedamzik, K., & Rehm, J. B. 2001, Phys. Rev. D, 64, 023510 [NASA ADS] [CrossRef] [Google Scholar]
 Kaplinghat, M., Steigman, G., Tkachev, I., & Walker, T. P. 1999, Phys. Rev. D, 59 [Google Scholar]
 Kaplinghat, M., Steigman, G., & Walker, T. P. 2000, Phys. Rev. D, 61, 103507 [NASA ADS] [CrossRef] [Google Scholar]
 Kelly, P. L., Hicken, M., Burke, D. L., Mandel, K. S., & Kirshner, R. P. 2010, ApJ, 715, 743 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kowalski, M., Rubin, D., Aldering, G., et al. 2008, ApJ, 686, 749 [NASA ADS] [CrossRef] [Google Scholar]
 KurkiSuonio, H., & Sihvola, E. 2000, Phys. Rev. D, 62, 103508 [NASA ADS] [CrossRef] [Google Scholar]
 Lohiya, D., Batra, A., Mahajan, S., & Mukherjee, A. 1998 [arXiv:grqc/9808031] [Google Scholar]
 Milne, E. A. 1933, Zeitschrift fur Astrophysik, 6, 1 [NASA ADS] [Google Scholar]
 Olive, K. A., & Skillman, E. D. 2004, ApJ, 617, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Omnès, R. 1972, Phys. Rep., 3, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Peebles, P. J. E. 1993, Principles of physical cosmology (Princeton Series in Physics, Princeton, NJ: Princeton University Press) [Google Scholar]
 Pérez, P., Liszkay, L., Barthe, M. F., et al. 2008, in Amer. Inst. Phys. Conf. Ser., ed. Y. Kanai & Y. Yamazaki, 1037, 35 [Google Scholar]
 Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Pettini, M., Zych, B. J., Murphy, M. T., Lewis, A., & Steidel, C. C. 2008, MNRAS, 391, 1499 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Pfenniger, D., Combes, F., & Martinet, L. 1994, A&A, 285, 79 [NASA ADS] [Google Scholar]
 Piran, T. 1997, General Relativity and Gravitation, 29, 1363 [NASA ADS] [CrossRef] [Google Scholar]
 Protheroe, R. J., Stanev, T., & Berezinsky, V. S. 1995, Phys. Rev. D, 51, 4134 [NASA ADS] [CrossRef] [Google Scholar]
 Ramani, A., & Puget, J. L. 1976, A&A, 51, 411 [NASA ADS] [Google Scholar]
 Rauscher, T., Applegate, J. H., Cowan, J. J., Thielemann, F.K., & Wiescher, M. 1994, ApJ, 429, 499 [NASA ADS] [CrossRef] [Google Scholar]
 Rehm, J. B., & Jedamzik, K. 2001, Phys. Rev. D, 63, 043509 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, ApJ, 116, 1009 [Google Scholar]
 Ryan, S. G., Beers, T. C., Olive, K. A., Fields, B. D., & Norris, J. E. 2000, ApJ, Lett., 530, L57 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Rybicki, G. B., & Lightman, A. P. 1979, Radiative processes in astrophysics (New York, WileyInterscience, 393) [Google Scholar]
 Schwarz, D. J. 2003, Annalen der Physik, 515, 220 [NASA ADS] [CrossRef] [Google Scholar]
 Seikel, M., & Schwarz, D. J. 2008, J. Cosmol. AstroPart. Phys., 2, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Sethi, G., Dev, A., & Jain, D. 2005, Phys. Lett. B, 624, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Sigl, G., Jedamzik, K., Schramm, D. N., & Berezinsky, V. S. 1995, Phys. Rev. D, 52, 6682 [NASA ADS] [CrossRef] [Google Scholar]
 Sihvola, E. 2001, Phys. Rev. D, 63, 103001 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJSS, 148, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Steigman, G. 2007, Ann. Rev. Nucl. Part. Sci., 57, 463 [NASA ADS] [CrossRef] [Google Scholar]
 Svensson, R., & Zdziarski, A. 1990, ApJ, 349, 415 [NASA ADS] [CrossRef] [Google Scholar]
 Tang, X., Azhari, A., Gagliardi, C. A., et al. 2003, Phys. Rev. C, 67, 015804 [NASA ADS] [CrossRef] [Google Scholar]
 Tsidil’kovskiĭ, I. M. 1975, Sov. Phys. Uspekhi, 18, 161 [NASA ADS] [CrossRef] [Google Scholar]
 Wagoner, R. V. 1969, ApJSS, 18, 247 [NASA ADS] [CrossRef] [Google Scholar]
 Witten, E. 1984, Phys. Rev. D, 30, 272 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Values of the different parameters for the fit using only the 71 highz SNe Ia of the SNLS.
Values of the different parameters for the fit using only the 71 highz SNe Ia of the SNLS.
Values of the different parameters for the fit using the complete data set of 115 SNe Ia.
Values of the different parameters for the fit using the complete data set of 115 SNe Ia.
All Figures
Fig. 1 Age of the universe as a function of temperature for the DiracMilne universe (full line, red) and for the fiducial ΛCDM cosmology (dashed line, blue). 

Open with DEXTER  
In the text 
Fig. 2 Comparison between expansion rates and weak reaction rates. The shortdashed green line represents the neutrontoproton conversion rate. The two solid lines represent the expansion rate of the DiracMilne universe (red) and the standard cosmology (blue). The long dashed lines represent the protontoneutron conversion rates. Weak interactions decouple when these rates become lower than the expansion rate. 

Open with DEXTER  
In the text 
Fig. 3 ^{4}He production rate for the DiracMilne (dashed red line). The rate drops abruptly when weak interactions decouple around T ~ 80keV. As a comparison, ^{4}He production rate for the standard BBN is also presented (dotted blue line). 

Open with DEXTER  
In the text 
Fig. 4 ^{4}He (top) and ^{7}Li (bottom) theoretical abundances for the DiracMilne universe as a function of the baryon density. Here η_{10} = 10^{10}η. Shaded areas correspond to the observational range. 

Open with DEXTER  
In the text 
Fig. 5 Abundances of light elements obtained in the Milne universe with a baryontophoton ratio η = 8.8 × 10^{9}. Both ^{4}He and ^{7}Li are produced at observationally compatible levels (see Fig. 4), but D and ^{3}He are almost totally destroyed by the very slow thermal nucleosynthesis. 

Open with DEXTER  
In the text 
Fig. 6 Comoving diffusion length (solid red line). The blue dashed line is the comoving thermalizing length of deuterium. Comoving Hubble distance is also shown in green. 

Open with DEXTER  
In the text 
Fig. 7 Quantity of D (solid green) and ^{3}He (dashed red) nuclei produced by ^{4}He nuclei photodisintegration for one annihilation. 

Open with DEXTER  
In the text 
Fig. 8 D and ^{3}He final abundance as a function of the size of the emulsion and the redshift of gravitational decoupling. Hatched region is forbidden because of CMB distortions constraints. 

Open with DEXTER  
In the text 
Fig. 9 Residuals of the Hubble diagram for the DiracMilne (a) and b)), flat ΛCDM (c) and d)), and Einstein deSitter (e) and f)) models. The left column represents the residuals obtained by the minimization of the χ^{2} defined by Eq. (25), in which the values of the intrinsic dispersion is adjusted so that χ^{2}/d.o.f. = 1. In the right column, this intrinsic dispersion is fixed to 0. 

Open with DEXTER  
In the text 