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.

We present a Bayesian analysis to constrain the equation of state of dense nucleonic matter by exploiting the available data from symmetric nuclear matter at saturation and from observations of compact X-ray sources and from the gravitational wave event GW170817. For the first time, such analysis is performed by using a class of models, the relativistic mean field models, which allow to consistently construct an equation of state in a wide range of densities, isospin asymmetries and temperatures. The selected class of models contains five nuclear physics empirical parameters at saturation for which we construct the joint posterior distributions. By exploring different types of priors, we find that the equations of state with the largest evidence are the ones featuring a strong reduction of the effective mass of the nucleons in dense matter which can be interpreted as an indication of a phase transition to a chiral symmetry restored phase. Those equations of state in turn predict $R_{1.4} \sim 12$ km. Finally, we present a preliminary investigation on the effect of including $\Lambda$ hyperons showing that they appear in stars more massive than about $1.6 M_{\odot}$ and lead to radii larger than about $R_{1.4} \sim 14$ km. Within the model here explored, the formation of such particles provide a poor agreement with the constraints from GW170817.

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.

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

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

The possible existence of two families of compact stars, neutron stars and quark stars, naturally leads to a scenario in which a conversion process between the two stellar objects occurs with a consequent release of energy of the order of $10^{53}$ erg. We discuss recent hydrodynamical simulations of the burning process and neutrino diffusion simulations of cooling of a newly formed strange star. We also briefly discuss this scenario in connection with recent measurements of masses and radii of compact stars.

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.

There are strong indications that the process of conversion of a neutron star into a strange quark star proceeds as a strong deflagration implying that in a few milliseconds almost the whole star is converted. Starting from the three-dimensional hydrodynamic simulations of the combustion process which provide the temperature profiles inside the newly born strange star, we calculate for the first time the neutrino signal that is to be expected if such a conversion process takes place. The neutrino emission is characterized by a luminosity and a duration that is typical for the signal expected from protoneutron stars and represents therefore a powerful source of neutrinos which could be possibly directly detected in case of events occurring close to our Galaxy. We discuss moreover possible connections between the birth of strange stars and explosive phenomena such as supernovae and gamma-ray-bursts.

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.

We discuss the possible impact of strange quark matter on the evolution of core-collapse supernovae with emphasis on low critical densities for the quark-hadron phase transition. For such cases the hot proto-neutron star can collapse to a more compact hybrid star configuration hundreds of milliseconds after core-bounce. The collapse triggers the formation of a second shock wave. The latter leads to a successful supernova explosion and leaves an imprint on the neutrino signal. These dynamical features are discussed with respect to their compatibility with recent neutron star mass measurements which indicate a stiff high density nuclear matter equation of state.

We compare two classes of hybrid equations of state with a hadron-to-quark matter phase transition in their application to core collapse supernova simulations. The first one uses the quark bag model and describes the transition to three-flavor quark matter at low critical densities. The second one employs a Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with parameters describing a phase transition to two-flavor quark matter at higher critical densities. These models possess a distinctly different temperature dependence of their transition densities which turns out to be crucial for the possible appearance of quark matter in supernova cores. During the early post bounce accretion phase quark matter is found only if the phase transition takes place at sufficiently low densities as in the study based on the bag model. The increase critical density with increasing temperature, as obtained for our PNJL parametrization, prevents the formation of quark matter. The further evolution of the core collapse supernova as obtained applying the quark bag model leads to a structural reconfiguration of the central proto-neutron star where, in addition to a massive pure quark matter core, a strong hydrodynamic shock wave forms and a second neutrino burst is released during the shock propagation across the neutrinospheres. We discuss the severe constraints in the freedom of choice of quark matter models and their parametrization due to the recently observed 2 solar mass pulsar and their implications for further studies of core collapse supernovae in the QCD phase diagram.

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.

The recent observation of the pulsar PSR J1614-2230 with a mass of 1.97 +/- 0.04 M_sun gives a strong constraint on the quark and nuclear matter equations of state (EoS). We explore the parameter ranges for a parameterized EoS for quark stars. We find that strange stars, made of absolutely stable strange quark matter, comply with the new constraint only if effects from the strong coupling constant and color-superconductivity are taken into account. Hybrid stars, compact stars with a quark matter core and an hadronic outer layer, can be as massive as 2 M_sun, but only for a significantly limited range of parameters. We demonstrate that the appearance of quark matter in massive stars depends crucially on the stiffness of the nuclear matter EoS. We show that the masses of hybrid stars stay below the ones of hadronic and pure quark stars, due to the softening of the EoS at the quark-hadron phase transition.

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.

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.

The implications of a QCD phase transition at high temperatures and densities for core-collapse supernovae are discussed. For a strong first order phase transition to quark matter, various scenarios have been put forward in the literature. Here, detailed numerical simulations including neutrino transport are presented, where it is found that a second shock wave due to the QCD phase transition emerges shortly after bounce. It is demonstrated that such a supernova banging twice results in a second peak in the antineutrino spectrum. This second peak is clearly detectable in present neutrino detectors for a galactic supernova.

