Within the covariant density functional theory of hypernuclear matter we build a series of equations of state for hypernuclear compact stars, by calibrating the coupling constants of the $\Xi$-hyperon to the experimental binding energy of the single-$\Xi$ hypernuclei $^{15}_{\Xi^-}$C and $^{12}_{\Xi^-}$Be. Coupling constants of the $\Lambda$-hyperon to nucleons have been calibrated on a vast collection of experimental data on single $\Lambda$-hypernuclei and we employ those values. Uncertainties on the couplings of the $\Sigma$-hyperon to nuclear matter, due to lack of experimental data, are accounted for by allowing for a wide variation of the well depth of $\Sigma$ at rest in symmetric saturated nuclear matter. To account for uncertainties in the nucleonic sector at densities much larger than the saturation density, a rich collection of parametrizations is employed, some of them in agreement with existing constraints from nuclear physics and astrophysics. Neutron star properties are investigated with all these calibrated equations of state. The effects of the presence of hyperons on the radius, on the tidal deformability, on the moment of inertia, and on the nucleonic direct Urca process are discussed. The sensitivity of the hyperonic direct Urca processes to uncertainties in the nucleonic and hyperonic sectors is also addressed. It is shown that the relative variations of the radius, tidal deformability and moment of inertia from the values that characterize purely nucleonic stars are linearly correlated with the strangeness fraction. The maximum radius deviation, obtained for most massive neutron stars, is $\approx 10\%$. The reduction of the maximum mass, triggered by nucleation of strangeness, is estimated at $\approx 15 - 20\%$, out of which 5\% comes from insufficient information on the $\Sigma$-hyperon interactions.

We search for possible correlations between neutron star observables and thermodynamic quantities that characterize high density nuclear matter. We generate a set of model-independent equations of state describing stellar matter from a Taylor expansion around saturation density. Each equation of state which is a functional of the nuclear matter parameters is thermodynamically consistent, causal and compatible with astrophysical observations. We find that the neutron star tidal deformability and radius are strongly correlated with the pressure, the energy density and the sound velocity at different densities. Similar correlations are also exhibited by a large set of mean-field models based on non-relativistic and relativistic nuclear energy density functionals. These model independent correlations can be employed to constrain the equation of state at different densities above saturation from measurements of NS properties with multi-messenger observations. In particular, precise constraints on the radius of PSR J0030+0451 thanks to NICER observations would allow to better infer the properties of matter around two times the nuclear saturation density.

