in the last few years a growing amount of high - quality data concerning white dwarfs ( wds ) , both from space as well as from ground based telescopes , have been produced .
in particular , the discovery of wd sequences in globular clusters ( paresce , de marchi , & romaniello 1995 , richer et al .
1997 ) and the observation of their faint ends in old open clusters ( von hippel , gilmore , & jones 1995 , von hippel & gilmore 2000 ) provide new tools for the study of the history of both the halo and the disk of the milky way .
the possible use of the observed wd sequences as age indicators for a variety of galactic components is particularly promising .
as firstly shown by mestel ( 1952 ) , the luminosity of a wd is largely supplied by its thermal energy content , so that the cooling time scale is inversely proportional to a certain power of the luminosity .
thus , comparisons between theoretical estimations of the age - luminosity relation for wds and the observed cooling sequences may be used to derive stellar ages . on the other side
, an accurate comparison between the predictions of the theoretical models with specific observations of wds could provide a better comprehension of the physical properties of high density matter . in this context the recent development of the astereoseismology applied to wds allows to check the internal structure of these compact objects ( bradley & winget 1991 , bradley , winget , & wood 1992 , clemens 1993 , metcalfe , winget , & charbonneau 2001 ) .
white dwarfs are the final destiny of the evolution of low and intermediate mass stars .
the majority of the presently observed wds are post - asymptotic giant branch stars ( agb ) .
the core of an agb star is made of ashes of the he - burning , namely a mixture of primary carbon and oxygen , with a minor amount of secondary @xmath2nen , during the h - burning , and , later on during the he - burning , into @xmath2ne via the chain @xmath3n@xmath4f@xmath5o@xmath6ne .
thus , the abundance ( in number ) of @xmath2ne in the co core of a wd roughly correspond to the abundance ( in number ) of the original cno . ] . during this evolutionary phase the core mass increases , due to the accretion of fresh material processed by the he - burning shell . in the meantime
, a huge mass loss erodes the h - rich envelope until it is reduced down to a critical value ( see e.g. castellani , limongi , & tornamb 1995 ) , and the star departs from the agb . about @xmath7 of the total mass of the resulting wd is the co core builded up during both the central and shell he - burning phases .
m@xmath8 is the minimum stellar mass for which a carbon ignition occurs before the onset of the agb and set the upper mass limit for the progenitors of co wds .
current stellar evolution theory predicts that the value of m@xmath8 varies between 6.5 and 7.5 m@xmath9 , for metallicity ranging between z=0 and 0.02 ( see cassisi , castellani , & tornamb 1996 and dominguez et al .
1999 for recent reviews of these calculations ) .
according to a well established theoretical scenario , single stars having mass lower than m@xmath8 end their life as co wds with masses ranging between about 0.55 and 1.05 m@xmath9 ( see e.g. dominguez et al .
the theoretical initial - final mass relation , as derived from stellar evolution models , roughly agrees with the few available measurements of wd masses in nearby open clusters ( weidemann & koester 1983 , weidemann 1987 , 2000 , herwig 1997 ) .
however , it should be recalled that stellar rotation could induce the formation of more massive co wd ( dominguez et al .
1996 ) .
few wds may have an he - rich core . in principle , stars having initial mass particularly low ( namely m@xmath100.5 m@xmath11 ) never attain temperatures high enough for the he - burning . after the main sequence , they develop a degenerate he core and their final fate should be a low mass he - rich compact remnant .
however , owing to the very long evolutionary time scale of these stars ( much longer than the hubble time ) , we can definitely exclude that the present generation of wds could include their remnants .
nevertheless , low mass he - rich wds may result from the deviated evolution of more massive progenitors when , under particular circumstances , a complete envelope removal takes place during the first ascent of the red giant branch .
for example , a stripping of the envelope may be caused by a close encounter with another star ( livne & tuchman 1988 ) or by roche lobe overflow in close binary systems ( kippenhahn , koll , & weigart 1967 , kippenhahn , thomas , & weigart 1968 , iben & tutukov 1986 ) .
it has been also supposed that massive wds could be generated by stars having an initial mass slightly larger than m@xmath8 ( dominguez , tornamb , & isern 1993 , ritossa , garcia - berro , & iben 1996 , 1999 ; garcia - berro , ritossa , & iben 1997 ) . in this case , in fact ,
an highly degenerate o - ne - mg core is left by the carbon burning .
then , after a rather normal agb phase , eventually characterized by an important mass loss , these stars could terminate their life as massive onemg wds ( m@xmath12 m@xmath9 ) .
in this paper , we confine our analysis to the co wds having an h - rich envelope , which are the most common type of these compact remnants ( the so called da ) .
let us note that the physical evolution of the internal structure of a white dwarf appears rather simple when compared to other stellar evolutionary phases .
the residual nuclear burnings provide only a negligible contribution to the total luminosity . the highly degenerate core is almost incompressible , so that only the contraction of the most external layers may still contribute to the release of gravitational energy .
only the brightest wds ( logl / l@xmath13 ) suffer a non - negligible energy loss through neutrino production .
the cooling time is then determined by the rate of temperature decrease , which depends on the efficiency of the energy transport from the hot core through the thin opaque envelope . a very efficient thermal conduction , due to the highly degenerate electrons , takes place in the whole co core that rapidly becomes almost isothermal . on the contrary , the more external layers of the wd ( i.e. an he - rich mantel and , eventually , an h - rich envelope , where electron are initially only partially or non degenerate ) are efficient insulators and regulate the energy loss rate
basically , the time scale of the wd evolution is controlled by the heat capacity of the co core and by the opacity of the external wrap .
thus , the main problem of this studies concerns the accurate evaluation of these two physical ingredients and no sophisticated algorithms for the computation of the stellar structure is , in principle , required .
nonetheless , we believe that the best approach should be that based on a full evolutionary code calculation , which allows a coherent description of both the progenitor and the wd evolution and may account for all the contributions to the energy production and transport .    in this paper we critically re - analyze the influence of the assumed theoretical scheme on the predicted cooling time .
in such a way , we try to evaluate the present uncertainty of the cosmo - chronology predictions based on wd models . a rough idea of the present level of uncertainty
may be obtained by comparing some recently published theoretical cooling sequences for h - rich atmosphere wd ( see figure 1 ) .
one may recognize sizeable differences , in particular near the faint end of the cooling sequence .
they are likely due to different evaluations of the basic ingredients of the theoretical cooking : the specific heat in the core , the radiative and the conductive opacities , the treatment of the external convection and the like ; but they could be also partially attributed to differences in the initial models .
we argue that the origin of these uncertainties should be identified before using cooling sequences as cosmic clocks .
owing to the complexity of the matter , our investigation is far from being exhaustive .
many previous works have addressed the question of the reliability of the input physics used in wd models computation . in some case
we do not repeat their analysis .
for example , we do not address the question of the chemical segregation occurring in the core when the crystallization of the oxygen component takes place , because this argument has been exhaustively discussed in the recent literature ( segretain et al .
1994 , salaris et al .
1997 , montgomery et al . 1999 , isern et al . 1997 , 2000 ) . in some case
we extent previous investigations up to the faint end of the cooling sequence , as in the case of the effects of microscopic diffusion .
in the following section we illustrate the input physics used to compute the various cooling sequences . in the third section
we discuss the uncertainties related to the progenitor evolution that determine the internal stratification of the core and the size of the thin external layers : the he - rich mantel and the h - rich envelope .
section four illustrates the dependence of the cooling time scale on some physical ingredients ( eos , opacity and the like ) .
a final discussion follows .
our wd models have been obtained by means of a full evolutionary code . in particular
we have used the franec in the version described by chieffi & straniero ( 1989 ) .
the input physics have been completely revised in order to account for the peculiar conditions developed in wd interiors .
we have adopted the radiative opacities of opal ( iglesias & rogers 1996 ) for high temperature ( log@xmath14 > 4.0 $ ] ) and the alexander & ferguson ( 1994 ) molecular opacities for the low temperatures ( log@xmath14\leq 4.0 $ ] ) .
the conductive opacities have been derived from the work made by itoh and coworkers ( itoh et al .
1983 , mitake , ichimaru , & itoh 1984 , itoh , hayashi , & kohyama 1993 ) .
additional models have been obtained by using the hubbard & lampe ( 1969 ) prescriptions and the table provided by a. potekhin ( see potekhin et al .
as it is well known ( see e.g. mazzitelli 1994 , dantona & mazzitelli 1990 ) , there is a region in the @xmath15 plane not covered by the opal radiative opacities and where the electron conductivity is not yet dominant .
thus one has to extrapolate the radiative opacities in order to fill this gap .
the problem is alleviated by the use , for the lower temperature , of the alexander & ferguson ( 1994 ) opacities , which extend at larger densities with respect to the opal .
present models have been computed adopting a linear extrapolation of the radiative opacity beyond the upper density provided by opal .
we have tested the reliability of this choice by changing the extrapolation method .
negligible effects on the cooling sequences have been found .
nuclear reaction are almost extinct during the wd evolution , but they play a relevant role in determine the chemical structure of the model at the beginning of the cooling sequence .
we have used the rates tabulated by the nacre collaboration ( angulo et al .
1999 ) , except for the @xmath16c@xmath17o . during the he - burning this reaction competes with the @xmath18 and regulates the final c / o ratio in the core ( iben 1972 ) .
the great influence on the resulting wd evolution has been deeply investigated ( dantona & mazzitelli 1990 ; salaris et al .
we have used two different rates , namely the one reported by caughlan et al .
( 1985 ) and that of caughlan & fowler ( 1988 ) .
note that the difference between these two compilation may be roughly considered as representative of the present experimental uncertainty level ( see e.g. buchmann 1996 ) .
nuclear reaction rates are corrected for the weak and intermediate electron screening by using the prescriptions of graboske et al .
1973 and de witt , graboske , & cooper 1973 , and for the strong screening by itoh , totsuji , & ichimaru 1977 and itoh et al . 1979 .
the efficiency of neutrino emission processes are taken from : haft , raffelt , & weiss 1994 ( plasma neutrinos ) , itoh et al .
1989 ( photo and pair neutrinos ) .      as recognized long ago ( salpeter 1961 ) coulomb interactions play a relevant role in wd interior , finally leading to crystallization of the stellar core ( lamb & van horn 1975 ) .
thus , a detailed treatment of the thermodynamic behavior of both liquid and solid matter is a necessary physical ingredient to study the evolution of cold wds .
the high pressure experienced in the whole core ensures the complete ionization of carbon and oxygen .
thus , we have updated and extended the eos for fully ionized matter described by straniero ( 1988 ) . in particular , we have revised the treatment of the electrostatic corrections up to the liquid - solid transition and beyond .
the free energy in the fluid phase can be written as @xmath19 where @xmath20 and @xmath21 are respectively the contribution of the ideal gas of ions and electrons . to compute these terms , as in straniero ( 1988 ) , we have assumed that ions follow the boltzmann distribution , while electrons are described by fermi - dirac integrals , for an arbitrary degree of degeneracy and relativistic state .
@xmath22 is the excess of ionic free energy due to ion - ion coulomb interactions . for this contribution
we adopted the analytical expression by potekhin & chabrier ( 2000 ) that in the region of high ionic coupling parameter ( @xmath23 ) fits the accurate results of recent monte carlo simulations of one component plasma ( de witt & slattery 1999 ) , while for @xmath24 it approaches to the cohen & murphy ( 1969 ) expansion and reproduces the debye - huckel limit for vanishing @xmath25 ( landau & lifshitz 1969 ) .
the quantum diffraction correction to ionic free energy has been neglected in the liquid phase .
@xmath26 represents the excess electron free energy due to electron - electron interactions . for this term we adopted the non relativistic expression by tanaka , mitake , & ichimaru 1985 .
we neglected the relativistic @xmath26 .
this is a reasonable approximation since at high density , when degenerate electrons become relativistic , the electron coupling parameter @xmath27 .
finally , for @xmath28 , that describes the excess due to ion - electron interactions , we adopted the analytical expression given by potekhin & chabrier ( 2000 ) .
as it is well known , the contribution of the ion - ion interactions increases with density , eventually leading to crystallization . according to potekhin & chabrier ( 2000 ) ,
the liquid - solid phase transition has been assumed at @xmath29 . in the solid phase ( @xmath30 ) the free energy can be written as @xmath31 where the electron terms @xmath21 and @xmath26 are the same as in the fluid phase . for the free energy of the ionic crystal @xmath32 we used the analytical expression by farouki & hamaguchi ( 1993 ) obtained by fitting numerical models of molecular dynamic .
the contribution of the ion - electron @xmath28 in the solid phase is from potekhin & chabrier ( 2000 ) .
finally , the term @xmath33 represents the quantum correction to the thermodynamic behavior of the ionic coulomb crystal .
this contribution is very important in late wd evolutionary phases , in fact , as the wd cools down , the crystallized co core reaches a quantum state ( diffraction parameter @xmath34 , where @xmath35 denotes the ion plasma frequency ) .
when this occurs , the ionic contribution to the specific heat decreases as @xmath36 , thus depleting the main thermal reservoir of the star : it is called debye cooling phase . for this term
we adopted the analytical expression described by stolzmann & blocker ( 2000 ) , which is based on the free energy of a coulomb crystal studied by chabrier , ashcroft , & de witt ( 1992 ) .
the various thermodynamic quantities for pure carbon and pure oxygen have been obtained by analytically deriving the corresponding free energy .
finally , an additive volume interpolation is used to calculate the thermodynamic quantities of a co mixturene . ] .    to describe the outer layers of partially ionized helium and hydrogen
, we adopted the eos of saumon , chabrier , & van horn ( 1995 , scvh ) .
this eos requires both high pressure and low pressure extensions . at high pressure ( p@xmath37 dyn/@xmath38 ) , which is in any case large enough to guarantee a full ionization of h and he
, we extended the scvh eos following the same procedure described above for the co core .
the match between the two eos is generally good . at low pressure
, we have used a perfect gas ( including h , h@xmath39 , h@xmath40 , h@xmath41 , he , he@xmath39 and he@xmath42 ) plus radiation . in such a case
the classical saha equation has been used to derive the population of the various species .
for each temperature , the precise value of the maximum pressure of the perfect gas ( or the minimum of the scvh ) has been varied in order to guarantee a smooth transition between the two eos .
tables of this eos are available on the web .      in order to fix the external boundary condition of a stellar model , appropriate model atmosphere
are needed . in the last decade
wd atmosphere theory has significantly improved ( bergeron , waesemael , & fontaine 1991 ; bergeron , saumon , & waesemel 1995 ; bergeron , waesemel , & beauchamp 1995 ; saumon & jacobson 1999 ; bergeron 2000 ) .
as it is well known , in the high density / low temperature atmospheres of an old wd the emerging electromagnetic flux significantly departs from the black body spectrum .
this is mainly due to the molecular hydrogen recombination and the consequent collision - induced absorption ( cia ) by h@xmath40 - h@xmath40 collisions .
cia is the main source of opacity in the infrared and the main cause of increasingly blue color indices for decreasing effective temperature ( t@xmath43 ) in cold wds ( hansen 1998 , 1999 ; saumon & jacobson 1999 ; bergeron 2000 ) .
detailed model atmospheres are also crucial ingredient to transform the theoretical quantities ( l , t@xmath43 ) into the observational ones ( magnitude , color ) . in the present computations
we have used a solar scaled t-@xmath44 relation ( see chieffi & straniero for details ) up to the onset of the cooling sequence . then , namely when log l / l@xmath45 , we adopt the model atmospheres of bergeron , saumon , & wesemael ( 1995 ) , plus a low temperature extension including the effects of the cia , which has been kindly provided us by p. bergeron . since no metals are included in these models , they are consistent with the usual hypothesis of complete sorting of the external layers ( see next subsection ) .
convective boundaries are fixed by using the schwarzschild criterion . concerning the progenitor evolution , the extension ( in mass ) of the convective core during the central he burning
have a great influence on the amount of c ( and o ) in the most internal layers of a wd ( imbriani et al .
2001 and dominguez et al . 2001 ) . during this phase
, we use the algorithm described by castellani et al .
( 1985 ) to take into account the growth of the convective core induced by the conversion of he ( low opacity ) into c and o ( large opacity ) and the resulting semiconvection . breathing pulses , which occur when the central he mass fraction decreases below about 0.1 , are usually neglected ( but see the discussion in section 3.3 ) .
note that any mechanism that increases the size of the well mixed region during the final part of the he - burning ( mechanical overshoot , semiconvection , breathing pulses or rotational induced mixing ) leads to a reduction of the resulting amount of c in the central region of the wd ( see e.g. imbriani et al .
2001 ) .
theoretical studies ( fontaine & michaud 1979 , iben & macdonald 1985 , 1986 , althaus et al .
2002 ) , supported by observational evidences ( see e.g. bergeron , ruiz , & legget 1997 ) , indicates that , as a consequence of the gravitational settling of heavy elements , da wds should have a practically pure h envelope and a pure he mantel .
in addition , as suggested by salaris et al .
( 1997 ) , the carbon and oxygen profiles left by the he - burning is smoothed by rayleigh - taylor instabilities . in order to account for the element diffusion , when log l / l@xmath45 , we adjust the composition of the envelope and the mantel : all the residual h is putted on the top of a pure he layer .
the total mass of h and he is conserved .
we also modify the internal c and o profile according the the prescription of salaris et al .
( 1997 ) .
owing to the larger @xmath25 , oxygen crystallization occurs when carbon is still liquid . for a certain time
, this occurrence produces an unstable stratification of the liquid phase and , in turn , an efficient mixing ( isern et al .
1997 , salaris et al . 1997 ) .
this phenomenon is presently not included in our calculations .
first of all one should know the physical and the chemical structure of the star at the beginning of the cooling sequence . since our investigation regards da white dwarfs , we need to know the internal stratification of the co core and the extension ( in mass ) of both the he and h - rich layers . the formation of the co core begins during the central he - burning phase and proceeds in the following agb phase . in principle ,
a lengthy computation through a series of recursive thermal pulses is required , until the mass loss erodes the envelope and the star becomes a wd ( see e.g. iben & renzini 1983 ) .    during the thermally pulsing agb phase ,
the mass of the co core is accreted with the ashes of the he - burning shell , while the mass between the two shells is progressively reduced ( see e.g. straniero et al .
thus , it is very important to determine the duration of the agb , which is evidently affected by mass loss : the longer the agb phase , the larger the final mass of the co core and the smaller the mass of the he - rich mantel .
agb mass loss rate ranges between 10@xmath46 and 10@xmath47 m@xmath0/yr , but it is still unclear how it is related to the changes of the chemical and physical structure occurring in the cool atmospheres of these giants ( see habing 1996 ) . in addition , if the wd progenitor is a low mass star ( m@xmath481.5 m@xmath0 ) , the mass loss during the first ascent to the tip of the red giant branch ( rgb ) must be also carefully taken into account .
the spread in color of the horizontal branch ( hb ) stars in globular clusters is commonly interpreted as an evidence of the different mass loss experienced by similar stars during the rgb phase : bluer hb stars should have a smaller envelope than red hb stars . since the agb phase of blue ( reduced envelope ) hb stars will be significantly shorter than that of a red hb star , the resulting wd will be different .
in other words , stars with similar initial mass and chemical composition might produce white dwarfs of different masses and/or different internal stratification , depending on the mass loss history .
not only the mass , but also the initial composition of a star determines its evolution and , possibly , its final fate . a detailed descriptions of the internal c and o profiles in wds generated by stars with mass ranging between 1 and 7 m@xmath0 and metallicity 0@xmath10z@xmath100.02 may be found in dominguez et al .
( 2001 ) .
this variety of possible progenitor evolutions could frighten the reader .
our scope is to quantify the uncertainty on the cooling time of wd , for different galactic components .
basing on some numerical experiments , we have evaluated the influence of the pre - wd history on the age calibration .
starting from zero age he - burning stellar structures ( i.e. an he core surrounded by an h - rich envelope , in which a suitable mixture of helium and metals is adopted ) , we have computed the evolutionary sequences up to the wd stage . at the beginning of the cooling phase ( @xmath49 )
the envelope has been completely sorted in a pure he buffer and a pure h outer envelope mimicking the effect of the very efficient gravitational settling .
table 1 reports the initial parameters of the models and the properties of the resulting white dwarfs , namely : identification label of the model , initial mass ( m@xmath50 ) , metallicity ( z ) , initial he - core mass ( m@xmath51 ) , final mass ( m@xmath52 ) , final amount of he ( @xmath53m@xmath51 ) , final amount of h ( @xmath53m@xmath54 ) , central carbon mass fraction ( x@xmath55 ) , mass of the homogeneous central region ( q = m@xmath56/m@xmath52 , see below for a definition of this parameter ) , wd age when log l / l@xmath0=-5.5 .
let us firstly discuss models b , b1 , b2 .
the total mass has been maintained constant ( m=0.6 m@xmath9 ) during the whole evolution for both case b and b1 .
they differ for the initial mass of the h - rich envelope , namely : @xmath57 and @xmath58 respectively for b and b1 . the third case ( b2 )
has been computed starting from a larger mass ( m=0.64 m@xmath9 ) and imposing a constant mass loss ( @xmath59 m@xmath60 ) since the onset of the thermally pulsing phase . in this case , the initial he - core mass is in agreement with recent evolutionary computation of h - burning models of low mass stars ( see e.g. straniero , chieffi , & limongi 1997 ) .
note that the initial total mass has been chosen to obtain a final ( wd ) mass similar to the ( constant ) one of case b and b1 .
a low metallicity ( z=@xmath61 ) , adequate for old halo wd progenitors , has been chosen for all these three models .
figure 2 illustrates the evolutionary tracks in the hr diagram .
the central he - burning phase is similar in the three cases . at the central he exhaustion
, model b1 evolves toward the red part of the hr diagram and develops the typical double shell structure of an agb star .
when the faster he - burning shell get closer to the h - burning one , it dies down and the star enters in the thermally pulsing phase ( for a more detailed description see e.g. iben & renzini 1983 ) .
after two major thermal pulses , the envelope mass is reduced down to about 0.001 m@xmath0 ; then , the star leaves the agb . during the post - agb ,
a further ( last ) thermal pulse takes place ; the star re - expand and comes back to the asymptotic giant branch .
finally , after a second horizontal crossing of the hr diagram , both shells definitely die down and the model settles on the wd cooling sequence .    owing to the very small envelope mass , after the central he exhaustion , model b remains in the blue side of the hr diagram , so skipping the classical agb phase .
when the outgoing he shell approaches the h - rich layer it dies down .
note that the h - burning shell has been practically inactive up to this moment .
the following contraction of the star induces some irregular and short h - ignitions .
finally this model attains its cooling sequence .    in model
b2 the duration of the agb and the final wd mass are mainly determined by the assumed mass loss rate , which is , in any case , quite typical for faint tp - agb stars .
this models experience 8 thermal pulses ; then , the envelope is reduced down to a certain critical value and the model moves toward the blue side of the hr diagram ( see figure 2 ) .
note that the mass loss has been stopped when the effective temperature becomes larger than 7000 k. at variance with case b1 , no more thermal pulses are found in the post - agb phase .
the resulting profiles of carbon and oxygen in the core of the three final wd models are shown in figure 3 .
the most internal region is produced during the central he - burning and its flatness is a consequence of the convective mixing occurring in the stellar core . as noted by salaris et al .
( 1997 ) , the region immediately outside the flat profile left by the central convection presents a positive gradient of mean molecular weight .
the consequent rayleigh - taylor instability induces a mixing and a further little reduction of the central carbon . in column 8 and 9 of table 1
we have reported the central c mass fraction and the mass of the homogeneous central region ( in unit of the total wd mass ) as they should result after the action of this rayleigh - taylor instability ( q = m@xmath56/m@xmath52 ) . the growing c profile in the more external layer is produced during the agb , when the he - burning shell advances in mass toward the stellar surface .
the carbon peak at the external border of the co core reflects the incomplete he - burning left by the last thermal pulse .
the three chemical structures are rather similar .
in particular the amount of carbon ( oxygen ) in the center of the wd and the extension of the inner flat region are practically the same .
the most relevant differences is the size ( in mass ) of the he - rich mantel ( see column 6 in table 1 ) . in model b1
the thickness of the mantel is about half of that found in case b. this is the memory of the last thermal pulse experienced by model b1 in the post - agb phase , when the h - burning shell was almost extinct .
the amount of hydrogen left in the envelope of the wd is practically the same in the three models .
this confirms the fact , already demonstrated in many papers , that for each metallicity and core mass , it exist a maximum mass of the h - rich envelope for a wd ( iben 1982 , fujimoto 1982a , b , castellani , deglinnocenti , & romaniello 1994 , piersanti et al .
when the h mass exceeds such a limit , the h - burning is efficiently activated at the base of the envelope .
it drives the expansion of the envelope leading to the formation of a giant star rather than to a white dwarf .
the cooling sequence obtained starting from the three pre - wd models illustrated in this subsection are shown in figure 4 .
they are very similar , with a little difference ( less than 0.5 gyr ) at the faint end .
in particular , the cooling of model b1 is a bit faster because the smaller he mantel .
note that these three models reach the wd stage through a very different path and they may be considered as representative of the possible changes occurring during the progenitor evolution . in summary , the knowledge of the cooling time scale at the faint end of the wd sequence ( namely at log l / l@xmath0=-5.5 ) is limited by a 3 - 4% error as a consequence of the uncertainty on the leading parameters characterizing the pre - wd history .
in many recent studies of cosmochronology ( salaris et al .
1997 , richer et al .
2000 ) wd models based on progenitors having solar chemical composition have been used .
this is certainly adequate to study the age of the disk and its components ( open clusters ) .
however , when the main goal are the oldest component of our galaxy , halo or globular cluster stars , low metallicity models are more appropriate . in this subsection
we would try to evaluate the influence on the cooling sequence of the progenitor chemical composition . for this reason we have computed an additional model producing a wd having m=0.6 m@xmath0 , but with a solar metallicity , namely z=0.02 .
the properties of this model are reported in the fourth row of table 1 ( model b3 ) .
a comparison with model b ( z=0.0001 ) reveals only small differences . as already discussed by dominguez et al .
( 2001 ) , the central carbon abundance is slightly lower at large metallicity ; they found a maximum 10% variation when the metallicity is varied between 0.0001 and 0.02 .
note that the extension of the innermost homogeneous region decreases by increasing the metallicity ( see the values reported in column 9 of table 1 ) .
as expected , since the h - burning is more efficient at large metallicity ( because the larger amount of cno ) , the mass of the h - rich envelope is significantly reduced when the star attains the cooling sequence . in spite of this difference , at the faint end model b3 is only 2% younger than model b. the cooling sequences obtained for the two metallicities are compared in figure 5 .
the amounts of carbon and oxygen left in the core of a wd have a great influence on its cooling rate ( dantona & mazzitelli 1990 , salaris et al .
1997 ) and plays an important role in determining the observable outcome of a type ia supernova ( dominguez et al .
the chemistry of a wd is determined by the competition of the two major nuclear reactions powering the he burning , namely the 3@xmath62 ( the carbon producer ) and @xmath16c@xmath17o ( the carbon destroyer ) . as deeply discussed by imbriani et al .
( 2001 ) , the final c / o not only depends on the rate of these two reactions , but it is significantly influenced by the efficiency of the convective mixing operating near the center of an he burning star .
the pre - wd models discussed so far , has been obtained by neglecting any source of extra - mixing , like mechanical overshoot , possibly occurring at the external edge of the convective core .
we have also suppressed the so called breathing pulses , which take place when the central mass fraction of he is reduced down to @xmath630.1 ( see castellani et al .
1985 and references therein ) . as shown by imbriani et al .
( 2001 ) , the amount of carbon left in the center is almost insensitive to the occurrence of the convective overshoot , but it is strongly influenced by any additional mixing occurring in the final part of the he - burning , when the carbon previously accumulated by the 3@xmath62 is efficiently converted into oxygen via the @xmath16c@xmath17o reaction .
this is exactly the job done by breathing pulses ( bps ) . in this subsection
we study the combined effects of a change of the @xmath16c@xmath17o rate and the inclusion of bps .
models b to b3 have been obtained by using a high rate for the @xmath16c@xmath17o reaction ( namely , the one reported by caughlan et al .
1985 , hereinafter c85 ) .
as it is well known , the adoption of a low rate leads to a significantly larger central abundance of carbon .
model b - low has been obtained by adopting the rate reported by caughlan & fowler ( 1988 , hereinafter cf88 ) , which is about 0.4 times the rate of c85 at the relevant temperatures ( around @xmath64 k ) .
recent experimental determinations of this rate ( buchmann 1996 and 1997 , kunz et al .
2001 , 2002 ) seem to confirm that the actual value should be intermediate between cf88 and c85 .
then , we have also repeated these two calculations by including the effect of the breathing pulses ( model b - bp and b - low - bp ) .
the internal chemical profiles of models b - low and b - low - bp are compared in figure 6 , while the influence on the cooling sequence is illustrated in figure 7 . in synthesis , the whole uncertainty ( convection plus nuclear reaction rate ) due to the he - burning phase induces a 9% uncertainty in the estimated cooling time at the faint end of the wd sequence ( see the values reported in the last column of table 1 ) .
it is worth noting that this is only a lower limit of the uncertainty caused by the possible variation of the internal composition .
in fact , as discussed by isern et al .
( 1997 ) and montgomery et al .
( 1999 ) , the delay of the cooling due to the chemical segregation upon crystallization is larger when @xmath65 .
therefore , the occurrence of this phenomenon should hampers the differences between case b - low and case b - bp .
a more quantitative estimation may be deduced from the results reported in table 2 of montgomery et al .. note that our b - low model have an internal stratification that roughly correspond to their 50:50 case , while our model b have a composition profile similar to the one of salaris et al .
( 1997 ) . then , according to montgomery et al .
, element segregation should increase the age difference between b - low and b - bp of at least 0.6 gyr , thus bringing the whole uncertainty on the age at the faint end of the cooling sequence up to 13 - 14% .
we have already discussed the existence of a maximum value for the mass of the h - rich envelope of a wd .
however , the actual amount of h at the top of a wd could be substantially lower .
in fact , stellar wind may continue during the post - agb phase so that the h - rich envelope could be strongly reduced and even completely lost .
it is well known that in an important fraction ( about 10 - 20 @xmath66 ) of the whole wd population the h - rich envelope is practically absent ( db - wds ) .
the minimum detectable amount of h is very low ( about @xmath67 m@xmath0 ) .
the post - agb mass loss is the easiest explanation of the lack of h in the atmosphere of the db wd .
an alternative possibility was advanced by iben & renzini ( 1983 ) .
they suggested that the progenitors of the db could have suffered a last thermal pulse during the post - agb , when the h - burning shell is almost extinguished . in this condition
, the convective shell generated by the thermal pulse could completely mix the he and the h - rich layers .
then , the h dredged down at high temperature is rapidly consumed .
we recall that a similar post - agb thermal pulse has been encountered in our model b1 , but no evidence of a deep h mixing has been found ( in this particular case ) .
dantona & mazzitelli ( 1989 ) and tassoul et al .
( 1990 ) have investigated the wd properties by changing the amount of h in the envelope . however , due to the lack of reliable opacity and eos , they did nt extent the computation down to the faint end of the cooling sequence .
then , we have computed two additional cooling sequences by artificially reducing the h - rich layer of model b , namely : @xmath53 m@xmath68 and @xmath69 for model b - lh and b - vlh , respectively ( row 5 and 6 of table 1 ) .
the comparison between these two models and model b is shown in figure 8 . as expected
, the cooling is faster in models with a smaller envelope .
the age reduction at log l / l@xmath70 would be as large as 1 gyr if the residual mass of the h - rich envelope was of the order of @xmath67 m@xmath0 .
in this section we address the uncertainty caused by some input physics directly affecting the wd evolution .      degenerate electrons dominates the heat conduction in white dwarf interiors .
two different tabulations of the thermal conductivity by electrons are commonly used in the computation of wd models , namely : the one presented by hubbard & lampe ( 1969 , hereinafter hl69 ) and those by itoh and coworkers ( itoh et al .
1983 , mitake et al .
1984 , itoh et al .
1993 , hereinafter i93 ) .
hl69 considered partially degenerate electrons scattered by ions .
they also included the effect of the electron - electron interaction , which is negligible in extreme degeneracy regime ( due to the pauli exclusion principle ) , but becomes important if degeneracy is not strong , as in the external region of a wd .
the main limitation of this calculation is due to the assumption of non - relativistic electrons . in practice ,
the hubbard and lampe conductive opacities are only adequate for densities lower than @xmath71 g/@xmath72 .
actually , a stronger limitation of this calculation is often ignored .
in fact , in order to smooth out the discontinuity in the conductive opacity occurring at the liquid / solid phase transition , the conductive opacities were artificially and progressively increased for @xmath7310 .
let us recall that at the end of the sixties it was believed that the phase transition from liquid to solid would occur for rather low values of the coulomb coupling parameter ( @xmath74 ) .
this explains the choice of a particularly low @xmath25 to manage the discontinuity occurring at the phase transition .
nevertheless , the most recent studies find that the crystallization begins at @xmath75 ( we adopt 175 ) , in any case well above @xmath76 .
the work made by itoh and coworkers is a substantial improvement of the calculation of electron conduction .
they treat inter - ionic correlations and electron - ion interactions by means of an improved theory , based on one - component monte carlo calculation .
however , since they assume fully degenerate electrons , their opacities are only accurate for large value of the degeneracy parameter . at lower density ( or at larger temperature ) , i93 probably underestimate the effects of the electron - electron interactions .    in figure 9
we illustrate the situation in the case of the core of our reference model ( b ) .
the empty area in the left panel shows the allowed region in the case of hl69 , while the right panel refers to i93 .
the two ( almost vertical ) series of arrows shows the evolution of the physical conditions at the center and at the external boundary of the core .
this figure demonstrate that hl69 is not adequate to describe the heat conduction in the major part of the wd structure , even if the density is somewhat lower than the relativistic limit . on the contrary
the condition of extreme degeneracy is well satisfied and the i93 assumption seems appropriated for the whole co core .
concerning the more external layers , figure 10 illustrate the case for the he - rich mantel . as in figure 9 ,
the left and the right panel show the region of validity for hl69 and i93 , respectively .
the two series of arrows marks the boundaries of the mantel .
hl69 should be preferred in the first part of the cooling sequence , since electrons may be partially degenerate while @xmath25 is generally rather small .
the situations change toward the faint end of the cooling sequence , during which the i93 should be more appropriate .
note that there exist a thin region in the t-@xmath77 plane in which both calculation are , in principle , valid .
it happens , for example in the case of an he - rich mixture , for @xmath78 . the thick line in figure 10 marks this condition .
however , a comparison along the @xmath79=0.1 line shows that the hl69 opacities are generally larger ( 20 - 40 % ) than the i93 ones .
such a mismatch reveals a more profound difference between the two calculations , which can not be easily understood .    in figure 11
we report the cooling sequences obtained by using the hl69 and i93 ( case b - hl69 and case b , respectively ) . in the former case
, we have generated the conductive opacities by means of the original code provided by hubbard .
as expected , the lower cooling rate obtained with hubbard & lampe reflects the larger value of the conductive opacities . note that the major discrepancies between the two cooling sequences occurs at about log(l / l@xmath80=-4 , exactly when electron conduction becomes the dominant mechanism of energy transportation at the base of the convective envelope .
as it is well known ( see e.g. tassoul et al .
1990 ) , this occurrence produces a temporary reduction of the cooling rate . since many authors use a combination of the hl69 and i93 , in figure 11 we have also reported an additional model obtained by using i93 for @xmath81 and hl69 elsewhere ( case b - comb ) .
this test shows that the discrepancy mainly concerns the conductive opacity in the weakly - degenerate regime . in conclusion ,
no ones of the prescriptions widely adopted in the computation of models of wds appear adequate to describe the whole structure of these cool stars .
unfortunately , the substantial disagreement , at temperature and density for which both computations should be correct , suggests extreme caution in using combined tables . since the differences of these two theoretical prescriptions imply a 17% variation in the estimated wd age ( at the faint end of the cooling sequence , see the last column in table 1 ) , we conclude that the calculation of the conductive opacities ( especially in partially degenerate regime ) deserves much attention in view of a reliable calibration of the wd cooling time . in this context
a recent paper by potekhin et al .
( 1999 , see also potekhin 1999 ) address this problem , extending the computation of electron conductivity to the weakly - degenerate regime .
other improvements , as , in particular , refined ion structure factors and multi - phonon scattering , have been included .
the conductive opacity obtained by potekhin and coworkers are intermediate between those of hl69 and i93 .
nevertheless , as shown in figure 11 , the resulting cooling sequence does not substantially differs from the one obtained by using the i93 . on the contrary ,
the cooling time is significantly smaller then those obtained in the model b - hl69 and b - comb . some properties of the model obtained by means of the electron conductivity calculated by potekhin et al .
( 1999 ) are reported in table 1 ( model b - pot ) .
chabrier et al .
2000 compare wd cooling sequences obtained by changing the prescriptions for the contribution of the ion - electron interactions to the free energy . in particular they found a significant decrease of the cooling time ( about 10% ) when the prescription of potekhin & chabrier ( 2000 ) are used in the solid phase instead of the result given by yakovlev & shalybkov ( 1989 ) .
as pointed out by potekhin and chabrier , the ion - electron interactions are generally negligible , except in the solid state at very high density . for carbon the ion - electron contribution to
the heat capacity becomes very important at t=@xmath82 k and @xmath83 g @xmath84 but at t=@xmath71 k , @xmath77 must be larger then @xmath85 g @xmath84 , a density never attained in co white dwarfs ( see figure 8 of potekhin & chabrier ) .
we have investigated the overall effect of this contribution to the heat capacity by computing an additional model ( model b - noie ) by fully neglecting the ion - electron interactions .
we recall that our `` standard '' computation are based on the potekhin & chabrier prescription .
the result is shown in figure 12 ( see also table 1 ) .
this test clearly demonstrates that the contribution of the ion - electron interaction derived by potekhin & chabrier produces a negligible effect on the cooling of co white dwarfs .
many studies have shown that the time - scale of diffusion is short enough to sort the external layers of a wd , thus explaining the observed mono - elemental feature of wd spectra ( fontaine & michaud 1979 , muchmore 1984 , paquette et al .
many authors have discussed the effects of diffusion on the wd structure and evolution ( see in particular iben & macdonald 1985 and 1986 , althaus et al .
all the models presented here have been obtained by assuming a complete sorting of the h and he layers at the beginning of the cooling sequence ( @xmath86 ) .
since we are mainly interested in quantify the overall variation of the cooling time , we have computed a new model in which the chemical stratification of both the envelope and the mantel is the one modeled by the nuclear burning and , eventually , by convection ( model b - nodif ) .
we note that the adopted model atmospheres does not include the effect of metals .
since the assumed metallicity is quite low ( z=0.0001 ) , the consequent inconsistency should not dramatically affect our result .
the comparison between the properties of this model and those of the _ standard _ one ( model b ) is illustrated in figure 13 and in table 1 .
in the standard case , the larger radiative opacity of the pure h induces the formation of a deeper convective envelope . as a consequence , the slow cooling phase , which takes place when the internal border of the convective layer reaches the region dominated by electron conduction ,
is anticipated . in conclusion , at the faint end of the cooling sequence , the model in which the effect of microscopic diffusion are neglected is a 11% younger than the standard model .
note that this is the younger model we found .
as it is well known the cooling time depends , among the other things , on the wd mass .
figure 14 shows models in the mass range 0.5 - 0.9 @xmath87 computed under the same prescriptions adopted for the case b. in table 2 we have summarized the properties of these models .
isochrones and luminosity functions based on these models are available on the web .
in the last decade cosmochronology based on wds has become a promising tool of scientific inquiry thanks to both the availability of a large amount of high quality data for cool and faint wds and the improvement of the understanding of the high density plasma physics .
however , as demonstrated by the rather large discrepancies among the recently published theoretical cooling sequences , a firm calibration of the age - luminosity relation is still not available .    we have analyzed some of the main sources of uncertainty affecting the theoretical cooling time .
physical and chemical parameters characterizing the white dwarfs and the progenitors evolutions have been revised .
concerning progenitors , we found that the larger uncertainty is due to the combined action of convective mixing a nuclear reactions operating during the central he - burning phase .
both these processes are largely affected by theoretical and experimental uncertainties .
they determine the amount of c ( and o ) left in the core of the wd and , in turn , have a great influence on the predicted cooling rate .
a conservative analysis allow us to conclude that the overall impact of the uncertainties due to the progenitors evolution on the estimated wd ages at logl / l@xmath88 is of the order of 2 gyr .
we hope that this uncertainty will be significantly reduced in the next future thanks to the renewed effort of nuclear physicists in measuring the @xmath16c@xmath17o reaction rate at low energy ( see gialanella et al .
2001 and kunz et al .
2001 ) .
concerning wd physics we emphasize the relevance of a reliable description of the electron conductivity which is the main mechanism of energy transport in wd interiors .
actually , a consistent part of the large discrepancies found in the comparisons of published cooling sequences may be due to the adopted conductive opacity .
the old hubbard and lampe ( 1969 ) conductive opacity are probably overestimated , even in the weakly degenerate regime .
the latest computations by itoh and coworkers and potekhin et al .
( 1999 ) imply a substantial reduction of the age . at log@xmath89 , we obtain models @xmath90 gyr younger than that obtained with hl69 .
let us conclude by noting that most of the uncertainties discussed in this paper mainly affect old wds , but have negligible influence on the age estimated for young wds , whose luminosity is larger than @xmath61 l@xmath0 .
this implies that the presently available theoretical scenario may be safely used to estimate the age of young stellar systems as , for example , the intermediate age open clusters ( t@xmath483 gyr ) .
we wish to thank p. bergeron for kindly providing us the models atmosphere and s. deglinnocenti , luciano piersanti and jordy isern for the many helpful comments .
we thank also s. shore for the careful reading of the manuscript and the several useful discussions .
we are deeply indebted to v. castellani for continuous encouragements and for a critical review of the manuscript .
this work has been supported by the miur italian grant cofin2000 .
alexander d.r . , & ferguson j.w .
1994 , , 437 , 879 althaus , l.g . , serenelli , a.m. , corsico , a.h . ,
& benvenuto , o.g .
2002 , , 330 , 685 angulo , c. et al . 1999 ,
a , 656 , 3 benvenuto , o.g . , & althaus , l.g .
1999 , , 303 , 30 bergeron ,
p. 2000 , private comunication bergeron , p. , ruiz , m.t .
, & leggett , s.k .
1997 , , 108 , 339 bergeron , p. , saumon , d. , & wesemael , f. 1995 , , 443 , 764 bergeron , p. , wesemael , f. , & beauchamp , a. 1995 , pasp , 107 , 1047 bergeron , p. , wesemael , f. , & fontaine , g. 1991 , , 367 , 253 bradley , p.a . , &
winget , d.e .
1991 , , 75 , 463 bradley , p.a .
, winget , d.e . , & wood , m.a .
1992 , , 391 , l33 buchmann , l. 1996 , , 468l , 127 buchmann , l. 1997 , , 479l , 153 cassisi , s. , castellani , v. , & tornamb , a. 1996 , , 459 , 298 castellani , m. , limongi , m. , & tornamb , a. 1995 , , 450 , 275 castellani , v. , chieffi , a. , tornamb , a. , & pulone , l. 1985 , , 296 , 204 castellani , v. , deglinnocenti , s. , & romaniello , m. 1994 , , 423 , 266 caughlan , g.r . , &
fowler , w.a .
1988 , atomic data nucl .
data tables , 40 , 283 caughlan , g.r . , fowler , w.a . , harris , m.j . , & zimmermann , b.a .
1985 , atomic data nucl .
data tables , 32 , 197 chabrier , g. , ashcroft , n.w .
, & de witt , h.e .
1992 , nat .
, 360 , 48 chabrier , g. , brassard , p. , fontaine , g. , & saumon , d. 2000 , , 543 , 216 chieffi a. , & straniero o. 1989 , , 71 , 47 clemens , j.c .
1993 , baltic .
astron . , 2 , 407 cohen , e.g.d . , & murphy , t.j . 1969 ,
fluids , 12 , 1404 dantona , f. , & mazzitelli , i. 1989 , , 347 , 934 dantona , f. , & mazzitelli , i. 1990 , ara&a , 28 , 139 de witt , h.e . , graboske , h.c .
, & cooper , m.s .
1973 , , 181 , 439 de witt , h. , & slattery , w. 1999 , contrib .
plasma phys . , 39 , 97 dominguez , i. , chieffi , a. , limongi , m. , & straniero , o. 1999 , , 524 , 226 dominguez , i. , hflich , p. , & straniero , o. 2001 , , 557 , 279 dominguez , i. , straniero , o. , tornamb , a. , & isern , j. 1996 , , 472 , 783 dominguez , i. , tornamb , a. , & isern , j. 1993 , , 419 , 268 farouki , r.t . , & hamaguchi , s. 1993 , phys .
e. , 47 , 4330 fontaine , g. , & michaud , g. 1979 , , 231 , 826 fujimoto , m.y .
1982a , , 257 , 767 fujimoto , m.y .
1982b , , 257 , 752 garcia - berro , e. , ritossa , c. , & iben , i.jr .
1997 , , 485 , 765 gialanella et al .
2001 , eur .
j.a . , 11 , 357 graboske , h. , de witt , h. , grossman , a. , & cooper m. 1973 , , 181 , 457 haft , m. , raffelt , g. , & weiss a. 1994 , , 425 , 222 habing , h.j .
1996 , a&a rev , 7 , 97 hansen , b.m .
1998 , nat . , 394 , 860 hansen , b.m . 1999 , , 520 , 680 herwig , f. 1997 , in proceedings of the 10th european workshop on white dwarfs , held in blanes , spain , 17 - 21 june 1996 .
j. isern , m. hernanz , and e. gracia - berro .
publisher : dordrecht : kluwer academic publishers , 1997 hubbard , w.b . , & lampe , m. 1969 , , 18 , 297 iben , i.j . , 1972 , , 178 , 433 iben , jr.i .
1982 , , 259 , 244 iben , jr.i . , & macdonald , j. 1985 , , 296 , 540 iben , jr.i . , & macdonald , j. 1986 , , 301 , 164 iben , jr.i . , & renzini , a. 1983 , ara&a , 21 , 271 iben , jr.i . , &
tutukov , a.v .
1986 , , 311 , 742 iglesias , a. , & rogers , f.j .
1996 , , 464 , 943 imbriani , g. , limongi , m. , gialanella , l. , terrasi , f. , straniero , o. , & chieffi , a. 2001 , , 558 , 903 isern , j. , garcia - berro , e. , hernanz , m. , & chabrier , g. 2000 , , 528 , 397 isern , j. , mochkovitch , r. , garcia - berro , e. , & hernanz , m. 1997 , , 485 , 308 itoh , n. , adachi , t. , nakagawa , m. , kohyama , y. , & munakata , h. 1989 , , 339 , 354 itoh n. , hayashi , h. , & kohyama , y. 1993 , , 418 , 405 itoh n. , mitake s. , iyetomi h. , & ichimaru s. , 1983 , , 273 , 774 itoh , n. , totsuji , h. , & ichimaru , s. 1977 , , 218 , 477 itoh , n. , totsuji , h. , ichimaru , s. , & de witt , h. 1979 , , 234 , 1079 kippenhahn , r. , kohl , k. , & weigert , a. 1967 , zeitschrift fur astrophysik , 66 , 58 kippenhahn , r. , thomas , h.c . , & weigert , a. 1968 , zeitschrift fur astrophysik , 69 , 265 kunz r. , jaeger f.m . , mayer a. , & hammer j.w .
2001 , phys .
lett . , 86 , 3244 kunz r. , jaeger f.m . , mayer a. , hammer j.w . , staudt , g. , harissopulos , s. , & paradellis , t. 2002 , , 567 , 643 lamb , d. q. , & van horn , h. m. 1975 , , 200 , 306 landau , l.d . , & lifshitz , e.m .
1969 , `` statistical physics '' , pergamon , oxford livne , e. , & tuchman , y. 1988 , , 332 , 271 mazzitelli , i. 1994 , `` the equation of state in astrophysics '' , proceedings of iau colloquium no .
147 , saint - malo , france , 14 - 18 june 1993 , cambridge : cambridge university press , 1994 , ed .
g. chabrier and e. schatzman , p.144 mestel , l. 1952 , , 112 , 583 metcalfe , t.s . ,
winget , d.e . , & charbonneau , p. 2001
, , 557 , 1021 mitake , s. , ichimaru , s. , & itoh , n. 1984 , , 277 , 375 montgomery , m.h . , klumpe , e.w .
, winget , d. e. , & wood , m.a .
1999 , , 525 , 482 muchmore , d. 1984 , , 278 , 769 paquette , c. , pelletier , c. , fontaine , g. , & michaud , g. 1986 , , 61 , 197 paresce , f. , de marchi , g. , & romaniello , m. 1995 , , 440 , 216 piersanti , l. , cassisi , s. , iben , i.jr . , &
tornamb , a. 2000 , , 535 , 932 potekhin , a.y . 1999 ,
a&a , 351 , 787 potekhin , a.y . , baiko , d.a . ,
haensel , p. , & yakovlev , d.g . 1999 ,
a&a , 346 , 345 potekhin , a.y . , & chabrier , g. 2000 , phys .
e. , 62 , 8554 richer , h.b .
1997 , , 484 , 741 richer , h.b .
2000 , , 529 , 318 ritossa , c. , garcia - berro , e. , & iben , jr.i 1996 , , 460 , 489 ritossa , c. , garcia - berro , e. , & iben , jr.i 1999 , , 515 , 381 salaris , m. , dominguez , i. , garcia - berro , e. , hernanz , m. , isern , j. , & mochkovitch , r. 1997 , , 486 , 413 salaris , m. , garcia - berro , e. , hernanz , m. , isern , j. , & saumon , d. 2000 , , 544 , 1036 salpeter , e.e .
1961 , , 134 , 669 saumon , d. , chabrier , g. , & van horn , h.m .
1995 , , 99 , 713 saumon , d. , & jacobson , s.b .
1999 , , 511 , l107 segretain , l. , chabrier , g. , hernanz , m. , garcia - berro , e. , isern , j. , & mochkovitch , r. 1994 , , 434 , 641 stolzmann , w. , & blocker , t. 2000 , a&a , 361 , 115 straniero , o. 1988 , a&a , 76 , 157 straniero , o. , chieffi , a. , & limongi , m. 1997 , , 490 , 425 straniero , o. , chieffi , a. , limongi , m. , busso , m. , gallino , r. , & arlandini , c. 1997 , , 478 , 332 tanaka , s. , mitake , s. , & ichimaru , s. 1985 , phys .
a , 32 , 1896 tassoul , m. , fontaine , g. , & winget , d.e .
1990 , , 72 , 335 von hippel , t. , gilmore , g. , & jones , d.h.p .
1995 , , 273l , 39 von hippel , t. , & gilmore , g. 2000 , , 120 , 1384 weidemann , v. 1987 , a&a , 188 , 74 weidemann , v. 2000 , a&a , 363 , 647 weidemann , v. , & koester , d. 1983 , a&a , 121 , 77 wood , m.a .
1995 , in `` white dwarfs '' , springer verlag , ed : d. koester & k. werner , 41 yakovlev , d. , & shalybkov , d. 1989 , astron .
space phys .
, 7 , 311    cccccccccc b & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & @xmath93 + b1 & 0.600 & 0.0001 & 0.580 & 0.600 & @xmath94 & @xmath95 & 0.258 & 0.542 & @xmath96 + b2 & 0.640 & 0.0001 & 0.515 & 0.601 & @xmath97 & @xmath98 & 0.249 & 0.516 & @xmath99 + b3 & 0.600 & 0.0200 & 0.599 & 0.600 & @xmath100 & @xmath101 & 0.248 & 0.509 & @xmath102 + b - lh & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath103 & 0.263 & 0.573 & @xmath104 + b - vlh & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath105 & 0.263 & 0.573 & @xmath106 +
b - bp & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath107 & @xmath108 & 0.151 & 0.316 & @xmath109 + b - low & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath110 & @xmath111 & 0.511 & 0.483 & @xmath112 + b - low - bp & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath113 & @xmath114 & 0.381 & 0.434 & @xmath115 + b - hl69 & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & @xmath116 +
b - comb & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & @xmath117 + b - pot & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & @xmath118 + b - noie & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & @xmath119 + b - nodif & 0.600 & 0.0001 & 0.599 & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & @xmath120 +    ccccccc a & 0.500 & @xmath121 & @xmath122 & 0.221 & 0.571 & 15.00 + b & 0.600 & @xmath91 & @xmath92 & 0.263 & 0.573 & 14.20 + c & 0.700 & @xmath123 & @xmath124 & 0.294 & 0.560 & 12.90 + d & 0.800 & @xmath125 & @xmath126 & 0.302 & 0.529 & 12.31 + e & 0.900 & @xmath127 & @xmath128 & 0.309 & 0.492 & 11.32 +