We discuss the formation of dark compact objects in a dark matter environment in view of the possible mass dependence of pulsars on the distribution of dark matter in the Galaxy. Our results indicate that the pulsar masses should decrease going towards the center of the Milky Way due to dark matter capture, thus becoming a probe for the existence and nature of dark matter. We thus propose that the evolution of the pulsar mass in a dark matter rich environment can be used to put constraints, when combined with future experiments, on the characteristics of our Galaxy halo dark matter profile, on the dark matter particle mass and on the dark matter self-interaction strength.

We show that the pulsar mass depends on the environment, and that it decreases going towards the center of the Milky Way. This is due to two combined effects, the capture and accumulation of self-interacting, non-annihilating dark matter by pulsars, and the increase of the dark matter density going towards the galactic center. We show that mass decrease depends both on the density profile of dark matter, steeper profiles producing a faster and larger decrease of the pulsar mass, and on the strength of self-interaction. Once future observations will provide the pulsar mass in a dark matter rich environment, close to the galactic center, the present result will be able to put constraints on the characteristics of our Galaxy halo dark matter profile, on the nature of dark matter, namely on its annihilating or non-annihilating nature, on its strength of self-interaction, and on the particle mass.

The theoretical status of the hyperon-nucleon and hyperon-hyperon interactions is reviewed, paying a special attention to chiral effective field theories. Results on hyperons in dense matter are presented and the consequences for the properties of neutron stars are analyzed.

We analyze the effects of including $\Delta(1232)$ isobars in an equation of state (EoS) for cold, $\beta$-stable neutron star matter, employing relativistic nuclear mean field theory. The selected EoS reproduces the properties of nuclear matter and finite nuclei and, in the astrophysical context, allows for the presence of hyperons in neutron stars having masses larger than 2$M_{\odot}$. We find that the composition and structure of neutron stars is critically influenced by the addition of the $\Delta$ isobars, which allows us to constrain their interaction with the meson fields taking into account astrophysical information. Imposing that the EoS is stable and ensures the existence of 2$M_{\odot}$ neutron stars, as well as requiring agreement with data of $\Delta$ excitation in nuclei, we find that, in the absence of other mechanisms stiffening the EoS at high densities, the interaction of the $\Delta$ isobars with the sigma and omega meson fields must be at least 10\% stronger than that of the nucleons. Moreover, the neutron star moment of inertia turns out to be sensitive to the presence of $\Delta$ isobars, whereas the inclusion of $\Delta$ isobars in the EoS allows for smaller stellar radii and for a lower value of the tidal deformability consistent with the analysis of the GW170817 merger event.

If a phase transition is allowed to take place in the core of a compact star, a new stable branch of equilibrium configurations can appear, providing solutions with the same mass as the purely hadronic branch and hence giving rise to twin-star configurations. We perform an extensive analysis of the features of the phase transition leading twin-star configurations and, at the same time, fulfilling the constraints coming from the maximum mass of $2M_\odot$ and the information following gravitational-wave event GW170817. In particular, we use a general equation of state for the neutron-star matter that parametrizes the hadron-quark phase transition between the model describing the hadronic phase and a constant speed of sound for the quark phase. We find that the largest number of twin-star solutions has masses in the neutron-star branch in the range $1-2M_\odot$ and twin-branch masses $\gtrsim 2M_\odot$. The analysis of the masses, radii and tidal deformabilities also reveals that when twin stars appear, the tidal deformability shows two distinct branches with the same mass, thus differing considerably from the behaviour expected for neutron stars. In addition, we find that the data from GW170817 is compatible with the existence of hybrid stars and could also be interpreted as produced by the merger of a binary system of hybrid stars or of a hybrid star with a neutron star. The presence of a hybrid star in the inspiral phase can be established clearly if future gravitational-wave detections measure chirp masses $\mathcal{M}\lesssim 1.2M_\odot$ and tidal deformabilities of $\Lambda_{1.4}\lesssim 400$ for $1.4M_\odot$ stars. Finally, combining all observational information available, we set constraints on the parameters that characterise the phase transition, the maximum masses, and the radii of $1.4M_\odot$ stars described by equations of state leading to twin-star configurations.

We present two recent parametrizations of the equation of state (FSU2R and FSU2H models) that reproduce the properties of nuclear matter and finite nuclei, fulfill constraints on high-density matter stemming from heavy-ion collisions, produce 2$M_{\odot}$ neutron stars, and generate neutron star radii below 13 km. Making use of these equations of state, cooling simulations for isolated neutron stars are performed. We find that two of the models studied, FSU2R (with nucleons) and, in particular, FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including nucleon pairing. This indicates that cooling observations are compatible with an equation of state that produces a soft nuclear symmetry energy and, thus, generates small neutron star radii. Nevertheless, both schemes produce cold isolated neutron stars with masses above $1.8 M_{\odot}$.

