New Dark Matter Generation Mechanism and Its Implications for the Cosmological Background

A new concept has been proposed and elaborated to account for recent observations deviating from ΛCDM and ΛWDM. Using an intermediate energy conversion mechanism in the expanding universe and assuming three neutrino families with identical mass, a neutrino mass of 1.19 0.19 ± eV/c has been found as well as a natural explanation for the difference in Hubble constant 0 H as measured by WMAP/Planck and obtained from LSS data. The value for the effective number of neutrinos eff N at the time of decoupling is compatible with the Planck result. The age of the universe is slightly younger at 13.5 0.5 ± Gyr. At late times, the closure parameter for the neutrino radiation dmr Ω drastically increases but still remains well below that of the baryons b Ω , among other energy contributions.


Introduction
The six parameter ΛCDM model [1], where CDM stands for Cold Dark Matter and Λ refers to the cosmological constant, has long been regarded as the standard model of modern cosmology.Despite its success in explaining many features of the Cosmic Microwave Background (CMB) power spectrum [2] among others, increasingly precise data from observations are more and more at odds with it.
In the ΛCDM scenario, an important role is played by cold dark matter.Many interesting scenarios have been proposed to explain the true nature of CDM, in-cluding extended theories of gravity [3], though some of them are experiencing difficulties explaining recent experiments and observations.Weakly Interacting Massive Particles (WIMPs) are very good candidates for cold dark matter.They are quite massive (about 100 GeV/c 2 ) and are thought to originate from thermal equilibrium processes in the very early universe.Despite a huge experimental effort, no convincing evidence for them could be found [4].
Another popular candidate is axions [4] which might be generated in the Coulomb field of the nucleon.In contrast to the WIMPs, their masses are expected to be very low (10 −6 -1 eV range) and to allow for structure formation they cannot be in thermal equilibrium with the other particles.Experiments indicate that the original axion scenario [5] does not describe how nature works.
With the currently accepted amount of CDM, N-body simulations [6] result in various phenomena including higher central densities in galaxies (cusp-core problem) [7] and more satellites (missing satellite problem) [8] than observed.One possible explanation comes from sterile neutrinos [9] in the keV range, referred to as Warm Dark Matter (WDM).They are assumed to not interact according to the known interactions.Boltzmann codes reveal, on the other hand, that there are a number of deviations in the CMB power spectrum such as the quadrupole anomaly, a lower Temperature-Temperature (TT) power spectrum at 20 -30 l = and a slightly enhanced first acoustic peak [2] [10].From oscillation experiments, we also know that at least two of the three neutrino families must be massive [11].
In this paper, dark matter is examined in a different cosmological model, using light dark matter particles in the eV range i.e.Hot Dark Matter (HDM), but still allowing for structure formation as observed.
First, a new concept is proposed for dark matter generation, starting from well-known physics and a single main additional assumption.Subsequently, the formulas used in the simulations to calculate the background are presented.In the third part, the simulations and results are dealt with, with the emphasis on the approach rather than the details.Finally, the predictions of these simulations are discussed, and three tests are proposed to compare theoretical results with observations.

Concept and Assumptions
The first Friedmann equation [12] describes how the scale factor evolves in a homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker (FLRW) universe with matter, radiation and dark energy (cosmological constant) and is essentially an energy equation [4].As energy is conserved locally in general relativity, one could imagine a mechanism converting kinetic and potential energy into mass.In the whole approach almost everything is based on well-known physics and there are very few assumptions.The main assumption is to take the proton with its characteristic radius of 0.88 fm p r = and typical string constant 1 R = GeV/fm [13] as representative for an underlying mechanism which is unknown.See Section 4 for a justification of the value of R.This mechanism is thought to be capable of following the expansion of the universe temporarily and in doing so building up some small amount of energy, which is eventually released by generating neutral weakly interacting massive particle/antiparticle pairs with their initial kinetic energies depending on cosmological time.
Note that their center of mass is that of the nucleon they originate from, an indispensable feature to allow for structure formation in the universe with light HDM particles which is otherwise not possible [4].This concept does not only allow for calculating the evolution of the characteristics of the universe but also results in specific predictions of, among others, the masses of the dark matter particles themselves.In addition to this main assumption, the universe is considered flat today with the possibility to vary towards the past and the Hubble constant, as obtained from measurements in the local universe, is taken as an input.

