Abstract

Possible modification in the velocity distribution in the nonresonant reaction rates leads to an extended reaction rate probability integral. The closed form representation for these thermonuclear functions is used to obtain the stellar luminosity and neutrino emission rates. The composite parameter that determines the standard nuclear reaction rate through the Maxwell-Boltzmann energy distribution is extended to by the extended reaction rates through a more general distribution than the Maxwell-Boltzmann distribution. The new distribution is obtained by the pathway model introduced by Mathai (2005). Simple analytic models considered by various authors are utilized for evaluating stellar luminosity and neutrino emission rates and are obtained in generalized special functions such as Meijer's G-function and Fox's H-function. The standard and extended nonresonant thermonuclear functions are compared by plotting them. Behaviour of the new energy distribution, which is more general than the Maxwell-Boltzmann, is also studied.

1. Introduction

The mystery behind the distant universe is explored so far by the understanding of the sun, the star near to us. It is the only star whose mass, radius, and luminosity are fairly accurately known. The structural change in the sun is due to the central thermonuclear reactor in it. Solar nuclear energy generation and solar neutrino emission are governed by chains of nuclear reactions in the gravitationally stabilized solar fusion reactor [1, 2]. Qualitative calculations of specific reaction rates require a large amount of experimental inputs and theoretical assumptions. By using the theories from nuclear physics and kinetic theory of gases one can determine the reaction rate for low-energy nonresonant thermonuclear reactions in nondegenerate plasma [3]. The formalization of the calculation of the reaction rate of interacting articles under cosmological or stellar conditions was presented by many authors [4, 5]. For the most common case, a nuclear reaction in which a particle of type strikes a particle of type producing a nucleus and a new particle is symbolically represented as where is the energy release given by , where denote the masses of the particles and denotes the velocity of light. The reaction rate of the interacting particles and is obtained by averaging the reaction cross section over the normalized density function of the relative velocity of the particles [57]. Let and denote the number densities of the particles and , respectively, and let be the reaction cross section where is the relative velocity of the particles and is the normalized velocity density; then the reaction rate is given by where is the Kronecker delta which is introduced to avoid double counting in the reaction if particles and are identical. is the thermally averaged product which is in fact the probability per unit time that two particles and confined to a unit volume will react with each other. is the reduced mass of the particles given by . is the kinetic energy of the particles in the centre of mass system. From literature [4, 5, 7] it may be noted that all the analytic expressions for astrophysically relevant nuclear reaction rates underline the hypothesis that the distribution of the relative velocities of the reacting particles always remains Maxwell-Boltzmann for a nonrelativistic nondegenerate plasma of nuclei in thermodynamic equilibrium. The Maxwell-Boltzmann relative kinetic energy distribution can be written as where is the Boltzmann constant and is the temperature. Substituting (3) in (2) we get The thermonuclear fusion depends on three physical variables, the temperature , the Gamow energy , and the nuclear fusion factor . If two nuclei of charges and collide at low energies below the Coulomb barrier, then the Gamow energy is given by [8, 9] where is the electromagnetic fine structure constant given by where is the quantum of electric charge, is Planck’s quantum of action, and is approximately [9] for our universe. Thus the Gamow factor, which is determined by the electromagnetic force, and the nuclear fusion factor set the nuclear reaction cross section at low energies for nonresonant charged particles as [7, 10] and is the cross section factor which is often found to be constant or a slowly varying function of energy over a limited range of energy given by [4, 5] Substituting (7) and (8) in (4) we obtain This is the nonresonant reaction rate probability integral in the Maxwell-Boltzmann case. The closed form evaluation of this integral can be seen in a series of papers by Mathai and Haubold; see, for example, Haubold and Mathai [6, 11], Mathai and Haubold [5], and so forth. The main aim of the present work is to extend the reaction rate probability integral given in (9) by replacing the Maxwell-Boltzmann energy distribution by a more general energy distribution called the pathway energy distribution obtained by using the pathway model of Mathai introduced in 2005.

The paper is organized as follows. In the next section we discuss a more general energy distribution than the Maxwell-Boltzmann distribution and obtain the extended reaction rate probability integral in the nonresonant case. We take advantage of the closed form representation of the extended thermonuclear reaction rate for finding the luminosity and the neutrino emission rate of the nonlinear stellar model under consideration in Section 3. Section 4 is devoted to finding the desired connection between stellar structure parameters and the neutrino emission of the stellar model by using the closed form analytic representation of the extended reaction rates. A comparison of the Maxwell-Boltzmann energy distribution with the pathway energy distribution is done with the help of graphs in Section 5. Also we try to discriminate the standard and extended reaction rates. Concluding remarks are included in Section 6.

2. Extended Nonresonant Thermonuclear Reaction Rate and Its Closed Forms

