Astrophysical thermonuclear functions

Stars are gravitationally stabilized fusion reactors changing their chemical composition while transforming light atomic nuclei into heavy ones. The atomic nuclei are supposed to be in thermal equilibrium with the ambient plasma. The majority of reactions among nuclei leading to a nuclear transformation are inhibited by the necessity for the charged participants to tunnel through their mutual Coulomb barrier. As theoretical knowledge and experimental verification of nuclear cross sections increases it becomes possible to refine analytic representations for nuclear reaction rates. Over the years various approaches have been made to derive closed-form representations of thermonuclear reaction rates (Critchfield, 1972; Haubold and John, 1978; Haubold, Mathai and Anderson, 1987). They show that the reaction rate contains the astrophysical cross section factor and its derivatives which has to be determined experimentally, and an integral part of the thermonuclear reaction rate independent from experimental results which can be treated by closed-form representation techniques in terms of generalized hypergeometric functions. In this paper mathematical/statistical techniques for deriving closed-form representations of thermonuclear functions, particularly the four integrals $$\begin{gathered} I_1 (z,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - zy^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ I_2 (z,d,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - zy^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ I_3 (z,t,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - z(y + 1)^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ I_4 (z,\delta ,b,v)\mathop = \limits^{def} \int\limits_0^\infty {y^v e^{ - y} e^{ - by^\delta } e^{ - zy^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } dy,} \hfill \\ \end{gathered} $$ will be summarized and numerical results for them will be given. The separation of thermonuclear functions from thermonuclear reaction rates is our preferred result. The purpose of the paper is also to compare numerical results for approximate and closed-form representations of thermonuclear functions. This paper completes the work of Haubold, Mathai, and Anderson (1987).