Research Article  Open Access
Riou Nakamura, Masaaki Hashimoto, Shinichiro Fujimoto, Katsuhiko Sato, "Constraint on Heavy Element Production in Inhomogeneous BigBang Nucleosynthesis from the Light Element Observations", Journal of Astrophysics, vol. 2013, Article ID 587294, 9 pages, 2013. https://doi.org/10.1155/2013/587294
Constraint on Heavy Element Production in Inhomogeneous BigBang Nucleosynthesis from the Light Element Observations
Abstract
We investigate the observational constraints on the inhomogeneous bigbang nucleosynthesis that Matsuura et al. (2005) suggested that states the possibility of the heavy element production beyond ^{7}Li in the early universe. From the observational constraints on light elements of ^{4}He and D, possible regions are found on the plane of the volume fraction of the highdensity region against the ratio between high and lowdensity regions. In these allowed regions, we have confirmed that the heavy elements beyond Ni can be produced appreciably, where p and/or rprocess elements are produced well simultaneously.
1. Introduction
Bigbang nucleosynthesis (BBN) has been investigated to explain the origin of the light elements, such as , D, , and , during the first few minutes [1â€“4]. Standard model of BBN (SBBN) can succeed in explaining the observation of those elements, [5â€“9], D [10â€“13], and [14, 15], except for . The study of SBBN has been done under the assumption of the homogeneous universe, where the model has only one parameter, the baryontophoton ratio . If the present value of is determined, SBBN can be calculated from the thermodynamical history with the use of the nuclear reaction network. We can obtain the reasonable value of by comparing the calculated abundances with observations. In the meanwhile, the value of is obtained as [1] from the observations of and D. These values agree well with the observation of the cosmic microwave background: [16].
On the other hand, BBN with the inhomogeneous baryon distribution also has been investigated. The model is called as inhomogeneous BBN (IBBN). IBBN relies on the inhomogeneity of baryon concentrations that could be induced by baryogenesis (e.g., [17]) or phase transitions such as QCD or electroweak phase transition [18â€“21] during the expansion of the universe. Although a largescale inhomogeneity is inhibited by many observations [16, 22â€“24], a small scale one has been advocated within the present accuracy of the observations. Therefore, it remains a possibility for IBBN to occur in some degree during the early era. In IBBN, the heavy element nucleosynthesis beyond the mass number has been proposed [17, 18, 25â€“35]. In addition, peculiar observations of abundances for heavy elements and/or could be understood in the way of IBBN. For example, the quasar metallicity of C, N, and Si could have been explained from IBBN [36]. Furthermore, from recent observations of globular clusters, a possibility of inhomogeneous helium distribution is pointed out [37], where some separate groups of different main sequences in blue band of low mass stars are assumed due to high primordial helium abundances compared to the standard value [38, 39]. Although baryogenesis could be the origin of the inhomogeneity, the mechanism of it has not been clarified due to unknown properties of the supersymmetric Grand Unified Theory [40].
Despite a negative opinion against IBBN due to insufficient consideration of the scale of inhomogeneity [41], Matsuura et al. have found that the heavy element synthesis for both  and processes is possible if [42], where they have also shown that the high regions are compatible with the observations of the light elements, and D [43]. However, their analysis is only limited to a parameter of a specific baryon number concentration. In this paper, we extend the investigations of Matsuura et al. [42, 43] to check the validity of their conclusion from a wide parameter space of the IBBN model.
In Section 2, we review and give the adopted model of IBBN which is the same one as that of Matsuura et al. [43]. Constraints on the critical parameters of IBBN due to light element observations are shown in Section 3, and the possible heavy elements of nucleosynthesis are presented in Section 4. Finally, Section 5 is devoted to the summary and discussion.
2. Model
In this section, we introduce the model of IBBN. We adopt the twozone model for the inhomogeneous BBN. In the IBBN model, we assume the existence of spherical highdensity region inside the horizon. For simplicity, we ignore in the present study the diffusion effects before and during the primordial nucleosynthesis , because the timescale of the neutron diffusion is longer than that of the cosmic expansion [25, 37].
To find the parameters compatible with the observations, we consider the average abundances between the high and lowdensity regions. We get at least the parameters for the extreme case by averaging the abundances in two regions. Let us define the notations, , , and as ,â€‰â€‰, and low baryon number densities. is the volume fraction of the high baryon density region. , , and are mass fractions of each element in average, and high and lowdensity regions, respectively. Then, basic relations are written as follows [43]: Here, we assume the baryon fluctuation to be isothermal as was done in previous studies (e.g., [18, 19, 30]). Under that assumption, since the baryontophoton ratio is defined by the number density of photon in standard BBN, (1) is rewritten as follows: where s with subscripts are the baryontophoton ratios in each region. In the present paper, we fix from the cosmic microwave background observation [16]. The values of and are obtained from both and the density ratio between high and lowdensity regions: .
To calculate the evolution of the universe, we solve the following Friedmann equation: where is the cosmic scale factor and is the gravitational constant. The total energy density in (4) is the sum of decomposed parts: Here, the subscripts , , and indicate photons, neutrino, and electrons/positrons, respectively. The final term is the baryon density obtained as .
We should note the energy density of baryon. To get the time evolution of the baryon density in both regions, the energy conservation law is used as follows: where is the pressure of the fluid. When we solve (6), initial values in both regions are obtained from (2) with and fixed. For , the baryon density in the highdensity region, , is larger than the radiation component at â€‰K. However, we note that the contribution to (5) is not , but . In our research, the ratio of to is about at BBN epoch. Therefore, we can neglect the final term of (5) in the same way as it has been done in SBBN during the calculation of (4).
3. Constraints from Light Element Observations
In this section, we calculate the nucleosynthesis in high and lowdensity regions with the use of the BBN code [44] which includes 24 nuclei from neutron to . We adopt the reaction rates of Descouvemont et al. [45], the neutron lifetime â€‰sec [1], and consider three massless neutrinos.
Let us consider the range of . For , the heavier elements can be synthesized in the highdensity regions as discussed in [33]. For , contribution of the lowdensity region to can be neglected, and therefore to be consistent with the observations of light elements, we need to impose the condition of .
Figure 1 illustrates the light element synthesis in the high and lowdensity regions with and = that corresponds to and . Light elements synthesized in these calculations are shown in Table 1. In the lowdensity region, the evolution of the elements is almost the same as the case of SBBN. In the highdensity region, while is more abundant than that in the lowdensity region, (or ) is much less produced. In this case, we can see that average values such as and D are overproduced as shown in Table 1. However, this overproduction can be saved by choosing the parameters carefully. We need to find the reasonable parameter ranges for both and by comparing with the observation of the light elements.

