# Big Bang nucleosynthesis with a non-Maxwellian distribution

###### Abstract

The abundances of light elements based on the big bang nucleosynthesis model are calculated using the Tsallis non-extensive statistics. The impact of the variation of the non-extensive parameter from the unity value is compared to observations and to the abundance yields from the standard big bang model. We find large differences between the reaction rates and the abundance of light elements calculated with the extensive and the non-extensive statistics. We found that the observations are consistent with a non-extensive parameter , indicating that a large deviation from the Boltzmann-Gibbs statistics () is highly unlikely.

## I Introduction

The cosmological big bang model is in agreement with many observations relevant for our understanding of the universe. However, comparison of calculations based on the model with observations is not straightforward because the data are subject to poorly known evolutionary effects and systematic errors. Nonetheless, the model is believed to be the only probe of physics in the early universe during the interval from min, after which the temperature and density of the universe fell below that which is required for nuclear fusion and prevented elements heavier than beryllium from forming. The model is inline with the cosmic microwave background (CMB) radiation temperature of 2.275 K Not11 (), and provides guidance to other areas of science, such as nuclear and particle physics. Big bang model calculations are also consistent with the number of light neutrino families . According to the numerous literature on the subject, the big bang model can accommodate values between (see, e.g., Ref. Oli02 ()). From the measurement of the width by LEP experiments at CERN one knows that LEP06 ().

In the big bang model nearly all neutrons end up in He, so that the relative abundance of He depends on the number of neutrino families and also on the neutron lifetime . The sensitivity to the neutron lifetime affects Big Bang Nucleosynthesis (BBN) in two ways. The neutron lifetime influences the weak reaction rates because of the relation between and the weak coupling constant. A shorter (longer) means that the reaction rates remain greater (smaller) than the Hubble expansion rate until a lower (larger) freeze-out temperature, having a strong impact on the equilibrium neutron-to-proton ratio at freeze-out. This -ratio is approximately given in thermal equilibrium by , where is the Boltzmann constant, the temperature at weak freeze-out, and is the neutron-proton mass difference. The other influence of is due to their decay in the interval between weak freeze-out ( s) and when nucleosynthesis starts ( s), reducing the ratio to . A shorter implies lower the predicted BBN helium abundance. In this work we will use the value of s, according to the most recent experiments Ser05 () (a recent review on the neutron lifetime is found in Ref. FG11 ()). Recently, the implications of a change in the neutron lifetime on BBN predictions have been assessed in Ref. Mat05 ().

The baryonic density of the universe deduced from the observations of the anisotropies of the CMB radiation, constrains the value of the number of baryons per photon, , which remains constant during the expansion of the universe. Big bang model calculations are compatible with the experimentally deduced value of from WMAP observations, Ko11 ().

Of our interest in this work is the abundances of light elements in big bang nucleosynthesis. At the very early stages (first 20 min) of the universe evolution, when it was dense and hot enough for nuclear reactions to take place, the temperature of the primordial plasma decreased from a few MeV down to about 10 keV, light nuclides as H, He, He and, to a smaller extent, Li were produced via a network of nuclear processes, resulting into abundances for these species which can be determined with several observational techniques and in different astrophysical environments. Apparent discrepancies for the Li abundances in metal poor stars, as measured observationally and as inferred by WMAP, have promoted a wealth of new inquiries on BBN and on stellar mixing processes destroying Li, whose results are not yet final. Further studies of light-element abundances in low metallicity stars and extragalactic H II regions, as well as better estimates from BBN models are required to tackle this issue, integrating high resolution spectroscopic studies of stellar and interstellar matter with nucleosynthesis models and nuclear physics experiments and theories Fie11 ().

The Maxwell-Boltzmann distribution of the kinetic energy of the ions in a plasma is one of the basic inputs for the calculation of nuclear reaction rates during the BBN. The distribution is based on several assumptions inherent to the Boltzmann-Gibbs statistics: (a) the collision time is much smaller than the mean time between collisions, (b) the interaction is local, (c) the velocities of two particles at the same point are not correlated, and (d) that energy is locally conserved when using only the degrees of freedom of the colliding particles (no significant amount of energy is transferred to and from collective variables and fields). If (a) and (b) are not valid the resulting effective two-body interaction is non-local and depends on the momentum and energy of the particles. Even when the one-particle energy distribution is Maxwellian, additional assumptions about correlations between particles are necessary to deduce that the relative-velocity distribution is also Maxwellian. Although the Boltzmann-Gibbs (BG) description of statistical mechanics is well established in a seemingly infinite number of situations, in recent years an increasing theoretical effort has concentrated on the development of alternative approaches to statistical mechanics which includes the BG statistics as a special limit of a more general theory Ts88 () (see also, Ren60 ()). Such theories aim to describe systems with long range interactions and with memory effects (or non-ergodic systems). A very popular alternative to the BG statistics was proposed by C. Tsallis Ts88 (); GT04 (), herewith denoted as non-extensive statistics (for more details on this subject, see the extensive reviews GT04 (); TGS05 (); Ts09 ()). Statistical mechanics assumes that energy is an “extensive” variable, meaning that the total energy of the system is proportional to the system size; similarly the entropy is also supposed to be extensive. This might be justified due to the short-range nature of the interactions which hold matter together. But if one deals with long-range interactions, most prominently gravity; one can then find that entropy is not extensive FL01 (); Lim02 (); TS03 (); TS04 (); CS05 ().

