We are going back to the roots of the original solar neutrino problem: analysis of data from solar neutrino experiments. The application of standard deviation analysis (SDA) and diffusion entropy analysis (DEA) to the SuperKamiokande I and II data reveals that they represent a non-Gaussian signal. The Hurst exponent is different from the scaling exponent of the probability density function and both Hurst exponent and scaling exponent of the probability density function of the SuperKamiokande data deviate considerably from the value of 0.5 which indicates that the statistics of the underlying phenomenon is anomalous. To develop a road to the possible interpretation of this finding we utilize Mathai's pathway model and consider fractional reaction and fractional diffusion as possible explanations of the non-Gaussian content of the SuperKamiokande data.

This paper deals with a general class of integrals, the particular cases of which are connected to outstanding problems in astronomy and physics. Reaction rate probability integrals in the theory of nuclear reaction rates, Kr\"atzel integrals in applied analysis, inverse Gaussian distribution, generalized type-1, type-2 and gamma families of distributions in statistical distribution theory, Tsallis statistics and Beck-Cohen superstatistics in statistical mechanics and the general pathway model are all shown to be connected to the integral under consideration. Representations of the integral in terms of generalized special functions such as Meijer's G-function and Fox's H-function are also pointed out.

The pathway model of Mathai (2005) mainly deals with the rectangular matrix-variate case. In this paper the scalar version is shown to be associated with a large number of probability models used in physics. Different families of densities are listed here, which are all connected through the pathway parameter 'alpha', generating a distributional pathway. The idea is to switch from one functional form to another through this parameter and it is shown that basically one can proceed from the generalized type-1 beta family to generalized type-2 beta family to generalized gamma family when the real variable is positive and a wider set of families when the variable can take negative values also. For simplicity, only the real scalar case is discussed here but corresponding families are available when the variable is in the complex domain. A large number of densities used in physics are shown to be special cases of or associated with the pathway model. It is also shown that the pathway model is available by maximizing a generalized measure of entropy, leading to an entropic pathway. Particular cases of the pathway model are shown to cover Tsallis statistics (Tsallis, 1988) and the superstatistics introduced by Beck and Cohen (2003).

The pathway model of Mathai (2005) is shown to be inferable from the maximization of a certain generalized entropy measure. This entropy is a variant of the generalized entropy of order 'alpha', considered in Mathai and Rathie (1975), and it is also associated with Shannon, Boltzmann-Gibbs, Renyi, Tsallis, and Havrda-Charvat entropies. The generalized entropy measure introduced here is also shown to haveinteresting statistical properties and it can be given probabilistic interpretations in terms of inaccuracy measure, expected value, and information content in a scheme. Particular cases of the pathway model are shown to be Tsallis statistics (Tsallis, 1988) and superstatistics introduced by Beck and Cohen (2003). The pathway model's connection to fractional calculus is illustrated by considering a fractional reaction equation.

An analytic study of the reaction probability integrals corresponding to the various forms of the slowly varying cross-section factor $S(E)$ is attempted. Exact expressions for reaction probability integrals are expressed in terms of the extended gamma functions.

The paper discusses the solution of a simple kinetic equation of the type used for the computation of the change of the chemical composition in stars like the Sun. Starting from the standard form of the kinetic equation it is generalized to a fractional kinetic equation and its solutions in terms of H-functions are obtained. The role of thermonuclear functions, which are also represented in terms of G- and H-functions, in such a fractional kinetic equation is emphasized. Results contained in this paper are related to recent investigations of possible astrophysical solutions of the solar neutrino problem.

An introductory account is given of the understanding of the structure of the universe. At present the most plausible theory of the origin of the universe is that it formed from the explosion of an extremely hot and dense fireball several billion years ago. During the first few seconds after the big bang, the energy density was so great that only fundamental particles (leptons, quarks, gauge bosons) existed. As the universe cooled and expanded after the big bang, nuclei and atoms formed and condensed into galaxies and stars and systems of them. The fundamental particles and a wide range of gravitational aggregates of them constitute the small-scale and large-scale structure of the present universe. The current knowledge of the elements of the structure of the universe based on the standard big bang model and the standard model of fundamental particles is considered.

Nuclear reactions govern major aspects of the chemical evolution od galaxies and stars. Analytic study of the reaction rates and reaction probability integrals is attempted here. Exact expressions for the reaction rates and reaction probability integrals for nuclear reactions in the case of nonresonant, modified nonresonant, screened nonresonant and resonant cases are given. These are expressed in terms of H-functions, G-functions and in computable series forms. Computational aspects are also discussed.

Nuclear reactions govern major aspects of the chemical evolution of galaxies and stars. Analytic study of the reaction rates and reaction probability integrals is attempted here. Exact expressions for the reaction rates and reaction probability integrals for nuclear reactions in the cases of nonresonant, modified nonresonant, screened nonresonant and resonant cases are given. These are expressed in terms of H-functions, G-functions and in computable series forms. Computational aspects are also discussed.

