We present several arguments which favor the scenario of two coexisting families of compact stars: hadronic stars and quark stars. Besides the well known hyperon puzzle of the physics of compact stars, a similar puzzle exists also when considering delta resonances. We show that these particles appear at densities close to twice saturation density and must be therefore included in the calculations of the hadronic equation of state. Such an early appearance is strictly related to the value of the L parameter of the symmetry energy that has been found, in recent phenomenological studies, to lie in the range $40<L<62$ MeV. We discuss also the threshold for the formation of deltas and hyperons for hot and lepton rich hadronic matter. Similarly to the case of hyperons, also delta resonances cause a softening of the equation of state which makes it difficult to obtain massive hadronic stars. Quark stars, on the other hand, can reach masses up to $2.75 M_{\odot}$ as predicted by perturbative QCD calculations. We then discuss the observational constraints on the masses and the radii of compact stars. The tension between the precise measurements of high masses and the indications of the existence of very compact stellar objects (with radii of the order of $10$ km) is relieved when assuming that very massive compact stars are quark stars and very compact stars are hadronic stars. Finally, we discuss recent interesting measurements of the eccentricities of the orbits of millisecond pulsars in low mass X-ray binaries. The high values of the eccentricities found in some cases could be explained by assuming that the hadronic star, initially present in the binary system, converts to a quark star due to the increase of its central density.

We propose a model for short duration gamma-ray bursts (sGRBs) based on the formation of a quark star after the merger of two neutron stars. We assume that the sGRB central engine is a proto-magnetar, which has been previously invoked to explain the plateau-like X-ray emission observed following both long and short GRBs. Here, we show that: i) a few milliseconds after the merger it is possible to form a stable and massive star made in part of quarks; ii) during the early cooling phase of the incompletely formed quark star, the flux of baryons ablated from the surface by neutrinos is large and it does not allow the outflow to achieve a bulk Lorentz factor high enough to produce a GRB; iii) after the quark burning front reaches the stellar surface, baryon ablation ceases and the jet becomes too baryon poor to produce a GRB; iv) however, between these two phases a GRB can be produced over the finite timescale required for the baryon pollution to cease; a characteristic timescale of the order of $\sim 0.1 $ s naturally results from the time the conversion front needs to cover the distance between the rotational pole and the latitude of the last closed magnetic field line; v) we predict a correlation between the luminosity of the sGRB and its duration, consistent with the data; vi) our model also predicts a delay of the order of ten seconds between the time of the merger event and the sGRB, allowing for the possibility of precursor emission and implying that the jet will encounter the dense cocoon formed immediately after the merger.

We discuss the formation of $\Delta$ isobars in neutron star matter. We show that their threshold density strictly correlates with the density derivative of the symmetry energy of nuclear matter, the $L$ parameter. By restricting $L$ to the range of values indicated by recent experimental and theoretical analysis, i.e. $40$ MeV $\lesssim L \lesssim 62$ MeV, we find that $\Delta$ isobars appear at a density of the order of 2$\div$3 times nuclear matter saturation density, i.e. the same range for the appearance of hyperons. The range of values of the couplings of the $\Delta$s with the mesons is restricted by the analysis of the data obtained from photoabsorption, electron and pion scattering on nuclei. If the potential of the $\Delta$ in nuclear matter is close to the one indicated by the experimental data then the equation of state becomes soft enough that a "$\Delta$ puzzle" exists, similar to the "hyperon puzzle" widely discussed in the literature.

The existence of stars with masses up to $2 M_{\odot}$ and the hints of the existence of stars with radii smaller than $\sim 11$ km seem to require, at the same time, a stiff and a soft hadronic equation of state at large densities. We argue that these two apparently contradicting constraints are actually an indication of the existence of two families of compact stars: hadronic stars which could be very compact and quark stars which could be very massive. In this respect, a crucial role is played, in the hadronic equation of state, by the delta isobars whose early appearance shifts to large densities the formation of hyperons. We also discuss how recent experimental information on the symmetry energy of nuclear matter at saturation indicate, indirectly, an early appearance of delta isobars in neutron star matter.