For finite chemical potential effective models of QCD predict a first order phase transition. In favour for the search of such a phase transition in nature, we construct an equation of state for strange quark matter based on the MIT bag model. We apply this equation of state to highly asymmetric core collapse supernova matter with finite temperatures and large baryon densities. The phase transition is constructed using the general Gibbs conditions, which results in an extended coexistence region between the pure hadronic and pure quark phases in the phase diagram, i.e. the mixed phase. The supernovae are simulated via general relativistic radiation hydrodynamics based on three flavor Boltzmann neutrino transport in spherical symmetry. During the dynamical evolution temperatures above 10 MeV, baryon densities above nuclear saturation density and a proton-to-baryon ratio below 0.2 are obtained. At these conditions the phase transition is triggered which leads to a significant softening of the EoS for matter in the mixed phase. As a direct consequence of the stiffening of the EoS again for matter in the pure quark phase, a shock wave forms at the boundary between the mixed and the pure hadronic phases. This shock is accelerated and propagates outward which releases a burst of neutrinos dominated by electron anti-neutrinos due to the lifted degeneracy of the shock-heated hadronic material. We discuss the radiation-hydrodynamic evolution of the phase transition at the example of several low and intermediate mass Fe-core progenitor stars and illustrate the expected neutrino signal from the phase transition.

The phase transition from hadronic to quark matter may take place already during the early post-bounce stage of core collapse supernovae when matter is still hot and lepton rich. If the phase transition is of first order and exhibits a barrier, the formation of the new phase occurs via the nucleation of droplets. We investigate the thermal nucleation of a quark phase in supernova matter and calculate its rate for a wide range of physical parameters. We show that the formation of the first droplet of a quark phase might be very fast and therefore the phase transition to quark matter could play an important role in the mechanism and dynamics of supernova explosions.

Explosive astrophysical systems, such as supernovae or compact star binary mergers, provide conditions where strange quark matter can appear. The high degree of isospin asymmetry and temperatures of several MeV in such systems may cause a transition to the quark phase already around saturation density. Observable signals from the appearance of quark matter can be predicted and studied in astrophysical simulations. As input in such simulations, an equation of state with an integrated quark matter phase transition for a large temperature, density and proton fraction range is required. Additionally, restrictions from heavy ion data and pulsar observation must be considered. In this work we present such an approach. We implement a quark matter phase transition in a hadronic equation of state widely used for astrophysical simulations and discuss its compatibility with heavy ion collisions and pulsar data. Furthermore, we review the recently studied implications of the QCD phase transition during the early post-bounce evolution of core-collapse supernovae and introduce the effects from strong interactions to increase the maximum mass of hybrid stars. In the MIT bag model, together with the strange quark mass and the bag constant, the strong coupling constant $\alpha_s$ provides a parameter to set the beginning and extension of the quark phase and with this the mass and radius of hybrid stars.

Some recent developments concerning the role of strange quark matter for astrophysical systems and the QCD phase transition in the early universe are addressed. Causality constraints of the soft nuclear equation of state as extracted from subthreshold kaon production in heavy-ion collisions are used to derive an upper mass limit for compact stars. The interplay between the viscosity of strange quark matter and the gravitational wave emission from rotation-powered pulsars are outlined. The flux of strange quark matter nuggets in cosmic rays is put in perspective with a detailed numerical investigation of the merger of two strange stars. Finally, we discuss a novel scenario for the QCD phase transition in the early universe, which allows for a small inflationary period due to a pronounced first order phase transition at large baryochemical potential.

We consider the possible formation of the quark hadron mixed phase in protoneutron stars. We discuss two cases: the first one, corresponding to a vanishingly small value of the surface tension of quark matter, is the well known mixed phase in which the global electric charge neutrality condition is imposed. In turn, this produces a non-constant pressure mixed phase. In the second case, corresponding to very large values of the surface tension, the charge neutrality condition holds only locally. However, the existence in protoneutron star matter of an additional globally conserved charge, the lepton number, allows for a new type of non-constant pressure mixed phase. We discuss the properties of the new mixed phase and the possible effects of its formation during the evolution of protoneutron stars.

The phase transition from hadronic matter to quark matter at high density might be a strong first order phase transition in presence of a large surface tension between the two phases. While this implies a constant-pressure mixed phase for cold and catalyzed matter this is not the case for the hot and lepton rich matter formed in a protoneutron star. We show that it is possible to obtain a mixed phase with non-constant pressure by considering the global conservation of lepton number during the stage of neutrino trapping. In turn, it allows for the appearance of a new kind of mixed phase as long as neutrinos are trapped and its gradual disappearance during deleptonization. This new mixed phase, being composed by two electric neutral phases, does not develop a Coulomb lattice and it is formed only by spherical structures, drops and bubbles, which can have macroscopic sizes. The disappearance of the mixed phase at the end of deleptonization might lead to a delayed collapse of the star into a more compact configuration containing a core of pure quark phase. In this scenario, a significant emission of neutrinos and, possibly, gravitational waves are expected.

