We study the possible occurrence of the hadron-quark phase transition (PT) during the merging of neutron star binaries by hydrodynamical simulations employing a set of temperature dependent hybrid equations of state (EoSs). Following previous work we describe an unambiguous and measurable signature of deconfined quark matter in the gravitational-wave (GW) signal of neutron star binary mergers including equal-mass and unequal-mass systems of different total binary mass. The softening of the EoS by the PT at higher densities, i.e. after merging, leads to a characteristic increase of the dominant postmerger GW frequency f_peak relative to the tidal deformability Lambda inferred during the premerger inspiral phase. Hence, measuring such an increase of the postmerger frequency provides evidence for the presence of a strong PT. If the postmerger frequency and the tidal deformability are compatible with results from purely baryonic EoS models yielding very tight relations between f_peak and Lambda, a strong PT can be excluded up to a certain density. We find tight correlations of f_peak and Lambda with the maximum density during the early postmerger remnant evolution. These GW observables thus inform about the density regime which is probed by the remnant and its GW emission. Exploiting such relations we devise a directly applicable, concrete procedure to constrain the onset density of the QCD PT from future GW measurements. We point out two interesting scenarios: if no indications for a PT are inferred from a GW detection, our procedure yields a lower limit on the onset density of the hadron quark PT. On the contrary, if a merger event reveals evidence for the occurrence of deconfined quark matter, the inferred GW parameters set an upper limit on the PT onset density. (abridged)

With the recent dawn of the multi-messenger astronomy era a new window has opened to explore the constituents of matter and their interactions under extreme conditions. One of the pending challenges of modern physics is to probe the microscopic equation of state (EoS) of cold and dense matter via macroscopic neutron star observations such as their masses and radii. Still unanswered issues concern the detailed composition of matter in the core of neutron stars at high pressure and the possible presence of e.g. hyperons or quarks. By means of a non-perturbative functional renormalization group approach the influence of quantum and density fluctuations on the quark matter EoS in $\beta$-equilibrium is investigated within two- and three-flavor quark-meson model truncations and compared to results obtained with common mean-field approximations where important fluctuations are usually ignored. We find that they strongly impact the quark matter EoS.

Neutrinos play an important role in compact star astrophysics: neutrino-heating is one of the main ingredients in core-collapse supernovae, neutrino-matter interactions determine the composition of matter in binary neutron star mergers and have among others a strong impact on conditions for heavy element nucleosynthesis and neutron star cooling is dominated by neutrino emission except for very old stars. Many works in the last decades have shown that in dense matter medium effects considerably change the neutrino-matter interaction rates, whereas many astrophysical simulations use analytic approximations which are often far from reproducing more complete calculations. In this work we present a scheme which allows to incorporate improved rates, for charged current interactions, into simulations and show as an example some results for core-collapse supernovae, where a noticeable difference is found in the location of the neutrinospheres of the low-energy neutrinos in the early post-bounce phase.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates on individual nuclei, the nuclear interaction for the EoS, the mass model for the nuclear statistical equilibrium distribution and the progenitor model. In agreement with previous works, we show that the individual EC rates are the most important source of uncertainty in the simulations, while the other inputs only marginally influence the results. A recently proposed analytic formula to extrapolate microscopic results on stable nuclei for EC rates to the neutron rich region, with a functional form motivated by nuclear-structure data and parameters fitted from large scale shell model calculations, is shown to lead to a sizable (16%) reduction of the electron fraction at bounce compared to more primitive prescriptions for the rates, leading to smaller inner core masses and slower shock propagation. We show that the EC process involves $\approx$ 130 different nuclear species around 86 Kr mainly in the N = 50 shell closure region, and establish a list of the most important nuclei to be studied in order to constrain the global rates.

We give an overview of the SImulation of X-ray TElescopes (SIXTE) software package, a generic, mission-independent Monte Carlo simulation toolkit for X-ray astronomical instrumentation. The package is based on a modular approach for the source definition, the description of the optics, and the detector type such that new missions can be easily implemented. The targets to be simulated are stored in a flexible input format called SIMPUT. Based on this source definition, a sample of photons is produced and then propagated through the optics. In order to model the detection process, the software toolkit contains modules for various detector types, ranging from proportional counter and Si-based detectors, to more complex descriptions like transition edge sensor (TES) devices. The implementation of characteristic detector effects and a detailed modeling of the read-out process allow for representative simulations and therefore enable the analysis of characteristic features, such as for example pile-up, and their impact on observations. We present an overview of the implementation of SIXTE from the input source, the imaging, and the detection process, highlighting the modular approach taken by the SIXTE software package. In order to demonstrate the capabilities of the simulation software, we present a selection of representative applications, including the all-sky survey of eROSITA and a study of pile-up effects comparing the currently operating XMM-Newton with the planned Athena-WFI instrument. A simulation of a galaxy cluster with the Athena- X-IFU shows the capability of SIXTE to predict the expected performance of an observation for a complex source with a spatially varying spectrum and our current knowledge of the future instrument.

