The observed thermal relaxation of transiently accreting neutron stars during quiescence periods in low-mass X-ray binaries suggests the existence of unknown heat sources in the shallow layers of neutron-star crust. Making use of existing experimental nuclear data, we estimate the maximum possible amount of heat that can be deposited in the outer crust of an accreting neutron star due to electron captures and pycnonuclear fusion reactions triggered by the burial of X-ray burst ashes.

Using a Wigner-Seitz approximation with spherical cells, we re-analyze a widely used single nucleus model of accreting neutron star crust. We calculate beta disequilibrium within the crust, which is sizable, and implies that neutron and baryon chemical potentials, $\mu_n$ and $\mu_{\rm b}$, are not equal. We include also non-equilibrium reactions, driven by matter compression, and proceeding in the reaction layers. The constancy of $e^{\Phi}\mu_n$, where the spacetime metric component $g_{00}=e^{2\Phi}$, in the shells between the reaction layers is not applicable, because single electron captures are blocked, so that the neutron fraction is fixed, and therefore neutrons are not an independent component of the crust matter. The absence of neutron diffusion in the shells between the reaction layers, stems from the constancy of the neutron fraction (concentration) in these shells. In the reaction layers, the outward force resulting from neutron fraction gradient is balanced by the inward gravitational force acting on unbound neutrons. Neglecting the thickness of the reaction layers compared to the shell thickness, we obtain condition $e^{\Phi(r)}f_Q(r)g(r)=$constant, where $g$ is Gibbs energy per nucleon, undergoing discontinuous drops on the reaction surfaces, and $f_Q(r)g(r)={\widetilde{g}}(r)$ is a continuous function, due to the factor $f_Q(r)$ canceling the discontinuities (drops) in $g(r)$. The function $f_Q(r)$ is calculated using the Tolman-Oppenheimer-Volkov equations from $f_Q(P)$ and $g(P)$ obtained from the equation of state (EOS) with discontinuites. The constancy of of $e^{\Phi(r)}{\widetilde{g}}(r)$ is an extension of the standard relation $e^{\Phi(r)}{{\mu}}_{\rm b}(r)=$constant, valid in hydrostatic equilibrium for catalyzed crust.

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.

We study the bulk viscosity in two (anti-spaghetti and Swiss cheese) phases of non- spherical nuclei in the mantle of an oscillating neutron star near the boundary with the stellar core. The bulk viscosity is produced by non-equilibrium Urca neutrino emis- sion processes. In the mantle, the direct Urca process may be open (Gusakov et al., 2004) if neutrons and protons move in a periodic potential created by a lattice of non- spherical nuclei (which allows the nucleons to have large quasi-momenta and satisfy direct Urca momentum-conservation). This bulk viscosity can dominate over that due to the modified Urca process in the outer stellar core and over the shear viscosity. The bulk viscosity depends strongly on temperature, oscillation frequency and nucleon superfluidity. The enhanced bulk viscosity in the mantle can control propagation and damping of neutron star oscillations.

We calculate modified Urca neutrino emission rates in the dense nuclear matter in neutron star cores. We find that these rates are strongly enhanced in the beta-stable matter in regions of the core close to the direct Urca process threshold. This enhancement can be tracked to the use of the in-medium nucleon spectrum in the virtual nucleon propagator. We describe the in-medium nucleon scattering in the non-relativistic Bruckner-Hartree-Fock framework taking into account two-body as well as the effective three-body forces, although the proposed enhancement does not rely on a particular way of the nucleon interaction treatment. Finally we suggest a simple approximate expression for the emissivity of the n-branch of the modified Urca process that can be used in the neutron stars cooling simulations with any nucleon equation of state of dense matter.

Observations of soft X-ray transients in quiescence suggest the existence of heat sources in the crust of accreting neutron stars. Heat is thought to be released by electroweak and nuclear processes triggered by the burying of ashes of X-ray bursts. The heating is studied using a fully quantum approach taking consistently into account nuclear shell effects. We have followed the evolution of ashes made of $^{56}$Fe employing the nuclear energy-density functional theory. Both the outer and inner crusts are described using the same functional, thus ensuring a unified and thermodynamically consistent treatment. To assess the role of the neutron-matter constraint, we have employed the set of accurately calibrated Brussels-Montreal functionals BSk19, BSk20, and BSk21 and for comparison the SLy4 functional. Due to nuclear shell effects, the fully accreted crust is found to be much less stratified than in previous studies. In particular, large regions of the inner crust contain clusters with the magic number $Z=14$. The heat deposited in the outer crust is tightly constrained by experimental atomic mass data. The shallow heating we obtain does not exceed $0.2$~MeV and is therefore not enough to explain the cooling of some soft X-ray transients. The total heat released in the crust is very sensitive to details of the nuclear structure and is predicted to lie in the range from $1.5$~MeV to $1.7$~MeV. The evolution of an accreted matter element and therefore the location of heat sources are governed to a large extent by the existence of nuclear shell closures. Ignoring these effects in the inner crust, the total heat falls to $\sim 0.6$~MeV. The neutron-matter constraint is also found to play a key role. The large amount of heat obtained by Steiner et al. (2012) could thus be traced back to unrealistic neutron-matter equations of state.

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.

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

Investigations of thermal evolution of neutron stars with hyperon cores require neutrino emissivities for many neutrino reactions involving strongly degenerate particles (nucleons, hyperons, electrons, muons). We calculate the angular integrals $I_n$ (over orientations of momenta of $n$ degenerate particles) for major neutrino reactions with $n$ =3, 4, 5 at all possible combinations of particle Fermi momenta. The integrals $I_n$ are necessary ingredients for constructing a uniform database of neutrino emissivities in dense nucleon-hyperon matter. The results can also be used in many problems of physical kinetics of strongly degenerate systems.

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.

We estimate the energy reservoir available in the deconfinement phase transition induced collapse of a neutron star to its hybrid star mass twin on the "third family" branch, using a recent equation of state of dense matter. The available energy corresponding to the mass-energy difference between configurations is comparable with energies of the most violent astrophysical burst processes. An observational outcome of such a dynamical transition might be fast radio bursts, specifically a recent example of a FRB with a double-peak structure in its light curve.

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.

The neutron-drip transition in the dense matter constituting the interior of neutron stars generally refers to the appearance of unbound neutrons as the matter density reaches some threshold density $\rho_\textrm{drip}$. This transition has been mainly studied under the cold catalyzed matter hypothesis. However, this assumption is unrealistic for accreting neutron stars. After examining the physical processes that are thought to be allowed in both accreting and nonaccreting neutron stars, suitable conditions for the onset of neutron drip are derived and general analytical expressions for the neutron drip density and pressure are obtained. Moreover, we show that the neutron-drip transition occurs at lower density and pressure than those predicted within the mean-nucleus approximation. This transition is studied numerically for various initial composition of the ashes from X-ray bursts and superbursts using microscopic nuclear mass models.

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.

The recent measurement of two solar mass pulsars has initiated an intense discussion on its impact on our understanding of the high-density matter in the cores of neutron stars. A task force meeting was held from October 7-10, 2013 at the Frankfurt Institute for Advanced Studies to address the presence of quark matter in these massive stars. During this meeting, the recent oservational astrophysical data and heavy-ion data was reviewed. The possibility of pure quark stars, hybrid stars and the nature of the QCD phase transition were discussed and their observational signals delineated.

Observations of massive ($M \approx 2.0~M_\odot$) neutron stars (NSs), PSRs J1614-2230 and J0348+0432, rule out most of the models of nucleon-hyperon matter employed in NS simulations. Here we construct three possible models of nucleon-hyperon matter consistent with the existence of $2~M_\odot$ pulsars as well as with semi-empirical nuclear matter parameters at saturation, and semi-empirical hypernuclear data. Our aim is to calculate for these models all the parameters necessary for modelling dynamics of hyperon stars (such as equation of state, adiabatic indices, thermodynamic derivatives, relativistic entrainment matrix, etc.), making them available for a potential user. To this aim a general non-linear hadronic Lagrangian involving $\sigma\omega\rho\phi\sigma^\ast$ meson fields, as well as quartic terms in vector-meson fields, is considered. A universal scheme for calculation of the $\ell=0,1$ Landau Fermi-liquid parameters and relativistic entrainment matrix is formulated in the mean-field approximation. Use of this scheme allow us to obtain numerical tables with the equation of state, Landau quasiparticle effective masses, adiabatic indices, the $\ell=0,1$ Landau Fermi-liquid parameters, and the relativistic entrainment matrix for the selected models of nucleon-hyperon matter. These data are available on-line and suitable for numerical implementation in computer codes modelling various dynamical processes in NSs, in particular, oscillations of superfluid NSs and their cooling.

One of the most intringuing questions about neutron stars concerns their maximum mass. The answer is intimately related to the properties of matter at densities far beyond that found in heavy atomic nuclei. The current view on the internal constitution of neutron stars and on their maximum mass, both from theoretical and observational studies, are briefly reviewed.

We calculate thermal conductivity and shear viscosity of nucleons in dense nuclear matter of neutron star cores in the non-relativistic Brueckner-Hartree-Fock framework. Nucleon-nucleon interaction is described by the Argonne v18 potential with addition of the Urbana IX effective three-body forces. We find that this three body force leads to decrease of the kinetic coefficients with respect to the two-body case. The results of calculations are compared with electron and muon transport coefficients as well as with the results of other authors.

The problem of interaction of the rotating magnetic field, frozen to a star, with a thin well conducting accretion disk is solved exactly. It is shown that a disk pushes the magnetic field lines towards a star, compressing the stellar dipole magnetic field. At the point of corotation, where the Keplerian rotation frequency coincides with the frequency of the stellar rotation, the loop of the electric current appears. The electric currents flow in the magnetosphere only along two particular magnetic surfaces, which connect the corotation region and the inner edge of a disk with the stellar surface. It is shown that the closed current surface encloses the magnetosphere. Rotation of a disk is stopped at some distance from the stellar surface, which is 0.55 of the corotation radius. Accretion from a disk spins up the stellar rotation. The angular momentum transferred to the star is determined.

