Experimental measurements in terrestrial laboratory, space and astrophysical observations of variation and fluctuation of nuclear decay constants, measurements of large enhancements in fusion reaction rate of deuterons implanted in metals and electron capture by nuclei in solar core indicate that these processes depend on the environment where take place and possibly also on the fluctuation of some extensive parameters and eventually on stellar energy production. Electron screening is the first important environment effect. We need to develop a treatment beyond the Debye-Huckel screening approach, commonly adopted within global thermodynamic equilibrium. Advances in the description of these processes can be obtained by means of q-thermostatistics and/or superstatistics for metastable states. This implies to handle without ambiguities the case q<1.

Stimulated by the recent debate on the physical relevance and on the predictivity of q-Gaussian formalism, we present specific analytical expressions for the parameters characterizing non-Gaussian distributions, such as the nonextensive parameter q, expressions that we have proposed for different physical systems, an important example being plasmas in the stellar cores.

Screened potential, modified by non standard electron cloud distributions responsible for the shielding effect on fusion of reacting nuclei in astrophysical plasmas, is derived. The case of clouds with depleted tails in space coordinates is discussed. The modified screened potential is obtained both from statistical mechanics arguments based on fluctuations of the inverse of the Debye-Huckel radius and from the solution of a Bernoulli equation used in generalized statistical mechanics. Plots and tables useful in evaluating penetration probability at any energy are provided.

This article illustrates how very small deviations from the Maxwellian exponential tail, while leaving unchanged bulk quantities, can yield dramatic effects on fusion reaction rates and discuss several mechanisms that can cause such deviations.

Density and temperature conditions in many stellar core (like the solar core) imply the presence of nonideal plasma effects with memory and long-range interactions between particles. This aspect suggests the possibility that the stellar core could not be in a global thermodynamical equilibrium but satisfies the conditions of a metastable state with a stationary (nonextensive) power law distribution function among ions. The order of magnitude of the deviation from the standard Maxwell-Boltzmann distribution can be derived microscopically by considering the presence of random electrical microfields in the stellar plasma. We show that such a nonextensive statistical effect can be very relevant in many nuclear astrophysical problems.

By minimizing free energy density, we show that the stellar core of a hydrogen burning star is not in a global thermodynamical equilibrium unless density, temperature, mass and composition assume given values. The core (as the solar interior) may be viewed more appropriately as a metastable state with very long lifetime. Slightly non-extensive distribution function could be the natural distribution for a weakly non-ideal plasma like a stellar core and represents a more appropriate approximation to this system than a Maxwellian distribution, without affecting bulk properties of stars.

We study two different physical scenarios of thermonuclear reactions in stellar plasmas proceeding through a narrow resonance at low energy or through the low energy wing of a wide resonance at high energy. Correspondingly, we derive two approximate analytical formulae in order to calculate thermonuclear resonant reaction rates inside very coupled and non ideal astrophysical plasmas in which non-extensive effects are likely to arise. Our results are presented as simple first order corrective factors that generalize the well known classical rates obtained in the framework of Maxwell-Boltzmann statistical mechanics. As a possible application of our results, we calculate the dependence of the total corrective factor with respect to the energy at which the resonance is located, in an extremely dense and non ideal carbon plasma.

We show that moderate deviations from the Maxwell-Boltzmann energy distribution can increase deuterium reaction rates enough to contribute to the heating of Jupiter. These deviations are compatible with the violation of extensivity expected from temperature and density conditions inside Jupiter.

