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.