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.

Using the records of the fluxes of solar neutrinos from the Homestake, GALLEX (GNO), and Super-Kamiokande experiments, their statistical analyses were performed to search for whether there existed a time variation of these fluxes. The results of the analysis for the three experiments indicate that these fluxes are varying quasi-biennially. This means that both efficiencies of the initial p-p and the pp-III reactions of the proton-proton chain are varying quasi-biennially together with a period of about 26 months. Since this time variation prospectively generated by these two reactions strongly suggests that the efficiency of the proton-proton chain as the main energy source of the Sun has a tendency to vary quasi-biennially due to some chaotic or non-linear process taking place inside the gravitationally stabilized solar fusion reactor. It should be, however, remarked that, at the present moment, we have no theoretical reasoning to resolve this mysterious result generally referred to as the quasi-biennial periodicity in the time variation of the fluxes of solar neutrinos. There is an urgent need to search for the reason why such a quasi-biennial periodicity is caused through some physical process as related to nuclear fusion deep inside the Sun.

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.

Possible modification in the velocity distribution in the non-resonant reaction rates leads to an extended reaction rate probability integral. The closed form representation for these thermonuclear functions are used to obtain the stellar luminosity and neutrino emission rates. The composite parameter {C} that determines the standard nuclear reaction rate through the Maxwell-Boltzmann energy distribution is extended to {C}^* 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 in 2005 [Linear Algebra and Its Applications, 396, 317-328]. 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 non-resonant thermonuclear functions are compared by plotting them. Behavior of the new energy distribution, more general than Maxwell-Boltzmann is also studied.

The method for the evaluation of the non-resonant thermonuclear function in the Maxwell-Boltzmann case with depleted tail is discussed. Closed forms of the analytical results are obtained in computational format, and written in terms of the H-function in two variables. The standard non-resonant cases are extended to Tsallis reaction rates through the pathway model. Behavior of the depleted non-resonant thermonuclear function is studied. A comparison of the Maxwell-Boltzmann energy distribution with a more general energy distribution called pathway energy distribution is also done.

In 2010, the United Nations Committee on the Peaceful Uses of Outer Space began consideration of a new agenda item under a three-year work plan on the International Space Weather Initiative (ISWI). The main objectives of ISWI are to contribute to the development of the scientific insight necessary to improve understanding and forecasting capabilities of space weather as well as to education and public outreach. The United Nations Programme on Space Applications, implemented by the Office for Outer Space Affairs, is implementing ISWI in the framework of its United Nations Basic Space Science Initiative (UNBSSI), a long-term effort, launched in 1991, for the development of basic space science and for international and regional cooperation in this field on a worldwide basis, particularly in developing countries. UNBSSI encompassed a series of workshops, held from 1991 to 2004, which addressed the status of basic space science in Africa, Asia and the Pacific, Latin America and the Caribbean, and Western Asia. As a result several small astronomical research facilities have been inaugurated and education programmes at the university level were established. Between 2005 and 2009, the UNBSSI activities were dedicated to promoting activities related to the International Heliophysical Year 2007 (IHY), which contributed to the establishment of a series of worldwide ground-based instrument networks, a node of which is also operated by the Office for Outer Space Affairs. Building on these accomplishments, UNBSSI is now focussing on the ISWI.

Since 1990, the UN Programme on Space Applications leads the United Nations Basic Space Science Initiative by contributing to the international and regional development of astronomy and space science through annual UN/ESA/NASA/JAXA workshops on basic space science, International Heliophysical Year 2007, and the International Space Weather Initiative. Space weather is the conditions on the Sun and in the solar wind, magnetosphere, ionosphere and thermosphere that can influence the performance and reliability of space-borne and ground-based technological systems and can endanger human life or health. The programme also coordinates the development of IHY/ISWI low-cost, ground-based, world-wide instrument arrays. To date, 14 world-wide instrument arrays comprising approximately 1000 instruments (GPS receivers, magnetometers, spectrometers, particle detectors) are operating in more than 71 countries. The most recent workshop was hosted by the Republic of Korea in 2009 for Asia and the Pacific. Annual workshops on the ISWI have been scheduled to be hosted by Egypt in 2010 for Western Asia, Nigeria in 2011 for Africa, and Ecuador in 2012 for Latin America and the Caribbean.

Geometric generalized Mittag-Leffler distributions having the Laplace transform $\frac{1}{1+\beta\log(1+t^\alpha)},0<\alpha\le 2,\beta>0$ is introduced and its properties are discussed. Autoregressive processes with Mittag-Leffler and geometric generalized Mittag-Leffler marginal distributions are developed. Haubold and Mathai (2000) derived a closed form representation of the fractional kinetic equation and thermonuclear function in terms of Mittag-Leffler function. Saxena et al (2002, 2004a,b) extended the result and derived the solutions of a number of fractional kinetic equations in terms of generalized Mittag-Leffler functions. These results are useful in explaining various fundamental laws of physics. Here we develop first-order autoregressive time series models and the properties are explored. The results have applications in various areas like astrophysics, space sciences, meteorology, financial modeling and reliability modeling.

