The tidal deformability probability distribution extracted from GW170817 alone, or including multi-messenger information, is confronted to astrophysical and nuclear physics additional constraints within a semi-agnostic approach for the dense matter equation of state. We use Bayesian statistics to combine together low density nuclear physics data, such as the ab-initio predictions based on $\chi$EFT interactions or the isoscalar giant monopole resonance, and astrophysical constraints from neutron stars, such as the maximum mass of neutron stars or the probability density function of the tidal deformability $\tilde{\Lambda}$ obtained from the GW170817 event. The posteriors probability distribution functions are marginalized over several nuclear empirical parameters ($L_\textrm{sym}$, $K_\textrm{sym}$, $Q_\textrm{sat}$ and $Q_\textrm{sym}$), as well as over observational quantities such as the $1.4M_\odot$ radius $R_{1.4}$ and the pressure at twice the saturation density $P(2n_\textrm{sat})$. The correlations between $L_\textrm{sym}$ and $K_\textrm{sym}$ and between $K_\textrm{sat}$ and $Q_\textrm{sat}$ are also further analyzed. Tension is found between the posteriors: the first one is localized in the tidal deformability probability distribution itself, depending whether multi-messenger analysis is included or not, and the second one is between the observational data and the nuclear physics inputs. These tensions impact the predictions for $L_\textrm{sym}$, $K_\textrm{sym}$ and $R_{1.4}$ with centroids which differ by 2-3$\sigma$. Implications for the nuclear equation of state are also discussed.

This paper presents a new analysis of the thermal emission from the neutron star surface to constrain the dense matter equation of state. It is based on the use of a Markov-Chain Monte Carlo algorithm combined with an empirical parametrization of the equation of state, as well as the consistent treatment of seven neutron star quiescent low-mass X-ray binaries in globular clusters with well-measured distances. Previous analyses have indicated that the thermal emission of these neutron stars tends to prefer low neutron star radii, questioning basic knowledge from nuclear physics. We show that it is possible to reconcile the thermal emission analyses with nuclear physics knowledge, with or without including a prior on the slope of the symmetry energy $L_{\rm sym}$. We obtain radii of the order of about 12~km without worsening the fit statistic. With an empirical parametrization of the equation of state, we obtain the following values for the slope of the symmetry energy, its curvature $K_{\rm sym}$, and the isoscalar skewness parameter $Q_{\rm sat}$: $L_{\rm sym}=37.2^{+9.2}_{-8.9}$ MeV, $K_{\rm sym}=-85^{+82}_{-70}$ MeV, and $Q_{\rm sat}=318^{+673}_{-366}$ MeV. For the first time, we measure the values of the empirical parameters $K_{\rm sym}$ and $Q_{\rm sat}$. These values are only weakly impacted by our assumptions, such as the distances or the number of free empirical parameters, provided they are taken within a reasonable range. We also study the weak sensitivity of our results to the set of sources analyzed, and we identify a group of sources that dominates the constraints. The resulting masses and radii obtained are also discussed in the context of the independent constraints from GW 170817 and its electromagnetic counterpart, AT 2017gfo.

Pulsating thermal X-ray emission from millisecond pulsars can be used to obtain constraints on the neutron star equation of state, but to date only five such sources have been identified. Of these five millisecond pulsars, only two have well constrained neutron star masses, which improve the determination of the radius via modelling of the X-ray waveform. We aim to find other millisecond pulsars that already have well constrained mass and distance measurements that show pulsed thermal X-ray emission in order to obtain tight constraints on the neutron star equation of state. The millisecond pulsar PSR~J1909--3744 has an accurately determined mass, M = 1.54$\pm$0.03 M$_\odot$ (1 $\sigma$ error) and distance, D = 1.07$\pm$0.04 kpc. We analysed {\em XMM-Newton} data of this 2.95 ms pulsar to identify the nature of the X-ray emission. We show that the X-ray emission from PSR~J1909--3744 appears to be dominated by thermal emission from the polar cap. Only a single component model is required to fit the data. The black-body temperature of this emission is kT=0.26\ud{0.03}{0.02} keV and we find a 0.2--10 keV un-absorbed flux of 1.1 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ or an un-absorbed luminosity of 1.5 $\times$ 10$^{30}$ erg s$^{-1}$. Thanks to the previously determined mass and distance constraints of the neutron star PSR~J1909--3744, and its predominantly thermal emission, deep observations of this object with future X-ray facilities should provide useful constraints on the neutron star equation of state.