In classical statistics, to calculate the average values of some quantities, such as the energy of the system, the number of molecules, the volume it occupies, etc, one searches for the probability distribution which maximizes the entropy, subject to the constraint that it gives the right average values of those quantities. As mentioned above, Tsallis proposed to replace the usual (BG) entropy with a new, non-extensive quantity, now commonly called the Tsallis entropy, and maximize that, subject to the same usual constraints. There is actually a whole infinite family of Tsallis entropies, indexed by a real-valued parameter , which quantifies the degree of departure from extensivity (one gets the usual entropy back again when ). It was shown in many circumstances that the classical results of statistical mechanics can be translated into the new theory Ts09 (). The importance of these families of entropies is that, when applied to ordinary statistical mechanics, they give rise to probabilities that follow power laws instead of the exponential laws of the standard case (for details on this see Ts09 ()). In most cases that Tsallis formalism is adopted, e.g. Ref. Pes01 (), the non-extensive parameter is taken to be constant and close to the value for which ordinary statistical mechanics is obtained (). Some works have also probed sizable deviations of the non-extensive parameter from the unity to explain a variety of phenomena in several areas of science GT04 ().

In the next sections, we shown that the Maxwell-Boltzmann distribution, a cornerstone of the big bang and stellar evolution nucleosynthesis, is strongly modified by the non-extensive statistics if strongly deviates from the unity. As a consequence, it also affects strongly the predictions of the BBN. There is no “a priori” justification for a large deviation of from the unity value during the BBN epoch. In particular, as radiation is assumed to be in equilibrium with matter during the BBN, a modification of the Maxwell distribution of velocities would also impact the Planck distribution of photons. Recent studies on the temperature fluctuations of the cosmic microwave background (CMB) radiation have shown that a modified Planck distribution based on Tsallis statistics adequately describes the CMB temperature fluctuations measured by WMAP with , which is close to unity but not quite Ber07 (). Perhaps more importantly, Gaussian temperature distributions based on the BG statistics, corresponding to the limit, do not properly represent the CMB temperature fluctuations Ber07 (). Such fluctuations, allowing for even larger variations of might occur during the BBN epoch, also leading to a change of the exponentially decaying tail of the Maxwell velocity distribution.

Based on the successes of the big bang model, it is fair to assume that it can set strong constraints on the limits of the parameter used in a non-extensive statistics description of the Maxwell-Boltzmann velocity distribution. In the literature, attempts to solve the lithium problem has assumed all sorts of “new physics” Fie11 (). The present work adds to the list of new attempts, although our results imply a much wider impact on BBN as expected for the solution of the lithium problem. If the Tsallis statistics appropriately describes the deviations of tails of statistical distributions, then the BNN would effectively probe such tails. The Gamow window (see figure 1) contains a small fraction of the total area under the velocity distribution. Thus, only a few particles in the tail of the distributions contribute to the fusion rates. In fact, the possibility of a deviation of the Maxwellian distribution and implications of the modification of the Maxwellian distribution tail for nuclear burning in stars have already been explored in the past MQ05 (); HK08 (); Deg98 (); Cor99 (). As we show in the next sections, a strong deviation from is very unlikely for the BNN predictions, based on comparison with observations. Moreover, if deviates from the unity value, the lithium problem gets even worse.

## Ii Maxwellian and non-Maxwellin distributions

Nuclear reaction rates in the BBN and in stelar evolution are strongly dependent of the particle velocity distributions. The fusion reaction rates for nuclear species 1 and 2 is given by , i.e., an average of the fusion cross section of with their relative velocity, described by a velocity distribution. It is thus worthwhile to study the modifications of the stellar reaction rates due to the modifications introduced by the non-extensive statistics.

### ii.1 Non-extensive Statistics

Statistical systems in equilibrium are described by the Boltzmann-Gibbs entropy,

(1) |

where is the Boltzmann constant, and is the probability of the i-th microstate. For two independent systems , , the probability of the system being in a state , with a microstate of and a microstate of , is

(2) |

Therefore, the Boltzmann-Gibbs entropy satisfies the relation

(3) |

Thus, the entropy based on the Boltzmann-Gibbs statistic is an extensive quantity.

In the non-extensive statistics Ts88 (), one replaces the traditional entropy by the following one:

(4) |

where is a real number. For , . Thus, the Tsallis statistics is a natural generalization of the Boltzmann-Gibbs entropy.

Now it follows that

(5) |

The variable is a measure of the non-extensivity. Tsallis has shown that a formalism of statistical mechanics can be consistently developed in terms of this generalized entropy Ts09 ().

A consequence of the non- extensive formalism is that the distribution function which maximizes is non- Maxwellian Sil98 (); Lim00 (); Mu06 (). For , the Maxwell distribution function is reproduced. But for , high energy states are more probable than in the extensive case. On the other hand, for high energy states are less probable than in the extensive case, and there is a cutoff beyond which no states exist.

### ii.2 Maxwellian Distribution

In stars, the thermonuclear reaction rate with a Maxwellian distribution is given by Fow67 ()

(6) | |||||

where is the fusion cross section, is the relative velocity of the -pair, is the number of nuclei of species , is the reduced mass of , is the temperature, is the astrophysical S-factor, and is the Sommerfeld parameter, with the i-th nuclide charge and is the relative energy of .

The energy dependence of the reaction cross sections is usually expressed in terms of the equation

(7) |

We write , where

(8) |

