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.

Using parametric equations of state (relativistic polytropes and a simple quark bag model) to model dense-matter phase transitions, we study global, measurable astrophysical parameters of compact stars such as their allowed radii and tidal deformabilities. We also investigate the influence of stiffness of matter before the onset of the phase transitions on the parameters of the possible exotic dense phase. The aim of our study is to compare the parameter space of the dense matter equation of state permitting phase transitions to a sub-space compatible with current observational constraints such as the maximum observable mass, tidal deformabilities of neutron star mergers, radii of configurations before the onset of the phase transition, and to give predictions for future observations. We studied solutions of the Tolman-Oppenheimer-Volkoff equations for a flexible set of parametric equations of state. We compare tidal deformabilities of stars with weak and strong phase transitions with the results of the GW170817 neutron star merger. Specifically, we study characteristic phase transition features in the $\Lambda_1-\Lambda_2$ relation, and estimate the deviations of our results from the approximate formul{\ae} for $\tilde{\Lambda}-R(M_1)$ and $\Lambda$-compactness proposed in the literature. We find constraints on the hybrid equations of state to produce stable neutron stars on the twin branch. For the exemplary equations of state most of the high-mass twins occur for the minimum values of the density jump $\lambda = 1.33-1.54$; corresponding values of the square of the speed of sound are $\alpha = 0.7-0.37$. We compare results with observations of gravitational waves and with the theoretical causal limit and find that the minimum radius of a twin branch is 9.5 - 10.5 km. For these solutions, the phase transition occurs below 0.56 [fm$^{-3}$].

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.

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

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.

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

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

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

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

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.

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.

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.

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.

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

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.

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

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.

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.

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.

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.

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.

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.

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.

The minimum value of the radius of strange star covered by the crust of nuclear matter is determined. The results for the maximum possible thickness of the crust (up to the neutron drip) as well as the possibility of thinner crust postulated by some authors are discussed. The minimum radius of the strange star with maximal crust is 5.5 km. The useful scaling formulae with respect to the main parameters describing strange matter and the density at the bottom of the crust are presented.

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.

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.

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.

It is shown that for very rapidly rotating low mass strange stars the marginally stable orbit is located above the stellar surface. This effect is explained by the very important role of the oblateness of the rotating strange star. The comparison with some ``academic'' examples is presented. This feature is purely Newtonian in its nature and has nothing to do with relativistic marginally stable orbit. The effect is very large and cannot be treated in a perturbative way. It seems that strange stars as a very dense self-bound objects are the only possibility in Nature to represent these toy models.

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.

It has been suggested that observations of quasi periodic oscillations (QPOs) in the X-ray fluxes from low mass X-ray binaries can be used to constrain the mass of the compact object and the equation of state of its matter. A specific suggestion that the kHz QPO frequency saturates at the maximum orbital frequency has been widely considered. We examine rotating compact stars described by a new model of strange quark matter (Dey et al. 1998). We calculate the maximum orbital frequencies for both normal and supramassive strange stars described by the Dey model, and present these frequencies for sequences of equilibrium models with constant baryon mass. The maximum orbital frequencies for these compact objects are always higher than the kHz QPO frequencies observed to date.

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.

We compute numerical models of uniformly rotating strange stars (SS) in general relativity for the recently proposed QCD-based equation of state (EOS) of strange quark matter (Dey et al. 1998). Static models based on this EOS are characterised by a larger surface redshift than strange stars within the MIT bag model. The frequencies of the fastest rotating configurations described by Dey model are much higher than these for neutron stars (NS) and for the simplest SS MIT bag model. We determine a number of physical parameters for such stars and compare them with those obtained for NS. We construct constant baryon mass equilibrium sequences both normal and supramassive. Similarly to the NS a supramassive SS, prior to collapse to a black hole, spins up as it loses angular momentum. We find the upper limits on maximal masses and maximal frequencies of the rotating configurations. We show that the maximal rotating frequency for each of considered evolutionary sequences is never the Keplerian one. A normal and low mass supramassive strange stars gaining angular momentum always slows down just before reaching the Keplerian limit. For a high mass supramassive SS sequence the Keplerian configuration is the one with the lowest rotational frequency in the sequence. The value of $T/W$ for rapidly rotating SS of any mass is significantly higher than those for ordinary NS. For Keplerian configurations it increases as mass decreases. The results are robust for all linear self-bound equations of state.

The equation of state of strange matter build of the u,d,s quarks in the framework of the MIT bag model is studied. The scaling relations with bag constant are discussed and applied to the determination of maximum frequency of a particle in stable circular orbit around strange star. The highest QPO frequency of 1.33 kHz observed so far is consistent with the strange stars models for which the maximum QPO frequency is 1.7-2.4 kHz depending on the strange quark mass and the QCD coupling constant. The linear approximation of the equation of state is found and the parameters of this EOS are determined as a functions of strange quark mass, QCD coupling constant and bag constant. This approximation reproduces exact results within an error of the order of 1% and can be used for the complete study of the properties of strange stars including microscopic stability of strange matter and determination of the total baryon number of the star.

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.

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.

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.

Radial pulsations of neutron stars and strange quark stars with nuclear crust are studied. The avoided crossing phenomenon occurring for the radial modes is found and discussed. Neutron star models are constructed using a realistic equation of state of dense matter and strange star models using a phenomenological bag model of quark matter. The eigenfrequencies of the three lowest modes of linear, adiabatic pulsations are calculated, using the relativistic equations for the radial oscillations.

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

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.

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.

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.

The role of direct URCA reactions in damping of the gravitational radiation driven instability is discussed. The temperature at which bulk viscosity suppresses completely this instability is calculated. The results are obtained analytically using recent calculations performed in the case of bulk viscosity due to the modified URCA processes (Lindblom 1995; Yoshida & Eriguchi 1995). The bulk viscosity caused by direct URCA reactions is found to reduce significantly the region of temperatures and rotation frequencies where a neutron star is subject to GRR instability.