We study the structure of compact objects that contain non-self annihilating, self-interacting dark matter admixed with ordinary matter made of neutron star and white dwarf materials. We extend the previous work Phys. Rev. D 92 123002 (2015) on these dark compact objects by analyzing the effect of weak and strongly interacting dark matter with particle masses in the range of 1-500 GeV, so as to set some constraints in the strength of the interaction and the mass of the dark matter particle. We find that the total mass of the compact objects increases with decreasing dark matter particle mass. In the strong interacting case and for dark matter particle masses in the range 1-10 GeV, the total mass of the compact objects largely exceeds the $2M_\odot$ constraint for neutron star masses and the nominal $1M_\odot$for white dwarfs, while for larger dark matter particle masses or in the weakly interacting case the compact objects show masses in agreement or smaller than these constraints, thus hinting at the exclusion of strongly self-interacting dark matter of masses 1-10 GeV in the interior of these compact objects. Moreover, we observe that the smaller the dark matter particle mass, the larger the quantity of dark matter captured is, putting constraints on the dark matter mass trapped in the compact objects so as to fullfill $\simeq 2M_\odot$ observations. Finally, the inhomogeneity of distribution of dark matter in the Galaxy implies a mass dependence of compact objects from the environment which can be used to put constraints on the characteristics of the Galaxy halo DM profile and on particle mass. In view of the these results, we discuss the formation of the dark compact objects in an homogeneous and non-homogeneous dark matter environment.

We construct parameter sets of the relativistic mean-field model fitted to the recent constraints on the asymmetry energy $J$ and the slope parameter $L$ for pure neutron matter. We find cases of unphysical behaviour, i.e.\ the appearance of negative pressures, for stiff parameter sets with low values of the effective mass $m^*/m$. In some cases the equation of state of pure neutron matter turns out to be outside the allowed band given by chiral effective field theory. The mass-radius relations of neutron stars for all acceptable parameter sets shows a maximum mass in excess of $2M_\odot$ being compatible with pulsar mass measurements. Given the constraints on the model in the low-density regime coming from chiral effective theory, we find that the radius of a $1.4M_\odot$ neutron star is nearly independent on the value of $L$. This is in contrast to some previous claims for a strong connection of the slope parameter with the radius of a neutron star. In fact, the mass-radius relation turns out to depend only on the isoscalar parameters of symmetric matter. The constraints of GW170817 on the tidal deformability and on the radius are also discussed.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

We perform cooling simulations for isolated neutron stars using recently developed equations of state for their core. The equations of state are obtained from new parametrizations of the FSU2 relativistic mean-field functional that reproduce the properties of nuclear matter and finite nuclei, while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions, measurements of massive 2$M_{\odot}$ neutron stars, and neutron star radii below 13 km. We find that two of the models studied, FSU2R (with nucleons) and in particular FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including extensive nucleon pairing. This suggests that the cooling observations are more compatible with an equation of state that produces a soft nuclear symmetry energy and, hence, generates small neutron star radii. However, both models favor large stellar masses, above $1.8 M_{\odot}$, to explain the colder isolated neutron stars that have been observed, even if nucleon pairing is present.

We reexamine the equation of state for the nucleonic and hyperonic inner core of neutron stars that satisfies the 2$M_{\odot}$ observations as well as the recent determinations of stellar radii below 13 km, while fulfilling the saturation properties of nuclear matter and finite nuclei together with the constraints on the high-density nuclear pressure coming from heavy-ion collisions. The recent nucleonic FSU2R and hyperonic FSU2H models are fine tuned improving the density dependence of pure neutron matter at subsaturation densities. The corresponding nuclear matter properties at saturation, the symmetry energy and its slope turn out to be compatible with recent experimental and theoretical determinations. We obtain the mass, radius and composition of neutron stars for the two updated models and study the impact on these properties of the uncertainties in the hyperon-nucleon couplings estimated from hypernuclear data. We find that the onset of appearance of each hyperon strongly depends on the hyperon-nuclear uncertainties, whereas the maximum masses for neutron stars differ by at most 0.1 $M_{\odot}$, although a larger deviation should be expected tied to the lack of knowledge of the hyperon potentials at the high densities present in the center of $2 M_\odot$ stars. For easier use, we provide tables with the results from the FSU2R and FSU2H models for the equation of state and the neutron star mass-radius relation.

