The total binding energy of compact stars is the sum of the gravitational binding energy $(BE)_g$ and the nuclear binding energy $(BE)_n$, the last being related to the microphysics of the interactions. While the first is positive (binding) both for hadronic stars and for strange quark stars, the second is large and negative for hadronic stars (anti-binding) and either small and negative (anti-binding) or positive (binding) for strange quark stars. A hadronic star can convert into a strange quark star with a larger radius because the consequent reduction of $(BE)_g$ is over-compensated by the large increase in $(BE)_n$. Thus, the total binding energy increases due to the conversion and the process is exothermic. Depending on the equations of state of hadronic matter and quark matter and on the baryonic mass of the star, the contrary is obviously also possible, namely the conversion of hadronic stars into strange quark stars having smaller radii, a situation more often discussed in the literature. We provide a condition that is sufficient and in most of the phenomenologically relevant cases also necessary in order to form strange quark stars with larger radii while satisfying the exothermicity request. Finally, we compare the two schemes in which quark stars are produced (one having large quark stars and the other having small quark stars) among themselves and with the third-family scenario and we discuss how present and future data can discriminate among them.

The detection of GW170817 and its electromagnetic counterparts allows to constrain the equation of state of dense matter in new and complementary ways. Very stiff equations of state are ruled out by the upper limit on the average tidal deformability, $\tilde\Lambda \lesssim 800$, imposed by the detected gravitational wave signal. A lower limit, $\tilde\Lambda \gtrsim 400$, can also be extracted by considering the large amount of ejected matter which powers the kilonova AT2017gfo. By using several microscopic nucleonic equations of state, we first confirm the existence of a monotonic relation between $R_{1.5}$ (the radius of the $1.5M_{\odot}$ configuration) and $\tilde\Lambda$. This translates the limits on $\tilde\Lambda$ into limits on the radius: $11.8\,\text{km} \lesssim R_{1.5} \lesssim 13.1\,$km. We then show that the monotonic relation is violated, if a second branch of compact stars composed of quark matter exists, as in the two-families or the twin-stars scenarios. In particular, it is possible to fulfill the limits on $\tilde\Lambda$ while having $R_{1.5}$ significantly smaller than $12\,$km. In both those scenarios the event GW170817/AT2017gfo originates from the merger of a hadronic star and a star containing quark matter.

We study the partial fragmentation of a strange quark star into strangelets during the process of the merger of two strange quark stars or of a strange quark star with a neutron star. We discuss the fate of the fragments considering their possible evaporation into nucleons. We show that only the small amount of large size strangelets produced by tidal forces can survive evaporation, while the strangelets formed during and after the merger evaporate completely into nucleons. In this way we show that: 1) the density of strangelets in the galaxy is too low to be able to affect significantly the evolution of stars and the probability of direct detection of a strangelet is negligible; 2) the matter ejected by the strange quark star during and after the merger evaporates into nucleons and therefore can contribute to the kilonova signal.

We shortly summarize the two-families scenario in which both hadronic stars and strange quark stars can exist and we describe the main predictions one can obtain from it. We then concentrate on the observables that most likely will be measured in the near future, i.e. masses, radii, tidal deformabilities and moments of inertia and we present a list of objects that are candidate strange quark stars in this scheme. We show that the estimates of the radii derived up to now from observations are all compatible with the two-families scenario and in particular all the objects having large radii can easily be interpreted as strange quark stars.