Using an axisymmetric numerical code, we perform an extensive study of the magnetic field configurations in non-rotating neutron stars, varying the mass, magnetic field strength and the equation of state. We find that the monopolar (spherically symmetric) part of the norm of the magnetic field can be described by a single profile, that we fit by a simple eighth-order polynomial, as a function of the star's radius. This new generic profile applies remarkably well to all magnetized neutron star configurations built on hadronic equations of state. We then apply this profile to build magnetized neutron stars in spherical symmetry, using a modified Tolman-Oppenheimer-Volkov (TOV) system of equations. This new formalism produces slightly better results in terms of mass-radius diagrams than previous attempts to add magnetic terms to these equations. However, we show that such approaches are less accurate than usual, non-magnetized TOV models, and that consistent models must depart from spherical symmetry. Thus, our `universal' magnetic field profile is intended to serve as a tool for nuclear physicists to obtain estimates of magnetic field inside neutron stars, as a function of radial depth, in order to deduce its influence on composition and related properties. It possesses the advantage of being based on magnetic field distributions from realistic self-consistent computations, which are solutions of Maxwell's equations.

We describe a multi-messenger interpretation of GW170817, which yields a robust lower limit on NS radii. This excludes NSs with radii smaller than about 10.7 km and thus rules out very soft nuclear matter. We stress the potential of this type of constraints when future detections become available. A very similar argumentation may yield an upper bound on the maximum mass of nonrotating NSs. We also discuss simulations of NS mergers, which undergo a first-order phase transition to quark matter. We point out a different dynamical behavior. Considering the gravitational-wave signal, we identify an unambiguous signature of the QCD phase transition in NS mergers. The occurrence of quark matter through a strong first-order phase transition during merging leads to a characteristic shift of the dominant postmerger frequency. The frequency shift is indicative for a phase transition if it is compared to the postmerger frequency which is expected for purely hadronic EoS models. A very strong deviation of several 100 Hz is observed for hybrid EoSs in an otherwise tight relation between the tidal deformability and the postmerger frequency. We address the potential impact of a first-order phase transition on the electromagnetic counterpart of NS mergers. Our simulations suggest that there would be no significant qualitative differences between a system undergoing a phase transition to quark matter and purely hadronic mergers. The quantitative differences are within the spread which is found between different hadronic EoS models. This implies on the one hand that GW170817 is compatible with a possible transition to quark matter. On the other hand these considerations show that it may not be easy to identify quantitative differences between purely hadronic mergers and events in which quark matter occurs considering solely their electromagnetic counterpart or their nucleosynthesis products. (abridged)

We identify an observable imprint of a first-order hadron-quark phase transition at supranuclear densities on the gravitational-wave (GW) emission of neutron star mergers. Specifically, we show that the dominant postmerger GW frequency f_peak may exhibit a significant deviation from an empirical relation between f_peak and the tidal deformability if a strong first-order phase transition leads to the formation of a gravitationally stable extended quark matter core in the postmerger remnant. A comparison of the GW signatures from a large, representative sample of microphysical, purely hadronic equations of state indicates that this imprint is only observed in those systems which undergo a strong first-order phase transition. Such a shift of the dominant postmerger GW frequency can be revealed by future GW observations, which would provide evidence for the existence of a strong first-order phase transition in the interior of neutron stars.

For core-collapse and neutron star merger simulations it is important to have at hand adequate equations of state, describing the underlying dense and hot matter as realistically as possible. Here, we present two newly constructed equation of state (EoS) including the entire baryon octet. Both EoS are compatible with the main constraints from nuclear physics, both experimental and theoretical. One of the EoS is equally describing maximum mass for cold $\beta$-equilibrated neutron stars of $2 M_\odot$ in agreement with recent observations. The predictions obtained with the new EoS are compared with the results obtained with DD2Y, the only presently existing EoS containing the baryon octet, that satisfies the same constraints within uncertainties. The main difference between our new EoS models and DD2Y is the harder symmetry energy of the latter. We show that the density dependence of the symmetry energy has a direct influence on the amount of strangeness inside hot and dense matter and consequently on thermodynamic quantities, e.g. the temperature for given entropy per baryon. We expect these differences to affect the evolution of a protoneutron star or binary neutron star mergers. We also propose several parametrizations calibrated to $\Lambda$ hypernuclei based on the DD2 and SFHo models that satisfy the two solar mass constraint.