We present new stable solutions of the Tolman Oppenheimer Volkoff equations for quark stars using a quark matter equation of state based on the SU(3) Quark-Meson model that exhibits the onset of the chiral phase transition. These new solutions appear as two stable branches in the mass-radius relation allowing for so called twin stars, i.e. two stable quark star solutions with the same mass, but distinctly different radii. We find solutions which are compatible with causality, the stability conditions of dense matter, the astrophysical constraints of the rotation of the millisecond pulsar PSR J1748-2446ad and the 2$M_{\odot}$ pulsar mass constraint.

We obtain a new equation of state for the nucleonic and hyperonic inner core of neutron stars that fulfills the 2$M_{\odot}$ observations as well as the recent determinations of stellar radii below 13 km. The nucleonic equation of state is obtained from a new parametrization of the FSU2 relativistic mean-field functional that satisfies these latest astrophysical constraints and, at the same time, reproduces the properties of nuclear matter and finite nuclei while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions. On the one hand, the equation of state of neutron star matter is softened around saturation density, which increases the compactness of canonical neutron stars leading to stellar radii below 13 km. On the other hand, the equation of state is stiff enough at higher densities to fulfill the 2$M_{\odot}$ limit. By a slight modification of the parametrization, we also find that the constraints of 2$M_{\odot}$ neutron stars with radii around 13 km are satisfied when hyperons are considered. The inclusion of the high magnetic fields present in magnetars further stiffens the equation of state. Hyperonic magnetars with magnetic fields in the surface of $ \sim 10^{15}$ G and with values of $\sim 10^{18}$ G in the interior can reach maximum masses of 2$M_{\odot}$ with radii in the 12-13 km range.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

We investigate compact objects formed by dark matter admixed with ordinary matter made of neutron star matter and white dwarf material. We consider non-self annihilating dark matter with an equation-of-state given by an interacting Fermi gas. We find new stable solutions, dark compact planets, with Earth-like masses and radii from few Km to few hundred Km for weakly interacting dark matter. For the strongly interacting dark matter case, we obtain dark compact planets with Jupiter-like masses and radii of few hundred Km. These objects could be formed primordially and accrete white dwarf material subsequently. They could be detected by observing exoplanets with unusually small radii. Moreover, we find that the recently observed 2 ${\rm M}_{\odot}$ pulsars set limits on the amount of dark matter inside neutron stars which is, at most, $10^{-6}{\rm M}_\odot$.

We study the shear and bulk viscosity coefficients as well as the thermal conductivity as arising from the collisions among phonons in superfluid neutron stars. We use effective field theory techniques to extract the allowed phonon collisional processes, written as a function of the equation of state and the gap of the system. The shear viscosity due to phonon scattering is compared to calculations of that coming from electron collisions. We also comment on the possible consequences for r-mode damping in superfluid neutron stars. Moreover, we find that phonon collisions give the leading contribution to the bulk viscosities in the core of the neutron stars. We finally obtain a temperature-independent thermal conductivity from phonon collisions and compare it with the electron-muon thermal conductivity in superfluid neutron stars.

We compute the contribution of phonons to the thermal conductivity in the core of superfluid neutron stars. We use effective field theory techniques to extract the phonon scattering rates, written as a function of the equation of state of the system. We also calculate the phonon dispersion law beyond linear order, which depends on the gap of superfluid neutron matter. With all these ingredients, we solve the Boltzmann equation numerically using a variational approach. We find that the thermal conductivity $\kappa$ is dominated by combined small and large angle binary collisions. As in the color-flavor-locked superfluid, we find that our result can be well approximated by $\kappa \propto 1/ \Delta^6$, where $\Delta$ is the neutron gap, the constant of proportionality depending on the density. We further comment on the possible relevance of electron and superfluid phonon collisions in obtaining the total contribution to the thermal conductivity in the core of superfluid neutron stars.

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

We analyze how recent computations of the shear viscosity $\eta$ in the core of superfluid neutron stars affect the r-mode instability window. We first analyze the contribution of superfluid phonons to the viscosity, both in their hydrodynamical and ballistic regime. We also consider the recent computation of $\eta$ arising from the collisions of electrons with electrons and protons by Shternin and Yakovlev, and discuss how the interactions among superfluid phonons and electrons might contribute to the shear viscosity. For assessing the r-mode instability window we compare the shear viscosity due to phonons in the hydrodynamical regime with respect to the shear viscosity due to electron collisions. Only at high temperatures the superfluid phonon contribution to $\eta$ starts to dominate the process of r-mode damping. While our results for the instability window are preliminary, as other dissipative processes should be taken into account as well, they differ from previous evaluations of the r-mode damping due to the shear viscosity in superfluid neutron stars.