It is well-known that r-mode oscillations of rotating neutron stars may be unstable with respect to the gravitational wave emission. It is highly unlikely to observe a neutron star with the parameters within the instability window, a domain where this instability is not suppressed. But if one adopts the `minimal' (nucleonic) composition of the stellar interior, a lot of observed stars appear to be within the r-mode instability window. One of the possible solutions to this problem is to account for hyperons in the neutron star core. The presence of hyperons allows for a set of powerful (lepton-free) non-equilibrium weak processes, which increase the bulk viscosity, and thus suppress the r-mode instability. Existing calculations of the instability windows for hyperon NSs generally use reaction rates calculated for the $\Sigma^-\Lambda$ hyperonic composition via the contact $W$ boson exchange interaction. In contrast, here we employ hyperonic equations of state where the $\Lambda$ and $\Xi^-$ are the first hyperons to appear (the $\Sigma^-$'s, if they are present, appear at much larger densities), and consider the meson exchange channel, which is more effective for the lepton-free weak processes. We calculate the bulk viscosity for the non-paired $npe\mu\Lambda\Xi^-$ matter using the meson exchange weak interaction. A number of viscosity-generating non-equilibrium processes is considered (some of them for the first time in the neutron-star context). The calculated reaction rates and bulk viscosity are approximated by simple analytic formulas, easy-to-use in applications. Applying our results to calculation of the instability window, we argue that accounting for hyperons may be a viable solution to the r-mode problem.

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.

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.

In the present study we analyse the effect of the density dependence of the symmetry energy on the hyperonic content of neutron stars within a relativistic mean field description of stellar matter. For the $\Lambda$-hyperon, we consider parametrizations calibrated to $\Lambda$-hypernuclei. For the $\Sigma$ and $\Xi$-hyperons uncertainties that reflect the present lack of experimental information on $\Sigma$ and $\Xi$-hypernuclei are taken into account. We perform our study considering nuclear equations of state that predict two solar mass stars, and satisfy other well settled nuclear matter properties. The effect of the presence of hyperons on the radius, the direct Urca processes, and the cooling of accreting neutron stars are discussed. We show that some star properties are affected in a similar way by the density dependence of the symmetry energy and the hyperon content of the star. To disentangle these two effects it is essential to have a good knowledge of the equation of state at supra-saturation densities. The density dependence of the symmetry energy affects the order of appearance of the different hyperons, which may have direct implications on the neutron star cooling as different hyperonic neutrino processes processes may operate at the center of massive stars. For models which allow for the direct Urca process to operate, hyperonic and purely nucleonic ones are shown to have a similar luminosity when hyperons are included in agreement with modern experimental data. It is shown that for a density dependent hadronic model constrained by experimental, theoretical and observational data, the low-luminosity of SAX J$1808.4-3658$ can only be modelled for a hyperonic NS, suggesting that hyperons could be present in its core.

Constraints set on key parameters of the nuclear matter equation of state (EoS) by the values of the tidal deformability, inferred from GW170817, are examined by using a diverse set of relativistic and non-relativistic mean field models. These models are consistent with bulk properties of finite nuclei as well as with the observed lower bound on the maximum mass of neutron star $\sim 2 ~ {\rm M}_\odot$. The tidal deformability shows a strong correlation with specific linear combinations of the isoscalar and isovector nuclear matter parameters associated with the EoS. Such correlations suggest that a precise value of the tidal deformability can put tight bounds on several EoS parameters, in particular, on the slope of the incompressibility and the curvature of the symmetry energy. The tidal deformability obtained from the GW170817 and its UV/optical/infrared counterpart sets the radius of a canonical $1.4~ {\rm M}_{\odot}$ neutron star to be $11.82\leqslant R_{1.4}\leqslant13.72$ km.

Observations of the spin distribution of rapidly rotating neutron stars show evidence for a lack of stars spinning at frequencies larger than $f\approx 700$ Hz, well below the predictions of theoretical equations of state. This has generally been taken as evidence of an additional spin-down torque operating in these systems and it has been suggested that gravitational wave torques may be operating and be linked to a potentially observable signal. In this paper we aim to determine whether additional spin-down torques are necessary, or whether the observed limit of $f\approx 700$ Hz could correspond to the mass-shedding frequency for the observed systems and is simply a consequence of the, currently unknown, state of matter at high densities. Given our ignorance with regard to the true equation of state of matter above nuclear saturation densities, we make minimal physical assumption and only demand causality in the core. We then connect our causally-limited equation of state to a realistic microphysical crustal equation of state for densities below nuclear saturation density. This produces a limiting model that will give the lowest possible maximum frequency, which we compare to observational constraints on neutron star masses and frequencies. We also compare our findings with the constraints on the tidal deformability obtained in the observations of the GW170817 event. We find that the lack of pulsars spinning faster than $f\approx 700$ Hz is not compatible with our causal limited `minimal' equation of state, for which the breakup frequency cannot be lower than $f_{\rm max}\approx 1200$ Hz. A low frequency cutoff, around $f\approx 800$ Hz could only be possible if we assume that these systems do not contain neutron stars with masses above $M\approx 2 M_\odot$. This would have to be due either to selection effects, or possibly to a phase transition in the interior of the neutron star.

Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray binaries and neutron stars that undergo thermonuclear bursts are sensitive to mass and radius. The mass-radius relation of neutron stars depends on the equation of state that governs their interior. Constraining this relation accurately is thus of fundamental importance to understand the nature of dense matter. In this context we introduce a pipeline to calculate realistic model spectra of rotating neutron stars with hydrogen and helium atmospheres. An arbitrarily fast rotating neutron star with a given equation of state generates the spacetime in which the atmosphere emits radiation. We use the Lorene/nrotstar code to compute the spacetime numerically and the ATM24 code to solve the radiative transfer equations self-consistently. Emerging specific intensity spectra are then ray-traced through the neutron star's spacetime from the atmosphere to a distant observer with the Gyoto code. Here, we present and test our fully relativistic numerical pipeline. To discuss and illustrate the importance of realistic atmosphere models we compare our model spectra to simpler models like the commonly used isotropic color-corrected blackbody emission. We highlight the importance of considering realistic model-atmosphere spectra together with relativistic ray tracing to obtain accurate predictions. We also insist on the crucial impact of the star's rotation on the observables. Finally, we close a controversy that has been appearing in the literature in the recent years regarding the validity of the ATM24 code.