We present a three-dimensional map of the hydrogen density distribution in the Galactic interstellar medium. The hydrogen equivalent column densities were obtained from the Exploring the X-ray Transient and variable Sky project ({\sc EXTraS}) which provides equivalent $N_{\rm H}$ values from X-ray spectral fits of observations within the {\it XMM-Newton} Data Release. {\sc EXTraS} include multiple fits for each source, allowing an accurate determination of the equivalent column densities, which depends on the continuum modeling of the spectra. A cross-correlation between the {\sc EXTraS} catalogue and the first {\it Gaia} Data Release was performed in order to obtain accurate parallax and distance measurements. We use a Bayesian method explained in \citet{rez17} in order to predict the most probable distribution of the density at any arbitrary point, even for lines of sight along which there are no initial observation. The resulting map shows small-scale density structures which can not been modeled by using analytic density profiles. In this paper we present a proof of concept of the kind of science possible with the synergy of these catalogs. However, given the systematic uncertainties connected to the source identification and to the dependence of $N_{\rm H}$ on the spectral model, the present maps should be considered qualitatively at this point.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

In this work we present a newly constructed equation of state (EoS) --applicable to stellar core collapse and neutron star mergers--, including the entire baryon octet. Our EoS is compatible with the main constraints from nuclear physics and, in particular, with a maximum mass for cold beta-equilibrated neutron stars of 2 solar masses in agreement with recent observations. As an application of our new EoS, we compute numerical stationary models for rapidly (rigidly) rotating hot neutron stars. We consider maximum masses of hot stars, such as proto-neutron stars or hypermassive neutron stars in the post-merger phase of binary neutron star coalescence. The universality of I-Q-relations at nonzero temperature for fast rotating models, comparing a purely nuclear EoS with its counterparts containing Lambda-hyperons or the entire baryon octet, respectively, is discussed, too. We find that the I-Q universality is broken when thermal effects become important, whatever the value of entropy gradients in our models. Thus, the use of I-Q relations for the analysis of proto-neutron stars or merger remnant data, including gravitational wave signals from the last stages of binary neutron star mergers, should be regarded with care.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We give an overview about equations of state (EOS) which are currently available for simulations of core-collapse supernovae and neutron star mergers. A few selected important aspects of the EOS, such as the symmetry energy, the maximum mass of neutron stars, and cluster formation, are confronted with constraints from experiments and astrophysical observations. There are just very few models which are compatible even with this very restricted set of constraints. These remaining models illustrate the uncertainty of the uniform nuclear matter EOS at high densities. In addition, at finite temperatures the medium modifications of nuclear clusters represent a conceptual challenge. In conclusion, there has been significant development in the recent years, but there is still need for further improved general purpose EOS tables.

We present a numerical model for uniformly rotating superfluid neutron stars, for the first time with realistic microphysics including entrainment, in a fully general relativistic framework. We compute stationary and axisymmetric configurations of neutron stars composed of two fluids, namely superfluid neutrons and charged particles (protons and electrons), rotating with different rates around a common axis. Both fluids are coupled by entrainment, a non-dissipative interaction which in case of a non-vanishing relative velocity between the fluids, causes the fluid momenta being not aligned with the respective fluid velocities. We extend the formalism by Comer and Joynt (2003) in order to calculate the equation of state (EoS) and entrainment parameters for an arbitrary relative velocity. The resulting entrainment matrix fulfills all necessary sum rules and in the limit of small relative velocity our results agree with Fermi liquid theory ones, derived to lowest order in the velocity. This formalism is applied to two new nuclear equations of state which are implemented in the numerical model. We are able to obtain precise equilibrium configurations. Resulting density profiles and moments of inertia are discussed employing both EoSs, showing the impact of entrainment and the dependence on the EoS.

In this paper, we study in detail the role of general relativity on the global dynamics of giant pulsar glitches as exemplified by Vela. For this purpose, we carry out numerical simulations of the spin up triggered by the sudden unpinning of superfluid vortices. In particular, we compute the exchange of angular momentum between the core neutron superfluid and the rest of the star within a two-fluid model including both (non-dissipative) entrainment effects and (dissipative) mutual friction forces. Our simulations are based on a quasi-stationary approach using realistic equations of state (EoSs). We show that the evolution of the angular velocities of both fluids can be accurately described by an exponential law. The associated characteristic rise time $\tau_{\text{r}}$, which can be precisely computed from stationary configurations only, has a form similar to that obtained in the Newtonian limit. However, general relativity changes the structure of the star and leads to additional couplings between the fluids due to frame-dragging effects. As a consequence, general relativity can have a large impact on the actual value of $\tau_{\text{r}}$: the errors incurred by using Newtonian gravity are thus found to be as large as $\sim 40 \%$ for the models considered. Values of the rise time are calculated for Vela and compared with current observational limits. Finally, we study the amount of gravitational waves emitted during a glitch. Simple expressions are obtained for the corresponding characteristic amplitudes and frequencies. The detectability of glitches through gravitational wave observatories is briefly discussed.