Recent measurement of mass of PSR J1614-2230 rules out most of existing models of equation of state (EOS) of dense matter with high-density softening due to hyperonization, based on the recent hyperon-nucleon and hyperon-hyperon interactions, leading to a "hyperon puzzle". We study a specific solution of "hyperon puzzle", consisting in replacing a too soft hyperon core by a sufficiently stiff quark core. We construct an analytic approximation fitting very well modern EOSs of 2SC and CFL color superconducting phases of quark matter. This allows us for simulating continua of sequences of first-order phase transitions from hadronic matter to the 2SC, and then to the CFL state of color superconducting quark matter. We obtain constraints in the parameter space of the EOS of superconducting quark cores, resulting from M_max> 2 M_sol. We also derive constraints that would result from significantly higher measured masses. For 2.4 M_sol required stiffness of the CFL quark core should have been close to the causality limit, the density jump at the phase transition being very small. Condition M_max > 2 M_sol puts strong constraints on the EOSs of the 2SC and CFL phases of quark matter. Density jumps at the phase transitions have to be sufficiently small and sound speeds in quark matter - sufficiently large. A strict condition of thermodynamic stability of quark phase results in the maximum mass of hybrid stars similar to that of purely baryon stars. Therefore, to get M_max>2 M_sol for stable hybrid stars, both sufficiently strong additional hyperon repulsion at high density baryon matter and a sufficiently stiff EOS of quark matter would be needed. However, it is likely that the high density instability of quark matter (reconfinement) indicates actually the inadequacy of the point-particle model of baryons in dense matter at very high densities.

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.

We review the history of neutron star physics in the 1930s that is related to L. Landau. According to recollections of Rosenfeld (1974, Proc. 16th Solvay Conference on Physics, p. 174), Landau improvised the concept of neutron stars in a discussion with Bohr and Rosenfeld just after the news of the discovery of the neutron reached Copenhagen in February 1932. We present arguments that the discussion took place in March 1931, before the discovery of the neutron, and that they in fact discussed the paper written by Landau in Zurich in February 1931 but not published until February 1932 (Phys. Z. Sowjetunion, 1, 285). In his paper Landau mentioned the possible existence of dense stars which look like one giant nucleus; this can be regarded as an early theoretical prediction or anticipation of neutron stars, prior to the discovery of the neutron. The coincidence of the dates of the neutron's discovery and the paper's publication has led to an erroneous association of the paper with the discovery of the neutron. In passing, we outline the contribution of Landau to the theory of white dwarfs and to the hypothesis of stars with neutron cores.

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.

Recent measurement of mass of PSR J1614-2230 rules out most of existing models of equation of state (EOS) of dense matter with high-density softening due to hyperonization or a phase transition to quark matter or a boson condensate. We look for a solution of an apparent contradiction between the consequences stemming from up-to-date hypernuclear data, indicating appearance of hyperons at 3 nuclear densities and existence of a two-solar-mass neutron star. We consider a non-linear relativistic mean field (RMF) model involving baryon octet coupled to meson fields. An effective lagrangian includes quartic terms in the meson fields. The values of the parameters of the model are obtained by fitting semi-empirical parameters of nuclear matter at the saturation point, as well as potential wells for hyperons in nuclear matter and the strength of the Lambda-Lambda attraction in double-Lambda hypernuclei. We propose a non-linear RMF model which is consistent with up-to-date semiempirical nuclear and hypernuclear data and allows for neutron stars with hyperon cores and M larger than 2 solar masses. The model involves hidden-strangenes scalar and vector mesons, coupled to hyperons only, and quartic terms involving vector meson fields. Our EOS involving hyperons is stiffer than the corresponding nucleonic EOS (with hyperons artificially suppressed) above five nuclear densities. Required stiffening is generated by the quartic terms involving hidden-strangeness vector meson.

Neutron star crust, formed via accretion of matter from a companion in a low-mass X-ray binary (LMXB), has an equation of state (EOS) stiffer than that of catalyzed matter. At a given neutron star mass, M, the radius of a star with an accreted crust is therefore larger, by DR(M), than for usually considered star built of catalyzed matter. Using a compressible liquid drop model of nuclei, we calculate, within the one-component plasma approximation, the EOSs corresponding to different nuclear compositions of ashes of X-ray bursts in LMXB. These EOSs are then applied for studying the effect of different formation scenarios on the neutron-star mass-radius relation. Assuming the SLy EOS for neutron star's liquid core, derived by Douchin & Haensel (2001), we find that at M=1.4 M_sun the star with accreted crust has a radius more than 100 m larger that for the crust of catalyzed matter. Using smallness of the crust mass compared to M, we derive a formula that relates DR(M) to the difference in the crust EOS. This very precise formula gives also analytic dependence of DR on M and R of the reference star built of catalyzed matter. The formula is valid for any EOS of the liquid core. Rotation of neutron star makes DR(M) larger. We derive an approximate but very precise formula that gives difference in equatorial radii, DR_eq(M), as a function of stellar rotation frequency.

We calculate stationary configurations of rapidly rotating compact stars in general relativity, to study the properties of circular orbits of test particles in the equatorial plane. We search for simple, but precise, analytical formulae for the orbital frequency, specific angular momentum and binding energy of a test particle, valid for any equation of state and for any rotation frequency of the rigidly rotating compact star, up to the mass-shedding limit. Numerical calculations are performed using precise 2-D codes based on multi-domain spectral methods. Models of rigidly rotating neutron stars and the space-time outside them are calculated for several equations of state of dense matter. Calculations are also performed for quark stars consisting of self-bound quark matter. At the mass-shedding limit, the rotational frequency converges to a Schwarzschildian orbital frequency at the equator. We show that orbital frequency for any orbit outside equator is also approximated by a Schwarzschildian formula. Using a simple approximation for the frame-dragging term, we obtain approximate expressions for the specific angular momentum and specific energy on the corotating circular orbits in the equatorial plane of neutron star, which are valid down to the stellar equator. The formulae recover reference numerical values with typically 1% of accuracy for neutron stars with M > 0.5 M_sun. They are less precise for quark stars consisting of self-bound quark matter.

We calculate the important quantity of superfluid hydrodynamics, the relativistic entrainment matrix for a nucleon-hyperon mixture at arbitrary temperature. In the nonrelativistic limit this matrix is also termed the Andreev-Bashkin or mass-density matrix. Our results can be useful for modeling the pulsations of massive neutron stars with superfluid nucleon-hyperon cores and for studies of the kinetic properties of superfluid baryon matter.

We calculate the relativistic entrainment matrix Y_ik at zero temperature for nucleon-hyperon mixture composed of neutrons, protons, Lambda and Sigma^- hyperons, as well as of electrons and muons. This matrix is analogous to the entrainment matrix (also termed mass-density matrix or Andreev-Bashkin matrix) of non-relativistic theory. It is an important ingredient for modelling the pulsations of massive neutron stars with superfluid nucleon-hyperon cores. The calculation is done in the frame of the relativistic Landau Fermi-liquid theory generalized to the case of superfluid mixtures; the matrix Y_ik is expressed through the Landau parameters of nucleon-hyperon matter. The results are illustrated with a particular example of the sigma-omega-rho mean-field model with scalar self-interactions. Using this model we calculate the matrix Y_ik and the Landau parameters. We also analyze stability of the ground state of nucleon-hyperon matter with respect to small perturbations.

Using numerical simulations based on solving the general relativistic hydrodynamic equations, we study the dynamics of a phase transition in the dense core of isolated rotating neutron stars, triggered by the back bending instability reached via angular momentum loss. In particular, we investigate the dynamics of a migration from an unstable configuration into a stable one, which leads to a mini-collapse of the neutron star and excites sizeable pulsations in its bulk until it acquires a new stable equilibrium state. We consider equations of state with softening at high densities, a simple analytic one with a mixed hadron-quark phase in an intermediate pressure interval and pure quark matter at very high densities, and a microphysical one that has a first-order phase transition, originating from kaon condensation. Although the marginally stable initial models are rigidly rotating, we observe that during the collapse (albeit little) differential rotation is created. We analyze the emission of gravitational radiation, which in some models is amplified by mode resonance effects, and assess its prospective detectability by interferometric detectors. We expect that the most favorable conditions for dynamic migration exist in very young magnetars. We find that the damping of the post-migration pulsations strongly depends on the character of the equation of state softening. The damping of pulsations in the models with the microphysical equation of state is caused by dissipation associated with matter flowing through the density jump at the edge of the dense core. If at work, this mechanism dominates over all other types of dissipation, like bulk viscosity in the exotic-phase core, gravitational radiation damping, or numerical viscosity.

We discuss whether the winding-up of the magnetic field by differential rotation in a new-born quark star can produce a sufficiently-high, energy, emission rate of sufficiently long duration to explain long gamma-ray bursts. In the context of magnetohydrodynamics, we study the torsional oscillations and energy extraction from a new-born, hot, differentially rotating quark star. The new-born compact star is a rapid rotator that produces a relativistic, leptonic wind. The star's torsional oscillation modulates this wind emission considerably when it is odd and of sufficient amplitude, which is relatively easy to reach. Odd oscillations may occur just after the formation of a quark star. Other asymmetries can cause similar effects. The buoyancy of wound-up magnetic fields is inhibited, or its effects are limited, by a variety of different mechanisms. Direct electromagnetic emission by the torsional oscillation in either an outside vacuum or the leptonic wind surrounding the compact object is found to be insignificant. In contrast, the twist given to the outer magnetic field by an odd torsional oscillation is generally sufficient to open the star's magnetosphere. The Poynting emission of the star in its leptonic environment is then radiated from all of its surface and is enhanced considerably during these open episodes, tapping at the bulk rotational energy of the star. This results in intense energy shedding in the first tens of minutes after the collapse of magnetized quark stars with an initial poloidal field of order of 10**14 Gauss, sufficient to explain long gamma-ray bursts.

We calculate Keplerian (mass shedding) configurations of rigidly rotating neutron stars and quark stars with crusts. We check the validity of empirical formula for Keplerian frequency, f_K, proposed by Lattimer & Prakash, f_K(M)=C (M/M_sun)^1/2 (R/10km)^-3/2, where M is the (gravitational) mass of Keplerian configuration, R is the (circumferential) radius of the non-rotating configuration of the same gravitational mass, and C = 1.04 kHz. Numerical calculations are performed using precise 2-D codes based on the multi-domain spectral methods. We use a representative set of equations of state (EOSs) of neutron stars and quark stars. We show that the empirical formula for f_K(M) holds within a few percent for neutron stars with realistic EOSs, provided 0.5 M_sun < M < 0.9 M_max,stat, where M_max,stat is the maximum allowable mass of non-rotating neutron stars for an EOS, and C=C_NS=1.08 kHz. Similar precision is obtained for quark stars with 0.5 M_sun < M < 0.9 M_max,stat. For maximal crust masses we obtain C_QS = 1.15 kHz, and the value of C_QS is not very sensitive to the crust mass. All our C's are significantly larger than the analytic value from the relativistic Roche model, C_Roche = 1.00 kHz. For 0.5 M_sun < M < 0.9 M_max,stat, the equatorial radius of Keplerian configuration of mass M, R_K(M), is, to a very good approximation, proportional to the radius of the non-rotating star of the same mass, R_K(M) = aR(M), with a_NS \approx a_QS \approx 1.44. The value of a_QS is very weakly dependent on the mass of the crust of the quark star. Both a's are smaller than the analytic value a_Roche = 1.5 from the relativistic Roche model.