We model the thermal states of both isolated neutron stars and accreting neutron stars in X-ray transients in quiescence and confront them with observations. We use an equation of state and superfluid baryon gaps, which are consistently calculated. We conclude that the direct Urca process is required to be consistent with low-luminous accreting neutron stars. In addition, proton superfluidity and sufficiently weak neutron superfluidity are necessary to explain the cooling of middle-aged neutron stars and to obtain a realistic distribution of neutron star masses.

We derive analytic approximations for the neutrino luminosities and the heat capacities of neutron stars with isothernal nucleon cores as functions of the mass and radius of stars. The neutrino luminosities are approximated for the three basic neutrino emission mechanisms, and the heat capacities for the five basic combinations of the partial heat capacities. The approximations are valid for for a wide class of equations of state of dense nucleonmatter. The results significantly simplify the theoretical interpretation of observations of cooling neutron stars as well as of quasistationary thermal states of neutron stars in X-ray transients. For illustration, we present an analysis of the neutrino cooling functions of nine isolated neutron stars taking into account the effects of their magnetic fields and of the presence of light elements in their heat blanketing envelopes. These results allow one to investigate the superfluid properties of neutron star cores.

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.

(Abridged) The recent measurement of the mass of two $2\, M_\odot$ pulsars has raised the question whether such large masses allow for the existence of exotic degrees of freedom, such as hyperons, inside neutron stars. In the present work we will investigate how the existing hypernuclei properties may constrain the neutron star equation of state and confront the neutron star maximum masses obtained with equations of state calibrated to hypernuclei properties with the astrophysical $2\,M_\odot$ constraint. The study is performed using a relativistic mean field approach to describe both the hypernuclei and the neutron star equations of state. A set of five models consistent with $2\,M_\odot$ for a purely nucleonic composition are employed. The $\Lambda$-meson couplings are determined for all the models considered. Hyperonic stars with the complete baryonic octet are studied, restricting the coupling of the $\Sigma$ and $\Xi$ hyperons to the $\omega-$, $\rho-$ and $\sigma-$mesons due to the lack of experimental data, and maximum star masses calculated for unified equations of state. We conclude that the currently available hypernuclei experimental data and the lack of constraints on the asymmetric equation of state of nuclear matter at high densities do not allow to further constrain the neutron star matter equation of state using the recent $2\, M_\odot$ observations. It is also shown that the $\Lambda$ potential in symmetric nuclear matter takes a value $\sim 30-32$ MeV at saturation for the $g_{\omega \Lambda}$ coupling given by the SU(6) symmetry, close to values generally used in the literature. However, the $\Lambda$ potential in $\Lambda$ matter varies between -16 and -8 MeV taking for vector mesons couplings the SU(6) values, at variance with generally employed values between $-1$ and $-5$ MeV.

We explore the implications of a strong first-order phase transition region in the dense matter equation of state in the interiors of rotating neutron stars, and the resulting creation of two disjoint families of neutron-star configurations (the so-called high-mass twins). We numerically obtained rotating, axisymmetric, and stationary stellar configurations in the framework of general relativity, and studied their global parameters and stability. The instability induced by the equation of state divides stable neutron star configurations into two disjoint families: neutron stars (second family) and hybrid stars (third family), with an overlapping region in mass, the high-mass twin-star region. These two regions are divided by an instability strip. Its existence has interesting astrophysical consequences for rotating neutron stars. We note that it provides a natural explanation for the rotational frequency cutoff in the observed distribution of neutron star spins, and for the apparent lack of back-bending in pulsar timing. It also straightforwardly enables a substantial energy release in a mini-collapse to another neutron-star configuration (core quake), or to a black hole.

We derive analytic approximations of neutrino luminosities and heat capacities of neutron stars with nucleon cores valid for a wide class of equations of state of dense nucleon matter. The neutrino luminosities are approximated for the three cases in which they are produced by (i) direct Urca or (ii) modified Urca processes in non-superfluid matter, or (iii) neutrino-pair bremsstrahlung in neutron-neutron collisions (when other neutrino reactions are suppressed by strong proton superfluidity). The heat capacity is approximated for the two cases of (i) non-superfluid cores and (ii) the cores with strong proton superfluidity. The results can greatly simplify numerical simulations of cooling neutron stars with isothermal interiors at the neutrino and photon cooling stages as well as simulations of quasi-stationary internal thermal states of neutron stars in X-ray transients. For illustration, a model-independent analysis of thermal states of the latter sources is outlined.