The paper provides a review of A.M. Mathai's applications of the theory of special functions, particularly generalized hypergeometric functions, to problems in stellar physics and formation of structure in the Universe and to questions related to reaction, diffusion, and reaction-diffusion models. The essay also highlights Mathai's recent work on entropic, distributional, and differential pathways to basic concepts in statistical mechanics, making use of his earlier research results in information and statistical distribution theory. The results presented in the essay cover a period of time in Mathai's research from 1982 to 2008 and are all related to the thematic area of the gravitationally stabilized solar fusion reactor and fractional reaction-diffusion, taking into account concepts of non-extensive statistical mechanics. The time period referred to above coincides also with Mathai's exceptional contributions to the establishment and operation of the Centre for Mathematical Sciences, India, as well as the holding of the United Nations (UN)/European Space Agency (ESA)/National Aeronautics and Space Administration (NASA) of the United States/ Japanese Aerospace Exploration Agency (JAXA) Workshops on basic space science and the International Heliophysical Year 2007, around the world. Professor Mathai's contributions to the latter, since 1991, are a testimony for his social conscience applied to international scientific activity.

The Maxwell-Boltzmannian approach to nuclear reaction rate theory is extended to cover Tsallis statistics (Tsallis, 1988) and more general cases of distribution functions. An analytical study of respective thermonuclear functions is being conducted with the help of statistical techniques. The pathway model, recently introduced by Mathai (2005), is utilized for thermonuclear functions and closed-form representations are obtained in terms of H-functions and G-functions. Maxwell-Boltzmannian thermonuclear functions become particular cases of the extended thermonuclear functions. A brief review on the development of the theory of analytic representations of nuclear reaction rates is given.

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.

The World Space Observatory is an unconventional space project proceeding via distributed studies. The present design, verified for feasibilty, consists of a 1.7-meter telescope operating at the second Largangian point of the Earth-Sun system. The focal plane instruments consist of three UV spectrometers covering the spectral band from Lyman alpha to the atmospheric cutoff with R~55,000 and offering long-slit capability over the same band with R~1,000. In addition, a number of UV and optical imagers view adjacent fields to that sampled by the spectrometers. Their performance compares well with that of HST/ACS and the spectral capabilities of WSO rival those of HST/COS.

We summarize the capabilities of the World Space Observatory (UV) Project (WSO/UV). An example of the importance of this project (with a planned launch date of 2007/8) for the study of Classical Novae is given.

Since 1990, the United Nations in cooperation with the European Space Agency is holding annually a workshop on basic space science for the benefit of the worldwide development of astronomy. These workshops have been held in countries of Asia and the Pacific (India, Sri Lanka), Latin America and the Caribbean (Costa Rica, Colombia, Honduras), Africa (Nigeria), Western Asia (Egypt, Jordan), and Europe (Germany, France). Additional to the scientific benefits of the workshops and the strengthening of international cooperation, the workshops lead to the establishment of astronomical telescope facilities in Colombia, Egypt, Honduras, Jordan, Morocco, Paraguay, Peru, Philippines, Sri Lanka, and Uruguay. The annual UN/ESA Workshops continue to pursue an agenda to network these astronomical telescope facilities through similar research and education programmes. Teaching material and hands-on astrophysics material has been developed for the operation of such astronomical telescope facilities in an university environment.

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 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.

In 1990, the United Nations, in cooperation with the European Space Agency, initiated the organization of a series of annual Workshops on Basic Space Science for the benefit of astronomers and space scientists in Asia and the Pacific, Latin American and the Caribbean, Africa, Western Asia, and Europe. This article summarizes accomplishments of these Workshops (1991-1998) and their follow-up projects with a view to enhance the worldwide development of astronomy and space science. The Workshops are being considered unique and a model for such an endeavor.

Wavelet analysis offers an alternative to Fourier based time-series analysis and is particularly useful when spectral features are time dependent. We analyze the solar neutrino capture rate of the radiochemical Homestake chlorine experiment with abbreviated Morlet wavelets, using Foster's (AJ, 111, 1709 (1996)) rescaled wavelet technique. We emphasize the complementarity of wavelet analysis to Fourier analysis. The Homestake data seem to contain a harmonic content with periodicities of 4.76 yr, 1.89 yr, 0.85 yr, and 0.51 yr. Wavelet analysis reveals that the 4.76 yr and 1.89 yr periods show an almost constant behavior over the 25 yr Homestake data record, while the 0.85 yr and 0.51 yr periods exhibit a transient phenomenon. The analysis does not show strong evidence for a period of the solar 11 yr cycle.

In 1990 the United Nations in cooperation with the European Space Agency initiated the organization of a series of annual Workshops on Basic Space Science for the benefit of astronomers and space scientists in Asia and the Pacific, Latin America and the Caribbean, Africa, Western Asia, and Europe. This article summarizes accomplishments of these Workshops and their follow-up projects.

The physical principles and the mathematical structure involved in deriving an analytical representation of the internal structure of the sun are discussed. For a two-parameter family of a non-linear matter density distribution, the run of mass, pressure, temperature, and luminosity throughout the sun are expressed in terms of Gauss' hypergeometric function. The system of differential equations governing hydrostatic equilibrium and energy conservation for the sun generates another field of application of special functions.

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.

The paper provides results of the Fourier spectrum analysis of the new Ar-37 production rate data of the Homestake solar neutrino experiment and compares them with results for earlier data, revealing the harmonic content in the Ar-37 production in the Homestake experiment.

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.

Physical principles and mathematical structure involved in deriving an analytical representation of the internal structure of the Sun is discussed. For a two-parameter family of a non-linear matter density distribution, the run of mass, pressure, temperature, and luminosity throughout the Sun is presented in terms of Gauss' hypergeometric function. The system of differential equations governing hydrostatic equilibrium and energy conservation for the spherical Sun is proved to be a laboratory for the application of special functions.

The Sun is forced into an orbit around the barycenter of the solar system because of the changing mass distribution of the planets. Solar-planetary-lunar dynamic relationships may form a new basis for understanding and predicting cyclic solar forcing functions on the Earth's climate.

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.

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.

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.

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.