Multi-messenger observations of neutron star (NS) mergers have the potential to revolutionize nuclear astrophysics. They will improve our understanding of nucleosynthesis, provide insights about the equation of state (EOS) of strongly-interacting matter at high densities, and enable tests of the theory of gravity and of dark matter. Here, we focus on the EOS, where both gravitational waves (GWs) from neutron-star mergers and X-ray observations from space-based detectors such as NICER will provide more stringent constraints on the structure of neutron stars. Furthermore, recent advances in nuclear theory have enabled reliable calculations of the EOS at low densities using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of GWs from GW170817 can be combined with modern calculations of the EOS to extract useful insights about the EOS of matter encountered inside neutron stars. We analyze the impact of various uncertainties, the role of phase transitions in the NS core, and discuss how future observations will improve our understanding of dense matter.

Gravitational waves from neutron-star mergers are expected to provide stringent constraints on the structure of neutron stars. At the same time, recent advances in nuclear theory have enabled reliable calculations of the low density equation of state using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of gravitational waves from GW170817 can be combined with modern calculations of the equation of state to extract useful insights about the equation of state of matter encountered inside neutron stars. We analyze the impact of various uncertainties and we show that the tidal deformability extracted from GW170817 is compatible, while less constraining, than modern nuclear physics knowledge.

The neutron star crustal EoS and transition point properties are computed within a unified meta-modeling approach. A Bayesian approach is employed including two types of filters: bulk nuclear properties are controlled from low density effective field theory (EFT) predictions as well as the present knowledge from nuclear experiments, while the surface energy is adjusted on experimental nuclear masses. Considering these constraints, a quantitative prediction of crustal properties can be reached with controlled confidence intervals and increased precision with respect to previous calculations: {$\approx 11\%$} dispersion on the crustal width and {$\approx 27\%$} dispersion on the fractional moment of inertia. The crust moment of inertia is also evaluated as a function of the neutron star mass, and predictions for mass and radii are given for different pulsars. The possible crustal origin of Vela pulsar glitches is discussed within the present estimations of crustal entrainment, disfavoring a large entrainment phenomenon if the Vela mass is above $1.4M_\odot$. Further refinement of the present predictions requires a better estimation of the high order isovector empirical parameters, e.g. $K_{sym}$ and $Q_{sym}$, and a better control of the surface properties of extremely neutron rich nuclei.

The correlation of the tidal polarizabilities $\Lambda_1$-$\Lambda_2$ for GW170817 is predicted by combining dense-matter equations of state (EOSs) that satisfy nuclear physics constraints with the chirp mass and mass asymmetry for this event. Our models are constrained by calculations of the neutron-matter EOS using chiral effective field theory Hamiltonians with reliable error estimates up to once or twice the nuclear saturation density. In the latter case, we find that GW170817 does not improve our understanding of the EOS. We contrast two distinct extrapolations to higher density: a minimal model (MM) which assumes that the EOS is a smooth function of density described by a Taylor expansion and a more general model parametrized by the speed of sound that admits phase transitions. This allows us to identify regions in the $\Lambda_1$-$\Lambda_2$ plots that could favor the existence of new phases of matter in neutron stars. We predict the combined tidal polarizability of the two neutron stars in GW170817 to be $80\le \tilde{\Lambda}\le 580$ ($280\le \tilde{\Lambda}\le 480$ for the MM), which is smaller than the range suggested by the LIGO-Virgo data analysis. Our analysis also shows that GW170817 requires a NS with $M=1.4M_\odot$ to have a radius $9.0<R_{1.4}<13.6$~km ($ 11.3 <R_{1.4}< 13.6$~km for the MM).

A quantitative knowledge of the nuclear Equation of State (EoS) requires an accurate estimation of the uncertainties on the EoS parameters and their mutual correlations. Such correlations are empirically observed in a large set of EoS models by different authors, but they are not always fully understood. We show that some of these correlations can be interpreted from basic physical constraints imposed on a simple Taylor expansion of the binding energy around saturation density. In particular, we investigate the correlations among the following empirical parameters: the symmetry energy $E_{sym}$, the slope and curvature of the symmetry energy $L_{sym}$ and $K_{sym}$, and the curvature and skewness of the binding energy in symmetric matter $K_{sat}$ and $Q_{sat}$. The uncertainties on these correlations is estimated through analytical modelling as well as a meta-modelling analysis of the EoS subject to physical constraints. We show that a huge dispersion of the correlations among low order empirical parameters is induced by the unknown higher order empirical parameters, such as $K_{sym}$ (second order) and $Q_{sat}$ and $Q_{sym}$ (third order). We also propose an explanation of the reason why $Q_{sat}$ is weakly constrained by present experimental data. We conclude that selective observables on high order parameters, such as $K_{sym}$ and $Q_{sat}$, should be better determined before the present uncertainties on the EoS can be further reduced.

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.

