neutrons stars are born in gravitational collapse of massive , degenerate stellar cores .
newly born neutron stars are hot and lepton rich objects , quite different from ordinary low temperature , lepton poor neutron stars . in view of these differences ,
newly born neutron stars are called _ protoneutron _ stars ; they transform into standard neutron stars on a timescale of the order of ten seconds , needed for the loss of a significant lepton number excess via emission of neutrinos trapped in the dense , hot interior .    in view of the fact that the typical evolution timescale of a protoneutron star ( seconds ) is some three orders of magnitude longer , than the dynamical timescale for this objects ( milliseconds )
, one can study its evolution in the quasistatic approximation ( burrows & lattimer 1986 ) . properties of static ( non - rotating ) protoneutron stars , under various assumptions concerning composition and equation of state ( eos ) of hot , dense stellar interior were studied by numerous authors ( burrows & lattimer 1986 , takatsuka 1995 , bombaci et al .
1995 , bombaci 1996 , bombaci et al .
1996 ) .
the scenario of transformation of a protoneutron star into a neutron star could be strongly influenced by a phase transition in the central region of the star .
brown and bethe ( 1994 ) suggested a phase transition implied by the @xmath0 condensation at supranuclear densities .
such a @xmath0 condensation could dramatically soften the equation of state of dense matter , leading to a low maximum allowable mass of neutron stars .
in such a case , the massive protoneutron stars could be stabilized by the effects of high temperature and of the presence of trapped neutrinos , and this would lead to maximum baryon mass of protoneutron star larger by some @xmath1 than that of cold neutron stars .
the deleptonization and cooling of protoneutron stars of baryon mass exceeding the maximum allowable baryon mass for neutron stars , would then inevitably lead to their collapse into black holes .
the dynamics of such a process was recently studied by baumgarte et al .
it should be mentioned , however , that the very possibility of existence of the kaon condensate ( or other exotic phases of matter , such as the pion condensate , or the quark matter ) at neutron star densities is far from being established .
recently , for instance , pandharipande et al . ( 1995 ) pointed out that kaon - nucleon and nucleon - nucleon correlations in dense matter raise significantly the threshold density for kaon condensation , possibly to the densities higher than those characteristic of stable neutron stars . in view of these uncertainties
, we will restrict in the present paper to a standard model of dense matter , composed of nucleons and leptons .
the calculations of the static models of protoneutron stars should be considered as a first step in the studies of these objects .
it is clear , in view of the dynamical scenario of their formation , that protoneutron stars are far from being static .
due to the nonzero initial angular momentum of the collapsing core , protoneutron stars are expected to rotate . on the other hand ,
the formation scenario involves compression ( with overshoot of central density ) and a hydrodynamical bounce , so that a newborn protoneutron star begins its life in a highly excited state , pulsating around its quasistatic equilibrium . in the present paper
we study the rotation of protoneutron stars ; pulsations of protoneutron stars will be discussed in a separate paper ( gondek , haensel & zdunik , in preparation ) .
some aspects of rapid uniform rotation of protoneutron stars have been recently studied in ( takatsuka 1995 , hashimoto et al .
however , the calculations reported by takatsuka ( 1995 ) were actually done for static ( non - rotating ) protoneutron stars , and were then used to estimate the maximum rotation frequency of uniformly rotating protoneutron stars , @xmath2 , via an `` empirical formula '' .
it should be stressed , that the validity of such an `` empirical formula '' , which expresses @xmath3 in terms of the mass and radius of the extremal _
static _ configuration with maximum allowable mass , had been checked only in the restricted case of _ cold _ neutron stars ( haensel & zdunik 1989 , friedman et al .
1989 , shapiro et al .
1989 , haensel et al .
1995 , nozawa et al .
. only isentropic equations of state were considered by takatsuka ( 1995 ) .
hashimoto et al . (
1995 ) calculated the structure of stationary configurations of uniformly rotating protoneutron stars , using a two - dimensional general relativistic code .
these authors restricted themselves to the case with zero trapped lepton number .
they assumed a constant temperature in the hot interior of the star , and used a zero temperature ( cold ) eos for @xmath4 .
it should be stressed , that the assumption of @xmath5 corresponds to an isothermal state in the newtonian ( flat space - time ) theory of gravitation . in general relativity
, we will define isothermal state by @xmath6 ( where @xmath7 is the lapse function and @xmath8 is the lorentz factor , see section 3.1 ) , and the effects of the space - time curvature will turn out to be rather important for massive neutron stars . also , their choice for the low density edge of the hot interior can be questioned .
finally , their criterion for finding maximally rotating configuration is actually valid only for cold ( @xmath9 ) or isentropic protoneutron stars : its use in the case of the @xmath5 hot interior is unjustified ( see section 3.2 for a correct statement of the stability criterion ) . in a recent paper , lai and shapiro ( 1995 ) have studied the secular evolution , secular `` bar instability '' , and the gravitational wave emission from the newly formed , rapidly rotating neutron stars .
however , these authors used unrealistic ( polytropic ) equations of state of neutron star matter .
moreover , the calculations were done within newtonian theory of gravitation . in view of this , the internal structure of their models of newly born neutron stars was quite different from that characteristic of the realistic models of protoneutron stars .
the problem of the secular `` bar instability '' in rapidly rotating neutron stars was also studied , using general relativity , by bonazzola et al ( 1995 ) . however , numerical calculations were done only for realistic equations of state of _ cold _ neutron star matter .    in the present paper
we study the properties of uniformly rotating protoneutron stars , using exact relativistic description of the rapid , stationary rotation , combined with realistic equations of state of hot dense matter , used in the whole range of temperatures and densities relevant for protoneutron stars .
in particular , we calculate the maximum frequency of uniform rotation of protoneutron stars and its dependence on their baryon mass , and on the thermal state and composition of stellar interior .
it is clear , that uniform rotation represents only an approximation to the actual rotational state of a newly born protoneutron star .
existing numerical simulations of gravitational collapse of rotating cores of massive stars produce differentially rotating protoneutron stars ( janka & moenchmeyer 1989a , b , moenchmeyer & mueller 1989 ) .
however , it should be stressed that the initial rotational state of collapsing core is unknown , and this implies uncertainty concerning the rotational state of resulting protoneutron star .
it is reasonable to say , that the actual degree of nonuniformity of rotation of a protoneutron star should be considered as unknown . in the present paper
we will not address the question of the physical mechanisms that could `` rigidify '' the rotational motion within the protoneutron star interior .
however , we will use the approximation of uniform rotation in order to limit the number the parameter space for our numerical calculation , and also because of the relative simplicity of the stability analysis in this specific , idealized case .    within our simplified model ,
the `` neutrinosphere '' ( which has actually a deformed , spheroidal shape ) will separate hot , neutrino - opaque interior of a protoneutron star ( hereafter referred to as `` hot interior '' ) from a significantly cooler , neutrino - transparent envelope .
the actual thermal state of the hot interior of protoneutron star is determined by its formation scenario , and is expected to be influenced by the dissipative processes ( damping of pulsations , viscous damping of differential rotation , neutrino diffusion ) . for simplicity , we will restrict ourselves to two limiting cases : an isothermal ( @xmath6 , see section 3 ) , and an isentropic ( entropy per baryon @xmath10 ) hot interior .
we will also consider two limiting cases of the lepton composition of the protoneutron star interior .
the first case , referring to the very initial state of protoneutron star , will correspond to a fixed trapped lepton number . in the second case
, neutrinos will not contribute to the lepton number of the matter , which will correspond to vanishing chemical potential of the electron neutrinos ; such a situation will take place after a deleptonization of a protononeutron star .
the position of the neutrinosphere will be located using a simple prescription based on specific properties of the neutrino opacity of hot dense matter . in all cases , the equation of state of hot dense matter
will be determined using one of the models of lattimer and swesty ( 1991 ) .
the plan of the paper is as follows . in section 2
we describe the physical state of the interior of protoneutron star , with particular emphasis on the eos of the hot interior at various stages of evolution .
we explain also our prescription for locating the `` neutrinosphere '' of a protoneutron star , and we give some details concerning the assumed temperature profile within a protoneutron star . using simple estimates of the timescales relevant for various transport processes , we justify the approximation of stationarity which is used throughout this paper . in section 3
we give a brief description of the exact equations , used for the calculation of stationary configuration of uniformly rotating protoneutron stars .
we discuss also stability of rotating configurations with respect to the axially - symmetric perturbations .
the numerical method , used for the calculation of rapidly rotating configurations of protoneutron star , is briefly described in section 4 , where we also discuss numerical precision of our solutions .
maximum rotation frequency , for various physical conditions prevailing in the hot stellar interior , calculated as a function of the baryon ( rest ) mass of protoneutron star , is presented in section 5 .
then , in section 6 we show the validity of an empirical formula , which enables one to express with a surprisingly high precision the maximum frequency of rotating protoneutron stars in terms of the mass and radius of the maximum mass configuration of static ( non - rotating ) protoneutron stars with same eos . in section 7
we study the evolutionary transformation of a rotating protoneutron star into a cold neutron star .
we show that , at fixed rest mass and angular momentum , maximum rotation frequency of protoneutron stars imposes severe constraints on the rotation frequency of solitary neutron stars .
finally , section 8 contains discussion of our results and conclusion .
we consider a protoneutron star ( pns ) just after its formation .
we assume it has a well defined `` neutrinosphere '' , which separates a hot , neutrino - opaque interior from colder , neutrino - transparent outer envelope .
important parameters , which determine the local state of the matter in the hot interior are : baryon ( nucleon ) number density @xmath11 , net electron fraction @xmath12 , and the net electron - neutrino fraction @xmath13 . the calculation of the composition of hot matter and of its eos is described below .
such a situation is characteristic of the very initial stage of existence of a pns .
matter is composed of nucleons ( both free and bound in nuclei ) and leptons ( electrons and neutrinos ; for simplicity , we do not include muons ) .
all constituents of the matter ( plus photons ) are in thermodynamic equilibrium at given values of @xmath11 , @xmath14 and @xmath15 .
the composition of the matter is calculated from the condition of beta equilibrium , combined with the condition of a fixed @xmath16 , @xmath17 where @xmath18 are the chemical potentials of matter constituents . at the very initial stage
we expect @xmath19 .
electron neutrinos are degenerate , with @xmath20 ( in what follows we measure @xmath14 in energy units ) .
the deleptonization , implying the decrease of @xmath16 , occurs due to diffusion of neutrinos outward ( driven by the @xmath21 gradient ) , on a timescale of seconds ( sawyer & soni 1979 , bombaci et al .
the diffusion of highly degenerate neutrinos from the central core is a dissipative process , resulting in a significant _ heating _ of the neutrino - opaque core ( burrows & lattimer 1986 ) .
this is the limiting case , reached after complete deleptonization .
there is no trapped lepton number , so that @xmath23 and @xmath24 , and therefore @xmath25 .
neutrinos trapped within the hot interior do not influence the beta equilibrium of nucleons , electrons and positrons , and for given @xmath11 and @xmath14 the equilibrium value of @xmath26 is determined from @xmath27 while @xmath28 . in practice , this approximation can be used as soon as electron neutrinos become non - degenerate within the opaque core , @xmath29 , which occurs after some @xmath30 seconds ( sawyer & soni 1979 , bombaci et al .
the neutrino diffusion is then driven by the temperature gradient , and the corresponding timescale of the heat transport ( pns cooling ) can be estimated as @xmath31  s ( sawyer & soni 1979 ) .      in principle ,
the temperature ( or entropy per nucleon ) profile within a pns has to be determined via evolutionary calculation , starting from some initial state , and taking into account relevant transport processes in the pns interior , as well as neutrino emission from pns .
transport processes within neutrino - opaque interior occur on timescales of seconds , some three orders of magnitude longer than dynamical timescales .
the very outer layer of pns becomes rapidly transparent to neutrinos , deleptonizes , and cools on a very short timescale via @xmath32 pair annihilation and plasmon decay .
it seemed thus natural to model the thermal structure of the pns interior by a hot core limited by a `` neutrinosphere '' , and an outer envelope of   @xmath33 mev .
the transition through the `` neutrinosphere '' is accompanied by a temperature drop , which takes place over some interval of density just above the `` edge''of the hot neutrino - opaque core , situated at some @xmath34 .    in view of the uncertainties in the actual temperature profiles within the hot interior of pns , we considered two extremal situations for @xmath35 , corresponding to an isentropic and an isothermal hot interior . in the first case ,
hot interior was characterized by a constant entropy per baryon @xmath10 . in the case of trapped lepton number
, this leads to the eos of the type : pressure @xmath36)$ ] , energy density @xmath37)$ ] , and temperature @xmath38)$ ] , with fixed @xmath39 and @xmath16 .
this eos will be denoted by eos[@xmath40 . in the case of an isentropic ,
zero trapped lepton number eos , we will have eos[@xmath41 .
the condition of isothermality , which corresponds to thermal equilibrium , is more complicated . due to the curvature of the space - time within pns
, the condition of isothermality corresponds to the constancy of @xmath42 ( see section 3.1 ) .
the significance of the @xmath43 condition will be discussed in section 3.1 . in the static case
, the isothermal state within the hot interior will be reached on a timescale corresponding to thermal equilibration , which is much longer than the lifetime of a pns . in the case of a rotating pns
, the situation can be expected to be even more complicated ( see , e.g. , chapter 8 of tassoul ( 1978 ) ) . nevertheless , we considered the @xmath43 models for several reasons .
first , as a limiting case so different from the @xmath10 one , it enables us to check the dependence of our results for rapidly rotating pns on the thermal state of the hot interior . moreover , for the isothermal pns we can apply the criterion of stability with respect to the axi - symmetric perturbations , and this will enable us to calculate the value of @xmath3 for the stable supramassive rotating pns models with an isothermal interior ( see section 3 ) .
our calculation of the `` neutrinosphere '' within the hot pns interior was done using a simple method , described below . for a given
pns model , the neutrinosphere radius , @xmath44 , has been located through the condition @xmath45 where @xmath46 is the neutrino mean free path which is calculated at the matter temperature @xmath14 , while @xmath47 is the mean energy of non - degenerate neutrinos at ( and above ) the neutrinosphere , @xmath48 .
we assumed that opacity above @xmath44 is dominated by the elastic scattering off nuclei and nucleons , so that @xmath49 . then , we determined the value of the density at the neutrinosphere , @xmath34 , for a given static pns model , combining eq .
( [ r_nu ] ) with that of hydrostatic equilibrium , and readjusting accordingly the temperature profile within the outer layers of pns .
this value of @xmath34 was then used in the calculations of rotating pns models . let us notice
, that similar approximation , in which the `` neutrinosphere '' in rapidly rotating protoneutron star was defined as a surface of constant density , was used by janka and moenchmeyer ( 1989a , b ) . some details concerning actual calculation of the temperature profile in the vicinity of the `` neutrinosphere '' will be given in subsection 3.3 .      the starting point for the construction of our eos for the pns models
was the model of hot dense matter of lattimer and swesty ( 1991 ) , hereafter referred to as ls .
actually , we used one specific ls model , corresponding to the incompressibility modulus at the saturation density of symmetric nuclear matter @xmath50mev .
for @xmath35 we supplemented the ls model with contributions resulting from the presence of trapped neutrinos of three flavours ( electronic , muonic and tauonic ) and of the corresponding antineutrinos .    in fig .
1a , b we show our eos in several cases , corresponding to various physical conditions in the hot , neutrino - opaque interior of pns . for the sake of comparison ,
we have shown also the eos for cold catalyzed matter , used for the calculation of the ( cold ) ns models . in fig .
1a we show eos at subnuclear densities . at these densities , both the temperature and the presence of trapped neutrinos stiffen the eos , as compared to the cold catalyzed matter one , and the stiffening is rather dramatic .
the constant @xmath14 eos stiffens considerably at lower densities , which is due to the weak dependence of the thermal contribution ( photons , neutrinos ) on the baryon density of the matter ( this effect will be to some extent moderated by the factor @xmath51 in the isothermal pns , see section 3.1 ) .
it is quite obvious , that @xmath5 eos becomes dominated by thermal effects below for @xmath52 . on the contrary , for isentropic eos
, the effect of the trapped lepton number ( @xmath53 ) turns out to be much more important than the thermal effects .
this can be see in fig .
1a , by comparing dash - dotted curve , @xmath54 $ ] , with the dotted line , which corresponds to an artificial ( unphysical ) case with small thermal effects , @xmath55 $ ] .
it is clear , that the correct location of the `` neutrinosphere '' , which separates hot interior from the colder outer envelope , should be important for the determination of the radius of pns , and in consequence , of the maximum rotation frequency of a given pns model
.    it may be useful to compare our subnuclear eos for pns with those used by other authors .
the subnuclear eos of pns , used in the papers of hashimoto et al .
( 1995 ) and bombaci et al .
( 1995 , 96 ) is very different from that used in the present paper .
in particular , hashimoto et al . (
1995 ) used a cold ( @xmath9 ) eos for the densities below @xmath56 . on the other hand , bombaci et al .
( 1995 , 96 ) stop their hot eos at the edge of the liquid interior , and use the @xmath9 ( cold catalyzed matter ) eos for the densities below @xmath57 ; in this way , they seriously underestimate thermal and neutrino trapping effects on the radius of pns .    our eos above nuclear density are plotted in fig .
the presence of a trapped lepton number softens the eos , while thermal effects always stiffen it with respect to that for cold catalyzed matter .
the softening of the supranuclear eos at fixed @xmath16 is due to the fact , that a significant trapped lepton number increases the proton fraction , which implies the softening of the nucleon contribution to the eos .
it should be stressed , that in contrast to hashimoto et al .
( 1995 ) and bombaci et al .
( 1995 , 96 ) we used a unified dense matter model , valid for both supranuclear and subnuclear densities .
also , the fact that we use various assumptions about the @xmath14 and @xmath39 profiles within pns , enables us to study the relative importance of the temperature profile and that of a trapped lepton number , for the pns models .
the mass - radius relation for the static pns models calculated using various versions of our eos for the hot interior is shown in fig .
we assumed @xmath58 , which was consistent with our definition of the `` neutrinosphere '' .
for the sake of comparison , we show also the mass - radius relation for the @xmath9 ( cold catalyzed matter ) eos , which corresponds to cold neutron star models . in the case of the isothermal hot interior with central temperature @xmath59  mev we note a very small increase of the maximum mass , as compared to the @xmath9 case ( c.f .
, bombaci et al . 1995 , 1996 ) .
however , the effect on the mass - radius relation is quite strong , and increases rapidly with decreasing stellar mass . in the case of the isentropic eos with a trapped lepton number , [ @xmath60
, the softening of the high - density eos due to the trapped @xmath16 leads to the decrease of @xmath61 compared to the @xmath9 case ; as far as the value of @xmath61 is concerned , the softening effect of @xmath16 prevails over that of finite @xmath39 ( this is consistent with results of takatsuka 1995 and bombaci et al .
1995 , 1996 ) .
however , the thermal effect on the radius is very important even in the case of @xmath53 .
this can be seen by comparing the @xmath54 $ ] curve with that corresponding to the unphysical , fictitious case of @xmath55 $ ] .
the very initial state of a pns corresponds to a significant trapped lepton number . with our assumption of a `` standard '' composition of dense matter ( i , e .
, excluding large amplitude @xmath0-condensate , or a very large percentage of hyperons in cold dense matter ) , the maximum baryon mass ( baryon number ) of pns is lower than that of cold ns .
therefore , in our case a stable pns transforms into a stable ns , and the scenario pns@xmath62black hole , considered by baumgarte et al .
( 1996 ) is excluded .    at a given mass ,
the radius of a pns is significantly larger than that of a cold ns . as remarked by hashimoto et al .
( 1995 ) , this should have important implications for rotating pns
. it should be stressed , however , that the value of radius , especially for pns which are not close to the @xmath61 configuration , turns out to be quite sensitive to the location of the edge of the hot neutrino - opaque interior ( i.e. , to the value of @xmath34 ) . the choice of bombaci et al .
( 1996 ) would lead to a much smaller effect on @xmath63 , while that of hashimoto et al .
( 1995 ) would result in larger values of the pns radii .
the eos of pns is evolving with time , due mainly to the deleptonization process , which changes the composition of matter , and also due to changes of the internal temperature of the star .
however , these changes occur on the timescales @xmath641 - 10  s , which are three or more orders of magnitude longer than the dynamical timescale , governing the readjustment of pressure and gravity forces .
this dynamical timescale @xmath65ms corresponds also to the characteristic periods of the pns pulsations and of their rapid rotation . in view of this , we are able to decouple pns evolution from its dynamics , and treat its rotation in the stationary approximation , with a well defined eos of the pns matter
.    one of the neglected dynamical processes , implied by the radiative processes and the evolution of the thermal structure of a rotating pns , is the meridional circulation of the matter .
strictly uniform rotation is incompatible with an assumption of a steady thermal state , resulting from the diffusive ( radiative ) equilibrium ( see , e.g. , tassoul 1978 ) .
the requirement of radiative equilibrium will necessarily imply the existence of a meridional circulation of the matter .
however , the velocity of this meridional circulation will be of the order of the stellar radius divided by the thermal timescale ( the timescale of changes of the entropy of the pns interior , which is of the order of neutrino diffusion timescale ) , which is much smaller than the rotational velocity ( see , eg . ,
section 8 of tassoul 1978 ) . in view of this
, we can neglect the effect of the meridional circulation when calculating the mechanical equilibrium of a rapidly rotating pns .
finally , let us notice that in a special , idealized case of @xmath66 , the radiative flux vanishes .
then , pure uniform rotation can be realized as a steady state of a pns .
the problem of the calculation of the stationary state of uniform rotation of cold neutron stars , within the framework of general relativity , was considered by numerous authors ( see the review article of friedman and ipser 1992 , and references therein ) .
extensive calculations for a broad set of realistic equations of state of cold dense matter were recently presented in ( cook et al . 1994 ) and ( salgado et al . 1994 ) . here , we will extend the methods used at @xmath9 to the case of hot pns .
we will use the notation and formalism developed in the paper of bonazzola et al .
( 1993 ) , hereafter referred to as bgsm .
one of the problems introduced by finite temperature is that , if one does not make any other assumption about the equilibrium , the equation of stationary motion does not have a first integral ( bardeen 1972 ) .    in the notation of bgsm , the equation of stationary motion ( eq
. 3.25 of bgsm ) reads , when the effects of temperature and rigid rotation are included , as @xmath67 where @xmath39 is the entropy per baryon , @xmath14 is the temperature , and @xmath68 $ ] is the so - called pseudo - enthalpy ( or log - enthalpy ) , @xmath7 is the lapse function appearing in the space - time metric , and @xmath8 is the lorentz factor due to rotation .
the equation of the stationary motion is to be supplemented with the equations determining the metric functions ( see bsgm for the derivation and the explicit form of the complete set of equations ) .    in the expression for @xmath69 , @xmath70 is the energy density ( which includes rest energy of matter constituents ) and @xmath71 is the nucleon rest mass .
this equation is the general relativistic equivalent of a well - known newtonian formula ( see e.g. tassoul 1978 ) .
it is straightforward to show that a _ sufficient _ condition for ( [ eqsta ] ) to be integrable is @xmath72 . in this case , one obtains a first integral of motion of the form : @xmath73 this enables us to calculate the density profile in the envelope of the pns , because we have assumed a specific @xmath74 profile within it ( see section 3.3 ) .
let us stress that , due to the integral term in ( [ firstint ] ) , the pseudo - enthalphy is no more an explicit function of the metric potentials , as it was in bgsm .
in fact , one must solve ( [ firstint ] ) to have @xmath69 within the star .    in the two particular cases , chosen by us for the hot interior , specific first integrals of eq .
( [ eqsta ] ) can be found .
first , in the case of general relativistic thermal equilibrium , @xmath75 , it is easy to show that @xmath76 where @xmath77 is the baryon chemical potential , is indeed the integral of ( [ eqsta ] ) .
the constancy of @xmath78 and @xmath79 corresponds to the general - relativistic thermodynamic equilibrium , if we neglect the time dependence of the eos of pns ( i.e. , if we freeze transport phenomena ) .
second , in the case of isentropic profile @xmath80 , the thermal term in ( [ eqsta ] ) vanishes and the first integral of eq .
( [ eqsta ] ) reads @xmath81      the instabilities that could develop in rapidly rotating pns , and which could limit the maximum angular frequency of these objects , are of secular and of dynamic type .
let us start with the problem of secular stability with respect to the axi - symmetric perturbation of rotating configurations .
we denote the total baryon number , total angular momentum , and total entropy of a pns by @xmath82 , @xmath83 , and @xmath84 , respectively .
we used the _ secular _ instability criterion of friedman et al .
( 1988 ) , extended to the case of finite temperature . for the purpose of completeness , we restate here their lemma ( and correct a misprint of their paper ) . +
* lemma :* consider a three - parameter family of uniformly rotating hot stellar models having an equation of state of the form @xmath85 .
suppose that along a continous sequence of models labelled by a parameter @xmath86 , there is a point @xmath87 at which @xmath88 , @xmath89 and @xmath90 vanish and where @xmath91 .
then the part of the sequence for which @xmath92 is unstable for @xmath86 near @xmath87 .
+ similarly as in the case of friedman et al . (
1988 ) , this lemma follows directly from theorem i of sorkin ( 1982 ) , with his function @xmath84 replaced by @xmath93 , the @xmath94 quantities replaced by @xmath82 , @xmath83 , @xmath84 and the @xmath95 ones by @xmath96 , @xmath97 , @xmath98 .
the conditions of the theorem are fulfilled because stellar models are configurations for which @xmath99 is minimized at fixed @xmath82 , @xmath83 and @xmath84 , _ and _ because difference in @xmath99 between two neighbouring equilibria can be expressed as ( bardeen 1970 ) @xmath100 let us stress here that this last equation requires either @xmath98 or @xmath39 to be constant through the whole star , which fixes , in each of these two cases , the temperature profile within the stellar model .
we investigated the stability of our models with a specific version of this criterion , in which we choose a continous sequence of equilibria to be a sequence at fixed @xmath82 and @xmath84 .
the point of the loss of stability is then simply the point of extremal @xmath83 , i.e. : @xmath101 where @xmath102 is the central density of the star .
the instability discussed above is a secular one ; it will develop on the timescale needed to transport the angular momentum within the perturbed model , in order to decrease the energy of the star while changing its shape and structure .
however , the timescale of the pns evolution is quite short ( seconds ) , and it is driven by the same transport processes involving neutrinos , as those which are needed to destabilize the star via the axi - symmetric perturbations . in view of this , one might expect that the secular instability described above is not efficient in disrupting the quasi - stationary configuration of rapidly rotating pns .
however , we should remember that the above considerations apply to the case of infinitesimal perturbations .
transport processes would be crucial for removing the energy barrier separating the initial , secularly unstable configuration , from the dynamically unstable one , which would eventually collapse into a black hole . however , newly born pns are expected to be in a highly excited state , in which various modes of stellar pulsations are excited .
one may expect that the energy contained in these pulsations is sufficient to overcome the energy barrier separating the actual metastable , secularly unstable state from the dynamically unstable , collapsing one . in view of this
, we expect that the secularly unstable configurations should be treated like unstable ones .
therefore , the critical configuration , given by eq.(9 ) , will be thus considered as the last stable one .
a rapidly rotating pns can be also susceptible to other types of secular instabilities .
the instability with respect to the non axisymmetric perturbations can be driven by the gravitational radiation reaction ( grr ) . however , detailed calculations performed for hot ns suggest , that at the temperature exceeding @xmath103k these instabilities can only slightly decrease the maximum rotation frequency of uniform rotation , due to the damping effect of the matter viscosity ( cutler et al .
1990 , ipser & lindblom 1991 , lindblom 1995 , yoshida & eriguchi 1995 , zdunik 1996 ) .
secular instabilities of rapidly rotating ns with respect to the non - axisymmetric `` bar '' mode were recently investigated by lai & shapiro ( 1995 ) and bonazzola et al .
while lai & shapiro ( 1995 ) addressed the problem of `` bar '' instability of newly formed ns , they assumed a rather unrealistic `` ellipsoidal model '' , involving only shear viscosity as a source of viscous dissipation , and performed their calculations within the newtonian theory of gravity . on the other hand , in their relativistic calculations bonazzola et al .
( 1995 ) considered only cold ns , and found that the secular `` bar '' instability can set in before the keplerian ( mass shedding ) limit is reached only for sufficiently stiff eos of ns matter .
the problem of the `` bar '' instability of pns , with realistic eos of the hot interior , will be investigated by us in the future . in any case ,
inclusion of possible additional secular instabilities can only decrease the maximum rotation frequency of pns below the values obtained using the simplified approach adapted in the present paper .
the fact of the existence of the maximum mass of rotating pns , @xmath104 ( and maximum baryon mass , @xmath105 ) ( see section 5 ) , puts a well defined and stringent limit on @xmath106 which can be reached by pns .
configurations with @xmath107 are _ dynamically unstable _ : no stationary solution exists above @xmath108 . at the same time , the angular frequency of rigid rotation can not exceed the keplerian value , @xmath106 .
these two conditions , combined with the instability criterion , expressed in eq .
( 9 ) , determine the maximum frequency ( or , strictly speaking , an upper bound on the frequency ) of rigid rotation of pns .
practical implementation of these criteria will be described in section 5 .
a consistent study of pns would require an exact treatment of the thermal transport in the frame of a non - stationary spacetime .
unfortunately , such a study is beyond the scope of the present work , and we have thus chosen to impose `` by hand '' the temperature profile within the star , following prescription described in the subsection 2.3 .
we divided the interior of pns into the hot interior ( @xmath109 ) , a layer corresponding to the temperature drop within the `` neutrinosphere '' ( @xmath110 ) , and the low temperature , neutrino - transparent outer envelope with @xmath111 .
we have chosen two types of temperature profiles within the hot , neutrino - opaque core @xmath112 :    * the isothermal profile : @xmath113 , with @xmath114 , * the isentropic profile : @xmath115 , with @xmath116 .
for the transition region and the low temperature envelope , we used a suitable profile @xmath74 : @xmath117 the function @xmath118 was chosen for computational convenience as a suitable combination of an exponential and a gaussian function , selected to lead to @xmath74 of class @xmath119 through the transition ( temperature drop ) region .
our typical choice was @xmath120 ; increasing this value up @xmath121 led to a very small increase of the stellar radius .
as the mass of a massive pns is almost entirely contained in the hot neutrino - opaque core ( a rough estimate for a typical star gives less than @xmath122 of the total mass for the cool envelope mass ) , we supposed that our stability criterion remained valid also in the case of the presence of the low temperature envelope .
we used a code based on the @xmath123 formulation of the einstein equations in stationary axisymmetric spacetimes ( bgsm ) .
the four elliptic equations obtained were solved by means of a spectral method , in which the functions are expanded in different polynomial bases ( chebyshev for @xmath124 , legendre for @xmath125 and fourier for @xmath126 ) .
we refer the reader to bgsm for a detailed description of the code , including the description of the `` virial '' indicator used for monitoring the convergence and the precision reached .
let us just say that the code was modified to take into account thermal effects and , contrary to salgado et al .
( 1994 ) , to converge to the solution with various quantities being held fixed ( for example @xmath127 or @xmath83 or @xmath84 ) .
here we briefly outline some of the numerical checks we made .
the two - parameters eos was interpolated using bicubic splines subroutine from the nag library . in this way , the thermodynamic functions are of class @xmath128 , but the thermodynamic consistency is not conserved .
the global relative error , evaluated by the means of the virial check ( see bgsm ) , is @xmath129 .
this relatively `` low '' precision is due to the thermodynamical inconsistencies .
we used two grids in @xmath124 for the star and one in @xmath130 for the exterior with @xmath131 in each , and one grid in @xmath125 with @xmath132 .
we checked that a greater @xmath133 or @xmath134 does not change the results by more than a few @xmath135 at most , which stays within the global precision reached .
let us stress that such a low number of grid points is sufficient to reach high accuracy within a spectral method ( that would not be the case with a finite difference scheme ) .
the temperature drop at the `` neutrinosphere '' is not always located on the border of the internal grid in @xmath124 , which could influence the precision ( remember that the spectral methods are very sensitive to discontinuities ) .
we checked that , in fact , even when the `` neutrinosphere '' is far from the border of the grid , the global physical quantities of the star did not change much ( we found also relative variations of a few @xmath135 at most ) .
for a `` simple '' model where @xmath136 were held fixed , convergence required no more than @xmath137 iterations and @xmath138 minute of cpu time on a silicon graphix indigo@xmath139 workstation . however , the number of iterations can reach @xmath140 if one needs to converge to , e.g. , a keplerian configuration at fixed supra - massive baryon mass .
finally , we compared our results at @xmath9 and @xmath141 with the results of another code used by one of us ( p. haensel ) , and found that the relative differences in the global properties of the stellar models are of the order of @xmath142 , which can be imputed to the different interpolation procedures .
for a given eos , rapidly rotating pns models can be divided into two families : normal and supramassive one .
normal models are those with baryon mass which does not exceed the maximum allowable baryon mass of static ( non - rotating ) configurations , @xmath143 . a _ normal _ rotating model can be transformed into a static configuration of the same baryon mass , through a continuous decrease of @xmath97 .
rotation of pns increases their maximum allowable mass with respect to the static case , up to @xmath144 .
rotating pns models with @xmath145 are called _
supramassive_. such a supramassive rotating model can not transform into a static model , because in the process of decreasing angular velocity it will collapse into a black hole .    the maximum angular velocity of the normal and supramassive rotating models results from different stability conditions . in the case of normal rotating configurations
the value of @xmath97 is bound by the mass shedding limit , which corresponds to the keplerian velocity at the stellar equator , @xmath146 .
the value of @xmath106 turns out to be rather sensitive to the location of the `` neutrinosphere '' within the pns .
this sensitivity is particularly large in the case of the isothermal profile of the hot neutrino - opaque core .
this is visualized in fig .
3 , where we show the dependence of the mass shedding limit @xmath106 for pns with baryon mass @xmath147 on the value of @xmath34 , for the central temperatures ranging between 10 and 20 mev . the dependence on the value of @xmath34 weakens with increasing @xmath127 . also , this effect is much less important in the case of isentropic pns .
llllll & & & & & + eos & @xmath148 & @xmath149 & @xmath104 & @xmath150 & @xmath151 + & @xmath152 $ ] & [ km ] & @xmath152 $ ] & [ km ] & [ hz ] + & & & & & + & & & & & + @xmath9 , @xmath153 & 2.048 & 10.59 & 2.430 & 14.34 & 1625 + @xmath154 , @xmath153 & 2.053 & 11.17 & 2.322 & 14.79 & 1521 + @xmath155 , @xmath53 & 1.957 & 10.85 & 2.180 & 14.40 & 1522 + @xmath156 , @xmath53 & 1.977 & 11.53 & 2.172 & 15.32 & 1388 + & & & & & + & & & & & +    for supramassive rotating pns , the limit on @xmath97 has been set by the condition of stability of rotating models with respect to the axi - symmetric perturbations , which was combined with the mass shedding stability condition .
our calculations have been performed for both isothermal ( @xmath43 ) and isentropic ( @xmath10 ) hot interiors of pns . in both cases ,
the absolute maximum of keplerian frequency for rotating models , which were stable with respect to the axi - symmetric perturbations , was obtained for a rotating configuration with a maximum baryonic mass ( and gravitational mass ) , @xmath157 [ @xmath104 ] . actually , rotating configuration with @xmath104 and that with @xmath2 do not generally coincide ( see , e.g. , cook et al .
1994 , stergioulas & friedman 1995 ) .
however , the difference is very small , and it could not be detected within the precision of our numerical code .
the value @xmath104 depends on the value of @xmath78 in the case of isothermal hot interior , but we preferred to parametrize it in terms of central temperature , @xmath159 .
the dependence of @xmath104 on the value of @xmath160 is displayed in fig .
4 a. the effect of temperature on @xmath104 is opposite to that seen for @xmath148 : thermal effects lower the value of @xmath104 as compared to that for cold neutron stars .
this is due to the thermal increase of the equatorial radius , which prevails over the thermal stiffening of the central core of pns .
the dependence of @xmath161 on the value of @xmath39 in the case of isentropic hot interior of pns is represented in fig .
4 b , in the case of @xmath162 ( trapped lepton number ) and for @xmath53 .
the decrease of @xmath161 is there smaller than in the case of the isothermal pns models .    in fig .
5 a , 5 b we show our results for the maximum rotation frequency of stable pns models , reached for supramassive configurations . in the case of the isothermal pns , thermal effects tend to decrease the value of @xmath2 , but even in the case of @xmath163mev the relative effect does not exceed ten percent .
irregularities at lower @xmath160 result from the limited accuracy of determination of the critical configuration with an isothermal core .
maximum rotation frequency in the case of the isentropic hot neutrino - trapped cores of pns is shown , for several values of the central entropy per baryon , in fig .
5 b. the value of @xmath2 is also decreasing with increasing @xmath39 .    in general , we find that the decrease of @xmath2 , as compared to cold neutron stars , is in our case smaller than that found by hashimoto et al .
( 1995 ) ; this may be due to the lower value of the density at the edge of the hot core , assumed by these authors .
some differences may also be due to different supranuclear eos , and to a different treatment of the thermal effects , in particular , to their assumption of @xmath5 . in order to visualize the importance of the relativistic effects on @xmath164 in the isothermal interior of a pns
, we show in fig .
6 the temperatures @xmath164 and @xmath78 , for @xmath165 @xmath166 , and a uniform keplerian rotation .
as one can see , in the central part of the pns , @xmath14 is about @xmath167 percent greater than @xmath168 .
a similar effect is seen also in the case of an isentropic hot interior , which in fig .
6 corresponds to @xmath169 .
the calculation of the rotating pns ( and ns ) models is incomparably more difficult than that of the static models . in the case of cold ns ,
one finds a surprisingly precise universal formula , which relates the maximum rotation frequency to the mass and radius of the maximum mass configuration for the static models ( haensel & zdunik 1989 , friedman at al .
1989 , shapiro et al .
1989 , haensel et al .
1995 , nozawa et al .
1996 ) , @xmath170 where the most recent value of @xmath171 , based on calculations performed for a broad set of realistic _
cold _ eos of dense matter , is @xmath172 ( haensel et al . 1995 ) .
the `` empirical formula '' , eq .
( 10 ) , reproduces the values of @xmath2 for cold ns models with typical precision better than @xmath173 ( although in some specific cases the deviations can reach nearly @xmath174 , see nozawa et al .
1996 ) .
the validity of the empirical formula for @xmath2 is of great practical importance . for cold eos of dense matter
, it enables one to get immediately a rather precise estimate of the value of @xmath2 , using easily calculated static neutron star models , and avoiding in this way incomparably more difficult 2-d calculations of rotating ns models and the analysis of their stability .
it is of interest to check , whether `` empirical formula '' is valid also in the case of hot pns . let us notice , that the subnuclear eos of pns , with possible lepton number trapping , is very different from that of cold catalyzed matter ( see fig
. 1 ) . in order to check the validity of the empirical formula
, we define the eos depending parameter @xmath175 , @xmath176 where the right - hand - side is calculated using exact results for a specific eos .
the values of @xmath175 for our models of pns are displayed in fig .
7 a , 7 b. let us notice , that for isentropic hot cores with trapped lepton number the value of @xmath175 shows a decreasing trend with increasing @xmath39 . in the case of isothermal hot cores with @xmath153 we see the opposite trend : the value of @xmath175 tends to increase with central temperature .
the value of @xmath177 will lead to a very good empirical formula for pns , with precison comparable to that used for cold ns . within the precision of the empirical formulae ,
the values of @xmath178 and @xmath179 can be considered as being equal .
maximum rotation frequencies and maximum baryon masses of pns with various eos , and corresponding frequencies of cold ns configurations , obtained via cooling at constant @xmath83 and @xmath127 ( @xmath83-losses due to neutrino emission are neglected ) . [ cols= " < , < ,
< , < " , ]     rapidly rotating hot , neutrino - opaque pns evolves eventually into a cold , rotating ns , which under favourable circumstances can be observed as a solitary pulsar .
let us assume , that such a transformation took place at constant baryon mass : this assumption is valid if mass accreted after the formation of a pns is negligibly small .
( notice , that we can _ define _ the moment of formation of a pns as that at which the accretion ends ) .
second assumption refers to the angular momentum of the star .
let us neglect for the time being the angular momentum loss due to emission of neutrinos ( mainly during deleptonization ) , as well as that due to gravitational radiation ( if rotating star deviated from the axial symmetry ) .
inclusion of both effects could only _ decrease _ the final angular momentum of the ns . under these two assumptions ,
the transformation of a pns into a solitary radio pulsar takes place at constant baryon mass , @xmath127 , and angular momentum , @xmath83 .
the maximum angular frequency of a solitary radio pulsar of a given baryon mass would be then determined by the parameters of a maximally rotating pns of the same baryon mass .    in the present paper
we restrict ourselves , similarly as hashimoto et al .
( 1995 ) , to the case of rigid rotation of pns .
we will first consider `` typical '' range of baryon masses expected for the ns born in a gravitational collapse of a massive stellar core , @xmath180 ( for our eos of dense matter this would correspond to @xmath181 ) .
all stars considered are normal , and the maximum initial frequency of pns will be then @xmath182 $ ] , and @xmath183 $ ] , for the limiting case of the zero trapped lepton number , and maximum trapped lepton number neutrino - opaque cores , respectively .
the values of these limiting angular frequencies , for the isothermal and isentropic pns cores , are plotted versus @xmath127 in fig . 8 ( thin dashed and dash - dotted lines , respectively ) .    our results displayed in fig .
8 lead to some conclusions , relevant to those neutron stars which remained solitary after their birth in the gravitational collapse of a massive stellar core .
if the starting configuration ( at the end of significant accretion , after the revival of a stagnated shock ) was that with trapped lepton number and isentropic , then the maximum frequency of such a neutron star , in the gravitational mass range of 1.3 - 1.8  @xmath184 , can not exceed 600 - 900  hz , the lower frequency limit corresponding to the lower mass limit .
the corresponding range of minimum periods would be from 1.6  ms to 1.1  ms , respectively .    if the `` initial isolated configuration '' ( i.e. , that with no subsequent accretion ) was a little older , and therefore deleptonized , and if we assume that it was isothermal , then the resulting constraints would be less stringent : for the range of masses 1.3 - 1.8  @xmath184 we get minimum rotation periods of 1.25 - 1.1  ms , respectively .    in both cases cooling resulted in speeding up of the rotation of the star :
such situation is characteristic of `` normal rotating models '' .
similar evolutionary considerations can be applied to the `` maximally rotating configurations '' : isentropic one with @xmath169 and isothermal ones with @xmath185mev ,  25  mev , respectively .
our results are shown in table 2 .
the initial `` hot '' configurations are supramassive , and evolve into more compact cool ones ( at constant @xmath186 and @xmath83 ) _ decreasing _ their angular frequency ( slowing down ) .
this purely relativistic effect was recently pointed out by hashimoto et al .
the `` relativistic slowing down '' during cooling of the maximally rotating isentropic configuration with @xmath169 corresponds to the increase of the period from @xmath187=0.72 $ ]  ms up to @xmath188  ms .
if such a scenario of formation of solitary pulsars is valid , an absolute limit on their period would be @xmath189 , and not @xmath190  ms , obtained for the @xmath9 eos .
the slowing down factor of @xmath191 coincides with that obtained by hashimoto et al .
( 1995 ) for a different eos and a different scenario of formation of solitary cold pulsars .
actually , they assumed that the initial configuration is that with @xmath192 .
their initial model would to some extent correspond to our isothermal models ( remember however the factor @xmath193 in our temperature profiles , which is absent in their calculations ) . using tables 1 , 2 we obtain @xmath194/p_{\rm min}[t=0]=0.92 $ ] , which is quite close to 1 .
in the present paper we studied rapid rotation of protoneutron stars , using a specific model of hot , dense matter .
our models of rapidly rotating protoneutron stars were based on several simplifications , which were necessary in order to make the problem tractable . in order to avoid difficulties and/or ambiguities of the largely unknown `` real situation '' , we restricted ourselves to studying idealized , limiting cases .
in many places we introduced approximations , which were crucial for making numerical calculations feasible .
we assumed rigid rotation within protoneutron star .
this simplified greatly our considerations , and reduced dramatically the number of stellar models .
actually , our formalism allows calculation of the quasi - stationary differentially rotating configurations , after introducing an additional @xmath195 term appearing on the right hand side of the equation of stationary motion , eq .
( 4 ) ( see section 5.1 of bgsm ) .
we plan to perform studies of differentially rotating protoneutron stars in the near future , as the next step in our investigations of dynamics of protoneutron stars .
while existing numerical simulations of gravitational collapse of rotating cores of massive stars yield a differentially rotating protoneutron star , the calculations stop too early after bounce .
in contrast to our models of protoneutron stars , the object which comes out from the simulations of moenchmeyer and mueller has an extended , very hot envelope , produced by the shock wave immediately after bounce ( janka & moenchmeyer 1989a , b , moenchmeyer & mueller 1989 ) .
also , the lack of knowledge of the initial angular velocity distribution within the collapsing core results in the uncertainty in the rotational state of produced protoneutron star . clearly , the study of quasistationary differential rotation of protoneutron stars should take into account all these uncertainties .
deleptonization of protoneutron star is connected with energy and angular momentum losses .
however , angular momentum taken away by neutrinos is expected to constitute at most a few percent of the total stellar angular momentum ( kazanas 1977 ) .
inclusion of this effect would slightly decrease some of our `` evolutionary limits '' of section 7 .
our treatment of the thermal state of the protoneutron star interior should be considered as very crude .
the temperature profile might be affected by convection . also , our method of locating the `` neutrinosphere '' was very approximate .
clearly , the treatment of thermal effects can be refined , but we do not think this will change our main results .
our calculations were performed for only one model of the nucleon component of dense hot matter .
the model was realistic , and enabled us to treat in a unified way the whole interior ( core as well as the envelope ) of the protoneutron star .
however , in view of the uncertainties in the eos of dense matter at supranuclear densities , one should of course study the whole range of theoretical possibilities , for a broad set - from soft to stiff - of supranuclear , high temperature eos .
an example of such an investigation , in the case of the _ static _ protoneutron stars , is the study of bombaci et al .
( 1996 ) . in view of the possible importance of the protoneutron star - neutron star connection for the properties of solitary pulsars , similar studies should also be done for rotating proto - neutron stars .
the calculations of the present paper were done under the assumption , that the frequency of uniform rotation of protoneutron stars is limited only by the mass shedding and the secular instability with respect to the axi - symmetric perturbations . in view of this , our results can be considered only as upper bounds to the maximum rotation frequency .
our main conclusion is that the minimum rotation period of solitary neutron stars , born as rapidly rotating protoneutron stars , is significantly larger , than the corresponding limit for cold neutron stars .
inclusion of additional secular instabilities in rapidly , uniformly rotating protoneutron stars can only strengthen this conclusion .
we are very grateful to e. gourgoulhon for his help at the initial stage of this project .
we are also very grateful to w. dziembowski for introducing us into the difficult subject of meridional circulation .
this research was partially supported by the jumelage program `` astronomie france - pologne '' of cnrs / pan , by the kbn grant no .
p304 014 07 , and by the mesr grant
94 - 3 - 1544 .
the numerical computations have been performed on the silicon graphics workstations , purchased thanks to the support of the spm department of the cnrs and the institut national des sciences de lunivers .
bardeen j.m .
, 1970 , apj , 162 , 71 bardeen j.m . , 1972 , in `` black holes - les astres occlus '' , les houches , ed .
de witt c. & de witt b.s . ,
gordon and breach sci . pub .
, new york baumgarte t.w .
shapiro , s.a .
teukolsky , 1996 , apj 458 , 680 bombaci i. , prakash m. , prakash m. , ellis p.j . , lattimer j.m . , brown g.e . , 1995 ,
a , 583 , 623 bombaci i. , 1996 , a&a 305 , 871 bombaci i. , prakash m. , prakash m. , ellis p.j . , lattimer j.m . , brown g.e . , 1996 ,
reports , in press bonazzola s. , gourgoulhon e. , salgado m. , marck j.a . , 1993 ,
a&a 278 , 421 bonazzola s. , frieben j. , gourgoulhon e. , 1995 , apj , 460 , 379 brown g.e . , bethe h.a . , 1994 ,
apj 423 , 659 burrows a. , lattimer j.m . , 1986 ,
apj , 307 , 178 cook g.b . ,
shapiro s.l . ,
teukolsky s.a . , 1994 ,
apj 424 , 823 cutler c. , lindblom l. , splinter r.j . , 1990 ,
apj 363 , 603 friedman j.l . ,
ipser j.r . ,
sorkin r.d . , 1988 ,
apj 325 , 722 friedman j.l . , ipser j.r .
, parker l. , 1989 , phys .
62 , 3015 friedman j.l .
, ipser j.r . , 1992 , phil . trans .
lond . a , 340 , 391 haensel p. , zdunik j.l . , 1989 ,
nature 340 , 617 haensel p. , salgado m. , bonazzola s. , 1995 , a&a 296 , 745 hashimoto m - a , oyamatsu k. , eriguchi y. , 1995 , apj 436 , 257 ipser j.r . ,
lindblom l. , 1991 , apj 373 , 213 janka h .- t .
, moenchmeyer r. , 1989a , a&a , 209 , l5 janka h .- t .
, moenchmeyer r. , 1989b , a&a , 226 , 69 kazanas d. , 1977 , nature 267 , 501 lai d. , shapiro s.l .
, 1995 , apj , 442 , 259 lattimer j.m . ,
swesty f.d , 1991 , nucl .
a535 , 331 lindblom l. , 1995 , apj 438 , 265 moenchmeyer r. , mueller e. , 1989 , in `` timing neutron stars '' , ed .
oegelman h. , van den heuvel e.p.j , kluwer , dordrecht nozawa t. , hashimoto m. , oyamatsu k. , eriguchi y. , 1996 , phys . rev .
d 53 , 1845 pandharipande v.r .
, pethick c.j . ,
thorsson v. , 1995 , phys.rev.lett .
75 , 4567 salgado m. , bonazzola s. , gourgoulhon e. , haensel p. , 1994
, a&a 291 , 155 sawyer r.f . , 1980 ,
apj 237 , 187 sawyer r.f . ,
soni a. , 1979 , apj 230 , 859 shapiro s.l . ,
teukolsky s.a . ,
wasserman i. , 1989 , nature 340 , 451 sorkin r.d . , 1982 , apj 257 , 847 stergioulas n. , friedman j.l . , 1995 , apj 444 , 306 takatsuka t. , 1995 , nucl .
a588 , 365 tassoul j.l . ,
1978 , `` theory of rotating stars '' , princeton university press , princeton yoshida s. , eriguchi y. , 1995 , apj 438 , 830 zdunik j.l .
, 1996 , a&a 308 , 828