(a)
(b)
Now, we put constraints on and by comparing the average values of and D obtained from (3) with the following observational values. First, we consider the primordial abundance reported in [8]: and [9]: We adopt abundances as follows:
Next, we take the primordial abundance from the D/H observation reported in [12]: and [13]: Considering those observations with errors, we adopt the primordial D/H abundance as follows:
Figure 2 illustrates the constraints on the plane from the above light element observations with contours of constant . The solid and dashed lines indicate the upper limits from (9) and (12), respectively. From the results, we can obtain approximately the following relations between and : The observation (9) gives the upper bound for , and the limit for is obtained from D observation (12). As shown in Figure 2, we can find the allowed regions which include the very highdensity region such as .
We should note that takes a larger value, nuclei which are heavier than are synthesized more and more. Then we can estimate the amount of total CNO elements in the allowed region. Figure 3 illustrates the contours of the summation of the average values of the heavier nuclei (), which correspond to Figure 2 and are drawn using the constraint from and D/H observations. As a consequence, we get the upper limit of total mass fractions for heavier nuclei as follows:
4. Heavy Element Production
In the previous section, we have obtained the amount of CNO elements produced in the twozone IBBN model. However, it is not enough to examine the nuclear production beyond because the baryon density in the highdensity region becomes so high that elements beyond CNO isotopes can be produced [17, 31, 32, 34, 42]. In this section, we investigate the heavy element nucleosynthesis in the highdensity region considering the constraints shown in Figure 2. Abundance change is calculated with a large nuclear reaction network, which includes 4463 nuclei from neutron and proton to Americium (Z = 95 and A = 292). Nuclear data, such as reaction rates, nuclear masses, and partition functions, are the same as the ones used in [46â€“49] except for the neutronproton interaction. We use the weak interaction of Kawano code [50], which is adequate for the hightemperature epoch of â€‰K.
As seen in Figure 3, heavy elements of are produced nearly along the upper limit of . Therefore, to examine the efficiency of the heavy element production, we select five models with the following parameters: = , , , , and corresponded to = , , , , and . Adopted parameters are indicated by filled squares in Figure 2.
First, we evaluate the validity of the nucleosynthesis code with nuclei. Table 2 shows the results of the light elements, , D, , , and . The results of the highdensity region are calculated by the extended nucleosynthesis code, and the abundances in the lowdensity region are obtained by BBN code. The average abundances are obtained by (3). Since the average values of and D are consistent with the observations, there is no difference between BBN code and the extended nucleosynthesis code in regard to the average abundances of light elements.
(a) For cases of and  
 
(b) For cases of and  

Figure 4 shows the results of nucleosynthesis in the highdensity regions with and . In Figure 4(a), we see the time evolution of the abundances of Gd and Eu for the mass number 159. First, (stable element) is synthesized and later and are synthesized through the neutron captures. After sec, decays to nuclei by way of â†’, where the halflifes of and are â€‰min and â€‰h, respectively.
(a)
(b)
For , the result is seen in Figure 4(b). , which is a protonrich nuclei is synthesized. After that, stable nuclei is synthesized by way of , where the halflifes of and are â€‰min and â€‰min, respectively. These results are qualitatively the same as Matsuura et al. [42].
In addition, we notice the production of radioactive nuclei of and , where is produced at early times, just after the formation of . Usually, nuclei such as and are produced in supernova explosions, which are assumed to be the events after the first star formation (e.g., [51]). In IBBN model, however, this production can be found to occur at an extremely highdensity region of as the primary elements without supernova events in the early universe.
Final results ( = â€‰K) of nucleosynthesis calculations are shown in Table 3. When we calculate the average values, we set the abundances of to be zero for lowdensity side. For , a lot of nuclei of are synthesized whose amounts are comparable to that of . Produced elements in this case include both element (i.e., ) and elements (for instance, and ). For , there are few elements while both elements (i.e., and ) and elements (i.e., and ) are synthesized such as the case of supernova explosions. For , the heavy elements are produced slightly more than the total mass fraction (shown in Figure 3) derived from the BBN code calculations. This is because our BBN code used in Section 3 includes the elements up to and the actual abundance flow proceeds to much heavier elements.
