it is well known that the cooling of a neutron star ( ns ) is strongly affected by the relationship between the internal temperature of the star , @xmath5 , and its effective surface temperature @xmath0 . throughout this paper
, @xmath5 denotes the temperature at the outer boundary of the isothermal internal region .
the @xmath5@xmath0 relationship is determined by equation of state ( eos ) and thermal conductivity of matter in the outer ns envelope . for non - magnetized nss with envelopes composed of iron
, this relationship has been thoroughly studied in a classical article of gudmundsson et al .
( 1983 , hereafter gpe ) .
several papers ( hernquist 1984 , 1985 , van riper 1988 , 1991 , schaaf 1988 , 1990 ) have considered strongly magnetized ns envelopes .    in this article , we will reconsider thermal transport in non - magnetized envelopes of nss and extend the results of the preceding authors in two respects .    first
, we analyze matter composed not only of iron , but of light elements as well .
the light elements can be provided by an accretion from a supernova remnant ( e.g. , chevalier 1989 , 1996 , brown & weingartner 1994 ) , from interstellar medium ( e.g. , miralda - escud et al .
1990 , blaes et al .
1992 , nelson et al .
1993 , morley 1996 ) , from a distant binary component , or by comets .
freshly accreted matter burns then into heavier elements ( he , c , o , fe ) while sinking within the ns .
recent multiwavelength observations of the geminga pulsar ( 1e 0630 + 17.8 ) suggest a possible h or he cyclotron feature in its spectrum ( bignami et al .
if confirmed , it may be a direct observational evidence of the presence of the light elements in the pulsar atmosphere .
chemical composition affects the eos and thermal conduction , and , therefore , the thermal structure and cooling of a ns . as a reference case
, we reconsider the outer ns envelopes composed of iron .
second , we implement new , advanced theoretical data on eos and thermal conductivity of dense matter . specifically , we employ the opacity library ( opal ) radiative opacities for h , he , and c , improved considerably with respect to the los - alamos opacities used in the previous studies . in the outermost
ns layers , we also implement the modern opal eos for fe and the saumon
chabrier eos for h and he .
we use improved thermal conductivities of degenerate electrons . for solidified matter
, we employ the thermal conductivity due to the electron - phonon scattering obtained with the inclusion of the debye
waller factor . for liquid matter
, we recalculate and implement the thermal conductivity due to coulomb electron - ion and electron - electron collisions .
the electron - ion scattering is described with the exact ( non - born ) coulomb cross sections and with the ion structure factors calculated when taking into account the response of the electron background .
the new physics input allows us to extend the results of previous studies to colder nss , with @xmath0 down to 50000  k.    the physics input is described in sec .  2 .
in sec .  3 , we calculate the @xmath5@xmath0 relationships for non - accreted and partly accreted nss and analyze the sensitivity of the results to the uncertainty in our knowledge of the electron thermal conductivity .
in addition , we simulate the cooling of nss with standard and enhanced neutrino energy losses .
we show that the cooling of a ns with the accreted envelope can be quite different from the cooling of a non - accreted ns .
this can change the conclusions on the internal structure of nss deduced from comparison of theoretical cooling curves with observations of ns thermal radiation .
useful analytical formulae for the physics input are given in the appendix .
consider thermal transport throughout the outer envelope of a ns that extends from the surface to the layers with the density @xmath6 g @xmath7 .
this envelope , of which study is the aim of the present paper , provides the main thermal insulation of the ns interior .
the envelope is thin ( @xmath8 m deep ) and contains very small fraction ( @xmath9 ) of the ns mass ( gpe ) . in these layers
, one can neglect the nonuniformity of the energy flux due to the neutrino emission and the variation of the gravitational acceleration .
then the temperature profile obeys the _ thermal structure equation _ ( gpe ) , which can be written as ( e.g. , van riper 1988 ) @xmath10 here , @xmath11 is the pressure , @xmath0 is the effective temperature , @xmath12 is the surface gravity , @xmath13 is the stellar mass , @xmath14 the stellar radius , @xmath15 the gravitational constant , and @xmath16 is the gravitational radius of the star .
furthermore , @xmath17 is the opacity , @xmath18 composed of the radiative opacity @xmath19 and the equivalent electron conduction opacity @xmath20 .
the latter is related to the electron thermal conductivity @xmath21 as @xmath22 where @xmath23 is the stefan  boltzmann constant .
the effective temperature @xmath0 is defined through the stefan s law , @xmath24 where @xmath25 is the local radiative luminosity .
the apparent luminosity measured by a distant observer is @xmath26 , and the apparent surface temperature inferred by the observer from the spectrum is @xmath27 ( e.g. , thorne 1977 ) .
we adopt the standard ( gpe ) outer boundary condition to eq .
( [ th - str ] ) by equating @xmath0 to the temperature @xmath28 at the stellar surface in the eddington approximation .
the radiative surface , determined in this way , lies at the optical depth @xmath29 , that is approximately at @xmath30 here and hereafter the subscript `` s '' denotes quantities taken at the surface , while the subscript `` b '' is assigned to those at the inner boundary @xmath31 .
previous calculations of the radiative transfer in the outer part of the atmosphere ( zavlin et al .
1996 ) have shown that , for the parameters of interest , there is no convective instability at @xmath32 which could invalidate the eddington approximation .
starting from the surface specified by eq .
( [ ps ] ) , we integrate eq.([th - str ] ) inward , using the classical 4-step 4th - order runge
kutta algorithm ( e.g. , fletcher 1988 ) . the integration step varies with depth and temperature to ensure the variations of @xmath33 and @xmath17 to be smaller than 5% at one step .
following previous studies ( e.g. , gpe , hernquist 1985 , van riper 1988 ) , we terminate the integration at @xmath34 g @xmath7 .
a comparison with analytically solvable models ( e.g. , hernquist & applegate 1984 ) shows that our numerical procedure determines @xmath5 at a given @xmath28 with an error @xmath35 1% .
although the density @xmath36 does not enter eq.([th - str ] ) explicitly , its value is required to determine the opacity @xmath17 at given @xmath37 .
thus , we need an eos of matter in the outer ns envelope . in this paper , we are specifically interested in h , he , c , o , and fe plasmas , which are likely to compose partly accreted ns envelopes . the parameters of interest cover the region which extends from the partially ionized non - degenerate non - relativistic plasma near the surface , through the domain of pressure ionization , to high - density layers of fully ionized relativistic plasma .
we assume a strongly stratified envelope , i.e. , that at each given density , the plasma is composed of electrons and of ions of one chemical element which can be in different ionization stages .
the plasma density can be characterized by the parameter @xmath38 where @xmath39 is the mean inter - electron distance , @xmath40 is the number density of electrons , and @xmath41 is the bohr radius . the parameter @xmath42 is directly related to the relativistic parameter of degenerate electrons @xmath43 ( e.g. , yakovlev & shalybkov 1989 ) : @xmath44 . here
, @xmath45 is the electron fermi momentum , @xmath46 is the electron mass , @xmath47 g @xmath7 , and @xmath48 and @xmath49 are the mean charge and mass number of ions , respectively . in our case , the averaging is over different ionization stages of ions of the same chemical element .
the electron degeneracy temperature is given by @xmath50 m_{\rm e}c^2     / k_{\rm b},\ ] ] where @xmath51 is the boltzmann constant .
let us also introduce the ion coupling parameter , @xmath52 , which measures typical energy of ion coulomb interaction relative to the thermal energy ( @xmath53 being the electron charge ) .
let us specify the high density regime as @xmath54 .
this corresponds to @xmath55 g @xmath7 , @xmath56 and @xmath57 being the charge and mass numbers of atomic nuclei . in this regime , @xmath58 is considerably smaller than the @xmath17-shell radius of an atom , and we have a fully ionized one - component plasma ( ocp ) of ions immersed in the `` rigid '' electron background .
then the pressure is @xmath59 where @xmath60 is the electron pressure , @xmath61 is the pressure of ideal gas of ions , and @xmath62 is the coulomb term .
the pressure of the ideal gas of electrons of any degeneracy is given by ( e.g. , landau & lifshitz 1986 ) , @xmath63     p^2{\rm d}p , \label{p_e}\ ] ] where @xmath64 is the energy of an electron with a momentum @xmath65 .
the electron chemical potential @xmath66 has to be determined from the equation @xmath67}. \label{ne}\ ] ] the pressure contribution from the ideal gas of ions in eq.([phigh ] ) is @xmath68 where @xmath69 is the ion number density .
the coulomb correction accounts for electrostatic interactions of plasma particles .
according to hansen & vieillefosse ( 1975 ) , the results of monte carlo simulations of the ocp of ions immersed in a rigid electron background can be fitted by @xmath70 , \label{p_c}\ ] ] where @xmath71 , @xmath72 , @xmath73 , and @xmath74 .
the fit error is smaller than 1% at @xmath75 , as confirmed by comparison with the recent monte carlo results of stringfellow et al .
( 1990 ) and dewitt et al .
equation ( [ p_c ] ) reproduces also the correct debye
hckel limit at @xmath76 .
equations ( [ phigh])([p_c ] ) can not be used at @xmath77 and @xmath78 because of the failure of the rigid electron background approximation .
they lead even to negative pressure at low temperatures and intermediate densities .
van riper ( 1988 ) modified the coulomb correction ( [ p_c ] ) in an _ ad hoc _ manner to avoid negative pressures .
we use another approach described below in sect .
[ sect - interdens ] .      in the _ low density regime _ , @xmath79 and @xmath76 ,
the plasma is nearly ideal and may be partially ionized , depending on @xmath33 and @xmath36
. theoretical description of partially ionized plasmas can be based either on the physical picture or on the chemical picture of the plasma ( e.g. , ebeling et al .
1977 ) . in the chemical picture
, the bound species ( atoms , molecules , ions ) are treated as elementary members of the thermodynamic ensemble , along with free electrons and nuclei . in the physical picture , nuclei and electrons ( free and bound ) are the only constituents of the ensemble . both pictures can be thermodynamically self - consistent , though the chemical picture has limited microscopic consistency ( see , e.g. , chabrier & schatzman 1994 , for review ) .
the most advanced results based on the physical picture are those obtained in the opal project ( rogers et al .
however , as an intrinsic limitation of the formalism , they are restricted to low densities and/or high temperatures .
we employ the opal eos for partially ionized iron , which composes non - accreted envelopes .    for the outermost layers of accreted envelopes , composed of h and he
, we use the eos derived by saumon et al .
( 1995 ) ( scvh ) within the framework of the chemical picture .
the second - order thermodynamic quantities provided by scvh suffer from some thermodynamic inconsistency in some regions of the phase diagram ( see scvh and potekhin 1996 , for a discussion ) , but the inconsistency is quite small at low densities . a comparison of the scvh and opal eoss for h and he reveals no discrepancies which could noticeably affect the temperature profiles in a ns envelope under the conditions of interest .
furthermore , the scvh tables cover a larger @xmath80 domain and provide also , along with the first- and second - order thermodynamic quantities , the fractions of h and he atoms , h@xmath81 molecules , and h@xmath82 , he@xmath83 , he@xmath84 ions , which enables us to determine an average ion charge , required for our calculations of the conductivities ( sect .
[ sect - conduct ] ) .
the intermediate density domain ( @xmath85 ) is most complicated because of the onset of recombination .
previous studies of the thermal structure of the ns envelopes ( gpe , van riper 1988 ) employed the model described by eqs.([phigh])([p_c ] ) , with the free electron density @xmath40 taken from the los alamos astrophysical opacity library ( huebner et al .
1977 , hereafter lao ) .
however , this approach fails at intermediate densities and relatively low temperatures ( see , e.g. , van riper 1988 ) , where it leads to unrealistic eos owing to too large values of the coulomb correction . in the afore - mentioned density region ,
the less dense and more strongly correlated electron gas is polarized by the external ionic field and eventually the electrons become bound by pressure recombination .
equation ( [ p_c ] ) becomes inadequate in this case , since it applies only to a _ fully ionized _ plasma , immersed in a _ rigid _ electron background .
fortunately , at temperatures of present interest ( @xmath86 } \ga 4.7 $ ] ) , the scvh tables extend up to @xmath87 ( @xmath88 g @xmath7 for he ) , i.e. , to the completely ionized region .
thus , they fully cover the intermediate density range for h and he .
we use these tables up to the highest tabulated densities , and we use the high - density eos given by eqs.([phigh])([p_c ] ) for still higher @xmath36 .    the high - density domain for c and o starts from @xmath89 g @xmath7 .
we do not need to consider c and o at lower densities in the present work .
the case of iron is more difficult .
the opal tables are bound from low-@xmath33 high-@xmath36 side approximately by the line @xmath90 g @xmath7 . in order to extend the eos beyond this boundary , we adopt the approach of fontaine et al .
( 1977 ) , who circumvented a similar difficulty by interpolating pressure logarithm along isotherms over the intermediate density region . however , the conventional cubic spline interpolation of @xmath91 vs  @xmath92 violates the condition @xmath93 although this `` normality condition '' is not required by the thermodynamic consistency , it is supposed to hold at intermediate densities ( see , e.g. , fontaine et al .  1977 ) .    in order to keep this condition , we introduce the following modification of the interpolation procedure .
( i ) pressure isotherms from @xmath94 to 8 are taken from the opal data at @xmath95 and calculated from eqs.([phigh])([p_c ] ) at @xmath96 .
( ii ) the differences @xmath97 between @xmath91 for neighbouring isotherms are interpolated across the gaps in @xmath92 between @xmath98 and @xmath99 by cubic splines .
( iii ) starting with the isotherm @xmath100 , for which eqs .
( [ phigh])([p_c ] ) are sufficiently accurate , the pressure is consecutively calculated between @xmath98 and @xmath99 along the isotherms @xmath101 , using the differences @xmath97 obtained at the preceding step : @xmath102 , where @xmath103 is the pressure at the preceding ( hotter ) isotherm .    in order to verify this procedure ,
we have applied it to he and found that it gives much better agreement with the scvh eos than the straightforward @xmath91@xmath92 interpolation across the intermediate densities .
we will treat the electron heat conduction in the mean ion approximation , that is , we consider the plasma as composed of electrons and one ionic species with an effective charge @xmath104 . in the case of h and he , we adopt the mean ion charge @xmath105 provided by the scvh tables , as discussed above . for carbon and oxygen , only the high - density ,
fully ionized regime is involved , for which @xmath106 and 8 , respectively .
finally , for iron we adjust the effective number density of free electrons @xmath40 for the pressure @xmath107 , calculated formally from eqs .
( [ p_e])([p_i ] ) , to reproduce the pressure @xmath11 given by the accurate eos at the same @xmath36 and @xmath33 .
the mean ion charge , that corresponds to @xmath40 , is then assumed to be equal to @xmath104 .
figure  [ fig - zeff ] shows the isotherms of @xmath104 corresponding to the interpolated fe eos .
the decreasing parts of the curves reflect the electron recombination due to the reduction of phase space per particle with increasing density . in this low - density region ,
our @xmath104 coincides with @xmath48 given by the opal data ( when available ) within @xmath108% .
the increase of @xmath104 at high densities corresponds to pressure ionization .
the non - monotonic behaviour at intermediate densities reflects consecutive pressure destruction of different electron shells of fe atoms .
in recent years , a significant progress has been made in calculations of the radiative opacities in dense atmospheres .
we use the most recent opal opacity library ( rogers et al .
1996 ) . for iron and hydrogen
, it displays the tables of the spectral opacities , which we convert into the rosseland opacities @xmath19 using the standard frequency - averaging procedure ( e.g. , schwarzschild 1958 ) . for helium , we implement the rosseland opacities readily available in the opal library .    for @xmath109 values within the table range ( @xmath95 )
, we have used bilinear interpolation of @xmath110 between the table entries in @xmath92 and @xmath111 . for @xmath112 ,
we employ similar bilinear extrapolation based on the nearest pairs of table entries in @xmath92 and @xmath111 .
the extrapolated values of @xmath19 have practically no effect on temperature profiles , because at high densities the heat transport is provided by electrons ( @xmath113 ) .
previous studies of the temperature profiles in ns envelopes have mainly used the electron thermal conductivity of completely ionized , strongly degenerate plasma . in particular
, one has often taken the conductive opacities from urpin & yakovlev ( 1980 ) and yakovlev & urpin ( 1980 ) ( hereafter yu ) . however , in a cold enough ns , the thermal conductivity of non - degenerate or partly degenerate electrons may become important .
a comparison with the tabular data of hubbard & lampe ( 1969 ) reveals that a straightforward extrapolation of the yu formulae from their validity domain ( fully ionized , degenerate plasma ) to the case of non - degenerate matter may lead to an underestimation of the electron thermal conductivity by orders of magnitude . as will be shown in sect .
[ sect - res ] , this may strongly affect the temperature profiles in the ns envelopes at @xmath114  k.    in order to determine the conductive opacities for any degeneracy , we employ the numerical code developed recently by potekhin & yakovlev ( 1996 ) ( hereafter py ) for calculating the electron transport coefficients along magnetic fields in the magnetized ns envelopes .
the code performs numerical thermal averaging of the effective energy - dependent electron relaxation time with the fermi
dirac distribution at any degeneracy .
this enables us to consider not only strongly degenerate , but also mildly degenerate matter . in addition
, py have improved the approach of yu by a more accurate calculation of the coulomb logarithm for the electron - ion scattering in the liquid phase of ns matter and by incorporating the debye
waller reduction factor ( e.g. , itoh et al .  1984 ) for the electron - phonon scattering in the solid phase .    in the present article
, we apply the above code for the particular case of zero magnetic field . for this purpose
, further modifications have been made .
first , in addition to the thermal conduction produced by the electron - ion scattering considered by py , we have incorporated also the contribution from the electron - electron scattering in the liquid regime ( see sect .
[ app - ee ] ) .
this process can be important for light elements ( h , he ) in a weakly degenerate matter .
second , we have improved the coulomb logarithm which determines the collision rate of strongly degenerate electrons with ions ( see sect .  [ sect - coulog ] ) .
third , in the solid regime , we have used the new expression for the electron thermal conductivity , derived recently by baiko & yakovlev ( 1995 ) . note that at temperatures much lower than the melting temperatures , the coulomb scattering of electrons by charged impurities may become important ( e.g. , yu ) .
it depends on the impurity charge and , most important , on its number density @xmath115 which is generally unknown in the ns envelopes .
we will present the results obtained for pure crystals ( @xmath116 ) .
our calculations including @xmath117ni impurities in iron envelopes show that the electron - impurity scattering may affect the temperature profiles in the crusts of cooling nss , whenever the impurity concentration is sufficiently high , @xmath118 .
finally , in the non - degenerate layers , we have performed the thermal averaging introducing modifications into the effective width of the landau levels which was used by py to account for collisional and inelastic broadenings of the landau levels in strongly degenerate matter .
the effective width proposed by py becomes inadequate in the non - degenerate regime , and we have switched it off by multiplying by the factor @xmath119 .
a comparison of our thermal conductivities with the tables of hubbard & lampe ( 1969 ) ( at @xmath120 g @xmath7 , where the tables are adequate ) shows maximum discrepancy up to 30% for carbon and up to 70% for lighter elements in the non - degenerate regime , and still smaller discrepancies in mildly and strongly degenerate matter .
more accurate values of the coulomb logarithm obtained in sect .  [ sect - coulog ]
modify @xmath20 by about 530% , which has only a little effect on the thermal structure of the envelope ( see sect .
[ sect - res ] ) .
the thermal conductivity @xmath21 of degenerate electrons ( @xmath121 ) in a liquid or a gas of ions with mass and charge numbers @xmath57 and @xmath56 can be written as @xmath122 where @xmath123 , @xmath124 is the effective electron collision frequency , @xmath125 is determined by the coulomb electron - ion collisions , and @xmath126 is determined by the electron - electron collisions .
new calculations and fits of @xmath126 are presented in sect .
[ app - ee ] . in this section ,
we examine @xmath125 .
it can be expressed as ( e.g. , yu ) @xmath127 where @xmath128 is the coulomb logarithm to be evaluated .
it is a slowly varying function of @xmath36 and @xmath33 .
first consider a strongly coupled coulomb liquid of ions , at @xmath129 , where @xmath130 is the critical value of the ion coupling parameter @xmath131 at which a coulomb crystal melts ( @xmath130 = 172 , for a classical ion system , nagara et al .
1987 ; the deviations from the classical melting curve due to the quantum corrections turn out to be small at @xmath132 g @xmath7 for carbon and heavier ions  see chabrier 1993 ) . in this case
, one usually assumes that the screening of the electron - ion coulomb interaction is static and produced by the ion - ion correlations .
the screening is described by the ion - ion structure factor @xmath133 .
extensive calculations of @xmath128 in the born approximation for different ion species , @xmath56 and @xmath57 , were performed by itoh et al .
the authors used the structure factors obtained in the approximation of rigid electron background .
they fitted their results by complicated expressions . later yakovlev ( 1987 ) evaluated non - born corrections to @xmath128 , using the weak screening approximation .
now we recalculate the coulomb logarithm for strongly degenerate electrons at @xmath129 from the expression @xmath134 where @xmath135 , @xmath136 , @xmath137 is a momentum transfer in an @xmath138 collision event , @xmath139 is the static longitudinal dielectric function of electrons , @xmath140 is the nuclear form factor to allow for finite nuclear size , and @xmath141 is the non - born correction factor specified below . the dielectric function has been calculated for the electron gas at any degeneracy .
however , we have verified that , for all the cases in the present study , this dielectric function yields the same results as the function obtained in the limit of strong electron degeneracy for the relativistic electron gas ( jancovici 1962 ) .
thus the bulk of the calculations has been done with this latter function .
we use the nuclear form factor @xmath140 appropriate to a uniform proton core of an atomic nucleus , with the same core radii as in itoh et al .
the finite sizes of the nuclei appear to have almost no effect on @xmath128 , for the conditions of study .
nevertheless we have kept @xmath140 in our calculations . following yakovlev ( 1987 ) , the non - born factor has been taken as @xmath142 , where @xmath143 is the exact differential coulomb scattering cross section for a momentum transfer @xmath144 , and @xmath145 is the cross section in the born approximation .
the exact cross sections are taken from doggett & spencer ( 1956 ) .    in our calculations from eq .
( [ l ] ) , we have used the ion structure factors evaluated for a _ responsive _ electron background including the local field corrections between electrons ( chabrier 1990 ) .
our calculations extend those of itoh et al .
( 1983 ) by including the non - born corrections and the effects of electron response on the ion structure factors .
on the other hand , we improve the results of yakovlev ( 1987 ) since we evaluate the non - born corrections beyond the weak screening approximation .
we have evaluated @xmath128 from eq .
( [ l ] ) for chemical elements with @xmath146 and for a dense grid of @xmath147 and @xmath148 g @xmath149 g @xmath7 provided @xmath150 . in the case of hydrogen , we have restricted ourselves to @xmath151 g @xmath7 since hydrogen should be completely burnt at higher densities . for light elements ,
the highest temperatures included in these data correspond to @xmath152 and appear to be much below @xmath153 .
accordingly , there exists a large temperature interval in the electron degeneracy domain where @xmath154 , the ion coupling ceases to be strong , the ion screening can not be treated in terms of ion - ion correlations , and eq .
( [ l ] ) is invalid . in principle , the coulomb logarithm for mildly - coupled ion plasma ( @xmath155 ) can be evaluated using the formalism developed by boerker et al .
( 1982 ) . in order to simplify our analysis
, we will restrict ourselves to the weak  coupling case ( @xmath76 ) , in which the ions constitute a boltzmann gas , the weak ion screening is dynamical and can be described by the dielectric function formalism , while the electron screening remains static ( e.g. , lampe 1968 ) . at @xmath76 and @xmath121 ,
we propose the expression @xmath156 where @xmath157 is the ion screening parameter in the weak coupling regime , @xmath158 is the debye radius of ions , @xmath159 , @xmath160 is the electron screening parameter , @xmath161 being the inverse screening length of a charge by degenerate plasma electrons and @xmath162 the fine - structure constant .
the first three terms represent the coulomb logarithm obtained for degenerate electrons and weakly coupled ions with the non - relativistic @xmath138 scattering cross section .
the expression presented is valid for @xmath163 ; it can be derived using the formalism of williams & dewitt ( 1969 ) to account for the dynamic ion screening and to integrate the @xmath138 collision rate over velocities of ions .
the fourth term is the relativistic correction in the weak  coupling limit ( e.g. , yu ) .
the fifth term is the second  born correction in the weak  coupling limit ( yakovlev 1987 ) .
thus , we have supplemented our table of the coulomb logarithms calculated for @xmath164 ( see above ) with new values evaluated from eq .
( [ dewitt ] ) for @xmath165 and @xmath121 .
analytic fits to the full set of about 31,000 values of @xmath128 are presented in sect .
[ app - l ] .
the tables of @xmath128 are freely available in the electronic form .    for illustration , in fig .
[ fig - l ] ( left panel ) we compare the coulomb logarithms in fully ionized fe and c plasmas with those calculated in the born approximation ( itoh et al .  1983 ) .
the non - born corrections enhance the small  angle coulomb scattering cross section , and increase @xmath128 .
the increase is especially pronounced for the heavier element , fe ( exceeds 20% ) , at densities @xmath166 g @xmath7 , where the electrons are relativistic .
the non - born corrections become less important with decreasing @xmath56 , and they are negligible for light elements ( h , he ) .
the inclusion of the electron response in the calculations of the ionic structure factors affects the coulomb logarithm @xmath128 for h and he at @xmath167 , as demonstrated on the right panel of fig .
[ fig - l ] .
one can observe an increase of @xmath128 up to 15% for h at @xmath168 g @xmath169 and @xmath170 with respect to the case of rigid electron background .
actually , the increase grows with @xmath131 and reaches about 40% at @xmath171 .
the effect is weaker for he , and small for heavier elements ( in the adopted parameter range ) .
thus , for the elements heavier than he , we have used the structure factors calculated in the rigid background approximation .
note that calculated values of @xmath128 determine also the electric conductivity of degenerate matter in a liquid or gaseous ionic plasma , @xmath172 .
in this section , we consider a non - accreted ns envelope composed solely of iron , which may be at various ionization stages .
we use the opal iron eos where available , the ideal - gas eos with the coulomb correction ( [ p_c ] ) at high densities , and the interpolated eos at intermediate densities as described in sect .
[ sect - interdens ] .
the conduction properties are described in the mean ion approximation with the effective ion charge @xmath104 determined consistently with the eos ( sect .
[ sect - zeff ] ) .
figure [ fig - prfe ] presents density dependence of the temperature in the ns envelope at various @xmath0 .
the integration of eq .
( [ th - str ] ) has been started with the surface density @xmath173 , shown in fig .
[ fig - rho - s ] , and terminated at @xmath34 g @xmath7 .
the value @xmath174 chosen in fig .
[ fig - prfe ] corresponds to a ns with a radius @xmath175 km and a mass @xmath176 .
circles on the curves are the points of equal radiative and conductive opacities .
the thermal flux is mainly carried by radiation at lower densities and by electrons at higher densities .
we also show the electron degeneracy curve and the melting curve .
in a wide range of @xmath0 , some layers ( fig .  [ fig - prfe ] ) turn out to be convectively unstable ( zavlin et al .
1996 , rajagopal & romani 1996 ) . in these layers ,
the energy transport is dominated by convection rather than by electron or radiative conduction .
we describe the convective energy flux in the adiabatic approximation ( schwarzschild 1958 ) , which assumes that the convective heat transport is much more efficient than other mechanisms .
accordingly , in the convective zones , we replace the temperature gradient given by eq .
( [ th - str ] ) by the adiabatic gradient provided by the eos tables .    in order to check the effect of convection , we have performed also calculations with `` frozen '' convection ( i.e. , neglecting the convection effect ) .
this is the extreme case opposite to the adiabatic one . in this approximation ,
we obtain ( fig .
[ fig - prfe ] ) slightly higher temperatures inside the convective part of the atmosphere .
the atmospheric temperature profiles derived by zavlin et al .
( 1996 ) , who solved numerically the radiative transfer equation at moderate optical depths and described the convection using the mixing - length theory , lie between our two extremes . in deeper layers ,
the two extreme profiles tend to merge , as clearly seen in fig .
[ fig - prfe ] .
this merging occurs because the factor @xmath177 on the r.h.s .  of eq .
( [ th - str ] ) rapidly decreases , with increasing @xmath33 , and thus ensures the isothermal profile at higher @xmath33 behind the convective zone ( at fixed @xmath0 ) .
thus , the thermal structure of the ns envelope at @xmath178 g @xmath7 is practically unaffected by the convection .    the dot - dashed curves in fig .  [ fig - prfe ]
are obtained with @xmath20 calculated according to py ( sect .  [ sect - conduct ] ) .
these curves almost coincide with the solid lines , calculated more accurately .
this confirms the validity of the py approximations for the non - accreted envelopes .
if @xmath179  k , the electron conduction becomes important not only in degenerate matter , but also in a non - degenerate region ( fig .
[ fig - prfe ] ) .
then the straight - line logarithmic extrapolation of the yu conductive opacities into the non - degenerate regions leads to a significant overestimation of the internal temperature ( dotted curves ) .
this effect is seen also in fig .
[ fig - link - fe ] , that shows @xmath180 as a function of @xmath181 for four values of the surface gravity .
if @xmath182>5.5 $ ] , the well - known simple fit of gpe , @xmath183 is rather accurate , but at lower @xmath28 the deviations from this fit are quite pronounced .
nevertheless , the scaling law @xmath184 holds with a good accuracy in the whole temperature  gravity range presented in this figure , except for the lowest @xmath5 values .
the fit which reproduces these new features is presented in sect .
[ app - link ] .
the chemical composition of an accreted envelope depends on the composition of the infalling matter and on the nuclear transformations ( thermo- or pycnonuclear burning , beta - captures ) which occur while the matter sinks within the ns .
the structure and stability of the accreted envelopes have been considered in number of works devoted to bursting nss ( see , e.g. , ayasli & joss 1982 , paczyski 1983 , ergma 1986 , miralda - escud et al .
1990 , blaes et al .
1992 , and references therein ) . however , most of the work focused either on high accretion rates or on slowly accreting old nss whose internal temperature is solely determined by the accretion . on the contrary , we are mainly interested in not too old , cooling nss with thin accreted envelopes , in which the heat release due to the nuclear transformations has no effect on the thermal structure .
let us assume that an accreted envelope consists of shells with different chemical compositions .
for an accretion of h / he , it is reasonable to expect that the outermost accreted layers are built of pure h , owing to the strong gravitational stratification ( alcock & illarionov 1980 ) , whereas he may exist in the deeper layers .
the accretion is accompanied by the nuclear transformations of h / he into heavier elements .
the details of these transformations are still uncertain , and we consider several models of accreted envelopes .    as a first example
, we take the typical temperatures and densities of the nuclear burning from ergma ( 1986 ) and assume the burning to be non - explosive . if @xmath0 is high and the accreted hydrogen layer reaches the region with @xmath185  k , then h burns efficiently into he in the thermonuclear regime . if temperature is lower and the hydrogen layer reaches the density @xmath186 g @xmath7 , a pycnonuclear burning of h into he proceeds . at higher @xmath187 g @xmath7 and/or @xmath188  k , helium burns into carbon ( schramm et al .
the pycnonuclear carbon burning occurs at @xmath189 g @xmath7 ( e.g. , yakovlev 1994 ) , which is already inside the isothermal region of the envelope .
figure [ fig - prc ] displays the thermal structure of a fully accreted envelope ( the accreted mass @xmath190 fills the layers up to @xmath191 ) of a ns with the surface gravity @xmath192 cm s@xmath193 .
the outer , intermediate , and inner shells of this envelope , separated by asterisks , are composed of h , he , and c , respectively .
one can observe significant differences from the thermal structure of a non - accreted fe envelope , as explained below .
for a not too cold ns ( @xmath194  k ) , the main temperature gradient occurs in a layer of degenerate electron gas with ions in the liquid state , where the thermal conduction is mostly provided by the electrons due to the electron - ion collisions .
the heavier the element , the smaller the thermal conductivity , and the steeper is the temperature growth inside the star . however , with decreasing @xmath0 , the width of the heat  insulating degenerate layer becomes smaller , and the effect is less pronounced . in a cooler ns ( @xmath195  k ) , the main temperature gradient shifts into the ns atmosphere ( to the optical depths @xmath196 ) . for heavier elements ,
the atmospheric layers appear generally denser ( at the same @xmath0 , see fig .
[ fig - rho - s ] ) . then the internal temperature gradient is weaker and the temperature grows slower inside the star .
thus , a not too cold accreted envelope is more heat  transparent than the iron envelope , whereas a cold accreted envelope is less heat - transparent . as seen from fig .
[ fig - rho - s ] , the effective surface temperature which separates these two regimes is almost independent of the surface gravity .
figure [ fig - py ] demonstrates the effect of the improvements in the conductive opacities .
the present results are compared with those in which the electron thermal conductivity is calculated with the py code .
although the difference is not too large , it is higher than for the non - accreted fe envelopes ( fig .
[ fig - prfe ] ) .
the difference comes mostly from the electron - electron collisions ( neglected in py ) .
dots in fig .
[ fig - py ] show the results obtained with the py code supplemented by the contribution from the electron - electron collisions to the thermal conduction ( sect .
[ app - ee ] ) .
the difference from the accurate results practically disappears .
now let us analyze the sensitivity of the temperature profiles to the positions of the h / he and he / c interfaces .
these interfaces , indicated by asterisks in fig .  [ fig - prc ] , are not well theoretically established .
for example , the densities of nuclear h and he burning presented by iben ( 1974 ) ( see also kippenhahn & weigert 1990 ) are much lower than those given by ergma ( 1986 ) .
we have significantly shifted the interfaces according to the results of iben ( 1974 ) , but the @xmath0@xmath5 relationship remains practically unchanged .
we have also calculated the temperature profiles for a purely helium ns envelope .
these profiles are slightly different from those in a burning accreted matter described above .
however , the @xmath0@xmath5 relationships remain practically the same .
recently the thermal structure of old , slowly accreting nss has been examined in several articles ( see miralda - escud et al .
1990 , blaes et al .
1992 , and references therein ) . in these studies ,
the internal stellar temperature is assumed to be determined by the accretion rate .
the composition of the accreted envelopes appears to be similar to the ones considered above except for the innermost layers . for instance , miralda - escud et al .
( 1990 ) assume that he transforms directly into fe while blaes et al .
( 1992 ) state that the burning of he into c at @xmath197 g @xmath7 is followed by the carbon transformation through the rapid @xmath198c@xmath199o process .
accordingly , we have recalculated the temperature profiles replacing @xmath198c by @xmath117fe or @xmath200o and found only very small differences .
the insensitivity of the temperature profiles to the composition of the innermost layers of the envelope is quite natural , because these layers are nearly isothermal .    summarizing the above discussion
we can conclude that the @xmath0@xmath5 relationship is remarkably insensitive to the parameters of nuclear burning and to the chemical composition of accreted matter ( h and/or he ) but differs from the one in a non - accreted ns .
this conclusion will not change if h burning is explosive ( see , e.g. , paczyski 1983 , miralda - escud et al .
1990 ) . during an outburst , the h - shell burns into he , without any noticeable effect on the @xmath0@xmath5 relation .
according to miralda - escud et al .
( 1990 ) helium burning in slowly accreting nss is unstable ; it produces nuclear outbursts with burning of all accreted material .
nevertheless , the helium layer should reach a very high density , @xmath187 g @xmath7 to trigger the instability in a not too hot ns . if the same is true for the conditions of our interest , our results are valid during a long period of time until the instability is reached . however , according to paczyski ( 1983 ) , the presence of the heat outflow from a not too cold , cooling ns essentially stabilizes nuclear burning in accreted envelope , and the he burning is mostly non - explosive for the parameters we are interested in .
figure  [ fig - link - a ] shows @xmath180 as a function of @xmath181 for various amounts of accreted matter . the dot - dashed line represents @xmath180 for a non - accreted envelope from fig .
[ fig - link - fe ] .
this function is strongly affected by the accretion .
even a thin hydrogen ( or he ) mantle of mass @xmath201 , which extends only to @xmath202 g @xmath7 , modifies strongly the @xmath5@xmath0 relation .
analytical fits of @xmath0 as a function of the internal temperature @xmath5 , surface gravity @xmath203 , and accreted mass @xmath204 are given in sect .
[ app - link ] .      consider briefly the effects of the possible presence of accreted matter in the ns envelopes on the thermal history of nss . for this purpose
, we simulate the cooling of nss using two ns models .
we use the cooling code of gnedin & yakovlev ( 1993 ) based on the approximation of isothermal ns interior .
the code follows accurately the cooling of a not too young ns whose thermal relaxation is already over ( @xmath205 yrs , see nomoto & tsuruta 1987 , umeda et al .
we adopt a moderately stiff eos of prakash et al .
( 1988 ) to describe matter in the ns cores .
this matter is assumed to consist of neutrons , protons , and electrons ( no exotic cooling agents such as quarks or kaon condensates ) .
the maximum ns mass , for this eos , is about 1.7 @xmath206 .
we examine two ns models : ( 1 ) @xmath207 , @xmath208 km , the central density @xmath209 g @xmath7 ; ( 2 ) @xmath210 , @xmath211 km , @xmath212 g @xmath7 .
the first model is a typical example of the _ standard _ cooling .
the central density is insufficient to switch on the most powerful neutrino cooling due to the direct urca process ( lattimer et al .
1991 ) . in this case
the main neutrino energy loss mechanisms in the ns cores are the modified urca reactions ( neutron and proton branches ) and the nucleon - nucleon bremsstrahlung ( friman & maxwell 1979 , yakovlev & levenfish 1995 ) .
our second ns model is an example of the _ rapid _ neutrino cooling . in this model , the central density slightly exceeds the threshold density @xmath213 g @xmath7 of the direct urca process ( for the given eos )
. then the ns has a small central kernel ( of radius 2.32 km and mass 0.035 @xmath206 ) where the direct urca process operates , and the neutrino luminosity exceeds the standard one by several orders of magnitude .    all the neutrino energy loss rates included in the calculations are described by levenfish & yakovlev ( 1996 ) . to simplify our analysis
, we assume that the ns cores are non - superfluid , and there is no internal reheating ( e.g. , energy dissipation due to differential rotation ) .
the results are presented in fig .
[ fig - cooling ] in the form of the cooling curves  effective temperatures @xmath214 as measured by a distant observer ( see sect .  [ sect - th - str - eq ] ) versus ns age @xmath215
.    figure [ fig - cooling ] presents the cooling curves under the assumption that all the accreted material is accumulated at the ns surface during the ns birth .
the accretion can occur , for instance , at the post - supernova stage ( e.g. , chevalier 1989 , 1996 , brown & weingartner 1994 ) especially if the pulsar activity is suppressed initially by burying the pulsar magnetic field under the fallen - back matter ( muslimov & page 1995 ) .
we show the standard and rapid neutrino cooling for different amount of accreted mass @xmath204 .
the fraction of accreted mass @xmath216 varies from 0 ( non - accreted nss ) to @xmath9 ( fully accreted ns envelope ) .
it is evident that further increase of @xmath204 does not affect the cooling .
the change of slopes of the cooling curves at @xmath217@xmath218 yrs reflects the change of the cooling regime .
initially , a ns is at the neutrino cooling stage .
it cools mainly via the neutrino emission ; the internal stellar temperature is ruled by this emission and is thus independent of the thermal insulation of the ns envelope .
the surface photon emission is determined by the @xmath5@xmath0 relation . since the accreted envelopes of not too cold nss are more thermally transparent ( sect .
[ sect - accr ] ) , the surface temperature of an accreted ns is noticeably higher than that of a non - accreted ns .
one can see that even a very small fraction of accreted matter , such as @xmath219 , may change appreciably the thermal history of the star .
the colder the star , the smaller is the fraction of accreted material which yields the same cooling curve as the fully accreted ns envelope . the effect is naturally more pronounced for the rapid cooling ( for colder nss ) .    if @xmath220@xmath218 yrs , the neutrino luminosity decreases , and the ns cools mainly via the surface emission of photons ( the photon cooling stage ) .
since the envelopes of the accreting ( and not too cold ) nss have lower thermal insulation , these stars cool now faster than the non - accreted ones ( fig .
[ fig - cooling ] ) .
the change of the thermal history of a ns by the accreted material should be taken into account in theoretical interpretation of observational data on thermal radiation from nss . for illustration , in fig .
[ fig - cooling ] we present the effective temperatures with @xmath221 confidence intervals of the pulsars 083345 vela ( gelman et al .
1993 , gelman 1995 ) , 0630 + 18 geminga ( halpern & ruderman 1993 ) , 0656 + 14 and 105552 ( greiveldinger et al .
1996 ) , rxj 0002 + 6246 ( hailey & craig 1995 ) , and the @xmath222 upper limit for the crab pulsar ( psr 0531 + 21 , becker & aschenbach 1995 ) .
these data have been deduced from the fits of spectral observations of the _ rosat _ x - ray observatory . for geminga ,
psr 105552 , and psr 0656 + 14 , we plot the values of @xmath0 obtained for `` soft '' components of the two- and three - component x - ray spectra models , since it is the soft component that is associated with the thermal flux from the stellar interior . for vela , we extend downward the confidence interval obtained by gelman et al .
( 1993 ) according to the more recent spectral fitting of gelman ( 1995 ) , that gives a more plausible ns radius than the previous fit .    the age of the crab pulsar in fig .
[ fig - cooling ] is its true age known from the chronicles , the age of vela is its possible true age determined recently by lyne et al .
( 1996 ) , the age of rxj 0002 + 6246 is an estimate for the related supernova remnant g117.7 + 06 , while for the three other pulsars we adopt their characteristic ages .
we plot the values of @xmath0 derived by the authors by fitting the assumed blackbody thermal spectrum to the measured x - ray spectrum .
the alternative often used is to derive @xmath0 from the apparent luminosity in the whole spectral range observed .
however , this latter way is less certain , because the inferred temperature strongly depends on the assumed ns distance ( usually poorly known ) and model - dependent radius .
note that the blackbody fits do not take into account the influence of the magnetized atmosphere .
the fits for psr 0656 + 14 ( anderson et al .
1993 ) , geminga ( meyer et al .  1994 ) , and vela ( page et al .
1996 ) , which employ the models of magnetized atmospheres , have resulted in a few times lower @xmath0 than the blackbody models .
our justification of using here the blackbody data is that the so far developed models of magnetized atmospheres do not take proper account of absorption by neutral atoms
. the neutral fraction can be significant at @xmath223 k , because the magnetic fields @xmath224  g increase the ground - state binding energy by about one order of magnitude .
the opacities of hydrogen in such strong fields are still under investigation ( see pavlov & potekhin ( 1995 ) and references therein ) .
the magnetic fields affect also the thermal insulation of the ns envelopes and , therefore , the cooling ( e.g. , yakovlev & kaminker 1994 ) .
we do not take this effect into account in the present article .
as suggested by hernquist ( 1985 ) and shown recently by shibanov & yakovlev ( 1996 ) , the dipole surface magnetic fields @xmath225  g affect the cooling much more weakly than predicted by the studies of simplified magnetic field geometries ( e.g. , van riper 1988 , 1991 ) .    as seen from fig .
[ fig - cooling ] , for the pulsars vela , geminga , and rxj 0002 + 6246 , the presence of accreted matter would simplify the explanation of the cooling by the standard neutrino emission without invoking superfluidity of ns matter , as favoured in previous work ( page 1994 , levenfish & yakovlev 1996 ) .
interestingly , the accretion rate of vela estimated by morley ( 1996 ) , multiplied by its age , yields @xmath226 , that just fits the middle of the error bar in fig .
[ fig - cooling ] . for psr 0656 + 14 ,
the scenarios with and without the accretion are equally consistent with the result of the fitting shown if fig .
[ fig - cooling ] , especially if we assume that its braking index is somewhat lower than 3 ( as happens for all 4 pulsars with known braking indices , see lyne et al .
1996 ) and therefore its true age is somewhat higher than its characteristic age . the relatively high temperature of psr 105552 can not be explained by the considered models of the cooling , that probably suggests some reheating process for this relatively old pulsar .
we have used , for illustration , the standard ns models and simplified blackbody fits to the observed spectra .
more sophisticated models and fits can also gain from the assumption of accreted envelopes .
for instance , page ( 1996 ) considers the cooling of nss with kaon condensed cores and reheating in the outer crusts , using our models of the envelopes .
we have considered thermal structure and evolution of nss whose envelopes are composed of non - accreted or accreted matter .
we have used new , state - of - the - art calculations of eos and opacities of ns envelopes , described in sect .
2 . in particular , we have recalculated ( sects .
[ sect - coulog ] , [ app - ee ] , and [ app - l ] ) the electron - electron and electron - ion collision frequencies , that determine the electron thermal conductivity , for a wide range of densities and temperatures of a degenerate electron gas and ionic liquid plasmas .    using this new physics input ,
we have calculated ( sects .
[ sect - non - accr ] , [ sect - accr ] ) the temperature profiles in the envelopes of non - accreted and partly accreted nss and obtained the relationships between the internal and effective temperatures of nss , @xmath5 and @xmath0 .
these relationships are important for simulating the ns cooling ; they are fitted by simple analytic expressions in sect .
[ app - link ] .
they appear to be very sensitive to the presence of accreted matter in the ns envelope ; even a small amount of the accreted matter , @xmath227 , can reduce substantially the thermal insulation of the envelopes . for a non - accreting ns
, our relationship between @xmath5 and @xmath0 extends the well  known result of gpe to lower temperatures ( down to @xmath228 k ) .    in sect .
[ sect - cooling ] , we have examined briefly the effect of the possible presence of accreted matter on the ns cooling .
we show that the accreted matter may increase the surface temperature ( photon thermal luminosity ) at the neutrino cooling stage , and decrease them at the subsequent photon cooling stage , as compared to the nss without accreted envelopes .
we have shown that these results can be important for a proper interpretation of observed thermal radiation from nss .
in particular , the presence of accreted matter facilitates the explanation of recent observational results concerning the pulsars vela and geminga , and rxj 0002 + 6246 , in the framework of the standard neutrino emission model ( without exotic matter , superfluidity , or direct urca processes ) .    in this paper , we have neglected effects of magnetic fields on the eos and the thermal conduction of matter , which can be significant ( e.g. , yakovlev & kaminker 1994 ) .
they do deserve further studies using improved eos and thermal conductivities of magnetized ns envelope ( e.g. , py ) and improved radiative opacities of magnetized ns atmosphere ( e.g. , pavlov & potekhin 1995 ) , with allowance for the possible presence of light elements in the surface layers .
finally , it is worth noting that the physics input used in the present calculations can be applied to a variety of other astrophysical problems concerning dense stellar matter , e.g. the thermal structure and bursting activity of x - ray bursters ( see , e.g. , miralda - escud et al .
1990 and references therein ) and the cooling of white dwarfs ( segretain et al .  1994 ) .
we are grateful to h.  dewitt for pointing out the article of boerker et al .
( 1982 ) , to dany page and yu.a .
shibanov for fruitful discussions , and to d.  page for reading the manuscript and making the most useful critical remarks .
this work was supported in part by the rbrf grant no .
96 - 02 - 16870a , intas grant no .
94 - 3834 , and dfg  rbrf grant no . 96 - 02 - 00177 g .
ayp gratefully acknowledges the hospitality of the theoretical astrophysics group of the ecole normale suprieure de lyon .
the effective collision frequency @xmath126 of non - relativistic degenerate electrons ( @xmath229 , @xmath121 ) was analyzed by lampe ( 1968 ) using the formalism of the dynamic screening of the electron
electron interaction .
the expression of @xmath126 for the relativistic degenerate electrons at @xmath230 was obtained by flowers & itoh ( 1976 ) . here , @xmath231 is the electron plasma temperature determined by the electron plasma frequency @xmath232 , @xmath233 @xmath234 .
the degeneracy temperature @xmath153 in the ns envelopes is typically higher than @xmath231 .
urpin & yakovlev ( 1980 ) extended the results of flowers & itoh ( 1976 ) to higher temperatures , @xmath121 . in the approximation of static electron screening of the coulomb interaction , urpin & yakovlev ( 1980 ) obtained @xmath235 where @xmath236 , @xmath237 , @xmath238 and @xmath239 ( see eq .
( [ dewitt ] ) ) .
now it is sufficient to calculate the function @xmath240 , presented by urpin and yakovlev ( 1980 ) as a 2d integral which depends parametrically on the relativistic parameter @xmath241 .
lampe ( 1968 ) analyzed this function in the static screening approximation at @xmath229 .
the asymptotes of @xmath242 were obtained by lampe ( 1968 ) for @xmath229 at @xmath243 and @xmath244 , by flowers & itoh ( 1976 ) for @xmath244 at any @xmath245 , and by urpin & yakovlev ( 1980 ) for @xmath243 and @xmath246 . timmes ( 1992 ) calculated @xmath240 numerically in the limit of @xmath246 .
however , the unified expression of @xmath240 at @xmath247 valid equally for relativistic and non - relativistic electrons has been absent .
note that the fit expression of @xmath240 obtained by timmes ( 1992 ) ( his eq .
( 10 ) ) is valid only at @xmath248 .
we have calculated @xmath242 numerically for a dense grid of @xmath241 and @xmath249 in the intervals @xmath250 and @xmath251 .
the results are fitted by the expression : @xmath252 \label{jfit}\end{aligned}\ ] ] which reproduces also all the asymptotic limits mentioned above .
the mean error of the fits is 3.7% , and the maximum error of 11% takes place at @xmath253 and @xmath254 .
the tables of the coulomb logarithms calculated from eqs .
( [ l ] ) and ( [ dewitt ] ) as described in sect .
[ sect - coulog ] are fitted by @xmath255 where @xmath256 and @xmath257 for @xmath258 , the fit parameters are given by : + @xmath259 , + @xmath260 , + @xmath261 , + @xmath262 , + @xmath263 , + @xmath264 , + @xmath265 .
if @xmath274 , the rms fit error over all our data set is about @xmath275 3.2% , and the maximum error of @xmath276 7.2% takes place at @xmath277 g @xmath7 and @xmath278 . if @xmath279 we have @xmath275 3.7% , @xmath276 7.8% at @xmath280 g @xmath7 and @xmath281 . for @xmath282 ,
we obtain @xmath283 2.4% , @xmath276 7.0% at @xmath284 ) \approx 2.5 $ ] and @xmath281 . for higher @xmath56 ,
the rms error remains about 1.51.7% , and the maximum error mainly decreases .
for instance , we have @xmath276 5.1% at @xmath285 and @xmath286 for @xmath287 ; @xmath276 3.7% at @xmath288 and @xmath289 for @xmath290 ; and @xmath276 3.4% at @xmath291 and @xmath292 for @xmath293 .
this fit accuracy is quite sufficient for studying the thermal structure of nss .
let us define @xmath298  k , @xmath299  k , and @xmath300 where @xmath13 is the ns mass and @xmath204 is the accreted mass . according to gpe , @xmath301 where @xmath302 is the pressure at the bottom of the accreted envelope .    for a purely iron ( non - accreted ) envelope , a very crude estimate ( with an error @xmath303 30% ) yields @xmath304
then @xmath305 is approximately an increase of the temperature through the iron envelope ( in @xmath306  k ) .
a more accurate fitting formula reads @xmath307 .
\label{t0}\ ] ] the typical fit error of @xmath308 is about 2% , with maximum 4.2% , over the @xmath309 domain indicated above .
finally , for the partially accreted envelopes at any temperatures within the indicated range , we have @xmath312 where @xmath313\,t_{\rm b9}^{5/3}.\ ] ] the typical fit error of eq.([ttfit ] ) for @xmath314 is about 3% , with maximum 5.2% , for all possible values of @xmath315 and any values of @xmath203 and @xmath0 within the indicated range .
the dependence ( [ t0 ] ) is recovered not only at sufficiently low accreted mass ( @xmath316 ) , but also at sufficiently high @xmath5 .
the latter result reflects the fact ( first demonstrated by gpe ) that at high @xmath0 the thermal insulation is mostly produced by the conductive opacities in the deep and hot layers of the envelope , in which light elements ( h , he ) burn into heavier ones . on the other hand , even at very low accreted mass ( @xmath317 ) , the approximation of fully accreted crust is good enough at sufficiently low temperature because in this case the thermal insulation is provided mainly by the low - density accreted surface layers .
page d. , shibanov yu.a .
, zavlin v.e .
1996 , temperature , distance and cooling of the vela pulsar . in :
zimmermann h .- u .
, trmper j. , yorke h. ( eds ) proc .
conf .  on x - ray astronomy and astrophysics , rntgenstrahlung from the universe .
mpe report 263 , garching , p.  173
yakovlev d.g .
, kaminker a.d .
, 1994 , neutron star crusts with magnetic fields . in : chabrier g. , schatzman e. ( eds ) proc .
147 , the equation of state in astrophysics .
cambridge univ .  press , cambridge , p.  214
zavlin v.e .
, pavlov g.g .
, shibanov yu.a . , rogers f.j . ,
iglesias c.a . , 1996 , x - ray spectra from convective photospheres of neutron stars .
in : zimmermann h .- u . , trmper j. , yorke h. ( eds ) proc .
conf .  on x - ray astronomy and astrophysics , rntgenstrahlung from the universe .
mpe report 263 , garching , p.  209