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

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.