Calculation of Background
The evolution of the scale factor ( ) a t of the expanding universe is described by the Friedmann equation [14] which can be derived from the 00-component of the Einstein field equations of general relativity.Accounting for curvature it reads [12] ( ) ( ) where H is the Hubble parameter, tot ρ stands for the total energy density, c is the speed of light, G denotes Newton's constant and κ is a parameter representing curvature, which is negative in case of an open universe with sub-critical density.
The curvature term can also be represented by a fictive density [12].The various energy contributions considered in this paper are the radiation from photons (subscript g), baryons (b), Standard Model (SM) neutrino radiation (nr) and neutrino mass (nm), generated dark matter radiation (dmr) and dark matter mass (dmm) as well as curvature (k) and dark energy (de).The SM neutrinos and the generated dark matter particles will turn out to be the same particles but their abundance and time and means of generation (non-thermal versus thermal) differ.
Taking into account the dependencies of the various species on the scale factor (see for example [4]) and introducing the closure parameter defined as ( ) The subscript 0 refers to the current cosmological time and all z dependent Journal of High Energy Physics, Gravitation and Cosmology quantities are calculated through to the present time.It should be noted already that ( ) ,0 k z Ω depends on the redshift and is only zero at the present cosmological time 0 t .Similarly, the contributions with dark matter have a z dependence as well.The amount of dark matter radiation can be inferred from the total amount of energy generated per decay minus the energy associated with the mass of the corresponding dark matter particles.In view of the conservation of the "center of energy" [15], it is considered as a matter contribution.
The photon (  ) abundances at 0 t follow directly from the equilibrium conditions in the early un- iverse (see for example [4]).Relating the corresponding photon density to the critical density for a Hubble constant 0 73.24 where m ν [eV] is the neutrino mass, h is the reduced Hubble constant ( 0 100 H ≡ , dimensionless) and the expression accounts for non-instantaneous decoupling.The fact that there are three families has been considered by inclusion of a factor of three.
For the generated dark matter and related curvature considered here, the situation is different and some care is required as the number densities and curvature value are dynamical and deviate from the common scaling.
( ) For the dark matter radiation, on the other hand,

( )
dmr z ρ has to be compared to the photon density ( ) Quantized packets of dark matter and radiation are calculated for logarithmically spaced intervals.At a specific redshift z the total amount of non-electromagnetically interacting matter and radiation consists of all previous contributions properly scaled to the redshift of interest as well as that of the new packet generated at that scale.These quantities can be converted then with Equation ( 4) and ( 5) in the desired z dependent closure parameters for the generated dark matter and radiation at 0 t .The closure parameter for the curvature is such that, when considered as a density, it makes the universe flat at that redshift, i.e. and dark energy de Ω .Note that the universe is assumed to be critical only at 0 t by requiring that the total amount of dark matter today is reached, thereby setting the curvature at the current cosmological time to zero.
The dark matter particles are gradually generated during the expansion history of the universe, starting at the moment of baryogenesis (taken here as 1 .In the initial stage enough energy has to be accumulated to account for the particle/antiparticle masses.After that, some additional energy build-up is required to get the system to decay.The average time to do so is crucial in order to calculate the evolution of the dark matter in the universe according to the concept presented here. In case of the weak interaction, combining the decay width of a (virtual) * Z boson with energy ) [19] with the correction for the fermion masses m ν and m ν as obtained from [20], the decay width The factor of three stems from the assumption of three neutrino families with degenerate masses and F G is Fermi's constant.The corresponding decay time directly follows from the relation τ ≡ Γ  with  the reduced Planck constant.
The peculiarity here is that * Z m , which represents the average energy available to a particular decay is not a constant but is steadily increasing.Its time dependence is given by where Δ m t is the time for mass generation of a particle and antiparticle.Note that the Hubble parameter H a a ≡  is used and not the derivative with respect to the scale factor a. See for example [12].
Finding the time for decay r t ∆ is a little more elaborate.Starting from the 1 Note already that the predictions and results obtained are quite robust against this value.Journal of High Energy Physics, Gravitation and Cosmology general expression for the probability of a decay ( ) , where τ stands for the decay constant, the likelihood that a decay occurs in an interval d t t + is given by the probability that it does not happen in the first n intervals dt multiplied by the probability that it occurs in the next interval dt .The probability ( ) where dt t n = has been substituted and t is the total elapsed time for decay.This approximation becomes quite accurate for n sufficiently large e.g.1000 n = .
The  It is well known that, as the universe expands, particles cool down and eventually become non-relativistic if they are massive.However, the type of distribution remains preserved.The key point is that the momentum is inversely proportional to the scale factor a [12] so that the total energy tot t ∆ of a particle with mass m scales as [1] ( ) ( ) where the subscript em stands for emission.Subtracting the mass energy The amount of radiation from relativistic particles is often expressed as an effective number of relativistic neutrinos.It is implicitly defined as [18] where ν ρ stands for all types of radiation, whether dmr ρ or nr ρ , from relati- vistic particles and g ρ is the photon density.For the SM neutrinos, the radia- tion density nr ρ is obtained from [14]   ( ) ( ) ( ) T T a ν ν = the momentum and neutrino temperature respectively and B k Boltzmann's constant.g ρ can be de- rived directly from Boltzmann thermodynamics and obeys the relation [4] ( ) are the spin degrees of freedom for the neutrinos and photons respectively.Note that the Fermi-Dirac distribution of the neutrinos is independent of their mass.Strictly speaking Equation ( 12) applies to particles with a thermal distribution but it is also used here to determine eff N for the generated neutrinos.Their total radiation density is obtained as described above.
Finally, a least squares method has been applied to astrophysical data to arrive at a justification of the value of R. To account for the magnitude of the error bars, the following fitting formula has been used [21] where