We analyse the phenomenological implications of the two-families scenario on the merger of compact stars. That scenario is based on the coexistence of both hadronic stars and strange quark stars. After discussing the classification of the possible mergers, we turn to detailed numerical simulations of the merger of two hadronic stars, i.e., "first family" stars in which delta resonances and hyperons are present, and we show results for the threshold mass of such binaries, for the mass dynamically ejected and the mass of the disk surrounding the post-merger object. We compare these results with those obtained within the one-family scenario and we conclude that relevant signatures of the two-families scenario can be suggested, in particular: the possibility of a rapid collapse to a black hole for masses even smaller than the ones associated to GW170817; during the first milliseconds, oscillations of the postmerger remnant at frequencies higher than the ones obtained in the one-family scenario; a large value of the mass dynamically ejected and a small mass of the disk, for binaries of low total mass. Finally, based on a population synthesis analysis, we present estimates of the number of mergers for: two hadronic stars; hadronic star - strange quark star; two strange quark stars. We show that for unequal mass systems and intermediate values of the total mass, the merger of a hadronic star and a strange quark star is very likely (GW170817 has a possible interpretation into this category of mergers). On the other hand, mergers of two strange quark stars are strongly suppressed.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

We discuss the different signals, in gravitational and electromagnetic waves, emitted during the merger of two compact stars. We will focus in particular on the possible contraints that those signals can provide on the equation of state of dense matter. Indeed, the stiffness of the equation of state and the particle composition of the merging compact stars, strongly affect e.g. the life time of the post-merger remnant and its gravitational wave signal, the emission of the short gamma-ray-burst, the amount of ejected mass and the related kilonova. The first detection of gravitational waves from the merger of two compact stars in August 2017, GW170817, and the subsequent detections of its electromagnetic counterparts, GRB170817A and AT2017gfo, is the first example of the era of "multi-messenger astronomy": we discuss what we have learned from this detection on the equation of state of compact stars and we provide a tentative interpretation of this event, within the two families scenario, as due to the merger of a hadronic star with a quark star.

The merger of two neutron stars is a very complex process. In order to disentangle the various steps through which it takes place it is mandatory to examine all the signals we can detect: gravitational waves and electromagnetic waves, in a huge spectrum ranging from X and $\gamma$-rays down to infrared and radio. Each of these signals provides a message and the totality of this information will allow us not only to understand the process of the merger but also the behavior of matter at those extreme conditions.

If only one family of "neutron stars" exists, their maximum mass must be equal or larger than $2 M_\odot$ and then only in less than about $18\%$ of cases the outcome of the merger of two neutron stars is a direct collapse to a black hole, since the newly formed system can avoid the collapse at least until differential rotation is present. In the so-called two-families scenario, stars made of hadrons are stable only up to about $1.5-1.6 M_\odot$, while the most massive compact stars are entirely made of strange quark matter. We show that in this scenario the outcome of the merger of two neutron stars is a prompt collapse in at least $34\%$ of the cases. It will therefore be easy to discriminate between the two scenarios once the gravitational waves emitted at the moment of the merger are detected.

The presence of exotic states of matter in neutron stars (NSs) is currently an open issue in physics. The appearance of muons, kaons, hyperons, and other exotic particles in the inner regions of the NS, favored by energetic considerations, is considered to be an effective mechanism to soften the equation of state (EoS). In the so-called two-families scenario, the softening of the EoS allows for NSs characterized by very small radii, which become unstable and convert into a quark stars (QSs). In the process of conversion of a NS into a QS material can be ablated by neutrinos from the surface of the star. Not only neutron-rich nuclei, but also more exotic material, such as hypernuclei or deconfined quarks, could be ejected into the atmosphere. In the NS atmosphere, atoms like H, He, and C should exist, and attempts to model the NS thermal emission taking into account their presence, with spectra modified by the extreme magnetic fields, have been done. However, exotic atoms, like muonic hydrogen $(p\,\mu^-)$ or the so-called Sigmium $(\Sigma^+\,e^-)$, could also be present during the conversion process or in its immediate aftermath. At present, analytical expressions of the wave functions and eigenvalues for these atoms have been calculated only for H. In this work, we extend the existing solutions and parametrizations to the exotic atoms $(p\,\mu^-)$ and $(\Sigma^+\,e^-)$, making some predictions on possible transitions. Their detection in the spectra of NS would provide experimental evidence for the existence of hyperons in the interior of these stars.