Few unified equations of state for neutron star matter where core and crust are described using the same nuclear model are available. However the use of non-unified equations of state with a simplified matching between the crust and the core has been shown to introduce uncertainties in the radius determination which can be larger than the expected precision of the next generation of X-ray satellites. We aim at eliminating the dependence of the radius and mass of neutron staron the detailed model for the crust and on the crust-core matching procedure. We solve the approximate equations of the hydrostatic equilibrium for the crust of neutron stars obtaining a precise formula for the radius which depends only on the core mass and radius, and on the baryon chemical potential at the core-crust interface and on the crust surface. For a fully accreted crust one needs additionally the value of the total deep crustal heating per one accreted nucleon. For typical neutron star masses the approximate approach allows to determine the neutron star radius with an error ~0.1% (~ 10 m, equivalent to a 1% inaccuracy in the crust thickness). The formalism applies to neutron stars with a catalyzed or a fully accreted crust. The difference in the neutron star radius between the two models is proportional to the total energy release due to deep crustal heating. For a given model of dense matter describing the neutron star core, the radius of a neutron star can be accurately determined independently of the crust model with a precision much better than the ~5% one expected from the next generation of X-ray satellites. This allows to circumvent the problem of the radius uncertainty which may arise when non-unified equations of state for the crust and the core are used.

We examine the correlations of neutron star radii with the nuclear matter incompressibility, symmetry energy, and their slopes, which are the key parameters of the equation of state (EoS) of asymmetric nuclear matter. The neutron star radii and the EoS parameters are evaluated using a representative set of 24 Skyrme-type effective forces and 18 relativistic mean field models, and two microscopic calculations, all describing 2$M_\odot$ neutron stars. Unified EoSs for the inner-crust-core region have been built for all the phenomenological models, both relativistic and non-relativistic. Our investigation shows the existence of a strong correlation of the neutron star radii with the linear combination of the slopes of the nuclear matter incompressibility and the symmetry energy coefficients at the saturation density. Such correlations are found to be almost independent of the neutron star mass in the range $0.6\text{-}1.8M_{\odot}$. This correlation can be linked to the empirical relation existing between the star radius and the pressure at a nucleonic density between one and two times saturation density, and the dependence of the pressure on the nuclear matter incompressibility, its slope and the symmetry energy slope. The slopes of the nuclear matter incompressibility and the symmetry energy coefficients as estimated from the finite nuclei data yield the radius of a $1.4M_{\odot}$ neutron star in the range $11.09\text{-}12.86$ km.

We study the effects of heat blanketing envelopes of neutron stars on their cooling. To this aim, we perform cooling simulations using newly constructed models of the envelopes composed of binary ion mixtures (H--He, He--C, C--Fe) varying the mass of lighter ions (H, He or C) in the envelope. The results are compared with those calculated using the standard models of the envelopes which contain the layers of lighter (accreted) elements (H, He and C) on top of the Fe layer, varying the mass of accreted elements. The main effect is that the chemical composition of the envelopes influences their thermal conductivity and, hence, thermal insulation of the star. For illustration, we apply these results to estimate the internal temperature of the Vela pulsar and to study the cooling of neutron stars of ages of 0.1 - 1 Myr at the photon cooling stage. The uncertainties of the cooling models associated with our poor knowledge of chemical composition of the heat insulating envelopes strongly complicate theoretical reconstruction of the internal structure of cooling neutron stars from observations of their thermal surface emission.