The physics of neutron star crusts is vast, involving many different research fields, from nuclear and condensed matter physics to general relativity. This review summarizes the progress, which has been achieved over the last few years, in modeling neutron star crusts, both at the microscopic and macroscopic levels. The confrontation of these theoretical models with observations is also briefly discussed.

We calculate parameters A and B of the Baym-Pines model of the hydro-elastic equilibrium of rotating neutron stars. Parameter A determines the energy increase of a non-rotating star due to a quadrupolar deformation of its shape. Parameter B determines residual quadrupolar deformation due to the crustal shear strain, in a neutron star that spun-down to a non-rotating state. The calculations of A are based on precise numerical 2-D calculations for rotating neutron stars with realistic equations of state (EOSs) of dense matter. An approximate, but quite precise, formula for B is used, which allows us to separate the contribution of the crust from the dependence on the stellar mass M and radius R. The elastic shear strain distribution within the crust is modeled following Cutler et al. (2003). A(M) and B(M) are calculated for 0.2Msun < M < 0.9M_max for seven EOSs of neutron star core, combined with several crust models. A standard formula based on the incompressible fluid model is shown to severely underestimate the value of A. For M<0.7Msun A(M) is nearly EOS independent and are given (within a few percent) by a universal formula A=3.87(M/Msun)^7/3 [10^53 erg]. We derive the scaling of B with respect to R and M, valid also for a thick crust. We show that B for accreted crust strongly depends on pycnonuclear fusions at rho>10^12 g/cm^3.

Fast rotation of compact stars (at submillisecond period) and, in particular, their stability, are sensitive to the equation of state (EOS) of dense matter. Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz (Kaaret et al. 2007). At such rotational frequency the effects of rotation on star's structure are significant. We study the interplay of fast rotation, EOS and gravitational mass of a submillisecond pulsar. We discuss the EOS dependence of spin-up to a submillisecond period, via mass accretion from a disk in a low-mass X-ray binary.

We analyze the thermal conductivity of ions (equivalent to the conductivity of phonons in crystalline matter) in a neutron star envelope. We calculate the ion/phonon thermal conductivity in a crystal of atomic nuclei using variational formalism and performing momentum-space integration by Monte Carlo method. We take into account phonon-phonon and phonon-electron scattering mechanisms and show that phonon-electron scattering dominates at not too low densities. We extract the ion thermal conductivity in ion liquid or gas from literature. Numerical values of the ion/phonon conductivity are approximated by analytical expressions, valid for T>10^5 K and 10^5 g cm^-3 < \rho < 10^14 g cm^-3. Typical magnetic fields B~10^12 G in neutron star envelopes do not affect this conductivity although they strongly reduce the electron thermal conductivity across the magnetic field. The ion thermal conductivity remains much smaller than the electron conductivity along the magnetic field. However, in the outer neutron star envelope it can be larger than the electron conductivity across the field, that is important for heat transport across magnetic field lines in cooling neutron stars. The ion conductivity can greatly reduce the anisotropy of heat conduction in outer envelopes of magnetized neutron stars.

Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz. Such rotational frequency would be the first for which the effects of rotation are significant. We study the consequences of very fast rotating neutron stars for the potentially observable quantities as stellar mass and pulsar period.

We present the general relativistic calculation of the energy release associated with a first order phase transition (PT) at the center of a rotating neutron star (NS). The energy release, E_rel, is equal to the difference in mass-energies between the initial (normal) phase configuration and the final configuration containing a superdense matter core, assuming constant total baryon number and the angular momentum. The calculations are performed with the use of precise pseudo-spectral 2-D numerical code; the polytropic equations of state (EOS) as well as realistic EOSs (Skyrme interactions, Mean Field Theory kaon condensate) are used. The results are obtained for a broad range of metastability of initial configuration and size of the new superdense phase core in the final configuration. For a fixed ``overpressure'', dP, defined as the relative excess of central pressure of a collapsing metastable star over the pressure of the equilibrium first-order PT, the energy release up to numerical accuracy does not depend on the stellar angular momentum and coincides with that for nonrotating stars with the same dP. When the equatorial radius of the superdense phase core is much smaller than the equatorial radius of the star, analytical expressions for the E_rel can be obtained: E_rel is proportional to dP^2.5 for small dP. At higher dP, the results of 1-D calculations of E_rel(dP) for non-rotating stars reproduce with very high precision exact 2-D results for fast-rotating stars. The energy release-angular momentum independence for a given overpressure holds also for the so-called ``strong'' PTs (that destabilise the star against the axi-symmetric perturbations), as well as for PTs with ``jumping'' over the energy barrier.

Heating associated with non-equilibrium nuclear reactions in accreting neutron-star crusts is reconsidered, taking into account suppression of neutrino losses demonstrated recently by Gupta et al. Two initial compositions of the nuclear burning ashes, A=56 and A=106, are considered. Dependence of the integrated crustal heating on uncertainties plaguing pycnonuclear reaction models is studied. One-component plasma approximation is used, with compressible liquid-drop model of Mackie and Baym to describe nuclei. Evolution of a crust shell is followed from 10^8 g/cm^3 to 10^(13.6) g/cm^3 The integrated heating in the outer crust agrees nicely with results of self-considtent multicomponent plasma simulations of Gupta et al.; their results fall between our curves obtained for A=56 and A=106. Total crustal heat per one accreted nucleon ranges between 1.5 MeV to 1.9 MeV for A=106 and A=56, respectively. The value of total crustal heat per nucleon depends weakly on the presence of pycnonuclear reactions at densities 10^(12)-10^(13) g/cm^3. Remarkable insensitivity of the total crustal heat on the details of the distribution of nuclear processes in accreted crust is explained.

We simulate the cooling of the neutron star in the X-ray transient KS 1731-260 after the source returned to quiescence in 2001 from a long (>~ 12.5 yr) outburst state. We show that the cooling can be explained assuming that the crust underwent deep heating during the outburst stage. In our best theoretical scenario the neutron star has no enhanced neutrino emission in the core, and its crust is thin, superfluid, and has the normal thermal conductivity. The thermal afterburst crust-core relaxation in the star may be not over.

We calculate the energy release associated with a strong first-order phase transition, from normal phase N to an "exotic" superdense phase S, in a rotating neutron star. Such a phase transition, accompanied by a density jump rho_N --> rho_S, is characterized by rho_S/rho_N > 3/2(1+P_0/rho_N c^2), where P_0 is the pressure, at which phase transition occurs. Configurations with small S-phase cores are then unstable and collapse into stars with large S-phase cores. The energy release is equal to the difference in mass-energies between the initial (normal) configuration and the final configuration containing an S-phase core, total stellar baryon mass and angular momentum being kept constant. The calculations of the energy release are based on precise numerical 2-D calculations. Polytropic equations of state (EOSs) as well as realistic EOS with strong first-order phase transition due to kaon condensation are used. For polytropic EOSs, a large parameter space is studied. For a fixed "overpressure", dP, defined as the relative excess of central pressure of collapsing metastable star over the pressure of equilibrium first-order phase transition, the energy release E_rel does not depend on the stellar angular momentum. It coincides with that for nonrotating stars with the same dP. Therefore, results of 1-D calculations of E_rel(dP) for non-rotating stars can be used to predict, with very high precision, the outcome of much harder to perform 2-D calculations for rotating stars with the same dP. This result holds also for dP_min < dP < 0, corresponding to phase transitions with climbing over the energy barrier separating metastable N-phase configurations from those with an S-phase core. Such phase transitions could be realized in the cores of newly born, hot, pulsating neutron stars.

The analysis of SN1987A led Brown and Bethe (1995) to conclusion, that the maximum mass of cold neutron stars is low, M_max ~ 1.5M_sun. Such a low M_max, due to a kaon condensation in the stellar core, implies collapse of a too massive deleptonized protoneutron star into a black hole. This would naturally explain the lack of a neutron star in the SN1987A remnant. On the other hand, recent evaluation of mass of PSR J0751+1807 gives M_max > 2M_sun. This contradicts the original Bethe-Brown model, but can be reconciled within scenarios proposed in the present Letter. We consider two types of dense matter models with high-density softening, due to a transition from a non-strange N-phase of matter to a strangeness carrying phase S: kaon condensation and deconfinement of quarks. Two scenarios of neutron star formation in stellar core collapse are considered. In the first scenario, realized in sufficiently hot and dense supernova cores, nucleation of an S-phase is sufficiently rapid so as to form an S-phase core, and implying M_max = M^S_max =~ 1.5M_sun. In the second scenario, nucleation of the S-phase at neutron star birth is too slow to materialize, and the star becomes cold without forming an S-phase core. Then, stellar mass can increase via accretion, until central density \rho_crit is reached, and the S phase forms. This N branch of neutron stars ends at M=M_crit. We select several models of N-phase satifying the necessary condition M^N_max > 2M_sun and combine them with models of kaon condensation and quark deconfinement. For kaon condensation, we get M_crit =~ M^S_max =~ 1.5M_sun, which is ruled out by PSR J0751+1807. On the contrary, for the EOSs with quark deconfinement we get M_crit =~ M^N_max > 2M_sun, which reconciles SN1987A and PSR J0751+1807.

Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz. Such rotation imposes bounds on the structure of neutron star in XTE J1739-285. These bounds may be used to constrain poorly known equation of state of dense matter. One-parameter families of stationary configurations rotating rigidly at 1122 Hz are constructed, using a precise 2-D code solving Einstein equations. Hydrostatic equilibrium solutions are tested for stability with respect to axi-symmetric perturbations. A set of ten diverse EOSs of neutron stars is considered. Hypothetical strange stars are also studied. For each EOS, the family of possible neutron star models is limited by the mass shedding limit, corresponding to maximum allowed equatorial radius, R_max, and by the instability with respect to the axi-symmetric perturbations, reached at the minimum allowed equatorial radius, R_min. We get R_min \simeq 10-13km, and R_max \simeq 16-18km, with allowed mass 1.4-2.3 M_\odot. Allowed stars with hyperonic or exotic-phase core are supramassive and have a very narrow mass range. Quark star with accreted crust might be allowed, provided such a model is able to reproduce X-ray bursts from XTE J1739-285.