In recent years, possible deviations of the velocity distribution of the plasma particles from the Maxwell-Boltzmann in connection with the production of neutrinos in the gravitationally stabilized solar fusion reactor have been pointed out [2, 10, 1215]. It was initiated by Tsallis, the originator of nonextensive statistical mechanics [1618], who has used -exponential function as the fundamental distribution instead of the Maxwell-Boltzmann distribution. An initial attempt to extend the standard theories of reaction rates to Tsallis statistics was done by many authors; see Mathai and Haubold [19] and Saxena et al. [20]. In 2005, Mathai introduced the pathway model by which even more general distributions can be incorporated in the theory of reaction rates [19, 21]. Initially, pathway model was introduced for the matrix variate case to cover many of the matrix variate statistical densities. The scalar case is a particular one there. Later, Mathai, his coworkers, and others found connection of pathway model with the information theory, the fractional calculus, the Mittag-Leffler functions, and so forth. The pathway model can be effectively used in any situation in which we need to switch between three different functional forms, namely, generalized type-1 beta form, generalized type-2 beta form, and generalized gamma form, using the pathway parameter . In practical purpose of fitting experimental data, pathway model can be utilized to switch between different parametric families with thicker or thinner tail. The pathway model for the real scalar case can be explained as follows: is the generalized type-1 beta form of the pathway model. This is a model with right tail cut-off for . The Tsallis statistics for can be obtained from this model by putting [1618]. Other cases available are the regular type-1 beta density, the Pareto density, the power function, and the triangular and related models [22]. The generalized type-2 beta form of the pathway model is given by Here also for we get the Tsallis statistics for [1618]. Other standard distributions coming from this model are regular type-2 beta density, -distribution, the Lévy model, and related models [22]. When , and will reduce to the generalized gamma form of the pathway model given by This model covers generalized gamma, gamma, exponential, chi-square, the Weibull, the Maxwell-Boltzmann, the Rayleigh, and related densities. , , and defined in (10), (11), and (12), respectively, are the normalizing constants if we consider statistical densities.

By a suitable modification of the Maxwell-Boltzmann distribution given in (3) through the pathway model, we get a more general energy distribution called the pathway energy distribution given by the density for , . The Maxwell-Boltzmann energy distribution can be retrieved from (13) by taking . Thus the reaction rate probability integral given in (4) can be modified by using (13) and we get the extended reaction rate as for , . Substituting and we obtain the above integral in a more convenient form as follows: for , . Here we consider the integral to be evaluated as The integral can be evaluated by the techniques in applied analysis and can be obtained in closed form via Meijer’s -function as [2, 23, 24] which yields the nonresonant reaction rate probability integral in the extended case as Meijer’s -function and its properties can be seen in Mathai and Saxena [25] and Mathai [26]. We can obtain series expansions of the -function given in (18) by combining the theories of residue calculus and generalized special functions; see Kumar and Haubold [24] for series expansions for all possible values of . In many cases the nuclear factor is approximately constant across the fusion window. Taking for and and taking , we obtain the extended reaction rate probability integral as The series representation for (19) can be obtained as where See the Appendix for detailed evaluation. As in (19), then, by using Stirling’s formula for gamma functions given by we get the reaction rate probability integral in the Maxwell-Boltzmann case as which is obtained in a series of papers by Mathai and Haubold; see, for example, Mathai and Haubold [5]. The integral in (23) is dominated by the minimum value of (say). The minimum value of the function , say , can be determined as and the function where . Now by using the Laplace method [27, 28] we can obtain an approximate value for (23) as In the next section we will obtain the mass, pressure, and temperature for the case of analytic stellar models characterized by density distribution and corresponding temperature distribution suggested by Haubold and Mathai [11].

3. Closed Forms of the Integral over the Stellar Nuclear Energy Generation Rate

Let us consider the density distribution considered by Haubold and Mathai [6, 11] and Mathai and Haubold [5] in the form where is the central density of the star, is an arbitrary distance from the center, and is the solar radius. This density function is capable of producing different density distributions by choosing the free parameter . Now we determine the quantities , and , the mass, the pressure, and the temperature at .

By the equation of the mass conservation we get From (30), we get the central density as If an element of a matter at a distance from the center of a spherical system is in hydrostatic equilibrium, then setting the sum of the radial forces acting on it to zero we obtain where is the gravitational constant. Assuming that the pressure at the center of the sun is and at the surface is zero, we get Using the boundary conditions , we get , where By the kinetic theory of gases, for a perfect gas, the pressure is given by For the temperature of interest for stellar models, we neglect the negligible radiation pressure from the total pressure and obtain from (35) the following: The central temperature , where is as defined in (34).

Thus we have obtained the mass, the pressure, and the temperature throughout the nonlinear stellar model with the density distribution defined in (28). Next our aim is to obtain analytical results for stellar luminosity and neutrino emission rates for various stellar models.

4. Stellar Luminosity and Neutrino Emission Rate