We show that in weakly non-ideal plasmas, like the solar interior, both non-extensivity and quantum uncertainty (a' la Galitskii and Yakimets) should be taken into account to derive equilibrium ion distribution functions and to estimate nuclear reaction rates and solar neutrino fluxes.

We describe the physical motivations why we are urged to look for an ionic equilibrium distribution function slightly different from the Maxwellian to calculate nuclear reaction rates in the solar core plasma answering, at the same time, to a recent paper by Bahcall et al. (astro-ph/0010055).

Nonextensive and quantum uncertainty effects (related to the quasiparticles composing the stellar core) have strong influence on the nuclear rates and, of course, affect solar neutrino fluxes. Both effects do coexist and are due to the frequent collisions among the ions. The weakly nonextensive nature of the solar core is confirmed. The range of predictions for the neutrino fluxes is enlarged and the solar neutrino problem becomes less dramatic.

The solar core, because of its density and temperature, is not a weakly-interacting or a high-temperature plasma. Collective effects have time scales comparable to the average time between collisions, and the microfield distribution influences the particle dynamics. In this conditions ion and electron diffusion is a non-Markovian process, memory effects are present and the equilibrium statistical distribution function differs from the Maxwellian one. We show that, even if the deviations from the standard velocity distribution that are compatible with our present knowledge of the solar interior are small, they are sufficient to sensibly modify the sub-barrier nuclear reaction rates. The consequent changes of the neutrino fluxes are comparable to the flux deficits that constitute the solar neutrino problem.

Density and temperature conditions in the solar core suggest that the microscopic diffusion of electrons and ions could be nonstandard: Diffusion and friction coefficients are energy dependent, collisions are not two-body processes and retain memory beyond the single scattering event. A direct consequence of nonstandard diffusion is that the equilibrium energy distribution of particles departs from the Maxwellian one (tails goes to zero more slowly or faster than exponentially) modifying the reaction rates. This effect is qualitatively different from temperature and/or composition modification: Small changes in the number of particles in the distribution tails can strongly modify the rates without affecting bulk properties, such as the sound speed or hydrostatic equilibrium, which depend on the mean values from the distribution. This mechanism can considerably increase the range of predictions for the neutrino fluxes allowed by the current experimental values (cross sections and solar properties) and can be used to reduce the discrepancy between these predictions and the solar neutrino experiments.

Nuclear reactions in stars occur between nuclei in the high-energy tail of the energy distribution and are sensitive to possible deviations from the standard equilibrium thermal-energy distribution. We are able to derive strong constraints on such deviations by using the detailed helioseismic information of the solar structure. If a small deviation is parameterized with a factor exp{-delta*(E/kT)^2}, we find that delta should lie between -0.005 and +0.002. However, even values of delta as small as 0.003 would still give important effects on the neutrino fluxes.

In this paper we will show that, because of the long-range microscopic memory of the random force, acting in the solar core, mainly on the electrons and the protons than on the light and heavy ions (or, equally, because of anomalous diffusion of solar core constituents of light mass and of normal diffusion of heavy ions), the equilibrium statistical distribution that these particles must obey, is that of generalized Boltzmann-Gibbs statistics (or the Tsallis non-extensive statistics), the distribution differing very slightly from the usual Maxwellian distribution. Due to the high-energy depleted tail of the distribution, the nuclear rates are reduced and, using earlier results on the standard solar model neutrino fluxes, calculated by Clayton and collaborators, we can evaluate fluxes in good agreement with the experimental data. While proton distribution is only very slightly different from Maxwellian there is a little more difference with electron distribution. We can define one central electron temperature as a few percent higher than the ion central temperature nearly equal to the standard solar model temperature. The difference is related to the different reductions with respect to the standard solar model values needed for $B$ and $CNO$ neutrinos and for $Be$ neutrinos.

We show that the observational data recently provided by Giovanelli et al. (1996 a, b) and discussed by Bahcall and Oh (1996) concerning the velocity distribution of clusters of galaxies can be naturally fitted by a statistical distribution which generalizes the Maxwell-Boltzmann one (herein recovered for the entropic index $q=1$). Indeed a recent generalization of the Boltzmann-Gibbs thermostatistics suggests for this problem that the probability function is given, within a simple phenomenological model, by \begin{eqnarray} P(>v) & = & \frac{ \int_{v}^{v_{max}} dv \left[ 1-(1-q)(v/v_0)^2 \right]^\frac{q}{1-q} } {\int_{0}^{v_{max}} dv \left[ 1-(1-q)(v/v_0)^2 \right]^\frac{q}{1-q}}, \nonumber \\ v_{max} & \equiv & \left\{ \begin{array}{ll} v_0 (1-q)^{-1/2} & \mbox{if $q<1$} \\ \infty & \mbox{if $q \geq 1$.} \end{array}\right. \nonumber \end{eqnarray} A remarkably good fitting with the data is obtained for $q=0.23^{+0.07}_{-0.05}$ and $v_0=490\pm 5$ km s$^{-1}$.

The generalized Tsallis statistics produces a distribution function appropriate to describe the interior solar plasma, thought as a stellar polytrope, showing a tail depleted respect to the Maxwell-Boltzmann distribution and reduces to zero at energies greater than about $20 \, k_{_B} T$. The Tsallis statistics can theoretically support the distribution suggested in the past by Clayton and collaborators, which shows also a depleted tail, to explain the solar neutrino counting rate.