If strange quark matter (SQM) is the true ground state of hadronic matter, then conversion of neutron stars (NS) into quark stars (QS) could release some 10^{53} erg. We describe a scenario of burning of a NS into a hot, differentially rotating QS. Emission of released non-baryonic energy through the QS surface is discussed. The role of magnetobuoyancy of SQM is mentioned. The outflow of gamma e+ e- lasting for up to ~1000 s could be at the origin of long GRBs. Advantages of hot, differentially rotating QS as an inner engine of long GRBs are reviewed.

The properties of superdense matter in neutron star (NS) cores control NS thermal states by affecting the efficiency of neutrino emission from NS interiors. To probe these properties we confront the theory of thermal evolution of NSs with observations of their thermal radiation. Our observational basis includes cooling isolated NSs (INSs) and NSs in quiescent states of soft X-ray transients (SXTs). We find that the data on SXTs support the conclusions obtained from the analysis of INSs: strong proton superfluidity with T_{cp,max} >= 10^9 K should be present, while mild neutron superfluidity with T_{cn,max} =(2*10^8 -- 2*10^9) K is ruled out in the outer NS core. Here T_{cn,max} and T_{cp,max} are the maximum values of the density dependent critical temperatures of neutrons and protons. The data on SXTs suggest also that: (i) cooling of massive NSs is enhanced by neutrino emission more powerful than the emission due to Cooper pairing of neutrons; (ii) mild neutron superfluidity, if available, might be present only in inner cores of massive NSs. In the latter case SXTs would exhibit dichotomy, i.e. very similar SXTs may evolve to very different thermal states.

We study possible impact of a softening of the equation of state by a phase transition, or appearance of hyperons, on the spin evolution of of isolated pulsars. Numerical simulations are performed using exact 2-D simulations in general relativity. The equation of state of dense matter at supranuclear densities is poorly known. Therefore, the accent is put on the general correlations between evolution and equation of state, and mathematical strictness. General conjectures referring to the structure of the one-parameter families of stationary configurations are formulated. The interplay of the back bending phenomenon and stability with respect to axisymmetric perturbations is described. Changes of pulsar parameters in a corequake following instability are discussed, for a broad choice of phase transitions predicted by different theories of dense matter. The energy release in a corequake, at a given initial pressure, is shown to be independent of the angular momentum of collapsing configuration. This result holds for various types of phases transition, with and without metastability. We critically review observations of pulsars that could be relevant for the detection of the signatures of the phase transition in neutron star cores.

We calculate energy release associated with a first order phase transition at the center of a rotating neutron star. The results are based on precise numerical 2-D calculations, in which both the polytropic equations of state (EOS) as well as realistic EOS of the normal phase are used. Presented results are obtained for a broad range of metastability of initial configuration and size of the new superdense phase core in the final configuration. For small radii of the superdense phase core analytical expressions for the energy release are obtained. For a fixed "overpressure" dP (the relative excess of central pressure of collapsing metastable star over the pressure of equilibrium first-order phase transition) the energy release remarkably does not depend on the stellar angular momentum and coincides with that for nonrotating stars with the same dP. The energy release is proportional to dP^2.5 for small dPs, when sufficiently precise brute force 2-D numerical calculations are out of question. For higher dPs, results of 1-D calculations of energy release for non-rotating stars are shown to reproduce, with very high precision, the exact 2-D results for rotating stars.

Hydrodynamical simulations of neutron star cores, based on a two fluid description in terms of a neutron-proton superfluid mixture, require the knowledge of the Andreev-Bashkin entrainment matrix which relates the momentum of one constituent to the currents of both constituents. This matrix is derived for arbitrary nuclear asymmetry at zero temperature and in the limits of small relative currents in the framework of the energy density functional theory. The Skyrme energy density functional is considered as a particular case. General analytic formulae for the entrainment parameters and various corresponding effective masses are obtained. These formulae are applied to the liquid core of a neutron star, composed of an homogeneous plasma of nucleons, electrons and possibly muons in beta equilibrium.

Phenomena such as pulsar frequency glitches are believed to be attributable to differential rotation of a current of ``free'' superfluid neutrons at densities above the ``drip'' threshold in the ionic crust of a neutron star. Such relative flow is shown to be locally describable by adaption of a canonical two fluid treatment that emphasizes the role of the momentum covectors constructed by differentiation of action with respect to the currents, with allowance for stratification whereby the ionic number current may be conserved even when the ionic charge number Z is altered by beta processes. It is demonstrated that the gauge freedom to make different choices of the chemical basis determining which neutrons are counted as ``free'' does not affect their ``superfluid'' momentum covector, which must locally have the form of a gradient (though it does affect the ``normal'' momentum covector characterising the protons and those neutrons that are considered to be ``confined'' in the nuclei). It is shown how the effect of ``entrainment'' (whereby the momentum directions deviate from those of the currents) is controlled by the (gauge independent) mobility coefficient K, estimated in recent microscopical quantum mechanical investigations, which suggest that the corresponding (gauge dependent) ``effective mass'' m* of the free neutrons can become very large in some layers. The relation between this treatment of the crust layers and related work (using different definitions of ``effective mass'') intended for the deeper core layers is discussed.

The back-bending phenomenon for compact stars is studied by means of analytical equations of state, for both constant-pressure phase transitions and the transitions through the mixed-phase region. We restrict ourselves to the case of normal rotating configurations, with baryon mass below the maximum allowable baryon mass for non-rotating stars. We use high-precision 2-D multi-domain spectral code LORENE to search the parameter space for possible instability regions, and possible changes in the stability character of rotating stars with phase transitions in their cores. Conditions on the density jump in constant-pressure phase transitions, leading to the existence of the unstable segments in the evolutionary sequences of spinning down isolated normal neutron stars, are derived. Conjectures concerning the existence of two disjoint families of non-rotating and rotating stationary configurations of neutron stars are formulated. Particular case of EOSs leading to marginal instability of static and rotating configurations is also studied: marginal instability point in non-rotating configurations continues to exist in all evolutionary spin-down tracks. The fate of rotating stars entering the region of instability is discussed. The change in radius, energy release, and spin-up associated with the corequake in rotating neutron star, triggered by the instability, are calculated. The energy release is found to be very weakly dependent on the angular momentum of collapsing star.

The influence of nucleons superfluidity on the beta relaxation time of degenerate neutron star cores, composed of neutrons, protons and electrons, is investigated. We numerically calculate the implied reduction factors for both direct and modified Urca reactions, with isotropic pairing of protons or anisotropic pairing of neutrons. We find that due to the non-zero value of the temperature and/or to the vanishing of anisotropic gaps in some directions of the phase-space, superfluidity does not always completely inhibit beta relaxation, allowing for some reactions if the superfluid gap amplitude is not too large in respect to both the typical thermal energy and the chemical potential mismatch. We even observe that if the ratio between the critical temperature and the actual temperature is very small, a suprathermal regime is reached for which superfluidity is almost irrelevant. On the contrary, if the gap is large enough, the composition of the nuclear matter can stay frozen for very long durations, unless the departure from beta equilibrium is at least as important as the gap amplitude. These results are crucial for precise estimation of the superfluidity effect on the cooling/slowing-down of pulsars and we provide online subroutines to be implemented in codes for simulating such evolutions.

The moment of inertia of the pulsar A in the neutron star binary J0737-3039 will soon be measurable through detailed measurements of the periastron advance. We present the calculation of the moment of inertia of neutron stars with the masses of the components of the binary J0737-3039 for a broad range of equations of state of dense matter and discuss the implications of such measurement for constraining the equation of state. An observational determination of the moment of inertia of the pulsar A in J0737-3039 with the accuracy of 10% shall narrow down considerably the range of viable equations of state. We also show that limits on maximal mass of a neutron star provide a complementary set of constraints on the properties of dense nuclear matter.

Long gamma-ray bursts (GRBs) are believed to be related to the explosion of type Ic supernovae, which have been stripped of their hydrogen and helium envelopes. There appear to be two types of these explosions: those which are approximately spherical (GRB980425/1998bw), and which are associated with weak bursts, and the classical GRBs which generate ultrarelativistic jets (GRB030329/SN2003dh). If this bimodality is real Swift will provide a clear evidence for it. We propose that classical powerful GRBs, which generate ultrarelativistic outflows, are a result of a formation of quark stars. Quark stars may provide an additional energy for the explosion of SN Ic, but far more important is a creation of a surface which acts as a membrane which cannot be penetrated by baryons. A surface of a quark star allows only ultrarelativistic matter to escape: photons, neutrinos, electron -- positron pairs and magnetic fields. The formation of a quark star follows the initial core collapse in several minutes. Possible evidence for this time delay is provided by BATSE precursors to GRBs, as analyzed by Lazzati (2005).

The relative current density $n^i$ of conduction neutrons in a neutron star crust beyond the neutron drip threshold can be expected to be related to the corresponding particle momentum covector $p_i$ by a linear relation of the form $n^i=K^{ij}p_j$ in terms of a physically well defined mobility tensor $K^{ij}$. This result is describable as an ``entrainment'' whose effect - wherever the crust lattice is isotropic - will simply be to change the ordinary neutron mass m to an effective mass $m_\star$ such that in terms of the relevant number density n of unconfined neutrons we shall have $K^{ij}=n/m_\star\gamma^{ij}$. In a preceding analysis based on an independent particle treatment using Bloch boundary conditions to obtain the distribution of energy $E_k$ and associated group velocity $v_k^i=\partial E_k/\partial\hbar k_i$ as a function of wavenumber $k_i$, it was shown that the mobility tensor would be given by $K^{ij}\propto\int d^3 k v_k^i v_k^j\delta\{E_k-\mu\}$, where $\mu$ is the Fermi energy. Using the approach due to Bogoliubov, it is shown here that the effect of BCS pairing with an energy gap \Delta_F and corresponding quasiparticle energy function \EE_k=\sqrt{(E_k-\mu)^2+\Delta_F^2}$ will just be to replace the Dirac distributional integrand by a smoother distribution : $K^{ij}\propto \int d^3 k v_k^i v_k^j \Delta_F^2 /\EE^3_k$. It is also shown how the pairing condensation gives rise to superfluidity in the technical sense of providing (meta) stability against resistive perturbations for a current that is not too strong (its momentum $p_i$ must be small enough to give $2|p_i v_k^i|<\EE_k^2/|E_k-\mu|$ for all modes). It is concluded that the prediction of a very large effective mass enhancement in the middle layers of the crust will not be significantly affected by the pairing mechanism.