The existence of strange quark stars has been proposed many years ago. More recently, the possible co-existence of a first family composed of "normal" neutron stars with a second family of strange quark stars has been proposed as a solution of problems related to the maximum mass and to the minimal radius of these compact stellar objects. In this paper we study the mass distribution of compact objects formed in binary systems and the relative fractions of quark and neutron stars in different subpopulations. We incorporate the strange quark star formation model provided by the two-families scenario and we perform a large-scale population synthesis study in order to obtain the population characteristics. In our model, below a critical gravitational mass $M_\mathrm{max}^H- \Delta M$ only normal (hadron) neutron stars exist. Then in the mass range $(M_\mathrm{max}^H- \Delta M) \leqslant M \leqslant M_\mathrm{max}^H$ strange quark stars and neutron stars coexist. Finally, above $M_\mathrm{max}^H$ all compact objects are strange quark stars. We argue that $M_\mathrm{max}^H$ is in the range $\sim 1.5-1.6 M_\odot$. According to our results, the main channel for strange quark star formation in binary systems is accretion from a secondary companion on a neutron star.This opens the possibility of having explosive GRB-like phenomena not related to supernovae and not due to the merger of two neutron stars. The enhancement in the number of compact objects in the co-existence mass range $(M_\mathrm{max}^H- \Delta M) \leqslant M \leqslant M_\mathrm{max}^H$ is not very pronounced. The number of double strange quark star's systems is rather small with only a tiny fraction which merge within a Hubble time. This drastically limits the flux of strangelets produced by the merger, which turns out to be compatible with all limits stemming from Earth and lunar experiments.

We discuss the possible existence of absolutely stable strange quark matter within three different types of chiral models. We will show that confinement plays a crucial role in determining the conditions for the Bodmer-Witten hypothesis to hold true. We discuss also which are the phenomenological signatures, related to measurements of masses and radii of compact stars, which would prove the existence of strange quark stars.

We investigate the possible implications of quark deconfinement on the phenomenology of Long Gamma-Ray Bursts focusing, in particular, on the possibility to describe multiple prompt emission phases in the context of the proto-magnetar model. Starting from numerical models of rotating Hadron Stars and Quark Stars in full general relativity we track the electromagnetic spin-down evolution in both the hadronic and quark phase, linking the two families through conservation of baryon number and angular momentum. We give estimates of the timescales and the energetics involved in the spin-down process deriving, in the relevant spin range, the relation between the initial and the final masses and rotational energies, whenever hadron-quark conversion is possible. We show how the results can be used in relevant astrophysical cases such as the double burst GRB 110709B.

We will follow the two-families scenario described in the accompanying paper, in which compact stars having a very small radius and masses not exceeding about 1.5$M_\odot$ are made of hadrons, while more massive compact stars are quark stars. In the present paper we discuss the dynamics of the transition of a hadronic star into a quark star. We will show that the transition takes place in two phases: a very rapid one, lasting a few milliseconds, during which the central region of the star converts into quark matter and the process of conversion is accelerated by the existence of strong hydrodynamical instabilities, and a second phase, lasting about ten seconds, during which the process of conversion proceeds till the surface of the star via production and diffusion of strangeness. We will show that these two steps play a crucial role in the phenomenological implications of the model. We will discuss the possible implications of this scenario both for long and for short Gamma Ray Bursts, using the proto-magnetar model as the reference frame of our discussion. We will show that the process of quark deconfinement can be connected to specific observed features of the GRBs. In the case of long GRBs we will discuss the possibility that quark deconfinement is at the origin of the second peak present in quite a large fraction of bursts. Also we will discuss the possibility that long GRBs can take place in binary systems without being associated with a SN explosion. Concerning short GRBs, quark deconfinement can play the crucial role in limiting their duration. Finally we will shortly revisit the possible relevance of quark deconfinement in some specific type of Supernova explosions, in particular in the case of very massive progenitors.

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 argue that the full conversion of a hadronic star into a quark or a hybrid star occurs within two different regimes separated by a critical value of the density of the hadronic phase $\overline{n_h}$. The first stage, occurring for $n_h>\overline{n_h}$, is characterized by turbulent combustion and lasts typically a few ms. During this short time-scale neutrino cooling is basically inactive and the star heats up thanks to the heat released in the conversion. In the second stage, occurring for $n_h<\overline{n_h}$, turbulence is not active anymore, and the conversion proceeds on a much longer time scale (of the order of tens of seconds), with a velocity regulated by the diffusion and the production of strange quarks. At the same time, neutrino cooling is also active. The interplay between the heating of the star due to the slow conversion of its outer layers (with densities smaller than $\overline{n_h}$) and the neutrino cooling of the forming quark star leads to a quasi-plateau in the neutrino luminosity which, if observed, would possibly represent a unique signature for the existence of quark matter inside compact stars. We will discuss the phenomenological implications of this scenario in particular in connection with the time structure of long gamma-ray-bursts.