The uncertainties in neutron star (NS) radii and crust properties due to our limited knowledge of the equation of state (EOS) are quantitatively analysed. We first demonstrate the importance of a unified microscopic description for the different baryonic densities of the star. If the pressure functional is obtained matching a crust and a core EOS based on models with different properties at nuclear matter saturation, the uncertainties can be as large as $\sim 30\%$ for the crust thickness and $4\%$ for the radius. Necessary conditions for causal and thermodynamically consistent matchings between the core and the crust are formulated and their consequences examined. A large set of unified EOS for purely nucleonic matter is obtained based on 24 Skyrme interactions and 9 relativistic mean-field nuclear parametrizations. In addition, for relativistic models 17 EOS including a transition to hyperonic matter at high density are presented. All these EOS have in common the property of describing a $2\;M_\odot$ star and of being causal within stable NS. A span of $\sim 3$ km and $\sim 4$ km is obtained for the radius of, respectively, $1.0\;M_\odot$ and $2.0\;M_\odot$ star. Applying a set of nine further constraints from experiment and ab-initio calculations the uncertainty is reduced to $\sim 1$ km and $2$ km, respectively. These residual uncertainties reflect lack of constraints at large densities and insufficient information on the density dependence of the EOS near the nuclear matter saturation point. The most important parameter to be constrained is shown to be the symmetry energy slope $L$ which exhibits a linear correlation with the stellar radius, particularly for masses $\sim 1.0\;M_\odot$. Potential constraints on $L$, the NS radius and the EOS from observations of thermal states of NS are also discussed. [Abriged]

A set of theoretical mass-radius relations for rigidly rotating neutron stars with exotic cores, obtained in various theories of dense matter, is reviewed. Two basic observational constraints are used: the largest measured rotation frequency (716 Hz) and the maximum measured mass ($2\;M_\odot$). Present status of measuring the radii of neutron stars is described. The theory of rigidly rotating stars in general relativity is reviewed and limitations of the slow rotation approximation are pointed out. Mass-radius relations for rotating neutron stars with hyperon and quark cores are illustrated using several models. Problems related to the non-uniqueness of the crust-core matching are mentioned. Limits on rigid rotation resulting from the mass-shedding instability and the instability with respect to the axisymmetric perturbations are summarized. The problem of instabilities and of the back-bending phenomenon are discussed in detail. Metastability and instability of a neutron star core in the case of a first-order phase transition, both between pure phases, and into a mixed-phase state, are reviewed. The case of two disjoint families (branches) of rotating neutron stars is discussed and generic features of neutron-star families and of core-quakes triggered by the instabilities are considered.

The recent mass measurements of two binary millisecond pulsars, PSR J1614-2230 and PSR J0751+1807 with a mass M=1.97+/-0.04 Msun and M= 1.26 +/- 0.14 Msun, respectively, indicate a wide range of masses for such objects and possibly also a broad spectrum of masses of neutron stars born in core-collapse supernovae. Starting from the zero-age main sequence binary stage, we aim at inferring the birth masses of PSR J1614-2230 and PSR J0751+1807 by taking the differences in the evolutionary stages preceding their formation into account. Using simulations for the evolution of binary stars, we reconstruct the evolutionary tracks leading to the formation of PSR J1614-2230 and PSR J0751+1807. We analyze in detail the spin evolution due to the accretion of matter from a disk in the intermediate-mass/low-mass X-ray binary. We consider two equations of state of dense matter, one for purely nucleonic matter and the other one including a high-density softening due to the appearance of hyperons. Stationary and axisymmetric stellar configurations in general relativity are used, together with a recent magnetic torque model and observationally-motivated laws for the decay of magnetic field. The estimated birth mass of the neutron stars PSR J0751+1807 and PSR J1614-2230 could be as low as 1.0 Msun and as high as 1.9 Msun, respectively. These values depend weakly on the equation of state and the assumed model for the magnetic field and its accretion-induced decay. The masses of progenitor neutron stars of recycled pulsars span a broad interval from 1.0 Msun to 1.9 Msun. Including the effect of a slow Roche-lobe detachment phase, which could be relevant for PSR J0751+1807, would make the lower mass limit even lower. A realistic theory for core-collapse supernovae should account for this wide range of mass.