where is the reduced mass in amu. The factor in the denominator of Eq. (6) corrects for the double-counting when . The S-factor has a relatively weak dependence on the energy , except when it is close to a resonance, where it is strongly peaked.

### ii.3 Non-Maxwellian Distribution

The non-extensive description of the Maxwell-Boltzmann distribution corresponds to the substitution , where Ts09 ()

(9) | |||||

If , Eq. (9) is real for any value of . However, if , is real only if the quantity in square brackets is positive. This means that

(10) |

Thus, in the interval one has and for one has .

With this new statistics, the reaction rate becomes

(11) |

and the rate integral, , is given by

(12) |

where the “modified” Gamow energy distribution is

(13) | |||||

is the non-extensive Maxwell velocity distribution, for and for , and in MeV units. is a normalization constant which depends on the temperature and the non-extensive parameter .

### ii.4 Non-Maxwellian distribution for relative velocities

It is worthwhile to notice that if the one-particle energy distribution is Maxwellian, it does not necessarily imply that the relative velocity distribution is also Maxwellian. Additional assumptions about correlations between particles are necessary to deduce that the relative-velocity distribution, which is the relevant quantity for rate calculations, is also Maxwellian. This has been discussed in details in Refs. MQ05 (); Ts88 (); Kan05 () where non-Maxwellian distributions, such as in Eq. (13), were shown to arise from non-extensive statistics.

Here we show that, if the non-Maxwellian particle velocity distribution is given by Eq. (9), then a two-particle relative can be modified to account for the center of mass recoil. Calling the kinetic energy of a particle , this distribution is given by

(14) |

The two-particle energy distribution is . We now exponentiate the Tsallis distribution.

(15) |

and the product reduces to

(16) |

Since , and thus, , where is the reduced mass of the two particles, , is the relative velocity, and the center of mass velocity, the product inside the natural logarithm can be reduced to

(17) |

Thus, the two-body distribution factorizes into a product of relative and center of mass parts

(18) |

where

(19) |

with the normalization constants obtained from the condition,

(20) |

Because the distribution factorizes, the unit normalization can be achieved for the relative and c.m. distributions separately. The distribution needed in the reaction rate formula is, therefore,

(21) |

which attains the same form for as the absolute distribution.

In the limit the two-particle distribution reduces to a Gaussian, with the last term in the left-hand-side of Eq. (17) dropping out,

(22) |

as expected.

### ii.5 Equilibrium with electrons, photons and neutrinos

One of the important questions regarding a plasma with particles (i.e., nuclei) described by the non-extensive statistics is how to generalize Fermi-Dirac, Bose-Einstein, and Tsallis statistics, to become more unified statistics with the distribution for the particles. This has been studied in Ref. BD00 (), where it was shown that a similar non-extensive statistic for the distribution can be obtained for fermions and bosons is given by

(23) |

where is the chemical potential. This reproduces the Fermi distribution, , for and the Bose-Einstein distribution for photons, , for and . Planck’s law for the distribution of radiation is obtained by multiplying in Eq. (23) (with ) by , where . The number density of electrons can be obtained from in Eq. (23) with the proper phase-factors depending if the electrons are non-relativistic or relativistic. Normalization factors also need to be introduced, as before.

The electron density during the early universe varies strongly with the temperature. At the electron density is about /cm, much larger than the electron number density at the center of the sun, /cm. The large electron density is due to the production by the abundant photons during the BBN. However, the large electron densities do not influence the nuclear reactions during the BBN. In fact, the enhancement of the nuclear reaction rates due to electron screening were shown to be very small WBB11 (). The electron Fermi energy for these densities are also much smaller than for the energy relevant for BBN, so that one can also use in Eq. (23) for .

In figure 2 we plot the the relative difference ratio between non-extensive, , and extensive, , statistics. For both distributions, we use . The solid curves are for Fermi-Dirac (FD) statistics, , and dashed curves are for Bose-Einstein (BE) statistics, . We show results for and , with . One sees that the non-extensive distributions are enhanced for and suppressed for , as compared to the respective FD and BE quantum distributions. The deviations for the FD and BE statistics also grow larger with the energy. For example, the non-extensive electron distribution is roughly a factor 2 larger than the usual FD distribution at MeV, at .

While we don’t obtain numerical results with modified Fermi-Dirac and Bose-Einstein distributions here, it can be expected that these generalizations will have a strong influence on the freezout temperature and the neutron to proton, n/p, ratio. A numerical study of this problem may be presented in another paper. The freezout temperature, occurs when the rate, , for weak reaction becomes slower than the expansion rate of the Universe. Because during the BBN the densities of all particles, including neutrinos, are low compared to , the chemical potential can be set to zero in the calculation of all reaction rates. The adoption of non-extensive quantum distributions such as in Eq. (23) will lead to the same powers of the temperature as those predicted by the FD distribution and the Bose-Einstein distribution. For example, Planck law for the total blackbody is , being form invariant with respect to non-extensivity entropic index which determines the the degree of non-extensivity BSD02 (). This result means that the weak decay reaction rates do not depend on the non-extensive parameter . The freezout temperature and n/p ratio remain the same as before.

More detailed studies have indicated that Planck’s law of blackbody radiation and other thermodynamical quantities arising from non-extensive quantum statistics can yield different powers of temperature than for the non-extensive case Ara03 (); Nau03 (); Tsa04 (); CDT04 (). If that is the case, then a study of the influence of non-extensive statistics on the weak-decay rates and electromagnetic processes during BBN is worth pursuing.

### ii.6 Thermodynamical equilibrium