The existence of neutron stars with masses of $\sim 2\,M_{\odot}$ requires a stiff equation of state at high densities. On the other hand, the necessary appearance also at high densities of new degrees of freedom, such as hyperons and $\Delta$ resonances, can lead to a strong softening of the equation of state with resulting maximum masses of $\sim 1.5\, M_{\odot}$ and radii smaller than $\sim 10$ km. Hints for the existence of compact stellar objects with very small radii have been found in recent statistical analysis of quiescent low-mass X-ray binaries in globular clusters. We propose an interpretation of these two apparently contradicting measurements, large masses and small radii, in terms of two separate families of compact stars: hadronic stars, whose equation of state is soft, can be very compact, while quark stars, whose equation of state is stiff, can be very massive. In this respect an early appearance of $\Delta$ resonances is crucial to guarantee the stability of the branch of hadronic stars. Our proposal could be tested by measurements of radii with an error of $\sim 1$ km, which is within reach of the planned LOFT satellite, and it would be further strengthened by the discovery of compact stars heavier than $\sim 2\,M_{\odot}$.

We investigate the bulk properties of protoneutron stars in the framework of a relativistic mean field theory based on nonextensive statistical mechanics, characterized by power-law quantum distributions. We study the relevance of nonextensive statistical effects on the beta-stable equation of state at fixed entropy per baryon, in presence and in absence of trapped neutrinos, for nucleonic and hyperonic matter. We show that nonextensive statistical effects could play a crucial role in the structure and in the evolution of the protoneutron stars also for small deviations from the standard Boltzmann-Gibbs statistics.

We discuss new limits on masses and radii of compact stars and we conclude that they can be interpreted as an indication of the existence of two classes of stars: "normal" compact stars and "ultra-compact" stars. We estimate the critical mass at which the first configuration collapses into the second.

We study the hydrodynamical transition from an hadronic star into a quark or a hybrid star. We discuss the possible mode of burning, using a fully relativistic formalism and realistic Equations of State in which hyperons can be present. We take into account the possibility that quarks form a diquark condensate. We also discuss the formation of a mixed phase of hadrons and quarks, and we indicate which region of the star can rapidly convert in various possible scenarios. An estimate of the final temperature of the system is provided. We find that the conversion process always corresponds to a deflagration and never to a detonation. Hydrodynamical instabilities can develop on the front. We estimate the increase in the conversion's velocity due to the formation of wrinkles and we find that, although the increase is significant, it is not sufficient to transform the deflagration into a detonation in essentially all realistic scenarios. Concerning convection, it does not always develop. In particular the system does not develop convection if hyperons are not present in the initial phase and if the newly formed quark phase is made of ungapped (or weakly gapped) quarks. At the contrary, the process of conversion from ungapped quark matter to gapped quarks always allows the formation of a convective layer. Finally, we discuss possible astrophysical implications of our results.

We study the transition from hadronic matter to a mixed phase of quarks and hadrons at high baryon and isospin densities reached in heavy ion collisions. We focus our attention on the role played by the nucleon symmetry energy at high density.In this respect the inclusion of a scalar isovector meson, the \delta-coupling, in the Hadron Lagrangian appears rather important. We study in detail the formation of a drop of quark matter in the mixed phase, and we discuss the effects on the quark drop nucleation probability of the finite size and finite time duration of the high density region. We find that, if the parameters of quark models are fixed so that the existence of quark stars is allowed, then the density at which a mixed phase starts forming drops dramatically in the range Z/A \sim 0.3--0.4. This opens the possibility to verify the Witten-Bodmer hypothesis on absolute stability of quark matter using ground-based experiments in which neutron-rich nuclei are employed. These experiments can also provide rather stringent constraints on the Equation of State (EoS) to be used for describing the pre-Supernova gravitational collapse. Consistent simulations of neutron rich heavy ion collisions are performed in order to show that even at relatively low energies, in the few AGeV range, the system can enter such unstable mixed phase. Some precursor observables are suggested, in particular a ``neutron trapping'' effect.

We discuss a model for long Gamma-Ray-Bursts in which the central engine is associated with the conversion process of a metastable hadronic star into a star containing quark matter. We analyze the observational signatures of the model, i.e. the Supernova-GRB temporal connection and the existence of long quiescent times in the temporal structure of Gamma-Ray-Bursts.

