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 1 strikes a particle of type 2 producing a nucleus 3 and a new particle 4 is symbolically represented as
(1)1+2⟶3+4+E12,
where E12 is the energy release given by E12=(m1+m2-m3-m4)c2, where mi,i=1,2,3,4 denote the masses of the particles and c denotes the velocity of light. The reaction rate r12 of the interacting particles 1 and 2 is obtained by averaging the reaction cross section over the normalized density function of the relative velocity of the particles [5–7]. Let n1 and n2 denote the number densities of the particles 1 and 2, respectively, and let σ(v) be the reaction cross section where v is the relative velocity of the particles and f(v) is the normalized velocity density; then the reaction rate r12 is given by
(2)r12=(1-12δ12)n1n2〈σv〉12=(1-12δ12)n1n2∫0∞vσ(v)f(v)dv=(1-12δ12)n1n2∫0∞σ(E)(2Eμ)1/2f(E)dE,
where δ12 is the Kronecker delta which is introduced to avoid double counting in the reaction if particles 1 and 2 are identical. 〈σv〉12 is the thermally averaged product which is in fact the probability per unit time that two particles 1 and 2 confined to a unit volume will react with each other. μ is the reduced mass of the particles given by μ=(m1m2)/(m1+m2). E=μv2/2 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
(3)fMBD(E)dE=2π(1πkT)3/2exp(-EkT)EdE,
where k is the Boltzmann constant and T is the temperature. Substituting (3) in (2) we get
(4)r12=(1-12δ12)n1n2(8πμ)1/2(1kT)3/2×∫0∞Eσ(E)exp(-EkT)dE.
The thermonuclear fusion depends on three physical variables, the temperature T, the Gamow energy EG, and the nuclear fusion factor S(E). If two nuclei of charges Z1e and Z2e collide at low energies below the Coulomb barrier, then the Gamow energy EG is given by [8, 9]
(5)EG=2μ(παZ1Z2c)2,
where α is the electromagnetic fine structure constant given by
(6)α=e2ℏc,
where e is the quantum of electric charge, ℏ is Planck’s quantum of action, and α is approximately 1/137 [9] for our universe. Thus the Gamow factor, which is determined by the electromagnetic force, and the nuclear fusion factor S(E) set the nuclear reaction cross section at low energies for nonresonant charged particles as [7, 10]
(7)σ(E)=S(E)Eexp[-(EGE)1/2],
and S(E) 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]
(8)S(E)≈S(0)+dS(0)dEE+12d2S(0)dE2E2=∑ν=02S(ν)(0)ν!Eν.
Substituting (7) and (8) in (4) we obtain
(9)r12=(1-12δ12)n1n2(8πμ)1/2×(1kT)3/2∑ν=02S(ν)(0)ν!×∫0∞Eνexp[-EkT-(EGE)1/2]dE.
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, 12–15]. It was initiated by Tsallis, the originator of nonextensive statistical mechanics [16–18], who has used q-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 q. 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:
(10)(1-q)nf1(x)=c1xγ-1[1-a(1-q)xδ]1/(1-q),a>0,δ>0,1-a(1-q)xδ>0,γ>0,q<1,
is the generalized type-1 beta form of the pathway model. This is a model with right tail cut-off for q<1. The Tsallis statistics for q<1 can be obtained from this model by putting γ=1 [16–18]. 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
(11)δ>0f2(x)=c2xγ-1[1+a(q-1)xδ]-1/(q-1),0<x<∞,q>1,a>0,γ>0,δ>0.
Here also for γ=1 we get the Tsallis statistics for q>1 [16–18]. Other standard distributions coming from this model are regular type-2 beta density, F-distribution, the Lévy model, and related models [22]. When q→1, f1(x) and f2(x) will reduce to the generalized gamma form of the pathway model given by
(12)f3(x)=c3xγ-1e-axδ,x>0.
This model covers generalized gamma, gamma, exponential, chi-square, the Weibull, the Maxwell-Boltzmann, the Rayleigh, and related densities. c1, c2, and c3 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
(13)fPD(E)dE=2π(q-1)3/2(πkT)3/2Γ(1/(q-1))Γ(1/(q-1)-3/2)×E[1+(q-1)EkT]-1/(q-1)dE,
for q>1, 1/(q-1)-3/2>0. The Maxwell-Boltzmann energy distribution can be retrieved from (13) by taking q→1. Thus the reaction rate probability integral given in (4) can be modified by using (13) and we get the extended reaction rate as
(14)r~12=(1-12δ12)n1n2(8πμ)1/2(q-1kT)3/2×Γ(1/(q-1))Γ(1/(q-1)-3/2)∑ν=02S(ν)(0)ν!×∫0∞Eν[1+(q-1)EkT]-1/(q-1)×exp[-(EGE)1/2]dE
for q>1, 1/(q-1)-3/2>0. Substituting y=E/kT and x=(EG/kT)1/2 we obtain the above integral in a more convenient form as follows:
(15)r~12=(1-12δ12)n1n2(8πμ)1/2(q-1)3/2×Γ(1/(q-1))Γ(1/(q-1)-3/2)∑ν=02(1kT)-ν+(1/2)S(ν)(0)ν!×∫0∞yν[1+(q-1)y]-1/(q-1)e-xy-1/2dy,
for q>1, 1/(q-1)-3/2>0. Here we consider the integral to be evaluated as
(16)I1q=∫0∞yν[1+(q-1)y]-1/(q-1)(q-1)×e-xy-1/2dy.
The integral can be evaluated by the techniques in applied analysis and can be obtained in closed form via Meijer’s G-function as [2, 23, 24]
(17)I1q=(π)-1/2(q-1)ν+1Γ(1/(q-1))×G1,33,1((q-1)x24|0,1/2,ν+12-1/(q-1)+ν)
which yields the nonresonant reaction rate probability integral in the extended case as
(18)r~12=(1-12δ12)n1n2(8μ)1/2π-1Γ(1/(q-1)-3/2)×∑ν=02(q-1kT)-ν+(1/2)S(ν)(0)ν!∑ν=02×G1,33,1[(q-1)EG4kT|0,1/2,ν+12-1/(q-1)+ν].
Meijer’s G-function and its properties can be seen in Mathai and Saxena [25] and Mathai [26]. We can obtain series expansions of the G-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 S(ν)(0) is approximately constant across the fusion window. Taking S(ν)(0)=0 for ν=1 and ν=2 and taking S0(0)=S(0), we obtain the extended reaction rate probability integral as
(19)r~12=(1-12δ12)n1n2[8(q-1)μkT]1/2π-1Γ(1/(q-1)-3/2)×S(0)G1,33,1[(q-1)EG4kT|0,1/2,12-1/(q-1)].
The series representation for (19) can be obtained as
(20)r~12=(1-12δ12)n1n2[8(q-1)μkT]1/2π-1Γ(1/(q-1)-3/2)×S(0){((q-1)EG4kT)rπΓ(1/(q-1)-1)×S(0)-2πΓ(1q-1-12)[(q-1)EG4kT]1/2×S(0)×F12(1q-1-12;32,12;-(q-1)EG4kT)×S(0)+(2π(q-1)EG4kT)×S(0)×∑r=0∞((q-1)EG4kT)r×S(0)(q-1)×[Ar-ln((q-1)EG4kT)]Br},
where
(21)Ar=Ψ(-12-r)+Ψ(1q-1+r)+Ψ(1+r)+Ψ(2+r),Br=(-1)rΓ(1/(q-1)+r)(3/2)rr!(1+r)!.
See the Appendix for detailed evaluation. As q→1 in (19), then, by using Stirling’s formula for gamma functions given by
(22)Γ(z+a)≈(2π)1/2zz+a-1/2e-z,|z|⟶∞,aisbounded,
we get the reaction rate probability integral in the Maxwell-Boltzmann case as
(23)r12=(1-12δ12)n1n2(8πμ)1/2(1kT)3/2S(0)×∫0∞exp[-EkT-(EGE)1/2]dE(24)=(1-12δ12)n1n2(8μkT)1/2S(0)π-1×G0,33,0[EG4kT|0,1/2,1-]
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 E/kT+(EG/E)1/2=g(E) (say). The minimum value of the function g(E), say E0, can be determined as
(25)ddE[EkT+(EGE)1/2]E=E0=1kT-12EG1/2E0-3/2=0⟹E0=EG1/2(kT2)3/2
and the function
(26)g(E0)=3(EG4kT)1/3=3Θ,
where Θ=(EG/4kT)1/3. Now by using the Laplace method [27, 28] we can obtain an approximate value for (23) as
(27)r12≈(1-12δ12)n1n28S(0)Θ2exp(-3Θ)3πμαZ1Z2c.
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 ϱ(r) considered by Haubold and Mathai [6, 11] and Mathai and Haubold [5] in the form
(28)ϱ(r)=ϱc[1-(rR)δ],δ>0,
where ϱc is the central density of the star, r is an arbitrary distance from the center, and R is the solar radius. This density function is capable of producing different density distributions by choosing the free parameter δ. Now we determine the quantities M(r),P(r), and T(r), the mass, the pressure, and the temperature at r.
By the equation of the mass conservation
(29)dM(r)dr=4πr2ϱ(r)
we get
(30)M(r)=4πϱc∫0rt2[1-(tR)δ]dt=4π3ϱcr3[1-3δ+3(rR)δ].
From (30), we get the central density as
(31)ϱc=3(δ+3)4πδM(R)R3.
If an element of a matter at a distance r 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
(32)dP(r)dr=-Gϱ(r)M(r)r2,
where G is the gravitational constant. Assuming that the pressure at the center of the sun is Pc and at the surface is zero, we get
(33)P(r)=Pc-G∫0rM(t)ϱ(t)t2dt=4π3Gϱc2R2[ξ-12(rR)2+δ+6(δ+2)(δ+3)(rR)δ+2=4π3Gϱc2R2-32(δ+1)(δ+3)(rR)2δ+2].
Using the boundary conditions P(R)=0, we get Pc=(4π/3)Gξϱc2R2, where
(34)ξ=12-δ+6(δ+2)(δ+3)+32(δ+1)(δ+3).
By the kinetic theory of gases, for a perfect gas, the pressure is given by
(35)P(r)=kNAμϱ(r)T(r).
For the temperature of interest for stellar models, we neglect the negligible radiation pressure from the total pressure and obtain from (35) the following:
(36)T(r)=μkNAP(r)ϱ(r)=4π3kNAGμϱcR2[1-(r/R)δ]×[(rR)2δ+2ξ-12(rR)2+δ+6(δ+2)(δ+3)(rR)δ+2-32(δ+1)(δ+3)(rR)2δ+2].
The central temperature Tc=(4/kNA)Gμξ(M(R)/R), 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 Lr=L(r) as the energy flux through the sphere of radius r, then we have
(37)dLrdr=4πr2ϱ(r)ε(r),
where ε(r) is the energy produced per second by nuclear reactions in each gram of stellar matter. ε(r) depends on the chemical composition in each gram of stellar matter. Here usually Lr is a constant but will be equal to L 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 ε(r) for one particular nuclear reaction. Now if we denote r~12(ϱ(r),T(r)) as the extended nonresonant thermonuclear reaction rate for the particles 1 and 2 defined by (19), then we will consider the energy generation rate ε12(r) and it can be written in terms of the extended reaction rates via
(38)ε12(r)=1ϱ(r)𝒞*ϱ2G1,33,1[(q-1)EG4kT|0,1/2,12-1/(q-1)],
where
(39)𝒞*=E12r~12(ϱ(r),T(r))ϱ2G1,33,1[(q-1)EG/4kT|0,1/2,12-1/(q-1)]
in which E12 is the amount of energy given off in a single reaction. It is to be noted that, by using the asymptotic behaviour of G1,33,1((q-1)EG/4kT) [26] and as q→1, 𝒞*→𝒞, the composite parameter considered by [9], which is defined as
(40)𝒞=E12r12ϱ2Θ2exp(3Θ)
for our universe 𝒞≈2×104 for proton-proton fusion under typical stellar conditions [9]. Then from (37) we have the total luminosity of the star by integration as follows:
(41)L(R)=∫0R4πr2ϱ(r)ε(r)dr.
If we are considering only one specific reaction defined as in (1), then we have
(42)L12(R)=∫0R4πr2ϱ(r)ε12(r)dr,
where the energy generation rate is defined in (38) and ϱ(r) is a suitable density distribution explaining the sun. Writing (42) in terms of r~12(ϱ(r),T(r)) we get
(43)L12(R)=∫0R4πr2𝒞*ϱ2G1,33,1[(q-1)EG4kT|0,1/2,12-1/(q-1)]dr=∫0R4πr2E12r~12(ϱ(r),T(r))dr.
The number density ni of a particle i, for a gas of mean density ϱ(r), can be expressed as
(44)ni(r)=ϱ(r)NAXiAi,
where NA stands for Avagadro’s constant, Ai is the atomic mass of particle i in atomic mass units, and Xi is the mass fraction of particle i such that ∑iXi=1. Substituting r~12(ϱ(r),T(r)) from (19) and using (44) we have
(45)L12(R)=∫0R4πr2E12(1-12δ12)NA2ϱ2(r)X1X2A1A2[8(q-1)μkT(r)]1/2×π-1Γ(1/(q-1)-3/2)S(0)G1,33,1×[(q-1)π2μ2kT(r)(Z1Z2e2ℏ)2|0,1/2,12-1/(q-1)]dr.
If we divide the “internal luminosity” L12(R⊙) by the amount of energy E12, then we get the total number of particles per second N12 liberated in the reaction given by (1) as follows:
(46)N12=L12(R)E12=4(1-12δ12)NA2X1X2A1A2[8(q-1)μk]1/2×1Γ(1/(q-1)-3/2)S(0)×∫0Rr2ϱ2(r)[T(r)]1/2G1,33,1[(Z1Z2e2ħ)2|0,(1/2),12-(1/(q-1))(q-1)π2μ2kT(r)×∫0Rr2ϱ2(r)[T(r)](1/2)2×(Z1Z2e2ℏ)2|0,1/2,12-1/(q-1)]dr=4(1-12δ12)NA2X1X2A1A2[8(q-1)μk]1/2×1Γ(1/(q-1)-3/2)S(0)12πi×∫LΓ(s)Γ(12+s)Γ(1+s)×Γ(1q-1-1-s)[(q-1)π2μ2k(Z1Z2e2ℏ)2]-s×∫0Rr2ϱ2(r)[T(r)]-1/2+sdrds.
For the density distribution defined in (28) introduced by Haubold and Mathai [5, 6] and the corresponding temperature distribution (36), we get
(47)∫0Rr2ϱ2(r)[T(r)]-1/2+sdr=[4πGμ3kNAϱcR2]s-1/2ϱc2×∫0Rr2[1-(rR)δ]5/2-s×[-12(rR)2+δ+6(δ+2)(δ+3)(rR)δ+2(δ+3)-32(δ+1)(δ+3)(rR)2δ+2]s-1/2dr,
where ξ is as defined in (34). If we put a substitution r=xR, then we get
(48)∫0Rr2ϱ2(r)[T(r)]-1/2+sdr=[4πGμ3kNAϱcR2]s-1/2ϱc2R3×∫01x2[1-xδ]5/2-s×[ξ-12x2+δ+6(δ+2)(δ+3)xδ+2(δ+2)-32(δ+1)(δ+3)x2δ+2]s-1/2dx.
Putting xδ=y and simplifying we obtain
(49)∫0Rr2ϱ2(r)[T(r)]-1/2+sdr=[4πGμξ3kNAϱcR2]s-1/2ϱc2R3δ×∫01y3/δ-1[1-y]5/2-s[1-u(y)]s-1/2dy,
where u(y) is defined as
(50)u(y)=y2/δξ[12-δ+6(δ+2)(δ+3)y+32(δ+1)(δ+3)y2].
As y→0, u(y)→0, and y→1, u(y)→1. If we take
(51)v(y)=12-δ+6(δ+2)(δ+3)y+32(δ+1)(δ+3)y2,
we have v(0)=1/2. The minimum value of v(y) is at y=(δ+1)(δ+6)/3(δ+2) and the value is 1/2-1/6((δ+6)2(δ+1)/(δ+3)(δ+2)2). Thus the minimum value is nonnegative since (δ+6)2(δ+1)/(δ+3)(δ+2)2 decreases steadily from 3 to 1 for all δ>0. Therefore v(y)≤0. Since ξ>0, for all δ>0, u(y)≤0 for all δ>0,ξ>0. Thus [1-u(y)]s-1/2≤0. Hence 0<u(y)<1 for 0<y<1 and for δ>0. Thus by using the binomial expansion we obtain
(52)[1-u(y)]s-1/2=∑m=0∞(1/2-s)mm![u(y)]m,
where (1/2-s)m is the Pochhammer symbol defined for a∈ℂ by
(53)(a)0=1,(a)m={a(a+1)⋯(a+m-1),Γ(a+m)Γ(a),m=1,2,…,a≠0
whenever Γ(a) exists. Taking δ=2 we get ξ=1/5 and
(54)[1-u(y)]s-1/2=(1-y)2s-1(1-12y)s-1/2.
Then from (49) we obtain
(55)∫0Rr2ϱ2(r)[T(r)]-1/2+sdr=[4πGμ15kNAϱcR2]s-1/2ϱc2R32×∫01y3/2-1[1-y]5/2+s-1(1-12y)s-1/2dy=[4πGμ15kNAϱcR2]s-1/2ϱc2R32×∑m=0∞(1/2-s)mm!12m=×∫01y3/2+m-1[1-y]5/2+s-1dy.
By using beta integral and using (53) we obtain
(56)∫0Rr2ϱ2(r)[T(r)]-1/2+sdr=[4πGμ15kNAϱcR2]s-1/2ϱc2R32×∑m=0∞Γ(1/2-s+m)m!Γ(1/2-s)12mΓ(3/2+m)Γ(5/2+s)Γ(4+s+m).
Now from (46) we obtain the total number of particles per second liberated in the reaction (1) as follows:
(57)N12=2NA2ϱc2R2μ(1-12δ12)X1X2A1A2[30(q-1)NAπGϱc]1/2×S(0)1Γ(1/(q-1)-3/2)×∑m=0∞12mΓ(3/2+m)m!12πi(3/2)×∫LΓ(s)Γ(12+s)Γ(1+s)Γ(52+s)(1+s)×Γ(1/(q-1)-1-s)Γ(1/2+m-s)Γ(1/2-s)Γ(4+m+s)(1+s)×[15π(q-1)NA8GϱcR2(Z1Z2e2ℏ)2]-sdsN12=2NA2ϱc2R2μ(1-12δ12)X1X2A1A2[30(q-1)NAπGϱc]1/2×S(0)1Γ(1/(q-1)-3/2)×∑m=0∞12mΓ(3/2+m)m!×G3,54,2[(Z1Z2e2ħ)2|0,1/2,1,5/2,1/22-(1/(q-1)),1/2-m,4+m15π(q-1)NA8GϱcR2(q-1)G3,54,2×(Z1Z2e2ℏ)2|0,1/2,1,5/2,1/22-1/(q-1),1/2-m,4+m].
For more details on G-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 E, they have to borrow an energy E from the thermal environment. The probability of such an energy is proportional to the Maxwell-Boltzmann energy exp[-E/kT]. The fusion will take place when the nuclei penerate the Coulomb barrier keeping them apart. The probability of penetration is given by the factor exp[-(EG/E)1/2]. The product of these two factors illustrates that fusion mostly occurs in the energy window given in the figure.
Schematic plot of the enrgy-dependent factors for the reaction rate probability integral: the Maxwell-Boltzmann energy density, nonresonant nuclear cross section, and the product of the Maxwell-Boltzmann density and the nonresonant cross section.
In Figure 2 the pathway energy density is plotted for q=0.7,0.9,1,1.2,1.4, respectively. For different values of q 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 q, namely, q=0.7,0.9,1,1.2,1.4, is also plotted. It is to be noted that as q→1 the pathway energy density coincides with the Maxwell-Boltzmann energy density and also the fusion window for the Maxwell-Boltzmann case in Figure 1.
Schematic plot of the energy-dependent factors for the extended reaction rate probability integral: pathway energy density, nonresonant nuclear cross section, and the product of the pathway density and the nonresonant cross section.
The curves in the figure represent pathway density [1+(q-1)(E/kT)]-1/(q-1) for (a) q=0.7, (b) q=0.9, (c) q=1, (d) q=1.2, and (e) q=1.4. (f) represents the nonresonant cross section exp[-(EG/E)1/2]. The product [1+(q-1)(E/kT)]-1/(q-1)exp[-(EG/E)1/2] is also represented for (g) q=0.7, (h) q=0.9, (i) q=1, (j) q=1.2, and (k) q=1.4.
Figure 3(a) shows the pathway energy density defined in (13) for q=1,1.2,1.3,1.4 and Figure 3(b) shows the Maxwell-Boltzmann energy density defined in (3). As q→1 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 q varies. If the Maxwell-Boltzmann density is the equilibrium situation, many other nonequilibrium situations are covered by the pathway energy density.
(a) Pathway energy density for q=1,1.2,1.3,1.4 and for k=1,T=100K. (b) The Maxwell-Boltzmann energy density for k=1,T=100K.
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 G-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 q, 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.
AppendixSeries 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 G-function given in (19). Consider the G-function as follows:
(A.1)G1,33,1((q-1)EG4kT|0,1/2,12-1/(q-1))=12πi∫c-i∞c+i∞Γ(s)Γ(12+s)Γ(1+s)=12πi×Γ(1q-1-1-s)((q-1)EG4kT)-sds.
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 Γ(s):s=0,-1,-2,…;
poles of Γ(1/2+s):s=-1/2,-3/2,-5/2,…;
poles of Γ(1+s):s=-1,-2,-3,….
Here the poles of Γ(s) and Γ(1+s) will coincide with each other at all points except at s=0. Note that the pole s=0 is a pole of order 1, s=-1/2,-3/2,-5/2,… are each of order 1, and s=-1,-2,-3,… are each of order 2. We know that
(A.2)lims→-r(s+r)Γ(s)=(-1)rr!,Γ(a-r)=(-1)rΓ(a)(1-a)r,Γ(a+m)=Γ(a)(a)m
when Γ(a) is defined, r=0,1,2,…;Γ(1/2)=π1/2,
(A.3)(a)r={a(a+1)⋯(a+r-1)ifr≥1,a≠01ifr=0.
The sum of the residues corresponding to the poles s=0 is given by
(A.4)R1=πΓ(1q-1-1).
The sum of the residues corresponding to the poles s=-1/2,-3/2,-5/2,… is
(A.5)R2=∑r=0∞(-1)rr!Γ(-12-r)Γ(12-r)×Γ(1q-1-12+r)[(q-1)EG4kT]-1/2+r=-2πΓ(1q-1-12)[(q-1)EG4kT]1/2×F21(1q-1-12;32,12;-(q-1)EG4kT),
where F21 is the hypergeometric function defined by
(A.6)F21(a;b,c;x)=∑r=0∞(a)r(b)r(c)rxrr!.
To obtain the sum of the residues corresponding to poles s=-1,-2,-3,… of order 2, we proceed as follows:
(A.7)R3=∑r=0∞lims→-1-r∂∂s[((q-1)EG4kT)-s(s+1+r)2Γ(1+s)Γ(s)Γ(12+s)=∑r=0∞lims→-1-r∂∂s×Γ(1q-1-1-s)((q-1)EG4kT)-s]=∑r=0∞lims→-1-r∂∂s[((q-1)EG4kT)-s(Γ2(2+s+r)Γ(1/2+s)=∑r=0∞lims→-1-r∂∂sn×Γ(1/(q-1)-1-s))=∑r=0∞lims→-1-r∂∂s×((s+r)2(s+r-1)2⋯(s+1)2s)-1=∑r=0∞lims→-1-r∂∂s×((q-1)EG4kT)-s]=∑r=0∞lims→-1-r∂∂sΦ(s),
where
(A.8)Φ(s)=Γ2(2+s+r)Γ(1/2+s)Γ(1/(q-1)-1-s)(s+r)2(s+r-1)2⋯(s+1)2s×((q-1)EG4kT)-s.
We have
(A.9)∂∂sΦ(s)=Φ(s)∂∂s[ln[Φ(s)]]lnΦ(s)=2ln[Γ(2+s+r)]+ln[Γ(12+s)]+ln[Γ(1q-1-1-s)]-sln((q-1)EG4kT)-2ln(s+r)-2ln(s+r-1)-⋯-2ln(s+1)-ln(s)∂∂s[ln[Φ(s)]]=2Ψ(2+s+r)+Ψ(12+s)+Ψ(1q-1-1-s)-ln((q-1)EG4kT)-2s+r-2s+r-1-⋯-2s+1-1slims→-1-r{∂∂sln[Φ(s)]}=Ψ(-12-r)+Ψ(1q-1+r)+Ψ(1+r)+Ψ(2+r)-ln((q-1)EG4kT),
where Ψ(z) is a Psi function or digamma function (see Mathai [26]) and Ψ(1)=-γ,γ=0.5772156649… is Euler’s constant. Now
(A.10)lims→-1-rΦ(s)=(-1)1+r2πΓ(1/(q-1)+r)(3/2)rr!(1+r)!×((q-1)EG4kT)1+r.
Then by using (A.7), (A.9), and (A.10) we get
(A.11)R3=∑r=0∞(-1)1+r2πΓ(1/(q-1)+r)(3/2)rr!(1+r)!((q-1)EG4kT)1+r×[Ψ(-12-r)+Ψ(1q-1+r)+Ψ(1+r)(2+r)RRR+Ψ(2+r)-ln((q-1)EG4kT)]=(2π(q-1)EG4kT)×∑r=0∞((q-1)EG4kT)r[Ar-ln((q-1)EG4kT)]Br,
where
(A.12)Ar=Ψ(-12-r)+Ψ(1q-1+r)+Ψ(1+r)+Ψ(2+r),Br=(-1)rΓ(1/(q-1)+r)(3/2)rr!(1+r)!.
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.
DavisR.Jr.Nobel Lecture: a half-century with solar neutrinos20037539859942-s2.0-014198774210.1103/RevModPhys.75.985HauboldH. J.KumarD.Extension of thermonuclear functions through the pathway model including Maxwell-Boltzmann and Tsallis distributions200829170762-s2.0-3864913934710.1016/j.astropartphys.2007.11.006HauboldH. J.JohnR. W.On the evaluation of an integral connected with the thermonuclear reaction rate in closed-form1978299522523210.1002/asna.19782990502FowlerW. A.CaughlanG. R.ZimmermanB. A.Thermonuclear rection rates1967552557010.1146/annurev.aa.05.090167.002521MathaiA. M.HauboldH. J.1988Berlin, GermanyAkademieHauboldH. J.MathaiA. M.Analytic representations of modified non-resonant thermonuclear reaction rates19863756856952-s2.0-002277869310.1007/BF00947916FowlerW. A.Experimental and theoretical nuclear astrophysics: the quest for the origin of the elements19845621491792-s2.0-3594902085210.1103/RevModPhys.56.149PhillipsA. C.19992ndChichester, UKJohn Wiley & SonsAdamsF. C.Stars in other universes: stellar structure with different fundamental constants200820088, article 102-s2.0-5634917079710.1088/1475-7516/2008/08/010CoradduM.KaniadakisG.LavagnoA.LissiaM.MezzoraniG.QuaratiP.Thermal distributions in stellar plasmas, nuclear reactions and solar neutrinos19992911531682-s2.0-0033095235HauboldH. J.MathaiA. M.On nuclear reaction rate theory198441380396CoradduM.LissiaM.MezzoraniG.QuaratiP.Super-Kamiokande hep neutrino best fit: a possible signal of non-Maxwellian solar plasma20033263-44734812-s2.0-004259450710.1016/S0378-4371(03)00359-5LavagnoA.QuaratiP.Classical and quantum non-extensive statistics effects in nuclear many-body problems20021335695802-s2.0-003648956210.1016/S0960-0779(01)00040-6LavagnoA.QuaratiP.Metastability of electron-nuclear astrophysical plasmas: motivations, signals and conditions200630532532592-s2.0-3384534075510.1007/s10509-006-9200-2LissiaM.QuaratiP.Nuclear astrophysical plasmas: ion distribution functions and fusion rates200536621121410.1051/epn:2005610TsallisC.Possible generalization of Boltzmann-Gibbs statistics1988521-24794872-s2.0-3364651648510.1007/BF01016429TsallisC.2009New York, NY, USASpringerGell-MannM.TsallisC.2004New York, NY, USAOxford University PressMathaiA. M.HauboldH. J.Pathway model, superstatistics, Tsallis statistics, and a generalized measure of entropy200737511101222-s2.0-3375154263310.1016/j.physa.2006.09.002SaxenaR. K.MathaiA. M.HauboldH. J.Astrophysical thermonuclear functions for Boltzmann-Gibbs statistics and Tsallis statistics20043443-46496562-s2.0-634425302710.1016/j.physa.2004.06.047MathaiA. M.A pathway to matrix-variate gamma and normal densities20053961–33173282-s2.0-1114429894410.1016/j.laa.2004.09.022MathaiA. M.HauboldH. J.On generalized distributions and pathways200837212210921132-s2.0-4304916728210.1016/j.physleta.2007.10.084HauboldH. J.KumarD.Fusion yield: Guderley model and Tsallis statistics20117711142-s2.0-7865041714110.1017/S0022377810000590KumarD.HauboldH. J.On extended thermonuclear functions through pathway model20104556987082-s2.0-7544912008410.1016/j.asr.2009.07.033MathaiA. M.SaxenaR. K.1973348New York, NY, USASpringerLecture Notes in MathematicsMathaiA. M.1993Oxford, UKClarendo PressErdélyiA.1956New York, NY, USADoverOlverF. W. J.1974New York, NY, USAAcademic PressHauboldH. J.MathaiA. M.Analytical results connecting stellar structure parameters and neutrino uxes1987442103116