We will review the approach used for studying the conversion of a hadronic star into a quark star based on the assumption of a infinitely thin combustion zone and we will discuss why, in this scheme, the combustion stops before the whole hadronic star is converted.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

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 recent measurement of two solar mass pulsars has initiated an intense discussion on its impact on our understanding of the high-density matter in the cores of neutron stars. A task force meeting was held from October 7-10, 2013 at the Frankfurt Institute for Advanced Studies to address the presence of quark matter in these massive stars. During this meeting, the recent oservational astrophysical data and heavy-ion data was reviewed. The possibility of pure quark stars, hybrid stars and the nature of the QCD phase transition were discussed and their observational signals delineated.

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 present ECHO-QGP, a numerical code for $(3+1)$-dimensional relativistic viscous hydrodynamics designed for the modeling of the space-time evolution of the matter created in high energy nuclear collisions. The code has been built on top of the \emph{Eulerian Conservative High-Order} astrophysical code for general relativistic magneto-hydrodynamics [\emph{Del Zanna et al., Astron. Astrophys. 473, 11, 2007}] and here it has been upgraded to handle the physics of the Quark-Gluon Plasma. ECHO-QGP features second-order treatment of causal relativistic viscosity effects in both Minkowskian or Bjorken coordinates; partial or complete chemical equilibrium of hadronic species before kinetic freeze-out; initial conditions based on the optical Glauber model, including a Monte-Carlo routine for event-by-event fluctuating initial conditions; a freeze-out procedure based on the Cooper-Frye prescription. The code is extensively validated against several test problems and results always appear accurate, as guaranteed by the combination of the conservative (shock-capturing) approach and the high-order methods employed. ECHO-QGP can be extended to include evolution of the electromagnetic fields coupled to the plasma.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

R-modes can generate strong magnetic fields in the core of accreting millisecond neutron stars (NSs). The diffusion of these fields outside the core causes the growth of the external magnetic field and thus it affects the evolution of the spin down rates $\dot{P}$ of the millisecond pulsars (MSPs). The diffusion of the internal magnetic field provides a new evolutionary path for the MSPs. This scenario could explain the large $\dot{P}$ of the pulsars J1823-3021A and J1824-2452A.

Differential rotation induced by the r-mode instability can generate very strong toroidal fields in the core of accreting, millisecond spinning neutron stars. We introduce explicitly the magnetic damping term in the evolution equations of the r-modes and solve them numerically in the Newtonian limit, to follow the development and growth of the internal magnetic field. We show that the strength of the latter can reach large values, $B \sim 10^{14} $ G, in the core of the fastest accreting neutron stars. This is strong enough to induce a significant quadrupole moment of the neutron star mass distribution, corresponding to an ellipticity $|\epsilon_B}| \sim 10^{-8}$. If the symmetry axis of the induced magnetic field is not aligned with the spin axis, the neutron star radiates gravitational waves. We suggest that this mechanism may explain the upper limit of the spin frequencies observed in accreting neutron stars in Low Mass X-Ray Binaries. We discuss the relevance of our results for the search of gravitational waves.