We find that the abundance of large clusters of nucleons in neutron-rich matter at sub-nuclear density is greatly reduced by finite temperature effects when matter is close to beta-equilibrium. Large nuclei and exotic non-spherical nuclear configurations called pasta, favored in the vicinity of the transition to uniform matter at $T=0$, dissolve at relatively low temperature. For matter close to beta-equilibrium we find that the pasta melting temperature is $T_m^\beta\simeq 4\pm 1$~MeV for realistic equations of state. The mechanism for pasta dissolution is discussed, and in general $T_m^\beta$ is shown to be sensitive to the proton fraction. We find that coherent neutrino scattering from nuclei and pasta makes a modest contribution to the opacity under the conditions encountered in supernovae and neutron star mergers. Implications for neutrino signals from galactic supernovae are briefly discussed.

A consistent Hartree-Fock study of the equation of state (EOS) of asymmetric nuclear matter at finite temperature has been performed using realistic choices of the effective, density dependent nucleon-nucleon (NN) interaction, which were successfully used in different nuclear structure and reaction studies. Given the importance of the nuclear symmetry energy in the neutron star formation, EOS's associated with different behaviors of the symmetry energy were used to study hot asymmetric nuclear matter. The slope of the symmetry energy and nucleon effective mass with increasing baryon density was found to affect the thermal properties of nuclear matter significantly. Different density dependent NN interactions were further used to study the EOS of hot protoneutron star (PNS) matter of the $npe\mu\nu$ composition in $\beta$-equilibrium. The hydrostatic configurations of PNS in terms of the maximal gravitational mass $M_{\rm max}$ and radius, central density, pressure and temperature at the total entropy per baryon $S/A= 1,2$ and 4 have been determined in both the neutrino-free and neutrino-trapped scenarios. The obtained results show consistently a strong impact of the symmetry energy and nucleon effective mass on thermal properties and composition of hot PNS matter. $M_{\rm max}$ values obtained for the (neutrino-free) $\beta$-stable PNS at $S/A=4$ were used to assess time $t_{\rm BH}$ of the collapse of 40 $M_\odot$ protoneutron progenitor to black hole, based on a correlation between $t_{\rm BH}$ and $M_{\rm max}$ found from the hydrodynamic simulation by Hempel {\it et al.}.

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.

The heat capacity of neutron matter is studied over the range of densities and temperatures prevailing in neutron-star crusts, allowing for the transition to a superfluid phase at temperatures below some critical temperature $T_{sf}$ and including the transition to the classical limit. Finite temperature Hartree-Fock-Bogoliubov equations (FTHFB) are solved and compared to existing approximate expressions. In particular, the formula given by Levenfish and Yakovlev is found to reproduce the numerical results with a high degree of accuracy for temperatures $T\leq T_{sf}$. In the non-superfluid phase, $T\geq T_{sf}$, the linear approximation is valid only at temperature $T\ll T_{{\rm F} n}$ ($T_{{\rm F} n}$ being the Fermi temperature of the neutron gas) which is rarely the case in the shallow layers of the neutron star's crust. A non-perturbative interpolation between the quantal and the classical regimes is proposed here. The heat capacity, conveniently parametrized solely in terms of $T_{sf}$, $T_{{\rm F} n}$, and the neutron number density $n_n$, can be easily implemented in neutron-star cooling simulations.

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 analyse how the structure of the inner curst is influenced by the pairing correlations. The inner-crust matter, formed by nuclear clusters immersed in a superfluid neutron gas and ultra-relativistic electrons, is treated in the Wigner-Seitz approximation. The properties of the Wigner-Seitz cells, i.e., their neutron to proton ratio and their radius at a given baryonic density, are obtained from the energy minimization at beta equilibrium. To obtain the binding energy of baryonic matter we perform Skyrme-HFB calculations with zero-range density-dependent pairing forces of various intensities. We find that the Wigner-Seitz cells have much smaller numbers of protons compared to previous calculations. For the dense cells the binding energy of the configurations with small proton numbers do not converge to a well-defined minimum value which precludes the determination of their structure. We show that for these cells there is a significant underestimation of the binding energy due to the boundary conditions at the border of the cells imposed through the Wigner-Seitz approximation.

The nuclear mean-field potentials obtained in the Hartree-Fock method with different choices of the in-medium nucleon-nucleon (NN) interaction have been used to study the equation of state (EOS) of the neutron star (NS) matter. The EOS of the uniform NS core has been calculated for the np$e\mu$ composition in the $\beta$-equilibrium at zero temperature, using version Sly4 of the Skyrme interaction as well as two density-dependent versions of the finite-range M3Y interaction (CDM3Y$n$ and M3Y-P$n$), and versions D1S and D1N of the Gogny interaction. Although the considered effective NN interactions were proven to be quite realistic in numerous nuclear structure and/or reaction studies, they give quite different behaviors of the symmetry energy of nuclear matter at supranuclear densities that lead to the \emph{soft} and \emph{stiff} scenarios discussed recently in the literature. Different EOS's of the NS core and the EOS of the NS crust given by the compressible liquid drop model have been used as input of the Tolman-Oppenheimer-Volkov equations to study how the nuclear symmetry energy affects the model prediction of different NS properties, like the cooling process as well as the gravitational mass, radius, and moment of inertia.