The physical appeal for non-extensivity is the role of long-range interactions, which also implies non-equilibrium. Accepting non-extensive entropy means abandoning the most important concept of thermodynamics, namely the tendency of any system to reach equilibrium. This also means that the concept of non-extensivity means renouncing to the second law of thermodynamics altogether!

The comments above, which seem to be shared by part of the community (see, e.g., Nau03 (); ZM03 (); ZM04 (); BDR06 (); Dau07 (); Tou13 ()), are worrisome when one has to consider a medium composed of particles obeying classical and quantum statistics. It is not clear for example if the non-extensive parameter has to be the same for all particle distributions, both classical and quantum. Even worse is the possibility that the temperatures are not the same for the different particle systems in the plasma.

In the present work, we will avoid a longer discussion on the validity of the Tsallis statistics for a plasma such as that existing during the BBN. We will only consider the effect of its use for calculating nuclear reaction rates in the plasma, assuming that it can be described by a classical distribution of velocities. This study will allow us to constrain the non-extensive parameter based on a comparison with observations.

## Iii Reaction rates during big bang nucleosynthesis

Based on the abundant literature on non-extensive statistics (see, e.g., MQ05 (); HK08 (); Deg98 (); Cor99 (); Ts88 ()), we do not expect that the non-extensive parameter differs appreciably from the unity value. However, in order to study the influence of a non-Maxwellian distribution on BBN we will explore values of rather different than the unity, namely, and . This will allow us to pursue a better understanding of the nature of the physics involved in the departure from the BG statistics. In figure 1 we plot the Gamow energy distributions of deuterons relevant for the reaction H(d,p)H at . The solid line, for , corresponds to the use of the Maxwell-Boltzmann distribution. Also shown are results for non-extensive distributions for (dotted line) and (dashed line). One observes that for , higher kinetic energies are more accessible than in the extensive case (). For high energies are less probable than in the extensive case, and there is a cutoff beyond which no kinetic energy is reached. In the example shown in the figure for , the cutoff occurs at MeV, or keV.

We will explore the modifications of the BBN elemental abundances due to a variation of the non-extensive statistics parameter . We will express our reaction rates in the form (in units of cm mol s), where is the Avogadro number and involves the integral in Eq. (6) with the Maxwell distribution replaced by Eq. (9). First we show how the reaction rates are modified for .

In figure 3 we show the S-factor for the reaction H(d,p)H as a function of the relative energy . Also shown is the dependence on (temperature in units of K) for the effective Gamow energy

(24) |

where is the reduced mass in amu. The data are from Refs. SC72 (); KR87 (); BR90 (); BO92 (); GR95 () and the solid curve is a chi-square polynomial function fit to the data.

Using the chi-square polynomial fit obtained to fit the data presented in figure 3, we show in figure 4 the reaction rates for H(d,p)H as a function of the temperature for two different values of the non-extensive parameter . The integrals in equation (12) are performed numerically. For charge particles, a good accuracy (witihin 0.1%) is reached using the integration limits between and , where is given by Eq. (25) below. The rates are expressed in terms of the natural logarithm of (in units of cm mol s). The solid curve corresponds to the usual Maxwell-Boltzmann distribution, i.e., . The dashed and dotted curves are obtained for and , respectively. In both cases, we see deviations from the Maxwellian rate. For the deviations are rather large and the tendency is an overall suppression of the reaction rates, specially at low temperatures. This effect arises from the non-Maxwellian energy cutoff which for this reaction occurs at MeV and which prevents a great number of reactions to occur at higher energies.

For the nearly similar result as with the Maxwell-Boltzmann distribution is due to a competition between suppression in reaction rates at low energies and their enhancement at high energies. The relevant range of energies is set by the Gamow energy which for a Maxwellian distribution is given by Eq. (24) and by the energy window,

(25) |

which for the reaction H(d,p)H amounts to MeV. This explains why, at , the range of relevant energies for the calculation of the reaction rate is shown by the solid curve in figure 1. For the Gamow window is larger and there is as much a contribution from the suppression of reaction rates at low energies compared to the Maxwell-Boltzmann distribution, as there is a corresponding enhancement at higher energies. This explains the nearly equal results shown in figure 4 for and . This finding applies to all charged particle reaction rates, except for those when the S-factor has a strong dependence on energy at, and around, . But no such behavior exists for the most important charged induced reactions in the BBN (neutron-induced reactions will be discussed separately).

The findings described above for the reaction H(d,p)H are not specific but apply to all charged particles of relevance to the BBN. We demonstrate this with one more example: the Li(p,)He reaction, responsible for Li destruction. In figure 5 we show the S-factor for this reaction as a function of the relative energy . One sees prominent resonances at higher energies. Also shown in the figure is the dependence of the reaction on . The data are from Refs. FR58 (); CA62a (); CA62b (); MA64 (); FI67 (); LE69 (); SP71 (); RO86 (); HA89 (); EN92a (); EN92b () and the solid curve is a chi-square function fit to the data using a sum of polynomials plus Breit-Wigner functions.

Using the chi-square function fit obtained to fit the data presented in figure 5, we show in figure 6 the reaction rates for Li(p,)He as a function of the temperature for two different values of the non-extensive parameter . The rates are given in terms of the natural logarithm of (in units of cm mol s). The solid curve corresponds to the usual Maxwell-Boltzmann distribution, i.e., . The dashed and dotted curves are obtained for and , respectively. As with the reaction presented in figure 4, in both cases we see deviations from the Maxwellian rate. But, as before, for the deviations are larger and the tendency is a strong suppression of the reaction rates as the temperature decreases. It is interesting to note that the non-Maxwellian rates for are more sensitive to the resonances than for . This is because, as seen in figure 1, for the velocity distribution is spread to considerably larger values of energies, being therefore more sensitive to the location of high energy resonances.