The existence of 2 Msun pulsars puts very strong constraints on the equation of state (EOS) of neutron stars (NSs) with hyperon cores, which can be satisfied only by special models of hadronic matter. The radius-mass relation for these models is sufficiently specific that it could be subjected to an observational test with future X-ray observatories. We want to study the impact of the presence of hyperon cores on the radius-mass relation for NS. We aim to find out how, and for which particular stellar mass range, a specific relation R(M), where M is the gravitational mass, and R is the circumferential radius, is associated with the presence of a hyperon core. We consider a set of 14 theoretical EOS of dense matter, based on the relativistic mean-field (RMF) approximation, allowing for the presence of hyperons in NSs. We seek correlations between R(M) and the stiffness of the EOS below the hyperon threshold needed to pass the 2 Msun test. For NS masses 1.0<M/Msun<1.6, we get R>13km, because of a very stiff pre-hyperon segment of the EOS. At nuclear density, the pressure is significantly higher than a robust upper bound obtained recently using chiral effective field theory. If massive NSs do have a sizable hyperon core, then according to current models the radii for M=1.0-1.6 Msun are necessarily >13km. If, on the contrary, a NS with a radius R<12 km is observed in this mass domain, then sizable hyperon cores in NSs, as we model them now, are ruled out. Future X-ray missions with <5% precision for a simultaneous M and R measurement will have the potential to solve the problem with observations of NSs. Irrespective of this observational test, present EOS allowing for hyperons that fulfill condition M_max>2 Msun yield a pressure at nuclear density that is too high relative to up-to-date microscopic calculations of this quantity.

Recent measurement of a high millisecond pulsar mass (PSR J1614-2230, 1.97+-0.04 Msun) compared with the low mass of PSR J0751+1807 (1.26+-0.14 Msun) indicates a large span of masses of recycled pulsars and suggests a broad range of neutron stars masses at birth. We aim at reconstructing the pre-accretion masses for these pulsars while taking into account interaction of the magnetic field with a thin accretion disk, magnetic field decay and relativistic 2D solutions for stellar configurations for a set of equations of state. We briefly discuss the evolutionary scenarios leading to the formation of these neutron stars and study the influence of the equation of state.

Using the intrinsic PSR J1903+0327 parameters evaluated from radio observations (mass, rotation period and dipole magnetic field deduced from the timing properties) we calculate the mass of its neutron star progenitor, M_i, at the onset of accretion. Simultaneously, we derive constraints on average accretion rate Mdot and the pre-accretion magnetic field B_i. Spin-up is modelled by accretion from a thin disk, using the magnetic-torque disk-pulsar coupling model proposed by Kluzniak and Rappaport (2007), improved for the existence of relativistic marginally-stable circular orbit. Orbital parameters in the disk are obtained using the space-time generated by a rotating neutron star in the framework of General Relativity. We employ an observationally motivated model of the surface magnetic field decay. We also seek for the imprint of the poorly known equation of state of dense matter on the spin-up tracks - three equations of state of dense matter, consistent with the existence of 2 Msun neutron star, are considered. We find that the minimum average accretion rate should be larger than 2-8 10^(-10) Msun/yr, the highest lower bound corresponding to the stiffest equation of state. We conclude that the influence of magnetic field in the "recycling" process is crucial - it leads to a significant decrease of spin-up rate and larger accreted masses, in comparison to the B=0 model. Allowed B_i-dependent values of M_i are within 1.0-1.4 Msun, i.e., much lower than an oversimplified but widely used B=0 result, where one gets M_i>1.55 Msun. Estimated initial neutron-star mass depends on the assumed dense-matter equation of state.

To estimate the feasibility of dense-matter phase transition, we studied the evolution of the central density as well as the baryon chemical potential of accreting neutron stars. We compared the thin-disk accretion with and without the magnetic field torque with the spin-down scenario for a selection of recent equations of state. We compared the prevalent (in the recycled-pulsar context) Keplerian thin-disk model, in which the matter is accreted from the marginally-stable circular orbit, with the recent magnetic-torque model that takes into account the influence of stellar magnetic field on the effective inner boundary of the disk. Calculations were performed using a multi-domain spectral methods code in the framework of General Relativity. We considered three equations of state consistent with the recently measured mass of PSR J1614-2230, 1.97 +- 0.04 M_sun (one of them softened by the appearance of hyperons). If there is no magnetic torque and efficient angular momentum transfer from the disk to the star, substantial central compression is limited to the region of initial stellar masses close to the maximum mass. Outside the maximum mass vicinity, accretion-induced central compression is significant only if the angular momentum transfer is inefficient. Accounting for the magnetic field effectively decreases the efficiency of angular momentum transfer and implies a significant central compression. An efficient angular momentum transfer from a thin disk onto a non-magnetized neutron star does not provide a good mechanism for the central compression and possible phase transition. Substantial central compression is possible for a broad range of masses of slowly-rotating initial configurations for magnetized neutron stars. Accretion-induced central compression is particularly strong for stiff equation of state with a high-density softening.

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.