Simulations and Results
As in the previous section, we restrict ourselves to the evolution of the background parameters.The main task is to solve the Friedmann equation-Equation (2) numerically with the approach and formulas as presented in Section 3.All (other) background quantities can then be derived from it.For this purpose numerical codes have been written in C++ and Python.What happens in the very early universe essentially is irrelevant for our purposes here.Therefore, the simulations start at In the Friedmann equation-Equation ( 2) the evolution of the densities of the standard components are implemented in the conventional way with their typical scaling.However, the dark matter fraction generated starts at zero and there is some curvature related to this deviation from "flatness".From an expansion point of view, the latter components are irrelevant as the universe is completely radiation dominated at that time.Integration is forward in time and as time evolves, dark matter and radiation are generated and curvature diminishes correspondingly.
After each integration loop, the z dependent dynamical contributions are updated in Equation ( 2) before a next cycle starts.This is iterated until 1 a = .The input parameters are ( ) The age of the universe follows directly from the simulations and turns out to be with its 13.5 0.5 ± Gyr slightly younger than in the SM of cosmology [2].
This value might still be compatible with the age of the oldest globular clusters and stars.The redshift of matter-radiation equality is now 2122 214 and dark energy takes over at 0.39 0.06 Figure 1 illustrates how the various contributions to the energy budget evolve for the concept presented here with the well-known epochs of radiation domination at early times (small scale factor), followed by matter domination and finally, very recently, a universe dominated by dark energy.A distinction is made between the neutrinos created from equilibrium reactions (SM) in the early universe (thermal) and those gradually generated during the expansion history (non-thermal).Note that, from an energy point of view, the deviation from flatness, represented by the curvature contribution, peaks around 0.25 a ≈ and never exceeds more than a few percent.
With the expansion history determined by the (dynamical) evolution of the different components, the energy contributions of dark radiation and matter as generated by the concept used here, can be plotted as a function of a. Figure 2 shows the result.The instantaneous situation is plotted; cooling due to expansion has not been considered in this figure.At early times the amount of radiation exceeds by far the matter contribution.However, the impact remains limited as the dark matter quantity is very low at this stage and radiation quickly decays as time evolves.A main part of the dark matter is generated in the last decade and in absolute value the total dark matter radiation also continuously increases.See also Figure 4.In the concept introduced, dark matter particles are thought to be representatively generated by an energy conversion mechanism where the nucleon potential plays a central role.Decays occur statistically and, even for a fixed cosmological time, the momenta of the decay products (neutrinos and antineutrinos) are distributions.peak values increase and distributions sharpen as time evolves, as might be expected.Equation ( 9) has been used with 800 n = as a good approximation everywhere.
The amount of radiation energy is often presented as an effective number of neutrinos eff N , as if they were massless.Figure 4 shows this evolution, once again as a function of the scale factor, as well as that of the neutrino number densities rescaled for today.For many decades in a, but a rather short period in time, the thermal neutrinos provide the main contribution to the total amount of dark matter.From Equation (3), the corresponding closure parameter is Additional structure formation power comes primarily from the fact that the generated dark matter originates from the baryons and inherits their instantaneous density perturbations due to the center of energy principle of the (relativistic) particle/antiparticle configuration.Consequently, baryons and dark matter are more strongly related than usual.The fact that dark matter is thought to originate from the baryon overdensities might also explain the mysterious conspiracy of both components as observed [22].
In ΛCDM the evolution of the universe is determined by one or the other species for long times and transitions are rather short.If dark matter is considered to be continuously generated and the neutrinos contribute considerably, the story is different.Figure 5 shows the evolution of the equation of state parameter for all species together ( tot w ) and the neutrinos (thermal and non-thermal) specifically ( dm n w + ).Equation ( 11) is used where the dynamical z dependent closure parameters are obtained from the simulations.Note the rather large transition regions.Up to this point, the amount of dark matter and baryons, constituting the total matter content, have been considered as fixed.It is a promising sign that the outcome is so good with a single additional assumption and without fine tuning.
Simulations allow us to vary the matter contribution and quantify the impact on the main quantities of interest here.Figure 6 illustrates how the simulated neutrino mass m ν , the effective number of neutrinos at the CMB and 0 t as well as the age of the universe 0 t evolve, if the total amount of matter ,0 m Ω is varied.
According to the simulations a large matter fraction above 0.35 is not possible in view of the resulting young age of the universe, the small value for A key parameter in the model, and the only new one here, is the string constant R which is monotonically related to the neutrino mass m ν .Higher values of R result in lower mass bounds.To investigate the influence of the value adopted here, the proposed model has been fitted to both ( ) H z measurements [23] and the Union2.1 compilation [24] simultaneously, using a least squares method.
Allowing the field strength to vary in the range 0.1 -2.2 GeV/fm in steps of 0.1 GeV/fm and considering the baryon content in between 0.042 and 0.049 (steps of 0.001), the total amount of matter in the range 0.20 -0.32 (steps of 0.01) and 0 H varying from 67.5 to 75 km/s/Mpc (steps of 0.5 km/s/Mpc) parameter space has been scanned to arrive at the best fit.A least squares method has been used supplemented by the criterion 3.84 eff N = [10].As might be expected, the result is fairly independent of b Ω .The best input parameters are then found to be Starting from the Planck results [2] The result of the best fit has been shown in Figure 7 for ΛCDM (red line) and the proposed concept (green line).The best fit input parameters used are