We investigate the qualitative properties of phase transitions in a general way, if not the single particle numbers of the system but only some particular charges like e.g. baryon number are conserved. In addition to globally conserved charges we analyze the implications of locally conserved charge fractions, like e.g. local electric charge neutrality or locally fixed proton or lepton fractions. The conditions for phase equilibrium are derived and it is shown, that the properties of the phase transition do not depend on the locally conserved fractions. Finally, the general formalism is applied to the liquid-gas phase transition and the hadron-quark phase transition for typical astrophysical environments like in supernovae, proto-neutron or a neutron stars. We demonstrate that the Maxwell construction known from cold-deleptonized neutron star matter with two locally charge neutral phases requires modifications and further assumptions concerning the applicability for hot lepton-rich matter. All possible combinations of local and global conservation laws are analyzed, and the physical meaningful cases are identified. Several new kinds of mixed phases are presented, as e.g. a locally charge neutral mixed phase in proto-neutron stars which will disappear during the cooling and deleptonization of the proto-neutron star.

We determine the Galactic production rate of strangelets as a canonical input to calculations of the measurable cosmic ray flux of strangelets by performing simulations of strange star mergers and combining the results with recent estimates of stellar binary populations. We find that the flux depends sensitively on the bag constant of the MIT bag model of QCD and disappears for high values of the bag constant and thus more compact strange stars. In the latter case strange stars could coexist with ordinary neutron stars as they are not converted by the capture of cosmic ray strangelets. An unambiguous detection of an ordinary neutron star would then not rule out the strange matter hypothesis.

The implications of the formation of strange quark matter in neutron stars and in core-collapse supernovae is discussed with special emphasis on the possibility of having a strong first order QCD phase transition at high baryon densities. If strange quark matter is formed in core-collapse supernovae shortly after the bounce, it causes the launch of a second outgoing shock which is energetic enough to lead to a explosion. A signal for the formation of strange quark matter can be read off from the neutrino spectrum, as a second peak in antineutrinos is released when the second shock runs over the neutrinosphere.

The possible role of a first order QCD phase transition at nonvanishing quark chemical potential and temperature for cold neutron stars and for supernovae is delineated. For cold neutron stars, we use the NJL model with nonvanishing color superconducting pairing gaps, which describes the phase transition to the 2SC and the CFL quark matter phases at high baryon densities. We demonstrate that these two phase transitions can both be present in the core of neutron stars and that they lead to the appearance of a third family of solution for compact stars. In particular, a core of CFL quark matter can be present in stable compact star configurations when slightly adjusting the vacuum pressure to the onset of the chiral phase transition from the hadronic model to the NJL model. We show that a strong first order phase transition can have strong impact on the dynamics of core collapse supernovae. If the QCD phase transition sets in shortly after the first bounce, a second outgoing shock wave can be generated which leads to an explosion. The presence of the QCD phase transition can be read off from the neutrino and antineutrino signal of the supernova.

We explore the implications of the QCD phase transition during the postbounce evolution of core-collapse supernovae. Using the MIT bag model for the description of quark matter and assuming small bag constants, we find that the phase transition occurs during the early postbounce accretion phase. This stage of the evolution can be simulated with general relativistic three-flavor Boltzmann neutrino transport. The phase transition produces a second shock wave that triggers a delayed supernova explosion. If such a phase transition happens in a future galactic supernova, its existence and properties should become observable as a second peak in the neutrino signal that is accompanied by significant changes in the energy of the emitted neutrinos. In contrast to the first neutronization burst, this second neutrino burst is dominated by the emission of anti-neutrinos because the electron-degeneracy is lifted when the second shock passes through the previously neutronized matter.

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 effect of the QCD phase transition is studied for the mass-radius relation of compact stars and for hot and dense matter at a given proton fraction used as input in core-collapse supernova simulations. The phase transitions to the 2SC and CFL color superconducting phases lead to stable hybrid star configurations with a pure quark matter core. In supernova explosions quark matter could be easily produced due to $\beta$-equilibrium, small proton fractions and nonvanishing temperatures. A low critical density for the phase transition to quark matter is compatible with present pulsar mass measurements.

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.

We study the equation state of strongly interacting quark matter within a NJL-like model in which the chiral condensates and the color superconducting gaps are computed self-consistently as a function of the baryon density. A vector interaction term is added to the Lagrangian in order to render the quark matter equation of state stiffer. For the low density hadronic phase we use a relativistic mean field model. The phase transition to quark matter is computed by a Maxwell construction. We show that stable CFL cores in hybrid stars are possible if the superconducting gap is sufficiently large. Moreover we find stable stellar configurations in which two phase transitions occur, a first transition from hadronic matter to 2SC quark matter and a second transition from 2SC quark matter to CFL quark matter.

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.

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

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.