The paper is based on lectures on the internal structure and evolution of the Sun. It contains an outline of observational and theoretical elements of the Standard Solar Model, emphasizing also recent research results on solar rotation, magnetic fields, thermonuclear energy generation, and solar activity.

This paper examines a number of analytic models which can describe the gravitationally stabilized fusion reactor of the Sun. An analytical multiparameter model is shown to reproduce various structure parameters such as density, mass, pressure and temperature throughout the solar core, which have been obtained by solving numerically the system of solar structure differential equations. For simplicity, numerical and analytical results are discussed in detail for a two-parameter model.

Hydrostatic equilibrium and energy conservation determine the conditions in the gravitationally stabilized solar fusion reactor. We assume a matter density distribution varying non-linearly through the central region of the Sun. The analytic solutions of the differential equations of mass conservation, hydrostatic equilibrium, and energy conservation, together with the equation of state of the perfect gas and a nuclear energy generation rate $\epsilon=\epsilon_0 \rho^nT^m$, are given in terms of Gauss' hypergeometric function. This model for the structure of the Sun gives the run of density, mass, pressure, temperature, and nuclear energy generation through the central region of the Sun. Because of the assumption of a matter density distribution, the conditions of hydrostatic equilibrium and energy conservation are separated from the mode of energy transport in the Sun.

Four operating neutrino observatories confirm the long standing discrepancy between detected and predicted solar neutrino flux. Among these four experiments the Homestake experiment is taking data for almost 25 years. The reliability of the radiochemical method for detecting solar neutrinos has been tested recently by the GALLEX experiment. All efforts to solve the solar neutrino problem by improving solar, nuclear, and neutrino physics have failed so far. This may also mean that the average solar neutrino flux extracted from the four experiments may not be the proper quantity to explain the production of neutrinos in the deep interior of the Sun. Occasionally it has been emphasized that the solar neutrino flux may vary over time. In this paper we do address relations among specific neutrino fluxes produced in the proton-proton chain that are imposed by the coupled systems of non-linear partial differential equations of solar structure and kinetic equations by focusing our attention on a statistical interpretation of selected kinetic equation of PPII/PPIII branch reactions of the proton-proton chain. A fresh look at the statistical implications for the outcome of kinetic equations for nuclear reactions may shed light on recent claims that the $^7$Be-neutrino flux of the Sun is suppressed in comparison to the pp- and $^8$B neutrino fluxes and may hint at that the solar neutrino flux is indeed varying over time as shown by the Homestake experiment.

Hydrostatic equilibrium and energy conservation determine the conditions in the gravitationally stabilized solar fusion reactor. We assume a matter density distribution varying non- linearly through the central region of the Sun. The analytic solutions of the differential equations of mass conservation, hydrostatic equilibrium, and energy conservation, together with the equation of state of the perfect gas and a nuclear energy generation rate $\epsilon=\epsilon_0\rho^nT^m$, are given in terms of Gauss' hypergeometric function. This model for the structure of the Sun gives the run of density, mass pressure, temperature, and nuclear energy generation through the central region of the Sun. Because of the assumption of a matter density distribution, the conditions of hydrostatic equilibrium and energy conservation are separated from the mode of energy transport in the Sun.

The study of solar neutrinos may provide important insights into the physics of the central region of the Sun. Four solar neutrino experiments have confirmed the solar neutrino problem but do not clearly indicate whether solar physics, nuclear physics, or neutrino physics have to be improved to solve it. Nonlinear relations among the different neutrino fluxes are imposed by two coupled systems of differential equations governing the internal structure and time evolution of the Sun. We assume that the results of the four neutrino experiments are correct and are concerned not with the discrepancy between the average rate and the predicted rate, but with a possible time dependence of the argon production rate as revealed in the Homestake experiment over a time period of 20 years. Based on the subtlety of the solar neutrino problem we review here quantitatively the physical laws employed for understanding the internal solar structure and conjecture that the interlink between specific nuclear reactions of the PPIII-branch of the proton-proton chain may allow the high-energy solar neutrino flux to vary over time.

The paper presents results for deriving closed-form analytic solutions of the non-relativistic linear perturbation equations, which govern the evolution of inhomogeneities in a homogeneous spatially flat multicomponent cosmological model. Mathematical methods to derive computable forms of the perturbations are outlined.

As theoretical knowledge and experimental verification of nuclear cross sections increases it becomes possible to refine analytic representations for nuclear reaction rates. In this paper mathematical/statistical techniques for deriving closed-form representations of thermonuclear functions are summarized and numerical results for them are given.The purpose of the paper is also to compare numerical results for approximate and closed-form representations of 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 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).

This article gives an brief outline of the applications of generalized special functions such as generalized hypergeometric functions, G-functions and H-functions into the general area of nuclear energy generation and reaction rate theory such as the energy generation in a simple stellar model and nuclear reaction rates in non-resonant and resonant as well as screened non-resonant

The paper presents results for deriving closed-form analytic solutions of the evolution of inhomogeneities in a homogeneous spatially flat multicomponent cosmological model. Mathematical methods to derive computable forms of the perturbations are outlined.