Predictions
A first firm prediction of the analytical work and simulations presented here is the mass of the dark matter particle ( 1.19 0.19 m ν = ± eV/c 2 ) and the fact that they are generated by the neutral current weak interaction and therefore are thought to be the standard model neutrinos.Note that the cosmology presented here is quite different from ΛCDM, which impacts the predictions considerably.
With a value of 1.19 m ν = eV/c 2 the predicted neutrino mass is just outside the range of experiments performed by Troitsk [25] and Mainz [26] some years ago but the Katrin experiment [27] starting up at this moment should be able to confirm or falsify this prediction in the very near future.
A second prediction is related to the effective number of neutrinos eff N and much effort is done to experimentally determine its value at various stages of the cosmological evolution.As can be inferred from Figure 4 eff N starts at the SM value of 3.046, becomes marginally larger at the CMB ( It is well known that there is some tension between the value of the Hubble constant as obtained from CMB data and that inferred from observations in the recent universe, see for example [16].In addition, as explained in [2], a positive correlation exists between

Summary
A new concept has been proposed for the generation of dark matter, based on the nucleon with its characteristic radius and potential, which is thought to be representative of an underlying mechanism converting energy from the expanding universe into dark matter and radiation.Energy build-up first accounted for the mass to be generated and was then subject to a decay which has been assumed to be mediated by the weak interaction.Furthermore, the universe has been taken to be flat today but allowing for a limited amount of curvature (open universe) towards the past, corresponding to the deficiency of dark matter at that redshift.
The background has been calculated/simulated from the first Friedmann equation with specific contributions for the standard thermal neutrinos and the generated dark matter.
Using commonly accepted input parameter values and the constraints required by the concept, the dark matter was found to be composed of the standard neutrinos but with a mass of momentum distribution.It will turn out to be fairly sharply peaked and can be replaced by its average momentum.This simplifies the computational effort considerably.
or antineutrino is readily obtained.To find the evolution of the radiation density dmr ρ of the neutrino population generated up to that time, all particle kinetic energies have to be rescaled correspondingly.A non-thermal distribution results, with a total radiation densioften used as an indicator for new physics beyond the ΛCDM model is the effective number of relativistic degrees of freedom eff N .It is positively correlated with the Hubble constant[2].As a rule of thumb,2 eff N h ∝has been adopted when applied to the CMB data.As we will see in the next sec-H.Defloor tion the results obtained here might give a natural explanation for the difference in Hubble constant as derived from CMB and Large Scale Structure (LSS) data.The content of the universe can be considered as a fluid obeying the equation of state equation of state (e.o.s.) parameter.It applies to a single component as well as to all species together.In the latter case it follows from the weighted average as is the resulting averaged value and the i w s are the e.o.s.parameters for the individual components.All types of radiation obey , i.e. a cosmological constant Λ.In view of conservation of the center of energy, both the generated dark matter and corresponding radiation are classified as 0 m w = , as are the baryons they originate from.
is sufficiently early as the amount of generated dark matter at that time is close to zero; the exact value of s t is non-critical and does not influence the results.

1 R
/Mpc as obtained from LSS data[16].Best results are found for a nucleon potential with = GeV/fm.These values apply throughout this paper.Requiring that the total amount of matter must be generated by 1 a = and that the universe must be flat today, the mass of the dark matter particles can be determined and the evolution of all background quantities and parameters is completely fixed.Error bars are found with a Monte Carlo method and starting from Gaussian distributions for the input parameters b Ω , out to be compatible with the observed energy contributions in the universe with the dark matter particles consisting of the standard model neutrinos (three families) with a (nearly) degenerate mass of 1.19 0), which is still within the matter dominated era.It should be noted however that , eff CMB N is composed differently here and consists of a contribution from the "SM neutrinos" ( , , 0.84 0.11 eff nr CMB N = ± ) and that of the generated dark matter neutrinos ( the corresponding closure parameter still Journal of High Energy Physics, Gravitation and Cosmology remains well below that of the baryons, among other energy contributions.

Figure 1 .
Figure 1.Evolution of the various energy components as a function of the scale factor a for the concept presented here.Contributions are as labelled and for the input values specified in the text.Total density is normalized to 1.

Figure 3 Figure 2 .
Figure 2. Relative importance of mass and radiation, instantaneously generated by the approach taken, as a function of a.The mass fraction is continuously increasing.

Figure 3 .
Figure 3. Momentum distributions (probability) as a function of the Hubble parameter as they are at their moment of creation.Note also the direction of time.
times, generated dark matter constitutes the main contribution.Simulations have been performed so that, at

Figure 4 .
Figure 4. Effective numbers of neutrinos eff N and neutrino number densities for thermal and generated neutrinos as a function of the scale factor.

Figure 5 .
Figure 5. Evolution of the equation of state parameters for all species together and the neutrinos (dark matter and standard neutrinos) specifically.

H
neutrino mass which would be in conflict with neutrino mass experiments.Finally, one could go a step further and try to fit the main model input parameters, the baryon content b Ω , the matter contribution ,0 m Ω , the local Hubble constant 0 as well as the potential R to astrophysical data.

Figure 6 .
Figure 6.Simulation results for the evolution of neutrino mass m ν , predicted age of the universe 0 t and effective number of neutrinos at the CMB and 0 t as a function of the total matter content ,0 m Ω .

1 R
= GeV/fm for the model proposed here.In addition, the model considered is assumed to be flat only today and open towards the past but with fixed de Ω .Equation (15) has been used.Note that the two fitted lines are indistinguishable in the lower plot.Least squares values are ( )

Figure 7 . 1 R
Figure 7. Best fit for ΛCDM (red) and proposed model (green) in case of ( ) H z and SN1a data respectively.Error bars are 1σ.Input parameter values are dramatically increases only in the last decade of the scale factor ( explained in Section 4 the value at the CMB consists of a contribution of the thermal and non-thermal neutrinos.Note that the kinetic H. Defloor energies of the neutrinos are very low (eV range at maximum) and therefore difficult to detect directly anyway.Comparison to results such as from Planck ( might indicate whether the presented model is potentially right.
Compatible values for the effective number of relativistic degrees of freedom were found around 3.9eff N ≈ which, starting from the Planck results[2], would imply a Hubble constant around 73 km/s/Mpc.This is a third prediction and consequence of the model presented here.In view of the progress in cosmological measurements, predictions from experiment might reveal the real nature of the Hubble constant and tell us more about the validity of our results.

2 χ
values of 21.0 and 562.4 respectively.Compared to those of ΛCDM, using the same parameter inputs, the proposed model has yielded a better fit to the ( ) H z data and an equivalently good one to the SN1a data.
value for i P t and almost coincides with its maximum.The resulting time for the decay process itself amounts then to .The mass energy generated is always the same but the amount of radiation for a specific decay varies and depends on the random decay character.Comparing the distribution for All these results are compatible with currently observed values or, at least, not excluded by them as far as the neutrino mass is concerned.