The OPERA Collaboration reported evidence for muonic neutrinos traveling slightly faster than light in vacuum. While waiting further checks from the experimental community, here we aim at exploring some theoretical consequences of the hypothesis that muonic neutrinos are superluminal, considering in particular the tachyonic and the Coleman-Glashow cases. We show that a tachyonic interpretation is not only hardly reconciled with OPERA data on energy dependence, but that it clashes with neutrino production from pion and with neutrino oscillations. A Coleman-Glashow superluminal neutrino beam would also have problems with pion decay kinematics for the OPERA setup; it could be easily reconciled with SN1987a data, but then it would be very problematic to account for neutrino oscillations.

We show that the r-mode instability can generate strong toroidal fields in the core of accreting millisecond quark stars by inducing differential rotation. We follow the spin frequency evolution on a long time scale taking into account the magnetic damping rate in the evolution equations of r-modes. The maximum spin frequency of the star is only marginally smaller than in the absence of the magnetic field. The late-time evolution of the stars which enter the r-mode instability region is instead rather different if the generated magnetic fields are taken into account: they leave the millisecond pulsar region and they become radio pulsars.

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.

In rotating neutron stars the existence of the Coriolis force allows the presence of the so-called Rossby oscillations (r-modes) which are know to be unstable to emission of gravitational waves. Here, for the first time, we introduce the magnetic damping rate in the evolution equations of r-modes. We show that r-modes can generate very strong toroidal fields in the core of accreting millisecond pulsars by inducing differential rotation. We shortly discuss the instabilities of the generated magnetic field and its long time-scale evolution in order to clarify how the generated magnetic field can stabilize the star.

In millisecond pulsars the existence of the Coriolis force allows the development of the so-called Rossby oscillations (r-modes) which are know to be unstable to emission of gravitational waves. These instabilities are mainly damped by the viscosity of the star or by the existence of a strong magnetic field. A fraction of the observed millisecond pulsars are known to be inside Low Mass X-ray Binaries (LMXBs), systems in which a neutron star (or a black hole) is accreting from a donor whose mass is smaller than 1 $M_\odot$. Here we show that the r-mode instabilities can generate strong toroidal magnetic fields by inducing differential rotation. In this way we also provide an alternative scenario for the origin of the magnetars.

We study matter at high density and temperature using a chiral lagrangian in which the breaking of scale invariance is regulated by the value of a scalar field, called dilaton \cite{Heide:1993yz,Carter:1995zi,Carter:1996rf,Carter:1997fn}. We provide a phase diagram describing the restoration of chiral and scale symmetries. We show that chiral symmetry is restored at large temperatures, but at low temperatures it remains broken at all densities. We also show that scale invariance is more easily restored at low rather than large baryon densities. The masses of vector mesons scale with the value of the dilaton and their values initially slightly decrease with the density but then they increase again for densities larger than $\sim 3 \rho_0$. The pion mass increases continuously with the density and at $\rho_0$ and T=0 its value is $\sim$ 30 MeV larger than in the vacuum. We show that the model is compatible with the bounds stemming from astrophysics, as e.g. the one associated with the maximum mass of a neutron star. The most striking feature of the model is a very significant softening at large densities, which manifests also as a strong reduction of the adiabatic index. While the softening has probably no consequence for Supernova explosion via the direct mechanism, it could modify the signal in gravitational waves associated with the merging of two neutron stars.

We describe possible scenarios of quark deconfinement in compact stars and we analyze their astrophysical implications. The quark deconfinement process can proceed rapidly, as a strong deflagration, releasing a huge amount of energy in a short time and generating an extra neutrino burst. If energy is transferred efficiently to the surface, like e.g. in the presence of convective instabilities, this burst could contribute to revitalize a partially failed SN explosion. We discuss how the neutrino observations from SN1987A would fit in this scenario. Finally, we focus on the fate of massive and rapidly rotating progenitors, discussing possible time separations between the moment of the core collapse and the moment of quark deconfinement. This mechanism can be at the basis of the interpretation of gamma ray bursts in which lines associated with heavy elements are present in the spectrum.