The entrainment matrix (also termed the Andreev-Bashkin matrix or the mass-density matrix) for a neutron-proton mixture is derived at a finite temperature in a neutron star core. The calculation is performed in the frame of the Landau Fermi-liquid theory generalized to account for superfluidity of nucleons. It is shown, that the temperature dependence of the entrainment matrix is described by a universal function independent on an actual model of nucleon-nucleon interaction employed. The results are presented in the form convenient for their practical use. The entrainment matrix is important, e.g., in kinetics of superfluid nucleon mixtures or in studies of the dynamical evolution of neutron stars (in particular, in the studies of star pulsations and pulsar glitches).

We study the spin-up of the accreting neutron stars with a realistic hyperon-softened equation of state. Using precise 2-D calculations we study the evolutionary tracks of accreting neutron stars in the angular-momentum - frequency plane. In contrast to the case of spinning-down solitary radio-pulsars, where a strong back-bending behavior has been observed, we do not see back-bending phenomenon in the accretion-powered spinning-up case. We conclude that in the case of accretion-driven spin-up the back-bending is strongly suppressed by the mass-increase effect accompanying the angular-momentum increase.

With (non-barotropic) equations of state valid even when the neutron, proton and electron content of neutron star cores is not in beta equilibrium, we study inertial and composition gravity modes of relativistic rotating neutron stars. We solve the relativistic Euler equations in the time domain with a three dimensional numerical code based on spectral methods, in the slow rotation, relativistic Cowling and anelastic approximations. Principally, after a short description of the gravity modes due to smooth composition gradients, we focus our analysis on the question of how the inertial modes are affected by non-barotropicity of the nuclear matter. In our study, the deviation with respect to barotropicity results from the frozen composition of non-superfluid matter composed of neutrons, protons and electrons, when beta equilibrium is broken by millisecond oscillations. We show that already for moderatly fast rotating stars the increasing coupling between polar and axial modes makes those two cases less different than for very slowly rotating stars. In addition, as we directly solve the Euler equations, without coupling only a few number of spherical harmonics, we always found, for the models that we use, a discrete spectrum for the $l = m = 2$ inertial mode. Finally, we find that, for non-barotropic stars, the frequency of this mode, which is our main focus, decreases in a non-negligible way, whereas the time dependence of the energy transfer between polar and axial modes is substantially different due to the existence of low-frequencies gravity modes.

We present approximate formulae describing the changes in neutron star parameters caused by the first-order phase transition to an ``exotic'' state (pion or kaon condensate, quark matter) resulting in formation of a mixed-phase core. The analytical formulae for the changes in radius, moment of inertia and the amount of energy released during the core-quake are derived using the theory of linear response of stellar structure to the core transformation. Numerical coefficients in these formulae are obtained for two realistic EOSs of dense matter. The problem of nucleation of the exotic phase as well as possible astrophysical scenarios leading to a core-quake phenomenon is also discussed.

In the inner crust of a neutron star, at densities above the ``drip'' threshold, unbound ``conduction'' neutrons can move freely past through the ionic lattice formed by the nuclei. The relative current density $n^i= n \bar v^i$ of such conduction neutrons will be related to the corresponding mean particle momentum $p_i$ by a proportionality relation of the form $n^i= {\cal K}p^i$ in terms of a physically well defined mobility coefficient $\cal K$ whose value in this context has not been calculated before. Using methods from ordinary solid state and nuclear physics, a simple quantum mechanical treatment based on the independent particle approximation, is used here to formulate $\cal K$ as the phase space integral of the relevant group velocity over the neutron Fermi surface. The result can be described as an ``entrainment'' that changes the ordinary neutron mass m to a macroscopic effective mass per neutron that will be given -- subject to adoption of a convention specifying the precise number density n of the neutrons that are considered to be ``free'' -- by $m_\star=n/{\cal K}$. The numerical evaluation of the mobility coefficient is carried out for nuclear configurations of the ``lasagna'' and ``spaghetti'' type that may be relevant at the base of the crust. Extrapolation to the middle layers of the inner crust leads to the unexpected prediction that $m_\star$ will become very large compared with m.

We present calculations of quasiequilibrium sequences of irrotational binary neutron stars based on realistic equations of state (EOS) for the whole neutron star interior. Three realistic nuclear EOSs of various softness and based on different microscopic models have been joined with a recent realistic EOS of the crust, giving in this way three different EOSs of neutron-star interior. Computations of quasiequilibrium sequences are performed within the Isenberg-Wilson-Mathews approximation to general relativity. For all evolutionary sequences, the innermost stable circular orbit (ISCO) is found to be given by mass-shedding limit (Roche lobe overflow). The EOS dependence on the last orbits is found to be quite important: for two 1.35 M_sol neutron stars, the gravitational wave frequency at the ISCO (which marks the end of the inspiral phase) ranges from 800 Hz to 1200 Hz, depending upon the EOS. Detailed comparisons with 3rd order post-Newtonian results for point-mass binaries reveals a very good agreement until hydrodynamical effects (dominated by high-order functions of frequency) become important, which occurs at a frequency ranging from 500 Hz to 1050 Hz, depending upon the EOS.

Analytical representations are derived for two equations of state (EOSs) of neutron-star matter: FPS and SLy. Each of these EOSs is unified, that is, it describes the crust and the core of a neutron star using the same physical model. Two versions of the EOS parametrization are considered. In the first one, pressure and mass density are given as functions of the baryon density. In the second version, pressure, mass density, and baryon density are given as functions of the pseudo-enthalpy, which makes this representation particularly useful for 2-D calculations of stationary rotating configurations of neutron stars.

We show that the direct Urca process of neutrino emission is allowed in two possible phases of nonspherical nuclei (inverse cylinders and inverse spheres) in the mantle of a neutron star near the crust-core interface. The process is open because neutrons and protons move in a periodic potential created by inhomogeneous nuclear structures. In this way the nucleons acquire large quasimomenta needed to satisfy momentum-conservation in the neutrino reaction. The appropriate neutrino emissivity in a nonsuperfluid matter is about 2--3 orders of magnitude higher than the emissivity of the modified Urca process in the stellar core. The process may noticeably accelerate the cooling of low-mass neutron stars.

The upper bound on the value of the surface gravity, g_s, for neutron stars with equations of state respecting v_sound <= c, is derived. This bound is inversely proportional to the maximum allowable mass M_max, and it reads g_s <= 1.411 x 10^15 (M_sun/M_max) cm/s^2. It implies an absolute upper bound 7.4 x 10^14 cm/s^2 if one uses the 2sigma lower bound on the neutron mass measured recently in 4U1700-37, 1.9 M_sun. A correlation between g_s and the compactness parameter 2GM/Rc^2 for baryonic stars is analyzed. The properties of g_s of strange quark stars and its upper bounds are discussed using the scaling properties of the strange-star models.

We use Landau's theory of a normal Fermi liquid to derive expressions for the static response of a system with a general tensor interaction that conserves the total spin and the total angular momentum of the quasiparticle-quasihole pair. The magnetic susceptibility is calculated in detail, with the inclusion of the center of mass tensor and cross vector terms in addition to the exchange tensor one. We also introduce a new parametrization of the tensor Landau parameters which significantly reduces the importance of high angular harmonic contributions. For nuclear matter and neutron matter we find that the two most important effects of the tensor interaction are to give a contribution from multipair states and to renormalize the magnetic moments. Response to a weak probe may be calculated using similar methods, replacing the magnetic moments with the matrix elements of the weak charges.

The effect of the hyperon softening of the equation of state (EOS) of dense matter on the spin evolution of isolated neutron stars is studied for a broad set of hyperonic EOSs. We use a multidomain 2-D code based on a spectral method, and show how important the precision of solving the equations of stationary motion is for the stability analysis. For some EOSs, the hyperon softening leads to spin-up by the angular momentum loss, in the form of the back bending phenomenon, for a rather broad range of stellar baryon mass. We show that large segments of the evolutionary tracks exhibiting the back bending behaviour in the moment-of-inertia - rotation-frequency plane, are actually unstable and therefore not astrophysically relevant. We show also that during the spin-up - angular-momentum-loss epoch, an isolated neutron star (e.g., a radio pulsar) can lose a sizable part of its initial angular momentum without significantly changing its rotation period. We propose also simple arguments and criteria allowing to connect the presence of a back bending epoch with the mass-radius relations and the stiffness and/or softness of the nucleon and hyperon EOSs of the neutron star core.

We study the viscosity driven (Jacobi-like) bar mode instability of rapidly rotating strange stars in general relativity. A triaxial, "bar shaped" compact star could be an efficient source of continuous wave gravitational radiation in the frequency range of the forthcoming interferometric detectors. We locate the secular instability point along several constant baryon mass sequences of uniformly rotating strange stars described by the MIT bag model. Contrary to neutron stars, strange stars with T/|W| (the ratio of the rotational kinetic energy to the absolute value of the gravitational potential energy) much lower than the corresponding value for the mass-shed limit can be secularly unstable to bar mode formation if shear viscosity is high enough to damp out any deviation from uniform rotation. The instability develops for a broad range of gravitational masses and rotational frequencies of strange quark stars. It imposes strong constraints on the lower limit of the frequency at the innermost stable circular orbit around rapidly rotating strange stars. The above results are robust for all linear self-bound equations of state assuming the growth time of the instability is faster than the damping timescale. We discuss astrophysical scenarios where triaxial instabilities (r-mode and viscosity driven instability) could be relevant in strange stars described by the standard MIT bag model of normal quark matter. Taking into account actual values of viscosities in strange quark matter and neglecting the magnetic field we show that Jacobi-like instability cannot develop in any astrophysicaly interesting temperature windows. The main result is that strange quark stars described by the MIT bag model can be accelerated to very high frequency in Low Mass X-ray binaries if the strange quark mass is ~ 200 MeV or higher.

A generic toy model of a cooling neutron star (NS) is used to analyze cooling of NSs with nucleon and exotic composition of the cores. The model contains the parameters which specify the levels of slow and enhanced neutrino emission as well as the lower and upper densities of the layer where the slow emission transforms into the enhanced one. The prospects to constrain these parameters from the present and future observations of isolated middle-aged NSs are discussed.

We study thermal states of transiently accreting neutron stars (with mean accretion rates $\dot{M} \sim 10^{-14}-10^{-9}$ M$_\odot$ yr$^{-1}$) determined by the deep crustal heating of accreted matter sinking into stellar interiors. We formalize a direct correspondence of this problem to the problem of cooling neutron stars. Using a simple toy model we analyze the most important factors which affect the thermal states of accreting stars: a strong superfluidity in the cores of low-mass stars and a fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark phases of dense matter) in the cores of high-mass stars. We briefly compare the results with the observations of soft X-ray transients in quiescence. If the upper limit on the quiescent thermal luminosity of the neutron star in SAX J1808.4-3658 (Campana et al. 2002) is associated with the deep crustal heating, it favors the model of nucleon neutron-star cores with switched-on direct Urca process.

