neutron star cooling is an important tool for the study of dense matter . by comparing cooling models with observations of thermal emission from these objects , we can gain insight into the equation of state ( eos ) of dense matter , the signatures of exotic particles , superfluid energy gaps , and magnetic field properties [ for a recent review , see @xcite ] .
a neutron star loses the thermal energy with which it was born initially through neutrino emission , and later through photon emission , the change between these two stages occurring about @xmath0 yr after birth .
several heating mechanisms may become important during late stages of the evolution of a neutron star , most of them related to the delayed adjustment to the progressively changing equilibrium state as the rotation slows down .
these include the dissipation of rotational energy due to interactions between superfluid and normal components of the star @xcite and release of strain energy stored by the solid crust due to spin - down deformation @xcite .
another of these mechanisms is _ rotochemical heating _
@xcite , which has its origin in deviations from beta equilibrium due to spin - down compression .
as a neutron star spins down , the centrifugal force acting on each fluid element diminishes , changing the local value of the pressure . since the composition of
neutron star matter in beta equilibrium is a one - parameter function , this compression results in a displacement of the equilibrium concentration of each particle species .
reactions which change the chemical composition are in charge of driving the system to the new equilibrium configuration .
but if the rate at which reactions do this task is slower than the change of the equilibrium concentrations due to spin - down compression , the system is permanently out of chemical equilibrium .
this implies an excess of energy , which is dissipated by enhanced neutrino emission and heat generation .
after its introduction , this heating mechanism was studied by several authors : @xcite applied it to quark stars , @xcite made an order - of - magnitude estimate of the effects of superfluidity , and @xcite studied the heating due to compositional transitions in the crust due to spin - down compression .
all studies agree that rotochemical heating is particularly important for old neutron stars with fast rotation and low magnetic fields , features that are characteristic of millisecond pulsars ( msps ) .
the most striking prediction associated with rotochemical heating is that , if the spin - down timescale is substantially longer than any other timescale involved ( with the likely exception of magnetic field decay ) , the star arrives at a _ quasi - equilibrium state _ , in which the temperature depends only on the current , slowly changing value of @xmath1 ( the product of the angular velocity and its time derivative ) , proportional to the spin - down power , and not on the star s previous history @xcite .
this provides a simple way to constrain the physics involved in theoretical models , once the spin parameters and observed surface temperature of a msp are known .    although the qualitative behavior of the thermal evolution of neutron stars with rotochemical heating is known , all previous studies made only order - of - magnitude estimates , ignoring the spatial structure of the star and thus not offering reliable predictions to be compared with observations .
also unknown is the dependence of the thermal evolution on eos and stellar mass .    in this work ,
we calculate the thermal evolution of msps with rotochemical heating , taking the structure of the star fully into account in the frame of general relativity , and using realistic eoss of dense matter . in order to do so
, we develop a general formalism to treat the evolution of the temperature and departures of chemical equilibrium . as a first approach
, we have chosen the simplest possible core composition , namely neutrons , protons , electrons , and muons ( @xmath2 matter ) , to see if we can explain observations with it , before invoking exotic particles .
the main difficulty of treating the spatial structure of the star is the lack of an expression for the spin - down compression in the frame of general relativity .
we develop such an expression with the aid of the slow - rotation approximation of @xcite .
the astrophysical situation to be modelled is the stage after the accretion - driven spin - up of the pulsar has finished .
although it is generally believed that crusts and cores of neutron stars have superfluid components ( e.g. , @xcite ) , in this work we ignore superfluidity .
our plan is to make an extensive study of rotochemical heating in non - superfluid stars and compare results with observations , in order to assess whether invoking superfluidity ( whose parameters are currently very uncertain ) is necessary to explain observations .
the structure of this paper is the following : in  2 , we develop our general formalism , as well as the expression for the spin - down compression rate . in  3 , we detail the input necessary for numerical calculations ( eos , reactions , heat capacities , etc . ) . in  4 , we describe our results , comparing our predictions with the recent detection of likely thermal ultraviolet emission from the nearest msp , psr j0437 - 4715 @xcite , and the lack of optical emission from the even closer `` classical '' pulsar psr j0108 - 1413 @xcite .
we also list some other msps for which the rotochemical heating luminosity might become detectable .
a short summary of our main conclusions is given in  5 .
it is conventional in neutron star cooling calculations to divide the star into two regions : a nearly isothermal _ interior _ , which ranges from the center out to density @xmath3 g @xmath4 , and a thin _ envelope _ , which ranges from @xmath5 to the surface and where a strong temperature gradient exists @xcite . since we are modelling the thermal evolution of a msp long after accretion has stopped , it is safe to assume that thermal relaxation from an initial non - uniform internal temperature profile has already occurred , so that the redshifted internal temperature , @xmath6 is uniform ( see , e.g. , @xcite ) . here
, @xmath7 is the time component of the metric of a non - rotating reference star , of which @xmath8 is the radial spherical coordinate .
although spherical symmetry is broken for a rotating star , we describe the latter as a perturbation to the corresponding non - rotating star ( with the same total baryon number ) .
we establish a lagrangian correspondence between the surfaces of constant @xmath8 in the non - rotating star with the constant - pressure surfaces of its rotating counterpart , on which all local thermodynamic quantities will be shown to be constant ( see   [ sec : perturbation ] ) .
the evolution of the internal temperature is given by the thermal balance equation @xcite , which for an isothermal interior is given by @xmath9,\ ] ] where @xmath10 is the total heat capacity of the star , @xmath11 is the total power released by heating mechanisms , @xmath12 is the total neutrino luminosity , and @xmath13 is the photon luminosity .
these quantities are calculated as @xmath14 respectively , where @xmath15 is the proper volume element , @xmath16 is the specific heat ( heat capacity per unit volume ) of each particle species , @xmath17 is the total neutrino emissivity contributed by reactions , @xmath18 is the total heating rate per unit volume , @xmath19 is the stefan - boltzmann constant , @xmath20 is the stellar coordinate radius , @xmath21 , @xmath22 is the effective radius as measured from infinity , and @xmath23 the redshifted effective temperature .
the surface temperature @xmath24 is obtained from the internal temperature by assuming an envelope model @xcite .
since the neutrino emissivity and heating rate are modified when the neutron star is out of chemical equilibrium @xcite , the evolution of the temperature depends on how strongly it departs from the beta equilibrium state . for @xmath2 matter
, this departure can be quantified by the chemical imbalances @xcite : @xmath25 where @xmath26 is the deviation from the equilibrium chemical potential of species @xmath27 , at a given pressure . since diffusion timescales are short compared with the evolutionary timescales to be considered @xcite , we assume uniform redshifted chemical potential deviations throughout the core : @xmath28 to obtain the time evolution of the chemical imbalances , we start by writing down the chemical potential of each particle species as a function of the number density of all particles : @xmath29 .
we assume small departures from chemical equilibrium , which can be expressed by requiring that @xmath30 . in this approximation ,
the departures from the equilibrium particle number densities @xmath31 are related to the @xmath32 by @xmath33 where the partial derivatives are evaluated at the beta equilibrium state . to eliminate the effect of particle diffusion between different regions of the star , we integrate equation  ( [ eq : dni_dmui ] ) over regions where free particles exist . after integration
, we obtain the deviation from the equilibrium number of particles @xmath34 as a function of the redshifted chemical potential lags : @xmath35 with @xmath36 where we have used equation  ( [ eq : delta_mu_infty ] ) .
since the @xmath37 do not depend on time , we invert and take the time derivative of equation  ( [ eq : dni_dmui ] ) , obtaining the evolution of the @xmath38 : @xmath39 the rate of change of @xmath34 is given by @xmath40 where @xmath41 is the change in the total number of particles of species @xmath27 due to reactions @xcite . here , @xmath42 is the net creation rate of particles of species @xmath27 per unit volume due to reaction @xmath43 .
it can be checked that the @xmath44 satisfy baryon number and charge conservation .
the quantity @xmath45 is the change in the equilibrium value of the total number of particles of species @xmath27 , which depends on the spin - down compression rate .
we defer its calculation to the next subsection .
the evolution of the redshifted chemical imbalances @xmath46 follows from equations  ( [ eq : eta1 ] ) , ( [ eq : eta2 ] ) , and ( [ eq : delta_mu_infty ] ) : @xmath47 equations  ( [ eq : dot_t ] ) , ( [ eq : dot_eta1 ] ) , and ( [ eq : dot_eta2 ] ) give a complete description of the thermal evolution of a neutron star with rotochemical heating and @xmath2 composition , given an expression for @xmath45 .      for slow enough rotation frequencies , the deviation of the star from the non - rotating configuration can be treated as a small perturbation , which we describe in terms of a lagrangian formalism . using the fact that the total baryon number @xmath48 of a star is conserved
, we can describe its interior in terms of surfaces of constant pressure @xmath49 enclosing a fixed number of baryons @xmath50 ( e.g. , the surface @xmath51 encloses @xmath48 baryons inside it ) . as the stellar rotation rate changes
, these surfaces will readjust their shapes , with a corresponding change in @xmath49 , but keeping @xmath50 constant . since there is a one - to - one relation between the enclosed number of baryons @xmath50 and the radial coordinate of the non - rotating configuration @xmath8 ( for fixed @xmath48 ) , we use them interchangeably to identify a given surface , relating them by @xmath52 .
the underlying assumption is that , on a surface of constant pressure , the other thermodynamical quantities are also constant .
this is straightforward for neutron - star matter in beta equilibrium ( described by a barotropic equation of state ) , but not so obvious for departures from this state ( which are not necessarily barotropic , as the pressure - density relation will also depend on the local particle abundances ) .
we show in appendix  [ sec : gradients ] that , in a uniformly rotating , perfect - fluid star in hydrostatic equilibrium , surfaces of constant pressure coincide with those of constant energy density and gravitational redshift , regardless of whether the equation of state is barotropic .. ] this , together with diffusive equilibrium ( uniform redshifted chemical potentials everywhere in the stellar interior ) , ensures that the chemical composition on each surface is also constant .
to quantify the spin - down compression rate and thus @xmath53 , we need an expression for the change in pressure with rotation frequency on each surface enclosing a constant @xmath50 , keeping @xmath48 constant .
@xcite developed a perturbative approach in the frame of general relativity , relating quantities on rotating and non - rotating stars with the same central energy density @xmath54 ( but , therefore , different @xmath48 ) , by expanding the metric of an axially symmetric non - rotating star in even powers of the angular velocity @xmath55 .
this perturbative approach is valid for stars rotating at frequencies much smaller than the kepler frequency @xmath56 , at which mass shedding from the equator occurs .
for realistic equations of state , this limiting frequency can be well approximated by means of empirical formulae @xcite .
the corresponding kepler periods , @xmath57 , are lower than the shortest msp periods measured to date [ see , e.g. , the atnf pulsar catalogue @xcite ] .
our task is to express our spin - down compression rate at constant @xmath48 in terms of hartle s results , which assume constant @xmath54
.    if the baryon number enclosed by a surface of constant pressure @xmath49 in a rotating star is @xmath58 , some partial - derivative manipulation yields @xmath59 where the subscripts denote quantities which are kept constant in carrying out the differentiation .
we note that @xmath60 and @xmath61 are just @xmath62 and @xmath63 , respectively , evaluated at the surface of the star ( @xmath51 ) . also , since we are perturbing around the non - rotating configuration , we can evaluate all the derivatives at @xmath64 .
@xcite found that the lowest - order change in a given stellar property with rotation is linear in the squared angular frequency @xmath65 , for a star with fixed central density .
we can thus write @xmath66 and consider it independent of @xmath65 . to find the number of baryons enclosed by a surface of constant pressure
, we recall the baryon number conservation law @xcite , and write @xmath67 where @xmath68 is the determinant of the metric , @xmath69 is the time component of the fluid s four - velocity , @xmath70 is the rest frame baryon number density , and @xmath71 is the spatial volume enclosed by the constant pressure surface ( throughout this section , we use @xmath72 ) . for the non - rotating case
, we have @xmath73 where @xmath74 is the radial component of the non - rotating metric . for the rotating case
, one has to express the perturbed metric of @xcite in the non - rotating coordinate system . keeping terms to order @xmath65 , we get @xmath75\\ \label{eq : ut_rot } u^t   & = & e^{-\phi}\left [ 1 - \xi\frac{d\phi}{dr } - h + \frac{1}{2}r^2\sin^2{\theta}{\bar{\omega}}^2e^{-2\phi}\right],\end{aligned}\ ] ] where @xmath76 , @xmath77 , and @xmath78 are metric - perturbing functions which are proportional to @xmath65 and can be separated into spherical and quadrupolar parts , @xmath79 is the difference between the star s rotation frequency and the frame - dragging frequency , and @xmath80 is a lagrangian displacement that relates surfaces of constant @xmath81 in rotating and non - rotating stars of the same @xmath54 @xcite .
the spherical polar angle is @xmath82 , and the mass enclosed inside the radius @xmath8 is @xmath83 . replacing equations ( [ eq : g_rot ] ) and ( [ eq : ut_rot ] ) in ( [ eq : baryon_law ] ) , keeping terms to order @xmath65 , taking the angular integral , and integrating by parts in @xmath8 , we get ) is not present , which is correct when integrating up to the surface , but not at intermediate layers . ]
@xmath84dy   + 4\pi r^2 e^{\lambda(\rho_c , r ) } \xi_0(\rho_c , r ) n(\rho_c , r).\end{aligned}\ ] ] the functions @xmath79 , @xmath85 , @xmath86 are obtained by solving first and second - order differential equations in the coordinate @xmath8 , as explained in @xcite , and which need a value for the rotation rate of the star @xmath55 and a non - rotating relativistic stellar structure as input .
the zero subscript on functions indicates their spherical part ( the quadrupolar part vanishes when performing the angular integral , and the function @xmath78 has no spherical part ) .
the remaining derivatives in equation ( [ eq : dpdomega2_final ] ) can be obtained from the non - rotating configuration .
we first have @xmath87 where @xmath88 can be calculated analytically from the relativistic stellar structure equations @xcite .
the other derivative has to be obtained by differentiating equation ( [ eq : n_norot ] ) : @xmath89 dy.\end{aligned}\ ] ] for numerical calculations , it is convenient to write the derivative in the first term of the right - hand side as @xmath90 in order to calculate the three derivatives with constant @xmath8 in the same way .    the resulting spin - down compression rate , equation ( [ eq : dpdomega2_final ] ) ,
is plotted in figure  [ fig : dpdomega2 ] , normalized by pressure and kepler frequency .
the left panel displays results obtained for different eoss , with fixed central pressure , while the right panel shows results for different stellar models and fixed eos .
as expected , the expression is negative , since spin - down compression means increase in @xmath49 with decrease of @xmath65 .
compression is strongest at the center of the star .
it can be shown that the central value is given by @xmath91 where @xmath92 is the adiabatic index , and the subscripts @xmath93 mean central values . since @xmath94 for the maximum mass non - rotating configuration , the divergence of the result for @xmath95 , as shown in the right panel of figure  [ fig : dpdomega2 ] , is easy to understand .
the graphs confirm the order - of - magnitude expression @xmath96 which implies that the perturbation is small as long as @xmath97 .    with our expression for the spin - down compression rate at hand , we can write the change in the equilibrium number of particles of species @xmath27 with time as @xmath98 where the derivative @xmath99 is calculated in beta equilibrium .
it is straightforward to check that the @xmath45 satisfy baryon number and charge conservation .
the variables @xmath50 and @xmath49 are clearly continuous functions of @xmath8 in the non - rotating star , regardless of the presence of a phase transition . in appendix
[ sec : continuity ] , we show explicitly that the same is true for @xmath100 . on the other hand
, the equilibrium particle concentrations may have a discontinuous jump @xmath101 at the transition , which would contribute a dirac delta function to @xmath102 .
therefore , there would be a finite contribution to @xmath45 from the transition , @xmath103_{trans}.\ ] ]    to order of magnitude , the fractional contribution of a phase transition to the total value of @xmath45 ( and to the total heating of the star , if the reactions are of the same kind everywhere ) is the ratio of the jump @xmath101 to the total variation in @xmath104 across the stellar core . in the present work , as described in the following sections , phase transitions are not important , as the only transition in one of our eoss is fairly weak in each integration step as the ratio of the ( not necessarily small ) increment in @xmath105 to the small increment in @xmath49 . ] .
in other cases , such as a star with a core of deconfined quarks of fairly uniform relative abundances , the phase transition may dominate the heating .    in order to maintain diffusive equilibrium , a finite number of particles per unit time , @xmath106 ,
must be brought to ( or , if negative , away from ) the phase transition much faster than the latter progresses due to the spin - down of the star .
therefore , a phase transition is the place where our assumption of diffusive equilibrium is most strongly put into question .
however , in the case of interest to us , the low temperatures ( implying long mean - free paths and therefore fast diffusion ) and the extremely slow evolution should still validate this approximation .
@xcite has shown that , due to the presence of two separate conserved charges ( baryon number and electric charge ) , phase transitions in neutron stars may not occur as a single discontinuity , but instead through a mixed phase , in which both phases coexist over a finite range in pressure ( or stellar radius ) . in this case , any discontinuities would be considerably weakened , if not fully eliminated , and we do not expect any special effects ( such as discussed above ) associated with the phase transitions .
to implement the formalism of  2 numerically for the non - superfluid case , we need as input an eos , neutrino emissivities , a specific heat , and an envelope model .
we then have to calculate the stellar structure and the corresponding integrals over the star .      since our formalism requires a detailed knowledge of the structure of the star , we need an eos to determine the spatial dependence of the thermodynamical quantities by solving the relativistic stellar structure equations @xcite .
in addition to the usual variables given in most published eoss ( pressure , energy density , chemical composition , and adiabatic index , all evaluated in chemical equilibrium ) , we need to know the partial derivatives of equation  ( [ eq : dni_dmui ] ) ( which are easy to compute if we know the energy density or energy per baryon of interacting particles as a function of baryon number density and proton fraction ) and the baryon effective masses for computing reaction rates and heat capacities
. we also want to cover a broad range of eoss , in order to assess the dependence of the quasi - equilibrium temperature on their properties .
given our requirements , we have used three sets of eoss :    1 .
the two most realistic eoss from @xcite ( apr ) for the core , supplemented with that of @xcite and @xcite , for the inner and outer crust , respectively ; 2 .   five representative eoss from @xcite ( bpal ) for the core , with the same eoss for the inner and outer crust as for the previous set ; and 3 .   a non - interacting fermi gas ( see , e.g. , @xcite ) for the whole star .
in table  [ tab : eos_par ] we list , for each eos , the mass @xmath83 , central energy density @xmath54 , coordinate radius @xmath20 , and effective radius as seen from infinity @xmath107 of the maximum - mass non - rotating configuration .
since in the apr set the thermodynamical quantities are obtained from an analytical fit to tabulated many - body calculations @xcite , care has to be taken when interpreting results obtained for densities higher than the highest tabulated value , since extrapolation can lead to significant errors .
this is the reason why we do not list the maximum - mass configuration for the a18 + @xmath108 eos .
another feature of this set is that the a18 + @xmath108 + uix * eos becomes non - causal at densities greater than @xmath109 , which is the central density of a star of @xmath110 , slightly below the maximum mass formally allowed for this eos , @xmath111 .
however , this eos is considered to be the state of the art among those derived from non - relativistic potentials , since it incorporates three - nucleon interactions and relativistic boost corrections @xcite .
for this reason , we will in what follows take it as our reference eos .
the eoss of the bpal set are labelled according to @xcite , who characterize them by the different magnitude of their nuclear incompressibility and density dependence of their symmetry energy .
the values for the fermi gas eos listed in table  [ tab : eos_par ] correspond to a central density slightly below the appearance of @xmath112 hyperons .
the a18 + @xmath113 + uix * eos has an additional feature : a phase transition associated with the appearance of a neutral pion condensate at a density @xmath114 g @xmath4 . in order to connect the normal phase with the pion - condensed one
, we assumed a maxwell transition , which implies a discontinuity in all thermodynamical quantities but the temperature , pressure and neutron chemical potential , resulting in an energy - density jump of 6.6% .
@xcite also considered the possibility of a mixed - phase transition ( e.g. , @xcite ) , concluding that ( for a standard , @xmath115 neutron star ) the shell containing the mixed phase is only @xmath116 m thick .
they also estimated the debye screening length of electrons to be smaller than the characteristic size of structures ( by a factor 6 - 12 ) , which would favor the maxwell transition as a better approximation to the real situation @xcite .
also listed in table  [ tab : eos_par ] are the kepler periods @xmath117 corresponding to most of the eoss , which are necessary in order to determine the extent to which the slow - rotation approximation holds for treating the spin - down compression .
the apr and pal values were computed with the empirical formula of @xcite .
this gives , for each eos , the maximum allowed rotation frequency in terms of the mass and radius of the maximum - mass non - rotating configuration .
the accuracy for non - causal eoss is somewhat lowered from its optimal value of 1.5% @xcite , which applies for a18 + @xmath118 + uix*. the reason why the entry corresponding to a18 + @xmath108 is empty is that its maximum mass configurations occur at a central density which requires substantial extrapolation , making the empirical formula an unreliable estimator of @xmath57 . for the fermi gas eos , we adopt the exact value calculated for a pure neutron gas by @xcite , since the empirical formula of @xcite does not work for this case .
so far , the fastest known msp is psr b1937 + 21 , with a period of 1.56 ms @xcite , which is about three times longer than kepler periods for all realistic eoss .      since we are not taking superfluidity into account , direct and modified urca reactions are the dominant neutrino emission mechanisms .
we have therefore neglected additional neutrino emission processes in the core and the crust of the star , and ignored possible modifications in urca rates due to the neutral pion condensate of the apr set ( see , e.g. , @xcite ) .
away from beta equilibrium , neutrino emissivities and net reaction rates per unit volume of urca - type reactions in non - superfluid matter are given by @xcite @xmath119 respectively , where @xmath120 is the neutrino emissivity in equilibrium due to reaction @xmath43 , @xmath121 is the chemical imbalance affected by reaction @xmath43 , @xmath122 is the local temperature , @xmath70 is the baryon number density , @xmath78 is boltzmann s constant , and @xmath123 and @xmath124 are dimensionless control functions which depend on whether the reaction is of direct of modified urca type .
these functions are given by @xcite @xmath125 where the subscripts @xmath126 and @xmath83 mean direct and modified urca , respectively .
a finite @xmath127 of either sign enhances neutrino emission due to the even nature of the functions @xmath123 .
the amount of energy released by each reaction of type @xmath43 is @xmath127 @xcite , thus the total energy dissipation rate per unit volume is @xmath128 with the signs defined by @xmath129 where @xmath130 is the rate per unit volume of the reaction that transforms the particle set @xmath48 to the set @xmath131 .
this sign convention ensures that reaction rates are enhanced in the direction which restores equilibrium , since the functions @xmath124 are odd .
the equilibrium emissivities of both direct and modified urca reactions can be written as @xmath132 where @xmath133 is a slowly varying function of @xmath70 , @xmath134 for direct urca , and @xmath135 for modified urca ( e.g. , @xcite ) .
this decoupling of the temperature from the spatial part introduces a great simplification . to obtain the non - equilibrium neutrino luminosities and net reaction rates
, we first define @xmath136 where we have used equations  ( [ eq : t_int ] ) , ( [ eq : l_nu ] ) and ( [ eq : delta_mu_infty ] ) , and where @xmath137 is the equilibrium neutrino luminosity due to reaction  @xmath43 . combining equations  ( [ eq : q_noneq ] ) , ( [ eq : dgamma_noneq ] ) , ( [ eq : xi_def ] ) and ( [ eq : int_qeq ] )
, we obtain @xmath138 we denote modified urca reactions with electrons and muons by @xmath43 = @xmath139 and @xmath140 , respectively , in each case adding the contributions of the neutron and proton branches @xcite .
direct urca with electrons or muons is denoted by @xmath43 = @xmath141 and @xmath142 , respectively . from the apr set , the only eos which allows direct urca is a18 + @xmath108 + uix * , exceeding the threshold for direct urca with electrons at @xmath143 g @xmath4 , which is the central density of a star of 2 @xmath144 .
direct urca with muons lies in the non - causal regime . with the exception of bpal31 , all the eoss of the pal set allow electron direct urca .
muon direct urca is forbidden for both bpal21 and bpal31 .    combining equations  ( [ eq : l_h ] ) , ( [ eq : delta_mu_infty ] ) , ( [ eq : dgamma_noneq ] ) , ( [ eq : q_h ] ) , ( [ eq : xi_def ] ) , and ( [ eq : int_dgamma ] ) , we obtain the heating luminosity corresponding to the reaction  @xmath43 @xmath145 we can combine the heating and cooling contributions of each reaction into a single expression , which will prove to be useful for physical insight . if we define @xmath146 we can write the difference of equations  ( [ eq : l_nu_j ] ) and ( [ eq : l_h_j ] ) as @xmath147 thus , the relative size of @xmath148 and @xmath149 will determine , through the functions @xmath150 , whether there is net cooling or heating due to the neutrino emitting reactions . since the functions @xmath150 are even , they do not depend on the sign of the @xmath151 .
we plot the functions @xmath152 for the direct and modified urca case in figure  [ fig : m_j ] . with @xmath153 , the constant term in @xmath123
gives the conventional cooling case .
as @xmath154 grows , the cooling is enhanced by the additional neutrino emission due to the chemical imbalance , reaching a maximum cooling at @xmath155 . for larger values of @xmath154 ,
the heating becomes important , completely balancing the cooling for @xmath156 . for very large values of @xmath154 , the heating dominates , growing as @xmath157 and @xmath158 in the direct and modified urca case , respectively . in this limit
, a fixed fraction of the energy released is emitted as neutrinos , 3/8 and 1/2 for modified and direct urca , respectively .
the remainder stays behind , heating the star .
since we want to simulate the evolution of a msp which has finished accretion , we use the fully accreted envelope model of @xcite to calculate the effective surface temperature @xmath159 where @xmath160 is the surface gravity in units of @xmath161 cm s@xmath162 , @xmath163 is the internal temperature in units of @xmath164 k , and @xmath165 .
the photon luminosity then follows from equation  ( [ eq : l_gamma ] ) .
an accreted envelope is more heat transparent than a pure iron one for @xmath166 k , leading to faster cooling , whereas for @xmath167 k the opposite happens @xcite .      to calculate the integral in equation  ( [ eq : capcal ] )
, we use the specific heat capacities at constant volume @xmath16 for degenerate , non - superfluid fermions ( see , e.g. , @xcite ) . in analogy with the neutrino emissivity
, we can separate the temperature dependence from the spatial one . from equation  ( [ eq : capcal ] )
, we can define @xmath168 where @xmath169 and @xmath170 are the fermi momentum and effective mass of species  @xmath27 , respectively .
the latter can be obtained analytically for the apr and pal eoss ( see , e.g. , @xcite ) , while for the noninteracting fermi gas we use @xmath171 , where @xmath172 is the chemical potential of species  @xmath27 .
the sum is performed over all particle species and the integral is done over the region where each particle species exists .
we take only free particles into account , neglecting the heat capacity due to the lattice of ions in the crust .      in this work
, we are neglecting processes in the crust which might also modify the total number of neutrons , protons , or electrons .
we may ask here if this choice is a good approximation for computing the evolution of the chemical imbalances , since the short diffusion timescale limit means that the chemical equilibrium everywhere in the star can be restored by reactions occurring wherever these are fastest .
@xcite calculated the heating of a neutron star due to spin - down compression in the crust .
they found that neutron absorptions in the inner crust are the dominant non - equilibrium process , producing a very small heating rate @xmath173 , where @xmath174 is the spin - down power .
although they did not consider neutron diffusion between different layers in the star , and modelled the spin - down compression as a one - dimensional process , we take their result to suggest a very small contribution to the total heating .
since the inner crust accounts for about 2% of the total baryons in the star and
the neutrino emissivities of the core overwhelm crustal emission processes , we integrated equations  ( [ eq : bij_def ] ) only over the core and neglected processes in the crust , under the assumption that the error introduced would be small .
our results confirm our assumption , since for the more realistic eoss we obtained heating rates due to core processes much bigger than iida & sato s results ( see equation [ eq : ratio_edot ] ) .
we have to warn that neglecting the crust may not be a good approximation if superfluidity is included , since non - equilibrium processes in the crust can even overtake suppressed urca reactions in the core .
it has also been suggested by @xcite that direct urca reactions could occur at the crust - core interface .
we have not explored this possibility .
however , we are aware that even a small amount of direct urca reactions in the inner crust would have a significant effect on the evolution of chemical imbalances .      finally , we write the coupled equations for the time evolution of the temperature and the chemical imbalances explicitly for the non - superfluid case with @xmath2 composition , taking direct or modified urca reactions into account . for the internal temperature , from equations  ( [ eq : dot_t ] ) , ( [ eq : capcal ] ) , ( [ eq : l_h ] ) , ( [ eq : l_nu ] ) , ( [ eq : l_gamma ] ) , ( [ eq : net_cooling ] ) , and ( [ eq : c_int ] ) , we get @xmath175.\end{aligned}\ ] ] for the chemical disequilibria , using equations  ( [ eq : dni_dmui ] ) , ( [ eq : dni_def ] ) , ( [ eq : dot_ni ] ) , ( [ eq : dot_eta1 ] ) , ( [ eq : dot_eta2 ] ) , and ( [ eq : nieq ] ) , we obtain @xmath176 where we have defined @xmath177 and @xmath178 typical values for the different constants are plotted in figure  [ fig : ctes ] . since @xmath179
, equations  ( [ eq : evolucion_eta1 ] ) and ( [ eq : evolucion_eta2 ] ) have a positive term , proportional to the spin - down power , which arises from the change in the equilibrium concentrations of each particle species with spin - down and makes the chemical imbalances grow .
the remaining terms , opposite in sign to @xmath180 and @xmath181 , account for the effect of reactions trying to restore beta equilibrium .
equation  ( [ eq : evolucion_ti ] ) has the photon luminosity as a negative term , and the net contribution of neutrino reactions , equation  ( [ eq : net_cooling ] ) , which may be positive or negative depending on the absolute value of @xmath182 and @xmath183 . for numerical calculations ,
the evolution of @xmath184 was computed assuming magnetic dipole braking with no field decay , relating magnetic field , rotation period , and period derivative by the conventional formula @xmath185 g , where @xmath49 is measured in seconds .
a typical solution of equations  ( [ eq : evolucion_ti ] ) , ( [ eq : evolucion_eta1 ] ) , and ( [ eq : evolucion_eta2 ] ) is shown in the left panel of figure  [ fig : evol_xi ] .
first , the temperature starts falling as in the conventional cooling case , while the chemical imbalances grow due to spin - down compression . at some point ,
the ratios @xmath182 and @xmath183 ( right panel of figure  [ fig : evol_xi ] ) are big enough , so that the net contributions of neutrino reactions , equation  ( [ eq : net_cooling ] ) , are positive , and their sum exceeds @xmath186 .
the temperature then starts rising .
the chemical imbalances continue to grow , more slowly than the temperature and hence reducing @xmath182 and @xmath183 .
finally , the star can arrive at a _ quasi - equilibrium _ state , where the rate at which spin - down modifies equilibrium concentrations is the same as the rate at which reactions drive the system towards the equilibrium configuration , with heating and cooling balancing each other @xcite .
the subsequent thermal evolution can then be approximated by the simultaneous solution of the equations @xmath187 and @xmath188 .
we discuss this _ quasi - equilibrium solution _ in detail in  4.2 .
the pure modified urca case ( @xmath189 ) is the simplest to analyze .
figure  [ fig : evol_tdmuvarios ] shows the thermal evolution of the same star as in figure  [ fig : evol_xi ] , this time with different initial conditions of internal temperature ( left ) and chemical imbalances ( right ) .
the spin - down parameters were chosen so that the arrival at the quasi - equilibrium state is clearly visible , the short - dashed lines being the quasi - equilibrium solutions .
we note that , in both cases , the value of the temperature at the quasi - equilibrium state and the time required to arrive do not depend on initial conditions .
if direct urca reactions are open , solutions change noticeably . in figure
[ fig : evol_urca ] we show the thermal evolution of three different stellar models built with the same eos , with identical initial conditions and spin - down parameters .
the 1.4@xmath144 star is the same as that of figure  [ fig : evol_tdmuvarios ] , with only modified urca reactions in operation .
the 2@xmath144 star has electron direct urca operating , with the corresponding muon reaction forbidden .
we see that the minimum temperature is higher and occurs at a much earlier time , after which the system reaches a  metastable " quasi - equilibrium state , corresponding to partial equilibration between @xmath190 and @xmath191 .
however , @xmath192 continues to grow , leading the system to a final quasi - equilibrium state , determined by muon modified urca .
the 2.17@xmath144 star has direct urca with electrons and muons .
this time the final quasi - equilibrium state occurs even earlier than the  metastable " quasi - equilibrium of the 2@xmath144 star , with the temperature at about the same level .      when only modified urca reactions operate , it is possible to solve analytically for the quasi - equilibrium values of the photon luminosity @xmath193 and chemical imbalances @xmath194 and @xmath195 , as function of stellar model and current value of @xmath184 .
what makes these approximations possible is the high value of @xmath182 and @xmath183 at quasi - equilibrium . from the right panel of figure
[ fig : evol_xi ] , we see that , at quasi - equilibrium ( @xmath196 yr ) , @xmath197 .
this enables us to ignore all but the greatest power in the functions @xmath198 and @xmath199 ( equations [ eq : f_m_def ] , [ eq : h_m_def ] , and [ eq : m_j_def ] ) : @xmath200 the error in this approximation , for @xmath201 , is less than one part in 100 .    to solve for @xmath193 , we set @xmath202 and use equations  ( [ eq : c_h_def ] ) and ( [ eq : c_m_def ] ) to rewrite ( [ eq : evolucion_ti ] ) , ( [ eq : evolucion_eta1 ] ) , and ( [ eq : evolucion_eta2 ] ) as @xmath203 solving for @xmath194 and @xmath195 in equations  ( [ eq : dot_eta1_qe ] ) and ( [ eq : dot_eta2_qe ] ) , we get @xmath204 where we have used equations  ( [ eq : w_npe_def ] ) and ( [ eq : w_npm_def ] ) and the fact that @xmath205
. both chemical imbalances at quasi - equilibrium are positive , since @xmath184 and @xmath206 are negative .
we finally replace equations ( [ eq : eta1_qe ] ) and ( [ eq : eta2_qe ] ) in ( [ eq : dot_ti_qe ] ) , to obtain the photon luminosity : @xmath207|\omega \dot{\omega}|^{8/7}.\ ] ] an interesting feature of equation  ( [ eq : l_gamma_qe ] ) is that it does not depend on envelope model .
the reason is that the limit @xmath208 implies that reactions are completely determined by chemical imbalances , eliminating any dependence on internal temperature .
the energy released is divided in fixed fractions between neutrino emission ( @xmath209 = 3/8 $ ] ) and heating ( @xmath210 ) , the latter being radiated entirely as photons .
note also that the ratio @xmath211 depends only on eos and stellar model , and is of order unity , thus we will make reference to either of them as @xmath212 .    for the entire range of eoss and stellar models ,
we can write @xmath213 where @xmath214 is the period derivative measured in units of @xmath215 , and @xmath216 the period in ms . values for different stellar models built with the a18 + @xmath108 + uix * eos are shown in figure  [ fig : ctes ] . using equation  ( [ eq : l_gamma ] ) , we can rewrite equation  ( [ eq : l_gamma_plnum ] ) in terms of the effective temperature : @xmath217 this time , the uncertainty due to eos and stellar model is smaller than 25% .
we can also write the ratio of quasi - equilibrium photon luminosity ( and thus total heating rate ) to the spin - down power @xmath174 @xmath218 it is also possible to calculate @xmath219 for the quasi - equilibrium state . its spin - down power dependence , @xmath220 where @xmath221 is the exponent of the envelope model of @xcite ( see  [ sec : envelope ] ) , shows that increasing @xmath222 reduces @xmath223 , decreasing the accuracy of the analytical approximation , although very slowly .      although we did not include hyperons in our calculations
, we may ask how their presence could modify the thermal evolution . the first particles to appear after muons are probably the @xmath224 and @xmath225 hyperons , which require the introduction of two additional chemical imbalances : @xmath226 once hyperons are present , the following non - leptonic reactions are open ( e.g. , @xcite ) : @xmath227 reaction ( [ eq : nnsp ] ) proceeds via weak interactions , since it does not conserve strangeness , while ( [ eq : nnln ] ) proceeds via strong interactions . both reactions
have timescales several orders of magnitude shorter than beta processes @xcite .
therefore , imbalances ( [ eq : eta_nnsp ] ) and ( [ eq : eta_nl ] ) remain small compared to @xmath228 and @xmath229 , and reactions ( [ eq : nnsp ] ) and ( [ eq : nnln ] ) contribute negligibly to the total heat generation . moreover , since chemical imbalances associated to direct or modified urca reactions with hyperons are linear combinations of @xmath230 , @xmath231 , @xmath228 , and @xmath229 , the latter two will determine the heat generation through both nucleon and hyperon reactions .    in order to assess the importance of including the urca processes involving hyperons in addition to the nucleon processes ,
consider , as an example , the @xmath112 direct urca reactions @xmath232 whose associated chemical imbalance is @xmath233 .
their net effect on equations  ( [ eq : evolucion_eta1 ] ) and ( [ eq : evolucion_eta2 ] ) is to enhance the electron direct urca rate according to @xmath234 , reducing the chemical imbalance and thus the stellar surface temperature .
this correction is small , since direct urca reactions with hyperons are at least a factor 5 weaker than their nucleon analogs @xcite . for the pure modified urca case
, it can be checked from equation  ( [ eq : l_gamma_qe ] ) that the correction to the surface temperature introduced by adding several reactions involving hyperons is a factor @xmath235^{1/28}$ ] , where @xmath236 and @xmath237 are the nucleon and hyperon urca luminosities , respectively .
thus , the corrections due to hyperons in a purely modified urca or purely direct urca scenario are fairly negligible .
the only important effect of hyperons in the context of rotochemical heating is that , in their presence , conditions for direct urca processes may become more easily satisfied , as is the case with reactions involving @xmath238 @xcite .
the steep increase of the proton concentration with the appearance of @xmath112 hyperons ( e.g. , @xcite ) can also cause the condition for nucleon direct urca to be satisfied at lower densities than if hyperons are excluded from the models .
we may now ask which of the known msps are likely to be in the quasi - equilibrium state . to answer this question , we need to know how long it takes to arrive at the quasi - equilibrium state , and if real pulsars are older than this time .
since the negative  equilibration " terms due to reactions in equations  ( [ eq : evolucion_eta1 ] ) and ( [ eq : evolucion_eta2 ] ) are important only for imbalances near their equilibrium values ( recall the limit @xmath239 ) , we can assume that the initial evolution of @xmath212 is due only to @xmath1 : @xmath240 integrating over time , we get @xmath241.\ ] ] for the true solution to have reached the quasi - equilibrium state , this approximate solution must exceed @xmath242 .
using equations  ( [ eq : eta1_qe ] ) or ( [ eq : eta2_qe ] ) , together with equation  ( [ eq : eta_ell_aprox ] ) , we may express this condition as an upper limit on the initial spin period @xmath243 with @xmath244 ( here and in the rest of this subsection , we assume that only modified urca processes are active . ) given the current spin parameters of a pulsar , condition  ( [ eq : p0_qe ] ) gives the highest initial period it can have had , so that enough time has elapsed for it to have reached the quasi - equilibrium state .
for msps , the constant @xmath48 is generally small . for a 1.4@xmath144 star built with the a18 + @xmath108 + uix * eos
, we get @xmath245 values for other stellar masses are shown in figure  [ fig : ctes ] . from equation  ( [ eq : etapto_aprox ] ) , we can estimate a characteristic timescale for equilibration as @xmath246    @xcite calculated the thermal evolution of neutron stars with rotochemical heating for different values of the magnetic field , using the magnetic dipole braking model with no field decay .
his conclusion was that a necessary condition for the quasi - equilibrium solution to be a good approximation to the exact one is that the spin - down timescale has to be longer than that for the growth of chemical imbalances . to quantify this , we estimate the spin - down timescale as the characteristic age , @xmath247 and take the ratio of @xmath248 to it , getting @xmath249 figure  [ fig : evol_chiteqvarios ] shows the internal temperature of our conventional star after the rise from the minimum temperature , for different initial values of @xmath48 , with fixed initial @xmath248 .
dotted lines are the quasi - equilibrium solutions at the instantaneous value of @xmath250 .
the figure suggests that , for high values of @xmath48 , exact and quasi - equilibrium solutions depart from each other , whereas for small @xmath48 , solutions remain `` together '' .
in fact , all solutions depart from each other , the departure being slower with time for smaller @xmath48 . for
more general spin - down laws with arbitrary braking index @xmath251 , it can be shown analytically that this is true for @xmath252 . the increase in @xmath253 , where @xmath254 is the exact solution and @xmath255 the quasi - equilibrium one ,
is roughly linear with time . for @xmath256 at @xmath257 , and assuming a standard braking index @xmath258 , @xmath259 is @xmath260 lower than @xmath254 at @xmath261 , the difference growing 2% every @xmath262 yr . rewriting equation  ( [ eq : a_pppto ] ) in terms of the magnetic field , @xmath263 where @xmath264 is the surface magnetic field in units of @xmath262 g , we can easily understand the results of @xcite , since increasing @xmath131 with fixed @xmath265 increases the initial value of @xmath48 .
relating @xmath48 at the present time and at @xmath257 for an arbitrary , but constant , braking index , @xmath266 shows that the former is greater than the latter for all reasonable values of @xmath70 .
thus , a small current value of @xmath48 is a good indicator that the star remains close to the quasi - equilibrium state , as long as it satisfies @xmath267 .
since constraints on the initial periods @xmath268 of a few msps have been obtained from the cooling ages of their white dwarf companions @xcite , we can check which objects are likely at the quasi - equilibrium state . in table
[ tab : predictions ] , further discussed in  [ sec : othermsps ] , we show the range in @xmath268 for some msps obtained by @xcite .
we made no distinction between the @xmath269 and @xmath258 case , and adopted the widest possible range in @xmath268 . in table
[ tab : predictions ] we also show the value of @xmath48 and the upper limit on initial spin period @xmath270 required for the msp to be currently in the quasi - equilibrium state . among msps with constraints on initial periods , psr j1012
+ 5307 is the only one not satisfying the conditions .
the small values of @xmath48 imply that any initial periods substantially shorter than the current ones will leave the msps in equilibrium by the present time .
recently , @xcite measured ultraviolet emission from psr j0437 - 4715 , the most probable explanation being that it corresponds to thermal radiation .
since this object satisfies @xmath271 with small @xmath48 , we assume that it is already at the quasi - equilibrium state . in figure
[ fig : ts_rinf ] , we plot the blackbody fit of @xcite as dashed lines , which correspond to 68% and 90% confidence contours .
we overplot our values of @xmath272 for the spin parameters of psr j0437 - 4715 @xcite , as function of @xmath107 , for the eoss listed in table  [ tab : eos_par ] .
the bold lines show the range in @xmath107 corresponding to the mass constraint of @xcite , @xmath273 .
the complete mass range for each eos goes from 1@xmath144 to 0.95@xmath274 .
the upper limit was chosen to avoid the divergence in the spin - down compression rate for masses near @xmath274 ( see figure  [ fig : dpdomega2 ] ) , which makes @xmath272 also diverge due to @xmath206 ( see equation  [ eq : l_gamma_qe ] ) .    the best agreement with observations
is reached if only modified urca reactions are present .
when direct urca reactions open , there are two abrupt drops in @xmath272 with increasing stellar mass , the first one ( @xmath275 ) due to electron direct urca , and the second ( @xmath276 ) due to muon direct urca .
this can be seen in figure  [ fig : ts_rinf ] for the curves calculated with bpal11 and a18 + @xmath108 + uix * , as examples of only electron direct urca , and in the curves corresponding to bpal32 and bpal33 , which have electron and muon direct urca .
perfect agreement could be obtained by allowing masses within @xmath277 of the maximum - mass non - rotating configuration , for the eoss with only modified urca .
however , we consider this scenario very unlikely .    using the eoss which have only modified urca reactions within the allowed mass for psr j0437 - 4715 , we can constrain the predicted effective temperature to the narrow range @xmath278 k , about 20% lower than the blackbody fit of @xcite .
there are three possible reasons why we are not matching their results :    1 .
we are not taking superfluidity into account .
this would reduce urca reaction rates , lengthening the equilibration timescale and raising the quasi - equilibrium temperature @xcite .
we are neglecting other heating mechanisms ( some of them directly related to superfluidity ) , which could further raise the temperature at any stage in the thermal evolution . nonetheless , in msps , all proposed mechanisms are less important than rotochemical heating @xcite .
3 .   the thermal spectrum could deviate from a blackbody , in the same way as the isolated neutron star rx j1856 - 3754 , which has a well - determined blackbody x - ray spectrum that underpredicts the optical flux @xcite , indicating a more complex spectral shape of its thermal emission .
kargaltsev and collaborators stress the fact that psr j0437 - 4715 , despite its much larger spin - down age , has a higher surface temperature than the upper limit for the younger , `` classical '' pulsar j0108 - 1431 , @xmath279 k , inferred from the optical non - detection by @xcite . in this regard , we note that the spin - down power ( @xmath280 ) of the latter pulsar is 680 times lower than that of j0437 , making rotochemical heating substantially less important .
its equilibration timescale , according to equation ( [ eq : tau_eq_def ] ) , is @xmath281 yr , longer than the age of the universe and certainly much longer than the spin - down age of the pulsar .
thus , it is not expected to be even close to reaching its quasi - equilibrium state ( although it is old enough to have lost its initial heat content ) .
its actual temperature should therefore be substantially smaller than the predicted quasi - equilibrium surface temperature , which is already @xmath282 times lower than that of j0437 , well below the observational upper limit .      in table
[ tab : predictions ] , we provide predictions of @xmath272 ( assuming modified urca reactions and no superfluidity ) for several msps . since the high - energy part of their expected thermal spectrum is highly absorbed by neutral hydrogen in the interstellar medium , we chose to order them by decreasing predicted quasi - equilibrium flux in the low - energy ( rayleigh - jeans ) regime , @xmath283 where @xmath284 is the distance to the object . each @xmath285 entry in table  [ tab : predictions ] is scaled by the value for psr j0437 - 4715 .
the first group of msps after psr j0437 - 4715 are nearby , single msps , for which @xmath285 is expected to be larger than for the binary msps , for which the initial period has been estimated from their white dwarf companion ( see  [ sec : conditions ] ) .
we consider the former as the primary targets for observations analogous to those of @xcite .
however , even for psr j0437 - 4715 , the nearest and brightest msp known to date , the detection of thermal emission is rather difficult , so the observation of other objects of table  [ tab : predictions ] will be a real challenge .
the second group of msps are those with estimates for @xmath268 .
however , their low fluxes make the detection of their thermal emission impossible for current instruments .
optical detections are even more daunting than those in the ultraviolet .
the extrapolated rayleigh - jeans spectrum of psr j0437 - 4715 gives magnitudes @xmath286 , @xmath287 , and @xmath288 , which are completely overwhelmed by the white dwarf companion . the most promising single msps have about three times lower expected fluxes at earth , i.  e. , are another 1.2 magnitudes fainter , beyond the reach of current telescopes .
in this work , we have made an extensive study of the effect of rotochemical heating in non - superfluid millisecond pulsars .
we have set up a general formalism for the thermal evolution in the framework of general relativity , which takes the stellar structure fully into account and can be extended to include both superfluidity and exotic particles . a key ingredient in this formalism is the spin - down compression rate based on the slow - rotation approximation of @xcite .
the main consequence of rotochemical heating is the arrival at a quasi - equilibrium state , in which the effective temperature of a neutron star depends only on the current value of the spin - down power .
we argue that most of the known msps are very likely in the quasi - equilibrium state .
if only modified urca reactions are allowed , the quasi - equilibrium bolometric photon luminosity in the non - superfluid case can be well approximated by @xmath289 independent of the neutron star envelope model .
the influence of eos and stellar model on the quasi - equilibrium state is very weak , the only significant factor being the occurrence of direct urca reactions .
if they are open for both muons and electrons , quasi - equilibrium temperatures are low , as in the conventional cooling case .
if they are open only for electrons , the system arrives at a  metastable " quasi - equilibrium , after which it proceeds as if only modified urca reactions with muons were present .
the highest temperatures are reached when all direct urca reactions are forbidden .
even our highest predicted quasi - equilibrium effective temperatures are lower than the blackbody fit of @xcite to the uv emission of psr j0437 - 4715 by about 20%
. the inclusion of superfluidity will likely raise our predicted temperatures , being the subject of future work .
we thank g.  pavlov and o.  kargaltsev for letting us know about their work in advance of publication and for kindly providing the data for figure  [ fig : ts_rinf ] .
the authors are also grateful to m.  taghizadeh , m.  van kerkwijk , r.  mignani , d.  page , m.  catelan , m.  a.  daz , c.  dib , and p.  jofr for discussions that benefited the present paper , and an anonymous referee for thoughtful comments that improved its final version .
this work made extensive use of nasa s astrophysics data system service , and received financial support from fondecyt through grant # 1020840 .
in this appendix , we show that , in uniformly rotating , relativistic stars in hydrostatic equilibrium , composed of a perfect fluid , the surfaces of constant pressure and those of constant energy density coincide , as in the newtonian case , even if the equation of state is not barotropic .
in addition , the gravitational redshift is constant on the same surfaces .
this result , together with the diffusive equilibrium condition , allows us to describe the spin - down compression in terms of lagrangian changes of thermodynamic variables on isobaric surfaces enclosing a fixed baryon number , since all thermodynamical quantities are constant on each isobar . as a by - product
, we also provide an alternative derivation of the equation of hydrostatic equilibrium for uniformly rotating , relativistic stars .    the assumption a perfect ( i.e. , non - dissipative ) fluid , whose stress - energy tensor depends only on the energy density @xmath81 and ( isotropic ) pressure @xmath49 , is of course not strictly true in the astrophysical situation of interest in this paper , in which entropy is generated by non - equilibrium weak interactions , particle diffusion , and heat conduction , and lost ( from the star ) through the emission of neutrinos and photons .
however , the total energy dissipated ( and eventually lost ) along the star s lifetime is much smaller than that associated with the mass , random motions , and interactions of its particles , as well as its rotation , and the time scale of the release of this energy is many orders of magnitude longer than the dynamical time and the rotation period of the star .
therefore , the contribution of these dissipative processes to the energy and momentum fluxes is extremely small and can be neglected in the stress - energy tensor .
this is in line with the usual approach ( also used in neutron star cooling calculations and in the rest of this paper ) of neglecting thermal effects when calculating the structure of the star , which is regarded as a fixed background on which thermal processes take place .
the metric of a stationary , axially symmetric system can be written in the form ( see , e.g. , @xcite ) : @xmath290 where the metric coefficients are functions of @xmath8 and @xmath82 only .
the components of the 4-velocity of any fluid element of the uniformly rotating star are @xmath291 , and @xmath292 , with @xmath293\right)^{-1/2}.\ ] ] we note that @xmath294 is the gravitational redshift factor , which corrects the energy of a freely moving particle , as it is measured by observers co - moving with fluid elements on different surfaces inside the star . , where @xmath295 are the covariant components of the particle s 4-momentum . due to the symmetries of the metric
, the expression within parenthesis is conserved along the particle s world line , and is its redshifted energy . ]
the energy - momentum conservation equation can then be written as @xmath296 where @xmath297 are the components of the stress - energy tensor of a perfect fluid . applying to this equation a projection operator orthogonal to the four - velocity , @xmath298 @xcite
, we get @xmath299 since the only non - vanishing terms are those with @xmath300 and @xmath301 taking the values @xmath302 , the relevant connection coefficients are @xmath303 hence , equation  ( [ eq : cons_enmom ] ) yields @xmath304 using equation  ( [ eq : ut_def ] ) , we arrive at the equation of hydrostatic equilibrium @xmath305 which is in agreement with the result of @xcite , for the special case of rigid rotation .
this shows that the redshift factor @xmath294 is also constant on the surfaces of constant @xmath49 .
thus , we can write @xmath306 which implies that @xmath81 is also constant on the same surfaces .
in this appendix , we show that each term of the spin - down compression rate , equation  ( [ eq : dpdomega2_final ] ) , is continuous across a phase transition , even if the energy density and baryon number density change discontinuously .
we use units with @xmath307 .
the discontinuous quantities in the derivative @xmath308 , equation  ( [ eq : dndp_def ] ) , are @xmath70 and @xmath88 .
the gibbs free energy per unit volume for @xmath2 matter in beta equilibrium at @xmath309 is @xmath310 since , at a phase transition , the pressure @xmath49 and neutron chemical potential @xmath311 are continuous , the fractional size of the number density jump can be written as @xmath312 where @xmath313 is the positive jump in number density . from the relativistic stellar structure equations , we have @xmath314 the fractional size of the corresponding jump across the phase transition
is @xmath315 thus , from equations  ( [ eq : dn_n ] ) and ( [ eq : ddpdr_dpdr ] ) , we see that the quantity @xmath316 is continuous across a phase transition , and so is @xmath317 .    regarding the derivative @xmath318 , equation  ( [ eq : dndomega2_n0 ] ) , since the jump in @xmath70 is integrated in equation  ( [ eq : n_norot ] ) , @xmath319 is continuous . for @xmath320 , the jump in @xmath70 due to the last term of equation  ( [ eq : n_rot ] )
is cancelled by the contribution from @xmath321 in the integral . using @xmath322 where @xmath301 is dirac s delta function and @xmath323 , with @xmath324 the pressure at which the phase transition occurs
, we can write @xmath325 where @xmath326 is the step function . equation  ( [ eq : disc_dndomega2 ] ) is equal to minus the jump ( increasing @xmath8 ) of the last term in equation  ( [ eq : n_rot ] ) .
finally , for the derivative @xmath327 , the contribution from @xmath328 to the integral in equation  ( [ eq : dndn0_p ] ) cancels out the discontinuity due to @xmath70 in the first term . using @xmath329 where @xmath330 , and @xmath331 the radial coordinate at which @xmath332 , we can write @xmath333 which is the negative of the jump ( increasing @xmath8 ) of the first term in equation  ( [ eq : dndn0_p ] ) , since the sign of @xmath334 is also @xmath335 .    the right panel of figure  [ fig : dpdomega2 ] shows the spin - down compression rate inside the core of several stellar models built with the a18 + @xmath108 + uix * eos , which has a phase transition . the occurrence of the transition can be seen as a ( fairly harmless ) discontinuity in the derivative of each curve .
lccccc a18 + @xmath108 & 1.55 & 1.86 & 9.81 & 13.42 & - + a18 + @xmath108 + uix@xmath336 & 2.19 & 2.78 & 9.97 & 16.79 & 0.51 + bpal 11 & 1.42 & 4.45 & 8.42 & 11.86 & 0.54 + bpal 21 & 1.70 & 3.46 & 9.33 & 13.69 & 0.56 + bpal 31 & 1.91 & 2.86 & 10.10 & 15.18 & 0.59 + bpal 32 & 1.95 & 2.66 & 10.56 & 15.64 & 0.63 + bpal 33 & 1.97 & 2.53 & 10.92 & 15.94 & 0.66 + fermi gas & 0.62 & 1.10 & 12.77 & 13.80 & 0.98    lccccccccc j0437 - 4715 & 5.76 & 1.86 & 0.14 & 0.72 & 1 & 0.17 & 5.33 & 2.45.3 & 1 , 2 + j1024 - 0719 & 5.16 & @xmath3371.85 & @xmath3370.25 & @xmath3370.79 & @xmath338 0.35 & @xmath3370.14 & @xmath3394.82 &  & 3 , 4 + j2124 - 3358 & 4.93 & 1.23 & 0.27 & 0.74 & 0.28 & 0.13 & 4.64 &  & 3 , 4 + j0030 + 0451 & 4.87 & @xmath3371.00 & 0.32 & @xmath3370.70 & @xmath3370.19 & @xmath3370.12 & @xmath3394.60 &  & 4 , 5 + j1744 - 1134 & 4.07 & 0.71 & 0.36 & 0.74 & 0.15 & 0.09 & 3.90 &  & 6 + b1257 + 12 & 6.22 & 4.26 & 0.45 & 0.86 & 0.12 & 0.22 & 5.63 &  & 3 , 4 + j0034 - 0534 & 1.88 & @xmath3370.67 & 0.53 & @xmath3371.41 & @xmath3370.14 & @xmath3370.03 & @xmath3391.85 & @xmath3371.4 & 2 , 4 , 7 + j1012 + 5307 & 5.26 & 1.34 & 0.41 & 0.71 & 0.11 & 0.14 & 4.93 & @xmath3395.1 & 2 , 4 , 8 + b1855 + 09 & 5.36 & 1.74 & 0.90 & 0.76 & 0.03 & 0.15 & 5.00 & @xmath3374.6 & 2 , 3 , 9 + j1713 + 0747 & 4.57 & 0.81 & 1.10 & 0.70 & 0.02 & 0.11 & 4.35 & 2.23.1 & 2 , 10 + j1640 + 2224 & 3.16 & 0.16 & 1.15 & 0.60 & 0.01 & 0.05 & 3.08 & @xmath3391.6 & 2 , 4 , 11 + j2019 + 2425 & 3.93 & @xmath3370.70 & 1.49 & @xmath3370.76 & @xmath3370.01 & @xmath3370.08 & @xmath3393.87 & 0.93.9 & 2 , 4 , 12