We calculate the three bulk viscosity coefficients as arising from the collisions among phonons in superfluid neutron stars. We use effective field theory techniques to extract the allowed phonon collisional processes, written as a function of the equation of state of the system. The solution of the dynamical evolution of the phonon number density allows us to calculate the bulk viscosity coefficients as function of the phonon collisional rate and the phonon dispersion law, which depends on the neutron pairing gap. Our method of computation is rather general, and could be used for different superfluid systems, provided they share the same underlying symmetries. We find that the behavior with temperature of the bulk viscosity coefficients is dominated by the contributions coming from the collinear regime of the $2\leftrightarrow3$ phonon processes. For typical star radial pulsation frequencies of $\omega \sim 10^{4} s^{-1}$, we obtain that the bulk viscosity coefficients at densities $n \gtrsim 4n_0$ are within 10% from its static value for $T \lesssim 10^9 K$ and for the case of strong neutron superfluidity in the core with a maximum value of the $^3P_2$ gap above 1 MeV, while, otherwise, the static solution is not a valid approximation to the bulk viscosity coefficients. Compared to previous results from Urca and modified Urca reactions, we conclude that at $T\sim 10^9$K phonon collisions give the leading contribution to the bulk viscosities in the core of the neutron stars, except for $n \sim 2n_0$ when the opening of the Urca processes takes place.

We study the implications on compact star properties of a soft nuclear equation of state determined from kaon production at subthreshold energies in heavy-ion collisions. On one hand, we apply these results to study radii and moments of inertia of light neutron stars. Heavy-ion data provides constraints on nuclear matter at densities relevant for those stars and, in particular, to the density dependence of the symmetry energy of nuclear matter. On the other hand, we derive a limit for the highest allowed neutron star mass of three solar masses. For that purpose, we use the information on the nucleon potential obtained from the analysis of the heavy-ion data combined with causality on the nuclear equation of state.

Measurements of kaon production at subthreshold energies in heavy-ion collisions point to a soft nuclear equation-of-state for densities up to 2-3 times nuclear matter saturation density. We apply these results to study the implications on compact star properties, especially in the context of the recent measurement of the two solar mass pulsar PSR J1614-2230. The implications are two-fold: Firstly, the heavy-ion results constrain nuclear matter at densities relevant to light neutron stars. Hence, a radius measurement could provide information about the density dependence of the symmetry energy which is a crucial quantity in nuclear physics. Secondly, the information on the nucleon potential obtained from the analysis of the heavy-ion data can be combined with restrictions from causality on the nuclear equation-of-state. From this we can derive a limit for the highest allowed compact star mass of three solar masses.

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

The phonon contribution to the shear viscosity $\eta$ in superfluid neutron stars is calculated by assuming neutron pairing in a $^1S_0$ channel. The shear viscosity is obtained by means of variational methods for the solution of the Boltzmann equation amended by a collision term which takes into account the binary collisions of phonons. Effective field theory techniques are used to extract the phonon scattering rates in terms of the equation of state (EoS) of the system. We find that $\eta \propto 1/T^5$, the proportionality factor depending on the EoS of the system. Our results indicate that the phonon contribution to $\eta$ might have important effects for the different oscillation modes of the star.

We compute the contribution of phonons to the shear viscosity $\eta$ in superfluid neutron stars, assuming neutron pairing in a $^1S_0$ channel. We use a Boltzmann equation amended by a collision term that takes into account the binary collisions of phonons. We use effective field theory techniques to extract the phonon scattering rates, written as a function of the equation of state (EoS) of the system. Our formulation is rather general, and can be used to extract the shear viscosity due to binary collisions of phonons for other superfluids, such as the cold Fermi gas in the unitarity limit. We find that $\eta \propto 1/T^5$, the proportionality factor depending on the EoS of the system. Our results indicate that the phonon contribution to $\eta$ cannot be ignored and might have relevant effects in the dynamics of the different oscillation modes of the star.

We calculate the neutron matter equation of state at finite temperature based on low-momentum two- and three-nucleon interactions. The free energy is obtained from a loop expansion around the Hartree-Fock energy, including contributions from normal and anomalous diagrams. We focus on densities below saturation density with temperatures T <= 10 MeV and compare our results to the model-independent virial equation of state and to variational calculations. Good agreement with the virial equation of state is found at low density. We provide simple estimates for the theoretical error, important for extrapolations to astrophysical conditions.

We calculate the neutron matter equation of state at finite temperature based on low-momentum nucleon-nucleon and three-nucleon interactions. Our results are compared to the model-independent virial equation of state and to variational calculations. We provide a simple estimate for the theoretical error, important for extrapolations to astrophysical conditions.