We develop a detailed and self-consistent numerical model for extremely-magnetised white dwarfs, which have been proposed as progenitors of overluminous Type Ia supernovae. This model can describe fully-consistent equilibria of magnetic stars in axial symmetry, with rotation, general-relativistic effects and realistic equations of state (including electron-ion interactions and taking into account Landau quantisation of electrons due to the magnetic field). We study the influence of each of these ingredients onto the white dwarf structure and, in particular, on their maximum mass. We perform an extensive stability analysis of such objects, with their highest surface magnetic fields reaching $\sim 10^{13}~G$ (at which point the star adopts a torus-like shape). We confirm previous speculations that although very massive strongly magnetised white dwarfs could potentially exist, the onset of electron captures and pycnonuclear reactions may severely limit their stability. Finally, the emission of gravitational waves by these objects is addressed, showing no possibility of detection by the currently planned space-based detector eLISA.

We review various theoretical approaches for the equation of state (EoS) of dense matter, relevant for the description of core-collapse supernovae, compact stars and compact star mergers. The emphasis is put on models that are applicable to all of these scenarios. Such EoS models have to cover large ranges in baryon number density, temperature and isospin asymmetry. The characteristics of matter change dramatically within these ranges, from a mixture of nucleons, nuclei, and electrons to uniform, strongly interacting matter containing nucleons, and possibly other particles such as hyperons or quarks. As the development of an EoS requires joint efforts from many directions we consider different theoretical approaches and discuss relevant experimental and observational constraints which provide insights for future research. Finally, results from applications of the discussed EoS models are summarized.

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analyzed. We also review its potential for public outreach in science museums.

With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.

The Square Kilometre Array (SKA) will make ground breaking discoveries in pulsar science. In this chapter we outline the SKA surveys for new pulsars, as well as how we will perform the necessary follow-up timing observations. The SKA's wide field-of-view, high sensitivity, multi-beaming and sub-arraying capabilities, coupled with advanced pulsar search backends, will result in the discovery of a large population of pulsars. These will enable the SKA's pulsar science goals (tests of General Relativity with pulsar binary systems, investigating black hole theorems with pulsar-black hole binaries, and direct detection of gravitational waves in a pulsar timing array). Using SKA1-MID and SKA1-LOW we will survey the Milky Way to unprecedented depth, increasing the number of known pulsars by more than an order of magnitude. SKA2 will potentially find all the Galactic radio-emitting pulsars in the SKA sky which are beamed in our direction. This will give a clear picture of the birth properties of pulsars and of the gravitational potential, magnetic field structure and interstellar matter content of the Galaxy. Targeted searches will enable detection of exotic systems, such as the ~1000 pulsars we infer to be closely orbiting Sgr A*, the supermassive black hole in the Galactic Centre. In addition, the SKA's sensitivity will be sufficient to detect pulsars in local group galaxies. To derive the spin characteristics of the discoveries we will perform live searches, and use sub-arraying and dynamic scheduling to time pulsars as soon as they are discovered, while simultaneously continuing survey observations. The large projected number of discoveries suggests that we will uncover currently unknown rare systems that can be exploited to push the boundaries of our understanding of astrophysics and provide tools for testing physics, as has been done by the pulsar community in the past.

We present a self-consistent model for the study of the structure of a neutron star in strong magnetic fields. Starting from a microscopic Lagrangian, this model includes the effect of the magnetic field on the equation of state, the interaction of the electromagnetic field with matter (magnetisation), and anisotropies in the energy-momentum tensor, as well as general relativistic aspects. We build numerical axisymmetric stationary models and show the applicability of the approach with one example quark matter equation of state (EoS) often employed in the recent literature for studies of strongly magnetised neutron stars. For this EoS, the effect of inclusion of magnetic field dependence or the magnetisation do not increase the maximum mass significantly in contrast to what has been claimed by previous studies.