We now turn to neutron induced reactions, which are only a few cases of high relevance for the BBN, notably the p(n,)d, He(n,p)t, and Be(n,p)Li reactions. For neutron induced reactions, the cross section at low energies is usually proportional to , where is the neutron velocity. Thus, it is sometimes appropriate to rewrite Eq. (7) as

(26) |

where is a slowly varying function of energy similar to an -factor. The distribution function within the reaction rate integral (12) is also rewritten as

The absence of the tunneling factor in Eq. (LABEL:MqET2) inhibts the dependence of the reaction rates on the non-extensive parameter .

In figure 7 we plot the kinetic energy distributions of nucleons relevant for the reaction p(n,)d at (upper panel) and (lower panel). The solid line, for , corresponds to the usual Boltzmann distribution. Also shown are results for the non-extensive distributions for (dotted line) and (dashed line). One observes that, as for the charged particles case, with higher kinetic energies are more probable than in the extensive case (). With high energies are less accessible than in the extensive case, and there is a cutoff beyond which no kinetic energy is reached. A noticeable difference form the case of charged particles is the absence of the Coulomb barrier and a correspondingly lack of suppression of the reaction rates at low energies. As the temperature increases, the relative difference between the Maxwell-Boltzmann and the non-Maxwellian distributions decrease appreciably. This will lead to a rather distinctive pattern of the reaction rates for charged compared to neutron induced reactions.

For neutron-induced reactions, a good accuracy (within 0.1%) for the numerical calculation of the reaction rates with Eq. (12) is reached using the integration limits between and . As an example we will now consider the reaction Be(n,p)Li. The energy dependence of for this reaction is shown in figure 8. The experimental data were collected from Refs. GM59 (); BP63 (); Sek76 (); Pop76 (); Koe88 ().

Using the chi-square fit with a sum of polynomials and Breit-Wigners obtained to reproduce the data in figure 8, we show in figure 9 the reaction rates for Be(n,p)Li as a function of the temperature for different values of the non-extensive parameter . The rates are given in terms of the natural logarithm of (in units of cm mol s). The solid curve corresponds to the usual Boltzmann distribution, i.e., . The dashed and dotted curves are obtained for and , respectively. In contrast to reactions induced by charged particles, we now see strong deviations from the Maxwellian rate both for and . For the deviations are larger at small temperatures and decrease as the energy increase, tending asymptotically to the Maxwellian rate at large temperatures. This behavior can be understood from figure 7 (for Be(n,p)Li the results are nearly the same as in Fig. 7). At small temperatures, e.g. , the distribution for is strongly enhanced at large energies and the tendency is that the reaction rates increase at low temperatures. This enhancement disappears as the temperature increase (lower panel of figure 7). For the reaction rate is suppressed, although not as much as for the charged-induced reactions, the reason being due a compensation by an increase because of normalization at low energies.

Having discussed the dependence of the reaction rates on the non-extensive parameter for a few standard reactions, we now consider the implications of the non-extensive statistics to the predictions of the BBN. It is clear from the results presented above that an appreciable impact on the abundances of light elements will arise.

## Iv BBN with non-extensive statistics

The BBN is sensitive to certain parameters, including the baryon-to-photon ratio, number of neutrino families, and the neutron decay lifetime. We use the values , , and s for the baryon-photon ratio, number of neutrino families, and neutron-day lifetime, respectively. Our BBN abundances were calculated with a modified version of the standard BBN code derived from Refs. WFH67 (); Kaw68 (); Kaw92 ().