Nuclear reactions in accreting neutron-star crusts and the heat release accompanying them are studied, under different assumptions concerning the composition of the outermost layer formed of the ashes of X-ray bursts. Particular examples of ashes containing nuclides with A ~ 90-110 are considered and compared with a standard A=56 case. In all cases, evolution of a crust shell is followed from 10^8 g/cm^3 to a few times 10^{13} g/cm^3. The total crustal heating produced in the non-equilibrium processes in the accreting crust is 1.1-1.5 MeV per one accreted nucleon. The composition of the accreted crust at densities exceeding the threshold for the pycnonuclear fusion is essentially independent of the assumed initial composition of the X-ray burst ashes.

A model of an accelerated expansion of the Crab Nebula powered by the spinning-down Crab pulsar is proposed, in which time dependence of the acceleration is connected with evolution of pulsar luminosity. Using recent observational data, we derive estimates of the Crab neutron-star moment of inertia. Correlations between the neutron star moment of inertia and its mass and radius allow for rough estimates of the Crab neutron-star radius and mass. In contrast to the previously used constant-acceleration approximation, even for the expanding nebula mass ~7 M_sun results obtained within our model do not stay in conflict with the modern stiff equations of state of dense matter.

Theoretical models of the equation of state (EOS) of neutron-star matter (starting with the crust and ending at the densest region of the stellar core) are reviewed. Apart from a broad set of baryonic EOSs, strange quark matter, and even more exotic (abnormal and Q-matter) EOSs are considered. Results of calculations of M_max for non-rotating neutron stars and exotic compact stars are reviewed, with particular emphasis on the dependence on the dense-matter EOS. Rapid rotation increases M_max, and this effect is studied for both neutron stars and exotic stars. Theoretical results are then confronted with measurements of masses of neutron stars in binaries, and the consequences of such a confrontation and their possible impact on the theory of dense matter are discussed.

Recent estimates of the properties of the Crab nebula are used to derive constraints on the moment of inertia, mass and radius of the pulsar. To this purpose, we employ an approximate formula combining these three parameters. Our "empirical formula" I = a(x)MR^2, where x = (M/M_sun)(km/R), is based on numerical results obtained for thirty theoretical equations of state of dense matter. The functions a(x) for neutron stars and strange stars are qualitatively different. For neutron stars a_NS(x) = x/(0.1+2x) for x <= 0.1 (valid for M > 0.2 M_sun) and a_NS(x) = (2/9)(1+5x) for x > 0.1. For strange stars a_SS(x) = (2/5)(1+x) (not valid for strange stars with crust and M < 0.1 M_sun). We obtain also an approximate expression for the maximum moment of inertia I_max,45 = (-0.37 + 7.12x_max) (M_max/M_sun)(R_(M_max)/(10 km))^2, valid for both neutron stars and strange stars. Applying our formulae to the evaluated values of I_Crab, we derive constraints on the mass and radius of the pulsar. A very conservative evaluation of the expanding nebula mass, M_neb = 2 M_sun, yields M_Crab > 1.2 M_sun and R_Crab = 10-14 km. Setting the most recent evaluation ("central value") M_neb = 4.6 M_sun rules out most of the existing equations of state, leaving only the stiffest ones: M_Crab > 1.9 M_sun, R_Crab = 14-15 km.

The adiabatic index Gamma_1 for perturbations of dense matter is studied under various physical conditions which can prevail in neutron star cores. The dependence of Gamma_1 on the composition of matter (in particular, on the presence of hyperons), on the stellar pulsation amplitude, and on the baryon superfluidity is analyzed. Timescales of damping of stellar pulsations are estimated at different compositions, temperatures, and pulsation amplitudes. Damping of pulsations by bulk viscosity in the neutron-star cores can prevent the stars to pulsate with relative amplitudes > (1-15)% (depending on the composition of matter).

The upper limit of the effective surface temperature of the neutron star (NS) PSR J0205+6449 in the supernova remnant 3C 58 obtained recently by Slane et al. (2002) is analyzed using a modern theory of NS cooling (Kaminker et al. 2002). The observational limit can be explained by cooling of a superfluid NS with the core composed of neutrons, protons and electrons, where direct Urca process is forbidden. However, combined with the data on the surface temperatures of other isolated NSs, it gives evidence (emphasized by Slane et al.) that direct Urca process is open in the inner cores of massive NSs. This evidence turns out to be less stringent than the evidence provided by the well known observations of Vela and Geminga.

Equilibrium configurations of cold neutron stars near the minimum mass are studied, using the recent equation of state SLy, which describes in a unified, physically consistent manner, both the solid crust and the liquid core of neutron stars. Results are compared with those obtained using an older FPS equation of state of cold catalyzed matter. The value of M_min\simeq 0.09M_sun depends very weakly on the equation of state of cold catalyzed matter: it is 0.094 M_sun for the SLy model, and 0.088 M_sun for the FPS one. Central density at M_min is significantly lower than the normal nuclear density: for the SLy equation of state we get central density 1.7 10^{14} g/cm^3, to be compared with 2.3 10^{14} g/cm^3 obtained for the FPS one. Even at M_min, neutron stars have a small liquid core of radius of about 4 km, containing some 2-3% of the stellar mass. Neutron stars with 0.09 M_sun <M<0.17 M_sun are bound with respect to dispersed configuration of the hydrogen gas, but are unbound with respect to dispersed Fe^56. The effect of uniform rotation on the minimum-mass configuration of cold neutron stars is studied. Rotation increases the value of M_min; at rotation period of 10 ms the minimum mass of neutron stars increases to 0.13 M_sun, and corresponds to the mass-shedding (Keplerian) configuration. In the case of the shortest observed rotation period of radio pulsars 1.56 ms, minimum mass of uniformly rotating cold neutron stars corresponds to the mass-shedding limit, and is found at 0.61 M_sun for the SLy EOS and 0.54 M_sun for the FPS EOS.

Bulk viscosity of neutron star cores containing hyperons is studied taking into account non-equilibrium weak process $n+n \rightleftharpoons p+\Sigma^-$. Rapid growth of the bulk viscosity within the neutron star core associated with switching on new reactions (modified Urca process, direct Urca process, hyperon reactions) is analyzed. The suppression of the bulk viscosity by superfluidity of baryons is considered and found out to be very important.

Apparent (radiation) radius of neutron star,R_\infty, depends on the star gravitational mass in quite a different way than the standard coordinate radius in the Schwarzschild metric, R. We show that, for a broad set of equations of state of dense matter, R_\infty(M_max) for the configurations with maximum allowable masses is very close to the absolute lower bound on R_\infty at fixed M, resulting from the very definition of R_\infty. Also, the value of R_infty at given M, corresponding to the maximum compactness (minimum R) of neutron star consistent with general relativity and condition v_sound<c, is only 0.6% higher than this absolute lower bound. Theoretical predictions for R_\infty are compared with existing observational estimates of the apparent radii of neutron stars.

An equation of state (EOS) of neutron star matter, describing both the neutron star crust and the liquid core, is calculated. It is based on the effective nuclear interaction SLy of the Skyrme type, which is particularly suitable for the application to the calculation of the properties of very neutron rich matter (Chabanat et al. 1997, 1998). The structure of the crust, and its EOS, is calculated in the T=0 approximation, and under the assumption of the ground state composition. The crust-core transition is a very weakly first-order phase transition, with relative density jump of about one percent. The EOS of the liquid core is calculated assuming (minimal) n-p-e-mu composition. Parameters of static neutron stars are calculated and compared with existing observational data on neutron stars. The minimum and maximum masses of static neutron stars are 0.094 M_sun and 2.05 M_sun, respectively. Effects of rotation on the minimum and the maximum mass of neutron stars are briefly discussed.

Recycling strange stars to millisecond periods is studied within the framework of general relativity. We employ equations of state of strange quark matter based on the MIT Bag Model, with massive strange quarks and lowest order QCD interactions. The presence of the crust of normal matter is taken into account, with a bottom density assumed to be equal to the neutron-drip one. The calculations are performed by solving the exact 2-D equations for rigidly rotating stationary configurations in general relativity. Evolutionary tracks of accreting strange stars are computed, and their dependence on the initial strange star mass and on the fraction of the angular momentum transferred to the star from the infalling matter, is studied. The differences between recycling strange stars and neutron stars are pointed out.

Cooling simulations of neutron stars (NSs) are performed assuming that stellar cores consist of neutrons, protons and electrons and using realistic density profiles of superfluid critical temperatures $T_{cn}(\rho)$ and $T_{cp}(\rho)$ of neutrons and protons. Taking a suitable profile of $T_{cp}(\rho)$ with maximum $\sim 5 \times 10^9$ K one can obtain smooth transition from slow to rapid cooling with increasing stellar mass. Adopting the same profile one can explain the majority of observations of thermal emission from isolated middle--aged NSs by cooling of NSs with different masses either with no neutron superfluidity in the cores or with a weak superfluidity, $T_{cn} < 10^8$ K. The required masses range from $\sim 1.2 M_\odot$ for (young and hot) RX J0822-43 and (old and warm) PSR 1055-52 and RX J1856-3754 to $\approx 1.45 M_\odot$ for the (colder) Geminga and Vela pulsars. Observations constrain the $T_{cn}(\rho)$ and $T_{cp}(\rho)$ profiles with respect to the threshold density of direct Urca process and maximum central density of NSs.

The diffusive thermal conductivity of neutrons in dense matter [$\rho \sim (1 - 8) \times 10^{14}$ g cm$^{-3}$] of neutron star cores is calculated. The contribution from neutron--neutron and neutron--proton collisions is taken into account. We use the transition probabilities calculated for symmetric dense nucleon matter on the basis of the Dirac--Brueckner approach to the in-medium effects and the Bonn model of bare nucleon--nucleon interaction. The diffusive thermal conductivity of neutrons in the presence of neutron and proton superfluidities is analyzed in a microscopic manner; the effects of superfluidity are shown to be significant. The low temperature behavior of the thermal conductivity appears to be extremely sensitive to the relation between critical temperatures of neutrons and protons. The results are fitted by simple analytic expressions. In combination with the formulae for the electron and muon thermal conductivities, obtained earlier, the present expressions provide a realistic description of the full diffusive thermal conductivity in the neutron star cores for normal and various superfluid phases.