The possibility to draw links between the isospin properties of nuclei and the structure of compact stars is a stimulating perspective. In order to pursue this objective on a sound basis, the correlations from which such links can be deduced have to be carefully checked against model dependence. Using a variety of nuclear effective models and a microscopic approach, we study the relation between the predictions of a given model and those of a Taylor density development of the corresponding equation of state: this establishes to what extent a limited set of phenomenological constraints can determine the core-crust transition properties. From a correlation analysis we show that a) the transition density $\rho_t$ is mainly correlated with the symmetry energy slope $L$, b) the proton fraction $Y_{p,t}$ with the symmetry energy and symmetry energy slope $(J,L)$ defined at saturation density, or, even better, with the same quantities defined at $\rho=0.1$ fm$^{-3}$, and c) the transition pressure $P_t$ with the symmetry energy slope and curvature $(J,K_{\rm sym})$ defined at $\rho=0.1$ fm$^{-3}$.

The slope of the nuclear symmetry energy at saturation density $L$ is pointed out as a crucial quantity to determine the mass and width of neutron-star crusts. This letter clarifies the relation between $L$ and the core-crust transition. We confirm that the transition density is soundly correlated with $L$ despite differences between models, and we propose a clear understanding of this correlation based on a generalised liquid drop model. Using a large number of nuclear models, we evaluate the dispersion affecting the correlation between the transition pressure $P_t$ and $L$. From a detailed analysis it is shown that this correlation is weak due to a cancellation between different terms. The correlation between the isovector coefficients $K_{\rm sym}$ and $L$ plays a crucial role in this discussion.

We discuss the thermalisation process of the neutron stars crust described by solving the heat transport equation with a microscopic input for the specific heat of baryonic matter. The heat equation is solved with initial conditions specific to a rapid cooling of the core. To calculate the specific heat of inner crust baryonic matter, i.e., nuclear clusters and unbound neutrons, we use the quasiparticle spectrum provided by the Hartree-Fock-Bogoliubov approach at finite temperature. In this framework we analyse the dependence of the crust thermalisation on pairing properties and on cluster structure of inner crust matter. It is shown that the pairing correlations reduce the crust thermalisation time by a very large fraction. The calculations show also that the nuclear clusters have a non-negligible influence on the time evolution of the surface temperature of the neutron star.

The inner crust of neutron stars, formed of a crystal lattice of uclear clusters immersed in a sea of unbound neutrons, may be the nique example of periodic nuclear systems. We have calculated the neutron specific heat in the shallow part of the crust using the band theory of solids with Skyrme nucleon-nucleon interactions. We have also tested the validity of various approximations. We have found that the neutron specific heat is well described by that of a Fermi gas, while the motion of the unbound neutrons is strongly affected by the nuclear lattice. These apparently contradictory results are explained by the particular properties of the neutron Fermi surface.

The inner crust of neutron stars formed of nuclear clusters immersed in a neutron sea has been widely studied in the framework of the Wigner-Seitz approximation since the seminal work of Negele and Vautherin. In this article, the validity of this approximation is discussed in the framework of the band theory of solids. For a typical cell of $^{200}$Zr, present in the external layers of the inner crust, it is shown that the ground state properties of the neutron gas are rather well reproduced by the Wigner-Seitz approximation, while its dynamical properties depend on the energy scale of the process of interest or on the temperature. It is concluded that the Wigner-Seitz approximation is well suited for describing the inner crust of young neutron stars and the collapsing core of massive stars during supernovae explosions. However the band theory is required for low temperature fluid dynamics.

The effects of the spin-orbit component of the particle-hole interaction on nuclear response functions and neutrino mean free path are examined. A complete treatment of the full Skyrme interaction in the case of symmetric nuclear matter and pure neutron matter is given. Numerical results for neutron matter are discussed. It is shown that the effects of the spin-orbit interaction remain small, even at momentum transfer larger than the Fermi momentum. The neutrino mean free paths are marginally affected.

It is shown that the large density fluctuations appearing at the onset of the first order nuclear liquid-gas phase transition can play an important role in the supernovae evolution. Due to these fluctuations, the neutrino gas may be trapped inside a thin layer of matter near the proto-neutron star surface. The resulting increase of pressure may induce strong particle ejection a few hundred milliseconds after the bounce of the collapse, contributing to the revival of the shock wave. The Hartree-Fock+RPA scheme, with a finite-range nucleon-nucleon effective interaction, is employed to estimate the effects of the neutrino trapping due to the strong density fluctuations, and to discuss qualitatively the consequences of the suggested new scenario.