Although BBN nucleosynthesis can involve reactions up to the CNO cycle Alan11 (), the most important reactions which can significantly affect the predictions of the abundances of the light elements [He, D, He, Li] are n-decay, p(n,)d, d(p,He, d(d,n)He, d(d,p)t, He(n,p)t, t(d,n)He, He(d,p)He, HeBe, tLi, Be(n,p)Li and Li(p,He. Except for these reactions, we have used the reaction rates needed for the remaining reactions from a compilation by NACRE Nacre () and that reported in Ref. Des04 (). For the 11 reactions mentioned above, we have collected data from Refs. SKM93 (); Nacre (); Des04 (), and references mentioned therein (data for n(p,,d) reaction was taken from the on-line ENDF database ENDF () - see also Cyb04 (); And06 ()), fitted the S-factors with a sum of polynomials and Breit-Wigner functions and calculated the reaction rates for Maxwellian and non-Maxwellian distributions.

### iv.1 Elemental abundances

In figure 10 we show the calculated deuterium abundance. The solid curve is the result with the standard Maxwell distributions for the reaction rates. Using the non-extensive distributions yields the dotted line for and the dashed line for . It is interesting to observe that the deuterium abundances are only moderately modified due to the use of the non-extensive statistics for . Up to temperatures of the order of , the abundance for D/H tends to agree for the extensive and non-extensive statistics. This is due to the fact that any deuterium that is formed is immediately destroyed (a situation known as the deuterium bottleneck). But, as the temperature decreases, the reaction rates for the p(n,)d reaction are considerably enhanced for (see figure 7), and perhaps more importantly, they are strongly suppressed for all other reactions involving deuterium destruction, as clearly seen in figure 4. This creates an unexpected over abundance of deuterons for the non-extensive statistics with . The deuterium, a very fragile isotope, is easily destroyed after the BBN and astrated. Its primordial abundance is determined from observations of interstellar clouds at high redshift, on the line of sight of distant quasars. These observations are scarce but allow to set an average value of Pet08 (). The predictions for the D/H ratio with the statistics (D/H = ) is about a factor 200 larger than those from the standard BBN model, clearly in disagreement with the observation.

Maxwell | Non-Max. | Non-Max. | Observation | |
---|---|---|---|---|

BBN | ||||

He/H | 0.249 | 0.243 | 0.141 | |

D/H | 2.62 | 3.31 | 570 | () |

He/H | 0.98 | 0.91 | 69.1 | () |

Li/H | 4.39 | 6.89 | 356. |

A much more stringent constraint for elemental abundances is given by He, which observations set at about Boe85 (); Ave10 (); IT10 (). The He abundance generated from our BBN calculation is plotted in figure 11. The solid curve is the result obtained with the standard Maxwell distribution for the reaction rates. Using the non-extensive distributions yields the dotted line for and the dashed line for . Again, the predicted abundances for deviate substantially from standard BBN results. This time only about half of He is produced with the use of a non-extensive statistics with . The reason for this is the suppression of the reaction rates for formation of He with through the charged particle reactions t(d,n)He, He(d,p)He.

A strong impact of using non-extensive statistics for both and values of the non-extensive parameter is seen in figure 12 for the He abundance. While for there is an overshooting in the production of He, for one finds a smaller value than the one predicted by the standard BBN. This is due to the distinct results for the destruction of He through the reaction He(n,p)t, which is enhanced for and suppressed for , in the same way as it happens for the reaction Be(n,p)Li, shown in figure 8. He is both produced and destroyed in stars and its abundance is still subject to large uncertainties, Ban02 (); Van03 ().

Non-extensive statistics for both and values also alter substantially the Li abundance, as shown in figure 13. For both values of the non-extensive parameter , and , there is an overshooting in the production of Li. The increase in the Li abundance is more accentuated for . The lithium problem is associated with a smaller value of the observed Li abundance as compared the predictions of BBN. There has been many attempts to solve this problem by testing all kinds of modifications of the parameters of the BBN or the physics behind it (a sample of this literature is found in Refs. Fie11 (); Cyb04 (); NB00 (); BC08 (); BP09 (); Boy10 (); WBB11 (); KD11 (); FP12 () and references therein). In the present case, the use of a non-Maxwellian velocity distribution seems to worsen this scenario. A recent analysis yields the observational value of Sb10 ().

### iv.2 Sensitivity study

We have calculated a window of opportunity for the non-extensive parameter with which one can reproduce the observed abundance of light elements. We chose the data for the abundances, , of He/H, D/H, He/H and Li as reference. We then applied the ordinary statistics, defined by the minimization of

(28) |

where are the abundances obtained with the non-extensive statistics with parameter , are the observed abundances, and the errors for each datum, and the sum is over all data mentioned in Table 1. From this chi-square fit we conclude that is compatible with observations.

No attempt has been made to determine which element dominates the constraint on . This might be important for a detailed study of the elemental abundance influence from nuclear physics inputs, namely, the uncertainty of the reaction cross sections. A study along these lines might be carried in a similar fashion as described in Ref. NB00 (). Weights on the reliability of observational data should also be considered for a more detailed analysis. For example, constraints arising from He/H abundance may not be considered trustworthy because of uncertain galactic chemical evolution. On the other hand, a constraint from the observation of He/D is more robust, and so on. Based on our discussion in Section II, it is more important to determine how the non-extensive statistics can modify more stringent conditions during the big bang, such as the modification of weak-decay rates and its influence on the n/p ratio which strongly affects the He abundance.

## V Conclusions

In table 1 we present results for the predictions of the BBN with Maxwellian and non-Maxwellian distributions. The predictions are compared with data from observations reported in the literature. It is evident that the results obtained with the non-extensive statistics strongly disagree with the data. The overabundance of Li compared to observation gets worse if . The three light elements D, He and Li constrain the primordial abundances rather well. For all these abundances, a non-extensive statistics with leads to a greater discrepancy with the experimental data.

Except for the case of He the use of non-extensive statistics with does not rule out its validity when the non-Maxwellian BBN results are compared to observations. He is at present only accessible in our Galaxy’s interstellar medium. This means that it cannot be measured at low metallicity, a requirement to make a fair comparison to the primordial generation of light elements. This also means that the primordial He abundance cannot be determined reliably. The result presented for the He abundance in table 1 is quoted from Ref. Ban02 (). Notice that our analysis does not include the changes that the non-extensive statistics would bring to the n/p conversion rates. The electron distributions would also be expected to change accordingly. This would change the freeze out temperature and a corresponding influence on the He abundance.

We conclude that it does not seem possible to change the Maxwell-Boltzmann statistics to reproduce the observed abundance of light elements in the universe without destroying many other successful predictions of big bang nucleosynthesis. A chi-square fit of our calculations with the observations of elemental abundance concludes that the non-extensive parameter is constrained to . This means that, should a non-Maxwellian distribution due to the use of the Tsallis non-extensive statistics be confirmed (with a sizable deviation from ), our understanding of the cosmic evolution of the universe would have to be significantly changed.

We would like to acknowledge beneficial discussions with Alain Coc and Richard Cyburt. This work was partially supported by the US-DOE grants DE-FG02-08ER41533 and DE-FG02-10ER41706, and by the Brazilian agencies, CNPq and FAPESP. One of the authors (C.B.) acknowledges the Helmholtz International Center for FAIR (HIC for FAIR) for supporting his visit to the GSI Helmholtzzentrum für Schwerionenforschung, where much of this work was done.

## References

- (1) P. Noterdaeme, P. Petitjean, R. Srianand, C. Ledoux, and S. Lopez, Astron. and Astrophys. 526, L7 (2011).
- (2) Keith A. Olive, Lectures at the First NCTS Workshop on Astroparticle Physics, December 6 - 9, 2001, Kenting, Taiwan, arXiv:astro-ph/0202486 (2002).
- (3) LEP Collaboration, Phys. Rep. 427, 257 (2006).
- (4) A. Serebrov et al., Phys. Lett. B 605, 72 (2005).
- (5) Fred E. Wietfeldt and Geoffrey L. Greene, Rev. Mod. Phys. 83, 1173 (2011).
- (6) G. J. Mathews, T. Kajino and T. Shima, Phys. Rev. D 71, 021302(R) (2005).
- (7) E. Komatsu, et al., ApJS 192, 18 (2011).
- (8) Brian D. Fields, Ann. Rev. Nucl. Part. Sci. 61, 47 (2011).
- (9) C. Tsallis, J. Stat. Phys. 52, 479 (1988).
- (10) A. Rényi, Proc. of the 4th Berkeley Symposium on Mathematics, Statistics and Probability, pp. 547Ð561 (1960).
- (11) M. Gell-Mann and C. Tsallis, eds., “Nonextensive Entropy-Interdisciplinary Applications, (Oxford University Press, New York, 2004). For a regularly updated bibliography see URL: http://tsallis.cat.cbpf.br/biblio.htm
- (12) C. Tsallis, M. Gell-Mann and Y. Sato, Proc. Natl. Acad. Sc. USA 102, 15377 (2005).
- (13) “Introduction to Nonextensive Statistical Mechanics”, C. Tsallis, Springer Verlag (2009), ISBN 978-0-387-85358-1.
- (14) K. Sau Fa and E.K. Lenzi, J. of Math. Phys. 42, 1148 (2001).
- (15) J.A.S. Lima, R. Silva, and J. Santos, Astron. and Astrophys., 396 309 (2002).
- (16) A. Taruya and M. Sakagami, Phys. Rev. Lett. 90, 181101 (2003).
- (17) A. Taruya and M. Sakagami, Cont.Mech. and Therm. 16 279 (2004).
- (18) P.H. Chavanis and C. Sire, Physica A 356 (2005) 419.
- (19) M.E. Pessah, D.F. Torres and H. Vucetich, Physica A 297, 164 (2001).
- (20) A. Bernui, C. Tsallis and T. Villela, Europhys. Lett. 78, 19001 ( 2007).
- (21) M. Lissia and P. Quarati, Europhysics News 36, 211 (2005).
- (22) J. Haubold and D. Kumar, Astroparticle Physics, 29 ,70 (2008).
- (23) S. Degl’Innocenti et al, Phys. Lett. B 441, 291 (1998).
- (24) M. Coraddu et al, Braz. J. Phys. 29, 153 (1999).
- (25) R. Silva Jr., A.R. Plastino, and J. Lima, Phys. Lett. A249, 401 (1998).
- (26) J.A. S. Lima, R. Silva Jr., and J. Santos, Phys. Rev. E61, 3260 (2000).
- (27) V. Muñoz, Nonlin. Processes Geophys. 13, 237 (2006).
- (28) W.A. Fowler, G.R. Caughlan, and B.A. Zimmerman, Ann. Rev. Astron. Astrophys. 5, 525 (1967).
- (29) G. Kaniadakis et al, Phys. Rev. E 71, 046128 (2005).
- (30) F. Buyukkilic, D. Demirhan, Eur. Phys. J B 14, 705 (2000).
- (31) Biao Wang, C.A. Bertulani and A.B. Balantekin, Phys. Rev. C 83, 018801 (2011).
- (32) F. Buyukkili, I. Sokmen and D. Demirhan, Chaos, Solitons and Fractals 13, 749 (2002).
- (33) H. Aragao-Rego et al, Physica A 317, 199 (2003).
- (34) M. Nauenberg, Phys. Rev. E 67, 036114 (2003).
- (35) C. Tsallis, Phys. Rev. E 69, 038101 (2004).
- (36) H. Chamati, A.Ts. Djankova, N.S. Tonchev, arXiv:cond-mat/0311234.
- (37) Damian H. Zanette, Marcelo M. Montemurro,Phys. Lett. A 316, 184 (2003).
- (38) Damian H. Zanette, Marcelo M. Montemurro, Phys. Lett. A 324, 383 (2004).
- (39) F. Bouchet, T. Dauxois, S. Ruffo, Europhysics News 37, 9 (2006).
- (40) T. Dauxois, Journal of Statistical Mechanics (2007) N08001.
- (41) H. Touchette, Phys. Lett A 377, 436 (2013).
- (42) R.L. Schulte, M. Cosack, A.W. Obst, J.L. Weil, Nucl. Phys. A192, 609 (1972).
- (43) A. Krauss, H.W. Becker, H.P. Trautvetter, C. Rolfs, K. Brand, Nucl. Phys. A465, 150 (1987). [
- (44) R.E. Brown, N. Jarmie, Phys. Rev. C41, 1391 (1990).
- (45) H.-S. Bosch, G.M. Hale, Nucl. Fusion 32, 611 (1992).
- (46) U. Greife, F. Gorris, M. Junker, C. Rolfs, D. Zahnow, Z. Phys. A351, 107 (1995).
- (47) J. M. Freeman, R. C. Hanna and J. H. Montague, Nucl. Phys. 5 (1958) 148.
- (48) Y. Cassagnou, J. M. Jeronymo,G. S. Mani, A. Sadeghi and P. D. Forsyth, Nucl. Phys. 33 (1962) 449.
- (49) Y. Cassagnou, J. M. Jeronymo,G. S. Mani, A. Sadeghi and P. D. Forsyth, Nucl. Phys. 41 (1963) 176
- (50) G. S. Mani, R. Freeman, F. Picard, A. Sadeghi and D. Redon, Nucl. Phys. 60 (1964) 588.
- (51) O. Fiedler and P. Kunze, Nucl. Phys. A96 (1967) 513.
- (52) G. M. Lerner and J. B. Marion, Nucl. Inst. Meth. 69 (1969) 115.
- (53) H. Spinka, T. Tombrello and H. Winkler, Nucl. Phys. A164 (1971) 1.
- (54) C. Rolfs and R. W. Kavanagh, Nucl. Phys. A455 (1986) 179.
- (55) J. F. Harmon, Phys. Rev. C8 (1973) 106.
- (56) S. Engstler, G. Raimann, C. Angulo, U. Greife, C. Rolfs, U. Schröder, E. Somorjai, B. Kirch and K. Langanke, Z. Phys. A342 (1992) 471.
- (57) S. Engstler, G. Raimann, C. Angulo, U. Greife, C. Rolfs, U. Schröder, E. Somorjai, B. Kirch and K. Langanke, Phys. Lett. B279 (1992) 20.
- (58) J.H. Gibbons, R.L. Macklin, Phys. Rev. 114 (1959) 571.
- (59) R.R. Borchers, C.H. Poppe, Phys. Rev. 129 (1963) 2679.
- (60) K.K. Sekharan, H. Laumer, B.D. Kern, F. Gabbard, Nucl. Instrum. Methods 133 (1976) 253.
- (61) C.H. Poppe, J.D. Anderson, J.C. Davis, S.M. Grimes, C. Wong, Phys. Rev. C 14 (1976) 438.
- (62) P.E. Koehler, C.D. Bowman, F.J. Steinkruger, D.C. Moody, G.H. Hale, J.W. Starner, S.A. Wender, R.C. Haight, P.W. Lisowski, W.L. Talbert, Phys. Rev. C 37 (1988) 917.
- (63) Robert Wagoner, William A. Fowler, and F. Hoyle, Astrophys. J. 148, 3 (1967).
- (64) L. Kawano (1988), FERMILAB Report No. PUB-88/34-A (unpublished).
- (65) L. Kawano, NASA STI/Recon Technical Report N 92, 25163 (1992).
- (66) Alain Coc, Stephane Goriely, Yi Xu, Matthias Saimpert, and Elisabeth Vangioni, arXiv:1107.1117.
- (67) C. Angulo et al., Nucl. Phys. A 656, 3 (1999).
- (68) P. Descouvemont, A. Adahchour, C. Angulo, A. Coc,, and E. Vangioni-Flam, At. Data Nucl. Data Tables 88, 203 (2004).
- (69) M.S. Smith, L.H. Kawano, R.A. Malaney, Astrophys. J. Suppl. 85 (1993) 219.
- (70) ENDF/B-VI, Online database at the NNDC Online Data Service, http://www.nndc.bnl.gov.
- (71) Richard Cyburt, Phys. Rev. D 70, 023505 (2004).
- (72) S. Ando et al., Phys. Rev. C 74, 025809 (2006).
- (73) M. Pettini, B.J. Zych, M. Murphy, A. Lewis, and C.C. Steidel, MNRAS 391, 1499 (2008).
- (74) A.M. Boesgaard, and G. Steigman, Ann. Rev. Astron. Astrophys. 23, 319 (1985).
- (75) E. Aver, K.A. Olive, and E.D. Skillman, J. Cosmol. Astropart. Phys. 5, 3 (2010).
- (76) Y.I. Izotov, and T.X. Thuan, ApJ, 710, L67 (2010).
- (77) T. Bania, R. Rood, and D. Balser, Nature 415, 54 (2002).
- (78) E. Vangioni-Flam, K.A. Olive, B.D. Fields, and M. Cass, ApJ, 585, 611 (2003).
- (79) Kenneth M. Nollett and Scott Burles, Phys. Rev. D 61, 123505 (2000).
- (80) R. H. Cyburt, B.D. Fields and K. Olive, J. Cos. Astropart. Phys. 11, 012 (2008).
- (81) R. H. Cyburt and M. Pospelov, arXiv:0906.4373 [astro-ph] (2009).
- (82) Richard N. Boyd, Carl R. Brune, George M. Fuller, and Christel J. Smith, Phys. Rev. D 82, 105005 (2010).
- (83) O. S. Kirsebom and B. Davids, Phys. Rev. C84, 058801 (2011).
- (84) Fabio Iocco and Miguel Pato, Phys. Rev. Lett. 109, 021102 (2012).
- (85) L. Sbordone, et al., Astron. Astrop. 522, 26 (2010).