Calculations of the properties of rotating strange stars with crusts are performed within the framework of general relativity. We employ an equation of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions. The crust is described by the BPS equation of state. A significant increase of the stellar radius is found close to the Keplerian (mass-shedding limit) configuration. This leads to the disappearance of the gap between the stellar surface and the innermost stable circular orbit (ISCO) at very high rotation rates, for a rather broad range of stellar masses. The Keplerian configuration for the strange star with crust corresponds to values of J, T/W which are about 20% smaller than in the case of bare strange stars. Because the Keplerian configuration is achieved due to the increase of the stellar oblateness, the Keplerian frequency (of the rotation) remains almost unaltered. The lack of the gap close to the Keplerian rotation could imply a more stringent limit on the frequency of the ISCO, if the existence of such a gap is supported by observations, as in the source 4U 1820-30 with the upper QPO frequency 1.07 kHz. During the spin-down of a magnetized strange quark star with crust, the crust matter is absorbed in the equatorial region by the strange matter core. The deconfinement of absorbed crust matter is a strongly exothermic process, which would influence the cooling curve of this compact object.

We study the bulk viscosity in neutron star cores due to modified Urca processes involving nucleons, electrons and muons and analize its reduction by singlet-state or triplet-state superfluidity of nucleons. In combination with the results of our previous paper on the bulk viscosity due to direct Urca processes, a realistic description of the bulk viscosity in superfluid neutron star cores is obtained. Switching off direct Urca processes with decreasing density in a nonsuperfluid matter lowers the bulk viscosity by 3-5 orders of magnitude. The presense of muons opens additional source of bulk viscosity due to muon Urca processes and lowers the threshold density of the electron direct Urca process. The superfluidity may strongly reduce the bulk viscosity and affect thus damping of neutron star vibrations.

We find constraints on minimum and maximum mass of ordinary neutron stars imposed by their early evolution (protoneutron star stage). We calculate models of protoneutron stars using a realistic standard equation of state of hot, dense matter valid for both supranuclear and subnuclear densities. Results for different values of the nuclear incompressibility are presented.

Constraints on minimum and maximum mass of ordinary neutron stars are imposed by the consideration of their early evolution (protoneutron star stage). Calculations are performed for a realistic standard model of hot, dense matter (Lattimer & Swesty 1991) valid for both supranuclear and subnuclear densities. Various assumptions concerning protoneutron star interior (large trapped lepton number, no trapped lepton number, isentropic, isothermal) are taken into account.

We review the main neutrino emission mechanisms in neutron star crusts and cores. Among them are the well-known reactions such as the electron-positron annihilation, plasmon decay, neutrino bremsstrahlung of electrons colliding with atomic nuclei in the crust, as well as the Urca processes and neutrino bremsstrahlung in nucleon-nucleon collisions in the core. We emphasize recent theoretical achievements, for instance, band structure effects in neutrino emission due to scattering of electrons in Coulomb crystals of atomic nuclei. We consider the standard composition of matter (neutrons, protons, electrons, muons, hyperons) in the core, and also the case of exotic constituents such as the pion or kaon condensates and quark matter. We discuss the reduction of the neutrino emissivities by nucleon superfluidity, as well as the specific neutrino emission produced by Cooper pairing of the superfluid particles. We also analyze the effects of strong magnetic fields on some reactions, such as the direct Urca process and the neutrino synchrotron emission of electrons. The results are presented in the form convenient for practical use. We illustrate the effects of various neutrino reactions on the cooling of neutron stars. In particular, the neutrino emission in the crust is critical in setting the initial thermal relaxation between the core and the crust. Finally, we discuss the prospects of exploring the properties of supernuclear matter by confronting cooling simulations with observations of the thermal radiation from isolated neutron stars.

We study effects of the strange quark mass and of the QCD coupling constant in the MIT bag model of quark matter with lowest order quark-gluon interactions on the rapid rotation of strange stars. The influence of rotation on global parameters of strange stars is greater than in the case of the neutron star models. We show that independently of the strange quark mass and the coupling constans the ratio of the rotational kinetic energy to the absolute value of the gravitational potential energy for rotating strange stars is significantly higher than for ordinary neutron stars. This might indicate that rapidly rotating strange stars could be important sources of gravitational waves.

The boundary between the solid crust, assumed to be in complete thermodynamic equilibrium (cold catalyzed matter), and the liquid core of a neutron star is studied using Skyrme SLy effective N-N interactions. An approximate value of the density at the inner edge of the crust is obtained from the threshold for the instability of homogeneous $npe$ matter with respect to small periodic density perturbations. Calculations of the bottom layer of the crust are performed within the Compressible Liquid Drop Model. Spherical nuclei are energetically preferred over exotic ones (cylinders, plates, tubes, bubbles), down to the inner edge of the crust, found at the density 0.08 fm^{-3}.

We compute the maximum mass of moderately rotating strange stars as a function of the strange quark mass, of the QCD coupling constant, $\alpha_c$, and of the bag constant (vacuum energy density), $B$, in the MIT bag model of quark matter with lowest order quark-gluon interactions. For a fixed value of $B$, the maximum stellar mass depends only weakly on $\alpha_c$, and is independent of this coupling in the limit of massless quarks. However, if it is the value of the chemical potential of quark matter at zero pressure which is held constant, for example at the value corresponding to the stability limit of nucleons against conversion to quark matter, the maximum mass of the strange star is higher by up to 25% for $\alpha_c=0.6$, than for non-interacting quarks, and this may be relevant in the discussion of kHz QPO sources. The maximum mass of a non-rotating strange star could be sufficiently high to allow an orbital frequency as low as 1.0 kHz in the marginally stable orbit. However, for all $\alpha_c<0.6$, the stellar mass cannot exceed $2.6M_\odot$ at any rotational period $\ge1.6 $ms.

The bulk viscosity of the neutron star matter due to the direct Urca processes involving nucleons, electrons and muons is studied taking into account possible superfluidity of nucleons in the neutron star cores. The cases of singlet-state pairing or triplet-state pairing (without and with nodes of the superfluid gap at the Fermi surface) of nucleons are considered. It is shown that the superfluidity may strongly reduce the bulk viscosity. The practical expressions for the superfluid reduction factors are obtained. For illustration, the bulk viscosity is calculated for two models of dense matter composed of neutrons, protons,electrons and muons. The presence of muons affects the bulk viscosity due to the direct Urca reactions involving electrons and produces additional comparable contribution due to the direct Urca reactions involving muons. The results can be useful for studying damping of vibrations of neutron stars with superfluid cores.

We analyze cooling of neutron stars, assuming the presence of localized protons in the densest region of their cores. Choosing a single threshold density for proton localization and adjusting neutron star mass, we reproduce the observational data on effective surface temperatures of Vela and PSR 0656+14, with or without an accreted hydrogen envelope. However, the presence of a tiny hydrogen envelope is mandatory, in this model, for reproducing the Geminga data.

Exact calculations of innermost stable circular orbit (ISCO) around rotating strange stars are performed within the framework of general relativity. Equations of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions, are used. The presence of a solid crust of normal matter on rotating, mass accreting strange stars in LMXBs is taken into account. It is found that, contrary to neutron stars, above some minimum mass a gap always separates the ISCO and stellar surface, independently of the strange star rotation rate. For a given baryon mass of strange star, we calculate the ISCO frequency as function of stellar rotation frequency, from static to Keplerian configuration. For masses close to the maximum mass of static configurations the ISCO frequencies for static and Keplerian configurations are similar. However, for masses significantly lower than the maximum mass of static configurations, the minimum value of the ISCO frequency is reached in the Keplerian limit. Presence of a solid crust increases the ISCO frequency for the Keplerian configuration by about ten percent compared to that for a bare strange star of the same mass.

We study neutrino energy emission rates (emissivities) due to electron bremsstrahlung produced by $ee$ and $ep$ collisions in the superfluid neutron star cores. The neutrino emission due to $ee$ collisions is shown to be the dominant neutrino reaction at not too high temperatures ($T \la 10^8$ K) in dense matter if all other neutrino reactions involving nucleons are strongly suppressed by neutron and proton superfluidity. Simple practical expressions for the $ee$ and $ep$ neutrino emissivities are obtained. The efficiency of various neutrino reactions in the superfluid neutron-star cores is discussed for the cases of standard neutrino energy losses and the losses enhanced by the direct Urca process.

Exact models of uniformly rotating strange stars, built of self bound quark matter, are calculated within the framework of general relativity. This is made possible thanks to a new numerical technique capable to handle the strong density discontinuity at the surface of these stars. Numerical calculations are done for a simple MIT bag model equation of state of strange quark matter. Evolutionary sequences of models of rotating strange stars at constant baryon mass are calculated. Maximally rotating configurations of strange stars are determined, assuming that the rotation frequency is limited by the mass shedding and the secular instability with respect to axisymmetric perturbations. Exact formulae which give the dependence of the maximum rotation frequency, and of the maximum mass and corresponding radius of rotating configurations, on the value of the bag constant, are obtained. The values of T/W for rapidly rotating massive strange stars are significantly higher than those for ordinary neutron stars. This might indicate particular susceptibility of rapidly rotating strange stars to triaxial instabilities.