We compute the bulk viscosity of a mixed quark-hadron phase. In the first scenario to be discussed, the mixed phase occurs at large densities and we assume that it is composed of a mixing of hyperonic matter and quarks in the Color Flavor Locked phase. In a second scenario, the mixed phase occurs at lower densities and it is composed of a mixing of nucleons and unpaired quark matter. We have also investigated the effect of a non-vanishing surface tension at the interface between hadronic and quark matter. In both scenarios, the bulk viscosity is large when the surface tension is absent, while the value of the viscosity reduces in the second scenario when a finite value for the surface tension is taken into account. In all cases, the r-mode instabilities of the corresponding hybrid star are suppressed.

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.

We study the effects of the gapless color-flavor locked (gCFL) phase on the equation of state of strongly interacting matter in the range of baryonic chemical potential involved in a compact star. We analyze the possibility of a phase transition from hadronic matter to gCFL quark matter and we discuss, for different values of the strange quark mass and diquark coupling strength, the existence of a gCFL phase in quark or hybrid stars. The mass-radius relation and the structure of compact stars containing the gCFL phase are shown and the physical relevance of this superconducting phase inside a stellar object is also discussed.

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 show that if color superconducting quark matter forms in hybrid or quark stars it is possible to satisfy most of recent observational boundaries on masses and radii of compact stellar objects. An energy of the order of $10^{53}$ erg is released in the conversion from a (metastable) hadronic star into a (stable) hybrid or quark star in presence of a color superconducting phase. If the conversion occurs immediately after the deleptonization of the proto-neutron star, the released energy can help Supernovae to explode. If the conversion is delayed the energy released can power a Gamma Ray Burst. A delay between the Supernova and the subsequent Gamma Ray Burst is possible, in agreement with the delay proposed in recent analysis of astrophysical data.

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 study the effects of color superconductivity on the structure and formation of compact stars. We show that it is possible to satisfy most of recent observational boundaries on masses and radii if a diquark condensate forms in a hybrid or a quark star. Moreover, we find that a huge amount of energy, of the order of $10^{53}$ erg, can be released in the conversion from a (metastable) hadronic star into a (stable) hybrid or quark star, if the presence of a color superconducting phase is taken into account. Accordingly to the scenario proposed in Astrophys.J.586(2003)1250, the energy released in this conversion can power a Gamma Ray Burst. This mechanism can explain the recent observations indicating a delay, of the order of days or years, between a few Supernova explosions and the subsequent Gamma Ray Burst.

We propose a model to explain how a Gamma Rays Burst can take place days or years after a supernova explosion. Our model is based on the conversion of a pure hadronic star (neutron star) into a star made at least in part of deconfined quark matter. The conversion process can be delayed if the surface tension at the interface between hadronic and deconfined-quark-matter phases is taken into account. The nucleation time (i.e. the time to form a critical-size drop of quark matter) can be extremely long if the mass of the star is small. Via mass accretion the nucleation time can be dramaticaly reduced and the star is finally converted into the stable configuration. A huge amount of energy, of the order of 10$^{52}$--10$^{53}$ erg, is released during the conversion process and can produce a powerful Gamma Ray Burst. The delay between the supernova explosion generating the metastable neutron star and the new collapse can explain the delay proposed in GRB990705 and in GRB011211.

We explore the dependence of the critical density, separating hadronic matter from a mixed phase of quarks and hadrons, on the ratio $Z/A$. We use both the MIT bag model and the Color Dielectric Model to describe the quark dynamics, while for the hadronic phase we employ various relativistic equations of state. We find that, if the parameters of quark models are fixed so that the existence of quark stars is allowed, then the critical density drops dramatically in the range $Z/A \sim $ 0.3--0.4. Moreover, for $Z/A \sim $ 0.3 the critical density is only slightly larger than the saturation density of symmetric nuclear-matter. This opens the possibility to verify the Witten-Bodmer hypothesis on absolute stability of quark matter using ground-based experiments in which neutron-rich nuclei are tested.

We show the possible existence of compact stars having a surface composed of a mixed phase of quarks and hadrons. This scenario can be realized both for self-bound stars, satisfying the so-called Witten-Bodmer hypothesis, and for gravitationally bound stars. This class of solutions of the Tolman-Oppenheimer-Volkoff equation can be obtained in all the models we discuss, within a physically acceptable range of values of the model parameters.

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.

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.