The formation of a strange or hybrid star from a neutron star progenitor is believed to occur when the central stellar density exceeds a critical value. If the transition from hadron to quark matter is of first order, the event has to release a huge amount of energy in a very short time and we would be able to observe the phenomenon even if it is at cosmological distance far from us; most likely, such violent quark deconfinement would be associated with at least a fraction of the observed gamma ray bursts. If we allow for temporal variations of fundamental constants like $\Lambda_{QCD}$ or $G_N$, we can expect that neutron stars with an initial central density just below the critical value can enter into the region where strange or hybrid stars are the true ground state. From the observed rate of long gamma ray bursts, we are able to deduce the constraint $\dot{G}_N/G_N \lesssim 10^{-17} {\rm yr^{-1}}$, which is about 5 orders of magnitude more stringent than the strongest previous bounds on a possible increasing $G_N$.

We show that r-mode instabilities severely constraint the composition of a compact star rotating at a submillisecond period. In particular, the only viable astrophysical scenario for such an object, present inside the Low Mass X-ray Binary associated with the x-ray transient XTE J1739-285, is that it has a strangeness content. Since previous analysis indicate that hyperonic stars or stars containing a kaon condensate are not good candidates, the only remaining possibility is that such an object is either a strange quark star or a hybrid quark-hadron star. We also discuss under which conditions sub-millisecond pulsars are rare.

Ozel, in a recent reanalysis of EXO 0748-676 observational data (astro-ph/0605106), concluded that quark matter probably does not exist in the center of compact stars. We show that the data is actually consistent with the presence of quark matter in compact stars.

We analyse the occurrence of quiescent times in the temporal structure of the Gamma-Ray-Bursts (GRBs) light curves. We show that if a long quiescent time is present, it is possible to divide the total duration of GRBs into three periods: the pre-quiescence emission, the quiescent time and the post-quiescence emission. We then discuss a model of the GRBs inner engine based on the formation of quark phases during the life of an hadronic star. Within this model the pre-quiescence emission is interpreted as due to the deconfinement of quark inside an hadronic star and the formation of 2SC quark matter. The post-quiescence emission is due to the conversion of 2SC into the Color-Flavor-Locking (CFL) phase. The temporal delay between these two processes is connected with the nucleation time of the CFL phase in the 2SC phase and it can be associated with the observed quiescent times in the GRBs light curves. The stability of CFL cores in compact stars is also discussed.

The energy released by the inner engine of GRBs can originate from structural readjustments inside a compact star. In particular, the formation of deconfined quark matter can liberate enough energy to power the burst. We show that the burning of a neutron star into a quark star likely proceeds as a deflagration and not as a detonation. In that way no strong baryon contamination is produced near the surface of the star. It is tempting to associate the temporal structures observed in the light curves with specific processes taking place inside the compact star. The so-called quiescent times, during which no signal is emitted in the highest energy band, correspond to pauses during the processes of readjustment. If the quark (or hybrid) star formed after these transformations is strongly magnetized and rotates rapidly, a prolonged gamma emission can be produced, as proposed by Usov years ago. This can explain the quasi-plateau observed by Swift in several GRBs.

We perform a statistical analysis of the temporal structure of long Gamma-Ray-Bursts (GRBs). First we consider a sample of bursts in which a long quiescent time is present. Comparing the pre-quiescence with the post-quiescence emission we show that they display similar temporal structures, hardness ratios and emitted powers, but, on the average, the post-quiescence emission is roughly twice as long as the pre-quiescence emission. We then consider a sample of long and bright GRBs. We show that the duration of each emission period is compatible with the duration of an active period computed in various inner engine models. At the contrary, if the inner engine is assumed to be always active, i.e. also during the quiescent times, in several cases the total duration of the burst largely exceeds the theoretical durations. Our analysis therefore does not support the interpretation of long quiescent times in terms of stochastic modulation of a continuous wind. Instead the quiescent times can be interpreted as dormancy periods of the inner engine. Before and after a dormancy period the inner engine produces similar emissions.

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.