This paper is the manual of the CompOSE service, see http://compose.obspm.fr, an equation of state (EoS) database for use in nuclear physics and in astrophysical applications, e.g. the simulation of core-collapse supernovae and compact stars. The main aims of the service are the following: to offer a robust, easily extendable data format to store and read EoS data tables and related physical quantities; to make these data available to the public -- respecting the original work by the contributors -- in a standard format designed for easy access and flexibility; and to provide routines to interpolate these data in the parameter space of density, temperature and hadronic charge fraction. The manual contains a detailed description of the service. It is divided into three parts, containing a general introduction as well as instructions for potential contributors and for users.

We present numerical simulations of stellar core-collapse with spherically symmetric, general relativistic hydrodynamics up to black hole formation. Using the CoCoNuT code, with a newly developed grey leakage scheme for the neutrino treatment, we investigate the effects of including pions and \Lambda-hyperons into the equation of state at high densities and temperatures on the black hole formation process. Results show non-negligible differences between the models with reference equation of state without any additional particles and models with the extended ones. For the latter, the maximum masses supported by the proto-neutron star are smaller and the collapse to a black hole occurs earlier. A phase transition to hyperonic matter is observed when the progenitor allows for a high enough accretion rate onto the proto-neutron star. Rough estimates of neutrino luminosity from these collapses are given, too.

In stellar core-collapse events matter is heated and compressed to densities above nuclear matter saturation density. For progenitors stars with masses above about 25 solar masses, which eventually form a black hole, the temperatures and densities reached during the collapse are so high that a traditional description in terms of electrons, nuclei, and nucleons is no longer adequate. We present here an improved equation of state which contains in addition pions and hyperons. They become abundant in the high temperature and density regime. We study the different constraints on such an equation of state, coming from both hyperonic data and observations of neutron star properties. In order to test the zero-temperature versions, we perform numerical simulations of the collapse of a neutron star with such additional particles to a black hole. We discuss the influence of the additional particles on the thermodynamic properties within the hot versions of the equation of state and we show that in regimes relevant to core-collapse and black hole formation, the effects of pions and hyperons on pressure, internal energy and sound speed are not negligible.

We study the spectrum of collective excitations in the inhomogeneous phases in the neutron star inner crust within a superfluid hydrodynamics approach. Our aim is to describe the whole range of wavelengths, from the long-wavelength limit which can be described by macroscopic approaches and which is crucial for the low-energy part of the spectrum, to wavelengths of the order of the dimensions of the Wigner-Seitz cells, corresponding to the modes usually described in microscopic calculations. As an application, we will discuss the contribution of these collective modes to the specific heat of the "lasagna" phase in comparison with other known contributions.

Surface effects in strange-quark matter play an important role for certain observables which have been proposed in order to identify strange stars, and color superconductivity can strongly modify these effects. We study the surface of color superconducting strange-quark matter by solving the Hartree-Fock-Bogoliubov equations for finite systems ("strangelets") within the MIT bag model, supplemented with a pairing interaction. Due to the bag-model boundary condition, the strange-quark density is suppressed at the surface. This leads to a positive surface charge, concentrated in a layer of ~1 fm below the surface, even in the color-flavor locked (CFL) phase. However, since in the CFL phase all quarks are paired, this positive charge is compensated by a negative charge, which turns out to be situated in a layer of a few tens of fm below the surface, and the total charge of CFL strangelets is zero. We also study the surface and curvature contributions to the total energy. Due to the strong pairing, the energy as a function of the mass number is very well reproduced by a liquid-drop type formula with curvature term.

We perform a study of the possible existence of hybrid stars with color superconducting quark cores using a specific hadronic model in a combination with an NJL-type quark model. It is shown that the constituent mass of the non-strange quarks in vacuum is a very important parameter that controls the beginning of the hadron-quark phase transition. At relatively small values of the mass, the first quark phase that appears is the two-flavor color superconducting (2SC) phase which, at larger densities, is replaced by the color-flavor locked (CFL) phase. At large values of the mass, on the other hand, the phase transition goes from the hadronic phase directly into the CFL phase avoiding the 2SC phase. It appears, however, that the only stable hybrid stars obtained are those with the 2SC quark cores.

We explore the relevance of color superconductivity inside a possible quark matter core for the bulk properties of neutron stars. For the quark phase we use an Nambu--Jona-Lasinio (NJL) type model, extended to include diquark condensates. For the hadronic phase, a microscopic many-body model is adopted, with and without strangeness content. In our calculations, a sharp boundary is assumed between the hadronic and the quark phases. For NJL model parameters fitted to vacuum properties we find that no star with a pure quark core does exist. Nevertheless the presence of color superconducting phases can lower the neutron star maximum mass substantially. In some cases, the transition to quark matter occurs only if color superconductivity is present. Once the quark phase is introduced, the value of the maximum mass stays in any case below the value of two solar masses.