As pointed out by Kutschera and W{\'o}jcik, very low concentration of protons combined with a specific density dependence of effective neutron-proton interaction could lead to a localization of ``proton impurities'' in neutron medium at densities exceeding four times normal nuclear matter density. We study consequences of the localization of protons for transport processes in dense neutron star cores, assuming random distribution of proton impurities. Kinetic equations, relevant for the transport of charge, heat and momentum, are solved using variational method. Localization of protons removes a T^{-2} factor from the transport coefficients, which leads, at lower temperatures, to a strong decrease of thermal conductivity, electrical conductivity and shear viscosity of neutron star matter, as compared to the standard case, where protons form a Fermi liquid. Due to the localization of protons a number of conventional neutrino emission processes (including modified URCA process) become inoperative in neutron star cores. On the other hand, the energy loss rate from neutrino-antineutrino pair bremsstrahlung due to electron and neutron scattering off (localized) protons, will have a specific T^6 dependence, which could modify the cooling of the neutron star core, as compared to the standard case. Possible astrophysical implications of the localization of protons for neutron star evolution and dynamics are discussed.

New calculations of the thermal and electrical electron conductivities are performed for a broad range of physical parameters typical for envelopes of neutron stars and cores of white dwarfs. We consider stellar matter composed of astrophysically important chemical elements from H to Fe in the density range from 10^2-10^4 g/ccm up to 10^7-10^{10} g/ccm, where atoms are fully ionized and electrons are strongly degenerate. We have used modified ion structure factors suggested in physics/9811052. In the ion liquid, these modifications take into account, in an approximate way, instantaneous electron-band structures that reduce the electron-ion scattering rate. In crystallized matter, the new structure factors include multi-phonon processes important at temperatures not very much lower than the melting temperature T_m. The transport coefficients obtained differ significantly from those derived earlier in the important temperature range T_m/5 < T < 5 T_m. The results of our numerical calculations are fitted by analytical expressions convenient for astrophysical applications.

The lower bound on the period of uniform rotation of neutron stars (NS) with causal equation of state (EOS), P_ min, is shown to be determined, with a high precision, by the upper bound on the surface redshift of static neutron star models, z_ max. We obtain z_max=0.851, which turns out to be about 5% lower than upper bound obtained by Lindblom (1984). Using then the empirical formula for the maximum rotation frequency of uniform rotation of neutron star models with realistic causal EOS, derived by Lasota et al. (1996), we establish an approximate relation between z_max and P_min, in which P_min is proportional to the maximum mass of static NS, M_s^ max. In this way we reproduce, within 2%, the exact formula of Koranda et al. (1997) based on extensive exact numerical calculations. Replacing M_s^max by the maximum measured NS mass, 1.442 M_sun, we get a lower bound P_min = 0.288 ms.

We calculate the neutrino energy emission rate due to singlet-state pairing of protons in the neutron star cores taking into account the relativistic correction to the non-relativistic rate. The non-relativistic rate is numerically small, and the relativistic correction appears to be about 10 -- 50 times larger. It plays thus the leading role, reducing great difference between the neutrino emissions due to pairing of protons and neutrons. The results are important for simulations of neutron star cooling.

We show that the neutron star minimum period of rotation is determined, for causal equations of state, by the maximum value of the relativistic (compactness) parameter 2GM/Rc^2, allowed by causality for static neutron stars, 0.7081, and by the largest measured mass of a neutron star. The relation between these three quantities, resulting from the extrapolation of the empirical formulae of Lasota et al. (1996), yields the minimum period of 0.288 ms, only 2% higher than an absolute lower bound obtained in extensive exact numerical calculations of Koranda et al. (1997).

We address the question concerning the maximum density, at which nuclei (and more generally -- nuclear structures) can exist in neutron star interiors. An absolute upper bound to such a density is obtained using the bulk approximation, in which surface and Coulomb effects are neglected. A very good approximation to the limiting density is given by the threshold for the instability of a uniform n-p-e plasma with respect to density modulations; this threshold is calculated using the Extended Thomas-Fermi approximation for the Skyrme energy functionals. For recent SLy Skyrme forces, which are particularly suitable for the description of very neutron rich nucleon systems, the maximum density at which nuclear structures can exist is 0.08 nucleons/fm^3 (i.e., half of the normal nuclear density); at this density protons constitute only 4\% of nucleons.

Models of differentially rotating protoneutron stars are calculated, using realistic equations of state of dense hot matter. Various conditions within the stellar interior, corresponding to different stages of protoneutron star evolution, are considered. Families of differentially rotating models of a given baryon mass are calculated, using a two-parameter formula describing the angular velocity profile within a rotating protoneutron star. Maximum angular momentum, which can be accommodated by a protoneutron star is calculated, for various thermal conditions in stellar interior, for a baryon mass of $1.5 M_\odot$. In the case of a thermally homogeneous (isentropic or isothermal) neutrino-opaque interior this maximum angular momentum turns out to be somewhat higher than that of a cold neutron star of the same baryon mass, rotating uniformly at the mass shedding angular velocity. However, if the protoneutron star has a thermal structure characteristic of initial state, with a low entropy (unshocked) core, and a high entropy (shocked) outer half of baryon mass, the maximum angular momentum is significantly lower. This leads to a minimum period of uniform rotation of cold neutron stars of baryon mass $\sim 1.5~M_\odot$, formed directly (i.e. without a subsequent significant accretion of mass) from protoneutron stars with shocked envelope, of about 1.7 ms and strengthens the hypothesis that millisecond pulsars are accretion accelerated neutron stars.

We study the synchrotron emission of neutrino pairs by relativistic, degenerate electrons in strong magnetic fields. Particular attention is paid to the case in which the dominant contribution comes from one or several lowest cyclotron harmonics. Calculations are performed using the exact quantum formalism and the quasiclassical approach. Simple analytic fits to the neutrino synchrotron emissivity are obtained in the domain of magnetized, degenerate and relativistic electron gas provided the electrons populate either many Landau levels or the ground level alone. The significance of the neutrino synchrotron energy losses in the interiors of cooling neutron stars is discussed.

Radial pulsations of newborn neutron stars (protoneutron stars) are studied for a range of internal temperatures and entropies per baryon predicted by the existing numerical simulations. Protoneutron star models are constructed using a realistic equation of state of hot dense matter, and under various assumptions concerning stellar interior (large trapped lepton number, zero trapped lepton number, isentropic, isothermal). Under prevailing conditions, linear oscillations of a protoneutron star can be treated as adiabatic, and evolutionary effects can be neglected on dynamic timescale. The dynamic behavior is governed by the adiabatic index, which in turn depends on the physical state of the stellar interior. The eigenfrequencies of the lowest radial modes of linear, adiabatic pulsations are calculated. Stability of protoneutron stars with respect to small radial perturbations is studied, and the validity of the static stability criteria is discussed.

We study the emission of neutrinos, resulting from the scattering of electrons off magnetic flux tubes (fluxoids) in the neutron star cores with superfluid (superconducting) protons. In the absence of proton superfluidity (T> T_{cp}), this process transforms into the well known electron synchrotron emission of neutrino pairs in a locally uniform magnetic field B, with the neutrino energy loss rate Q proportional to B^2 T^5. For temperatures T not much below T_{cp}, the synchrotron regime (Q \propto T^5) persists and the emissivity Q can be amplified by several orders of magnitude due to the appearance of the fluxoids and associated enhancement of the field within them. For lower T, the synchrotron regime transforms into the bremsstrahlung regime (Q \propto T^6) similar to the ordinary neutrino-pair bremsstrahlung of electrons which scatter off atomic nuclei. We calculate Q numerically and represent our results through a suitable analytic fit. In addition, we estimate the emissivities of two other neutrino-production mechanisms which are usually neglected -- neutrino-pair bremsstrahlung processes due to electron-proton and electron-electron collisions. We show that the electron-fluxoid and electron-electron scattering can provide the main neutrino production mechanisms in the neutron star cores with highly superfluid protons and neutrons at T < 5 10^8 K. The electron-fluxoid scattering is significant if the initial, locally uniform magnetic field B > 10^{13} G.

Rapid uniform rotation of newborn neutron stars (protoneutron stars) is studied for a range of internal temperatures and entropies per baryon predicted by the existing numerical simulations. Calculations are performed using general relativistic equations of hydrostatic equilibrium of rotating, axially symmetric stars. Stability of rotating configurations with respect to mass shedding and the axially symmetric perturbations is studied. Numerical calculations are performed for a realistic dense matter equation of state, under various assumptions concerning neutron star interior (large trapped lepton number, no trapped lepton number, isentropic, isothermal). For configurations with baryon mass well below the maximum one for the non-rotating models, the mass shedding limit depends quite sensitively on the position of the ``neutrinosphere''. The absolute upper limit on rotation frequency is, to a good approximation, obtained for the maximum baryon mass of rotating configurations. Empirical formula for the maximum rotation frequency of uniformly rotating protoneutron stars is shown to be quite precise; it actually coincides with that used for cold neutron stars. Evolutionary sequences at fixed baryon mass and angular momentum, which correspond to evolution of protoneutron stars into cold neutron stars are studied, and resulting constraints on the maximum rotation frequency of solitary pulsars are discussed.

We derive exact expressions for the cross section of neutrino scattering on electrons in dense, hot stellar matter, in the presence of strong magnetic fields. Numerical calculations of the scattering cross sections at various densities, temperatures and magnetic fields, are performed. Strong, quantizing magnetic fields modify significantly the angular and energy dependence of the scattering cross section.

Contents: 1.Introduction 2. Formation and composition of the solid crust 3. The bottom layers of the crust: exotic nuclear shapes and topologies 4. Crust contribution to stellar mass and moment of inertia 5. Solid cores 6. Deformation, elastic strain and mountains 7. Shape and energy of rotating, partially solid neutron star 8. Oblateness and precession 9. Gravitational radiation from precessing pulsars 10. Gravitational radiation from precessing pulsars and pulsar timing 11. Detectability

Possible consequences of ferromagnetic transition in dense matter suggested recently by Kutschera and W{\'o}jcik, for the magnetic properties of neutron stars, are studied. Specific model of dense matter, in which a small admixture of protons is completely polarized due to their interaction with neutrons, is considered. Magnetic field of neutron stars with a ferromagnetic core is calculated within the framework of general relativity. Two types of boundary conditions at the ferromagnetic core edge are considered, corresponding to normal and superconducting liquid envelope, respectively. Numerical results for the neutron star magnetic dipole moment are confronted with pulsar timing. To be consistent with observations, ferromagnetic cores surrounded by a non-superconducting envelope, should consist of weakly ordered ferromagnetic domains. If domains are highly ordered, ferromagnetic core should be screened by a superconducting envelope.

Neutrino emissivity from the electron $\nu\bar\nu$ bremsstrahlung in the liquid layers of the neutron star crusts is studied. Nuclear composition of matter in neutron star crusts is considered for various scenarios of neutron star evolution. For the deep layers of the crust, the compositions of cold catalyzed matter, accreted matter and hot matter ($T \ga 5 \times 10^9$~K) are shown to be very different, and this implies differences in the neutrino emissivity at given density and temperature. Neutrino-pair bremsstrahlung, due to collisions of relativistic degenerate electrons in a Coulomb liquid of atomic nuclei, is considered. The neutrino energy loss rate is expressed in a simple form through a Coulomb logarithm $L$ -- a slowly varying function of density, temperature, and nuclear composition. Non-Born corrections, thermal width of electron Fermi sphere, and finite sizes of nuclei are taken into account. Implications for cooling of neutron stars are discussed.

We show that for realistic equations of state of dense matter, the universal proportionality factor relating the maximum rotation rate of neutron stars due to mass-shedding limit to the mass and radius of maximum allowable mass configuration of non-rotating models results from a universal proportionality between masses and radii of static maximum-mass neutron stars and those of maximally rotating configurations. These empirical relations cannot be obtained in the slow rotation approximation.

We study the consequences of recent progress in the experimental determination of masses of neutron rich nuclei for our knowledge of the ground state of cold dense matter. The most recent experimental data determine the ground state of cold dense matter up to $ \rho \simeq 10^{11}~{\rm g~cm^{-3}} $. The composition and the equation of state of the ground state of matter, in this density interval, are calculated.