The temporal structure of Gamma Ray Bursts can be interpreted assuming as a inner engine a neutron star which undergoes a progressive compactification via production of strangeness (hyperons and kaons) and quarks. We will propose a tentative identification of various emission periods of the burst with specific structural changes of the star. Each of these modifications of the composition of the compact star takes place as a deflagration and not as a detonation, so the energy released in the transition goes mainly into heat and not into a mechanical wave. This is important in order to avoid an excessive baryonic contamination of the region surrounding the compact star. In this way a ultrarelativistic plasma of electron-positron pairs and of photons can be obtained, powering the Gamma Ray Burst.

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 show that sudden variations in the composition and structure of an hybrid star can be triggered by its rapid spin-down, induced by r-mode instabilities. The discontinuity of this process is due to the surface tension between hadronic and quark matter and in particular to the overpressure needed to nucleate new structures of quark matter in the mixed phase. The consequent mini-collapses in the star can produce highly energetic gravitational wave bursts. The possible connection between the predictions of this model and the burst signal found by EXPLORER and NAUTILUS detectors during the year 2001 is also investigated.

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.

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

The Pseudo-Goldstone Boson (PGB) emission could provide a very efficient mechanism for explaining the cosmic Gamma Ray Bursts (GRBs). The PGBs could be produced during the merging of two compact objects like two neutron stars or neutron star - black hole, and then decay into electron-positron pairs and photons at distances ~ 100-1000 km. In this way, a huge energy (up to more than 10^{54} erg) can be deposited into the outer space with low baryon density in the form of ultrarelativistic e^+e^- plasma, the so called fireball, which originates the observed gamma-ray bursts. The needed ranges for the PGB parameters are: mass of order MeV, coupling to nucleons g_{aN} ~ few x 10^{-6} and to electrons g_{ae} ~ few x 10^{-9}. Interestingly enough, the range for coupling constants correspond to that of the invisible axion with the Peccei-Quinn symmetry breaking scale V ~ few x 10^5 GeV, but the mass of the PGB is many orders of magnitude larger than what such a scale would demand to an axion. Neither present experimental data nor astrophysical and cosmological arguments can exclude such an ultramassive axion, however the relevant parameters' window is within the reach of future experiments. Another exciting point is that our mechanism could explain the association of some GRBs with supernovae type Ib/c, as far as their progenitor stars have a radius ~ 10^4 km. And finally, it also could help the supernova type II explosion: PGBs emitted from the core of the collapsing star and decaying in the outer shells would deposit a kinetic energy of the order of 10^{51} erg. In this way, emission of such an axion-like particle could provide an unique theoretical base for understanding the gamma ray bursts and supernova explosions.

The transition to quark matter can take place in neutron stars. The structure of a hybrid star, containing a core made of quark matter is discussed. The maximum mass of the non-rotating hybrid star turns out to be 1.6 M_s. Possible signatures of the quark phase as pulsar's timing, cooling rate etc. are briefly analyzed. The deconfinement transition can also take place during the pre-supernova collapse. This possibility is studied by introducing a finite temperature EOS. The dependence of the latter on the proton fraction is shown to be crucial. The softening of the EOS at densities just above nuclear matter saturation density for Z/A = 0.3 helps in obtaining an explosion. At the same time, at larger densities the EOS is stiff enough to support a neutron star compatible with observations.

We study the equation of state of quark matter at finite temperature, using a confinement model in which chiral symmetry remains broken in the deconfined phase. Implications for type II supernova explosion and for the structure and evolution of the proto-neutron star are discussed.

We study the equation of state (EOS) of quark matter at zero temperature, using the Color Dielectric Model (CDM) to describe confinement. Sensible results are obtained in the version of the CDM for which confinement is imposed {\it smoothly}. The two--phases version of the model turns out to give unrealistic results for the EOS. Chiral symmetry plays a marginal r\^ole and the quarks are massive till high densities. The deconfinement phase transition is smooth and unlikely to be first order. Instabilities of the quark matter and the gap equation are discussed.