The energy conservation equation states that the net increase in the rate of energy flux coming out of a spherical shell from the inside is the same as the energy produced within the shell [29]. If we denote as the energy flux through the sphere of radius , then we have where is the energy produced per second by nuclear reactions in each gram of stellar matter. depends on the chemical composition in each gram of stellar matter. Here usually is a constant but will be equal to at the surface of the star. We assume here that the star is chemically homogeneous (that is a star where chemical composition throughout is a constant). Also we assume the energy generation rate for one particular nuclear reaction. Now if we denote as the extended nonresonant thermonuclear reaction rate for the particles and defined by (19), then we will consider the energy generation rate and it can be written in terms of the extended reaction rates via where in which is the amount of energy given off in a single reaction. It is to be noted that, by using the asymptotic behaviour of [26] and as , , the composite parameter considered by [9], which is defined as for our universe for proton-proton fusion under typical stellar conditions [9]. Then from (37) we have the total luminosity of the star by integration as follows: If we are considering only one specific reaction defined as in (1), then we have where the energy generation rate is defined in (38) and is a suitable density distribution explaining the sun. Writing (42) in terms of we get The number density of a particle , for a gas of mean density , can be expressed as where stands for Avagadro’s constant, is the atomic mass of particle in atomic mass units, and is the mass fraction of particle such that . Substituting from (19) and using (44) we have If we divide the “internal luminosity” by the amount of energy , then we get the total number of particles per second liberated in the reaction given by (1) as follows: For the density distribution defined in (28) introduced by Haubold and Mathai [5, 6] and the corresponding temperature distribution (36), we get where is as defined in (34). If we put a substitution , then we get Putting and simplifying we obtain where is defined as As , , and , . If we take we have . The minimum value of is at and the value is . Thus the minimum value is nonnegative since decreases steadily from to for all . Therefore . Since , for all , for all . Thus . Hence for and for . Thus by using the binomial expansion we obtain where is the Pochhammer symbol defined for by whenever exists. Taking we get and Then from (49) we obtain By using beta integral and using (53) we obtain Now from (46) we obtain the total number of particles per second liberated in the reaction (1) as follows:

For more details on -function and its properties see [25, 26]. Thus we have obtained the the total number of particles per second liberated in the reaction (1) in terms of the density distribution considered by Haubold and Mathai [6, 29].

5. Comparison of Pathway Energy Density and the Maxwell-Boltzmann Energy Density

In Figure 1 it can be obtserved that, for the nuclei to react at energy , they have to borrow an energy from the thermal environment. The probability of such an energy is proportional to the Maxwell-Boltzmann energy . The fusion will take place when the nuclei penerate the Coulomb barrier keeping them apart. The probability of penetration is given by the factor . The product of these two factors illustrates that fusion mostly occurs in the energy window given in the figure.

In Figure 2 the pathway energy density is plotted for , respectively. For different values of we get different energy densities (curves (a), (b), (c), (d), and (e)). The nonresonant cross section is also plotted. The product of the pathway energy density and the nonresonant cross section for different values of , namely, , is also plotted. It is to be noted that as the pathway energy density coincides with the Maxwell-Boltzmann energy density and also the fusion window for the Maxwell-Boltzmann case in Figure 1.

The curves in the figure represent pathway density for (a) , (b) , (c) , (d) , and (e) . (f) represents the nonresonant cross section . The product is also represented for (g) , (h) , (i) , (j) , and (k) .

Figure 3(a) shows the pathway energy density defined in (13) for and Figure 3(b) shows the Maxwell-Boltzmann energy density defined in (3). As the pathway energy density reduces to the Maxwell-Boltzmann energy density. The pathway energy density covers many stable and unstable situations as the value of varies. If the Maxwell-Boltzmann density is the equilibrium situation, many other nonequilibrium situations are covered by the pathway energy density.

6. Concluding Remarks

In this paper we have modified the energy distribution for a nonresonant reaction rate probability integral. The composition parameter considered by [9] is extended to by the pathway energy density. Considering the analytic density distributions developed by Haubold and Mathai [6, 29], they are used to obtain the stellar luminosity and the neutrino emission rates and are obtained in generalized special functions such as Meijer’s -function. The pathway energy density considered here covers many density functions and hence the extended reaction rate integral covers a wider class of integral. Pathway energy density helps us to obtain various fusion windows by giving different values to , the pathway parameter, which in turn leads to a new opening in the fusion research. The graphs plotted here are by using Maple 14 in Windows XP platform.

Appendix

Series Representation

The series representation for the right-hand side of (19) can be obtained through the following procedure. Here we apply residue calculus on the -function given in (19). Consider the -function as follows: The right-hand side is the sum of the residues of the integrand. The poles of the gammas in the integral representation in (19) are as follows:poles of ;poles of ;poles of .Here the poles of and will coincide with each other at all points except at . Note that the pole is a pole of order 1, are each of order 1, and are each of order 2. We know that when is defined, , The sum of the residues corresponding to the poles is given by The sum of the residues corresponding to the poles is where is the hypergeometric function defined by

To obtain the sum of the residues corresponding to poles of order 2, we proceed as follows: where We have where is a Psi function or digamma function (see Mathai [26]) and is Euler’s constant. Now Then by using (A.7), (A.9), and (A.10) we get where Thus from (A.4), (A.5), and (A.11) we get (20).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the Department of Science and Technology, Government of India, New Delhi, for the financial assistance for this work under Project no. SR/S4/MS:287/05 and the Centre for Mathematical Sciences for providing all facilities.