arXiv:astro-ph/0001295v1 17 Jan 2000 oa ooecm oso os.I esynthesize we stars If of noise. sample Poisson al. large overcome et a statistically to evolve Eck synthesis as to Van so population necessary 1997; is make al. it to et meaningful order Boer In evoluti de stellar 1998). 1997; of aspects (Baglin many theory of re-evaluation can a techniques an initiate observing is improved how Catalogue of (Hoffleit Hipparcos example 1997), excellent The al. Catalogue synthesis. et Star large-scale (Perryman population Catalogue Bright of using Hipparcos findings the the and the as 1983) reproduce such do to to surveys, way attempt star One to observations. is against this tested be must lution in (1995). problems al. of et overview Noels results recent see the evolution a update For stellar to need calculations. pulses, ongoing these an of (the thermal is improved and there being example) overshooting continually for convective is a of theory, it, evolution treatment model stellar to As ability out collisions. our the stellar simulating and of clusters comes star Ex- of galaxie astrophysics. ages of the evolution of determining chemical areas the modelling many include are in amples calculations applications evolution for stellar quired detailed of results The INTRODUCTION 1 lso atclrtp fsa ehv nerrof error an have we star of type particular a of ples 8My2018 May 28 o.Nt .Ato.Soc. Astron. R. Not. Mon. 1 E-mail: 2 ardR Hurley R. Jarrod metallicity a and as mass evolution of stellar function for formulae analytic Comprehensive oee,dtie vlto oe a aesvrlhusto hours several take possible can codes sample of accurate evolution millions sufficently detailed a often However, get stars to required rarer are for progenitors that means which c nttt fAtooy aige od abig B 0HA, CB3 Cambridge Road, Madingley Astronomy, of Institute nttt eAto´sc eCnra,c i ´ce /,E L´actea s/n, Via c/ Canarias, Astrof´ısica de de Instituto 99RAS 1999 swt l hoisorudrtnigo tla evo- stellar of understanding our theories all with As [email protected],on@lice,[email protected] [email protected], [email protected], 000 1 noR Pols R. Onno , –9(99 rne 8My21 M L (MN 2018 May 28 Printed (1999) 1–29 , tr:Pplto I–glxe:selrcontent stellar galaxies: – II Population stars: asls rsrpinta a eue ihteefrua n inv and formulae these with words: used Key be can distribution. that remnant stellar prescription resulting loss mass a rpris n npriua o s ncmiainwith useful combination are in use formulae accurate for but These particular rapid very in models. co require and detailed that find properties, synthesis we of population part 5% as most such the within For to stages. remnant accurate the including, and to, safnto fage, of function a as ABSTRACT ag fmass of range epeetaayi omleta prxmt h vlto fsta of evolution the approximate that formulae analytic present We M ehd:aayi tr:eouin–sas udmna param fundamental stars: – evolution stars: – analytic methods: 2 n n metallicity and n hitpe .Tout A. Christopher and exam- ± 320L aua eeie Spain Tenerife, Laguna, La -38200 √ M re- nd on .uk s, n - . and UK o-nfr oeeto tri h Hertzsprung-Russell the in highly star the ap- a for analytic of nature movement find simple non-uniform to conveniently practice et a in Tout of difficult of proximations more lines is the It of along (1996). work (1989) al. single the Tout re- of & expanding set Fitchett the thus Eggleton, a on formulae, of (SSE) report construction evolution we metallicity approach, star paper and second present the mass the of of evolu- sults In stellar range form. of wide tabular results pos- a in the paper of for presented range previous calculations we full a tion 1998) the In al. cover masses. et that stellar (Pols models and of compositions idea set the sible a so have Stel- metallicity of 1998). to range al. is et wide et a (Schaller Mowlavi cover 1993; years populations al. many lar et for have Charbonnel form 1992; models al. tabular Stellar in simultaneously. available both we been so on procedures 1996), worked Both (Eggleton metallicity. have disadvantages and and mass advantages have age, formulae of interpolation functions of t number as is a incor- by other data The be the tables. to approximate these within is ra interpolate parameter a and overshoot porated) if and/or especially as con- metallicities large, runs to of rather is stellar-evolution approach (necessarily of One tables them. series struct require a that projects of for output data the using of t stage. utilise some to later in simple a them relatively at present is results and it gener- which models to in detailed desirable form of is convenient it set Thus large star. a one ate just of model a evolve Z o l hssfo h eoaemi-euneup main-sequence zero-age the from phases all for , Z hr r w lentv prahst h problem the to approaches alternative two are There tla uioiy aisadcr asaegiven are mass core and radius luminosity, Stellar . A T E tl l v1.4) file style X 1 N bd oe.W describe We codes. -body vlaino stellar of evaluation tnosformulae ntinuous sfrawide a for rs o purposes for siaethe estigate tr – eters nge he o l 2 J. R. Hurley, O. R. Pols and C. A. Tout diagram (HRD) than it is to interpolate in tables, but the ing as they evolve off the MS but for higher mass stars with resulting code is very much more compact and adaptable to convective cores the transition is not so smooth. Owing to the requirements of, for example, an N-body code (Aarseth mixing in the core there is a sudden depletion of fuel over 1996) or variable mass loss. This is reinforced in the circum- a large region which leads to a rapid contraction over the stance where one wishes to include binary-star interactions, inner region at core hydrogen-exhaustion. This causes the such as Roche-lobe overflow, common-envelope evolution, hydrogen-exhausted phase gap, or MS hook, which occurs and magnetic braking with tidal friction, for example (Tout on a thermal timescale. The different features of MS evo- et al. 1997). lution are illustrated by comparing the evolution tracks for In Section 2 we provide a brief overview of how stars the 1.0 M⊙ and 1.6 M⊙ stars in Figure 1. behave as they evolve in time which introduces some of the The immediate post-MS evolution towards the right in terminology that we use and will hopefully facilitate the un- the HRD occurs at nearly constant luminosity and is very derstanding of this paper. Section 3 describes the detailed rapid. For this reason very few stars are seen in this phase, models from which the formulae are derived and justifies and this region of the HRD is called the Hertzsprung gap the inclusion of enhanced mixing processes. In Section 4 we (HG), or the sub-giant branch. During this HG phase the ra- outline the procedure to be used for generating the SSE dius of the star increases greatly causing a decrease in Teff . package. The evolution formulae are presented in Section 5 For cool envelope temperatures the opacity increases causing for all nuclear burning phases from the main-sequence to the a convective envelope to develop. As the convective envelope asymptotic giant branch. Our formulae are a vast improve- grows in extent the star will reach the giant branch (GB) ment on the work of Eggleton, Fitchett & Tout (1989) not which is the nearly vertical line corresponding to a fully only due to the inclusion of metallicity as a free parameter convective star, also known as the Hayashi track. All stars but also because we have taken a great deal of effort to pro- ascend the GB with the hydrogen-exhausted core contract- vide a more detailed and accurate treatment of all phases of ing slowly in radius and heating while the hydrogen-burning the evolution. Features such as main-sequence formulae that shell is eating its way outwards in mass and leaving behind are continuous over the entire mass range and the modelling helium to add to the growing core. As the stars move up of second dredge-up and thermal pulses will be discussed. the GB convection extends over an increasing portion of the Section 6 discusses the behaviour of a star as the stellar en- star. The convective envelope may even reach into the previ- velope becomes small in mass and outlines what happens ously burnt (or processed) regions so that burning products when the nuclear evolution is terminated. We also provide are mixed to the surface in a process called dredge-up. formulae which model the subsequent remnant phases of Eventually a point is reached on the GB where the core evolution. In Section 7 we describe a comprehensive mass temperature is high enough for stars to ignite their central loss algorithm which can be used in conjunction with the helium supply. For massive stars, M > 2.0 M⊙, this takes evolution formulae, as well as a method for modelling stellar place gently. When core helium burning∼ (CHeB) begins the rotation. Various uses for the formulae and future improve- star descends along the GB until contraction moves the star ments are discussed in Section 8 along with details of how away from the fully convective region of the HRD and back to obtain the formulae in convenient subroutine form. towards the MS in what is called a blue loop. During CHeB, carbon and oxygen are produced in the core. Eventually core helium is exhausted and the star moves back to the right in the HRD. The size of the blue loop generally increase with 2 STELLAR EVOLUTION OVERVIEW mass, as can be seen by comparing the 4.0 M⊙ and 10.0 M⊙ A fundamental tool in understanding stellar evolution is the tracks in Figure 1. Lower mass stars have degenerate helium Hertzsprung-Russell diagram (HRD) which provides a cor- cores on the GB leading to an abrupt core-helium flash at relation between the observable stellar properties of lumi- helium ignition (HeI). The star then moves down to the nosity, L, and effective surface temperature, Teff . Figure 1 zero-age horizontal branch (ZAHB) very quickly. The initial shows the evolution of a selection of stars in the HRD from position of a star along the ZAHB depends on the mass of the zero-age main-sequence (ZAMS), where a star adjusts the hydrogen-exhausted core at the time of ignition and on itself to nuclear burning equilibrium, until the end of their the mass in the overlying envelope. Those stars with lower nuclear burning lifetimes. As stars take a relatively short mass, ie. shallower envelopes, appear bluer because there is time to reach the ZAMS all ages are measured from this less mass to shield the hot hydrogen burning shell. It is also point. The length of a stars life, its path on the HRD and possible for stars of very high mass, M > 12.0 M⊙, to reach its ultimate fate depend critically on its mass. high enough central temperatures on the∼ HG for helium to Stars spend most of their time on or near the main- ignite before reaching the GB. The 16.0 M⊙ star in Figure 1 sequence (MS) burning hydrogen to produce helium in their is such an example. As a result these stars by-pass the GB cores. To first order, the behaviour of a star on the MS can phase of evolution. be linked to whether it has a radiative or convective core. Evolution after the exhaustion of core-helium is very Stars with M < 1.1 M⊙ have radiative cores while in higher similar to evolution after core-hydrogen exhaustion at the mass stars a convective∼ core develops as a result of the steep end of the MS. The convective envelope deepens again so temperature gradient in the interior. During core hydrogen that the star once more moves across towards the Hayashi burning on the MS, low-mass stars will move upwards in L track to begin what is called the asymptotic giant branch and to higher Teff on the HRD while higher mass stars will (AGB). On the AGB the star consists of a dense core com- also move upwards in L but to a region of lower Teff . The MS posed of carbon and oxygen surrounded by a helium burn- evolution will end when the star has exhausted its supply of ing shell which adds carbon to the degenerate core. Initially hydrogen in the core. Low-mass stars will continue expand- the H-burning shell is extinguished so that the luminosity is
c 1999 RAS, MNRAS 000, 1–29 Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity 3 supplied exclusively by the He-burning shell; characterizing large enough that complete collapse to a black hole (BH) the early AGB (EAGB) phase. If the star is massive enough occurs. the convective envelope can reach into the H-exhausted re- Stars with M > 15 M⊙ are severely affected by mass gion again (second dredge-up). When the He-burning shell loss during their entire∼ evolution and may lose their en- catches up with the H-rich envelope the H-shell reignites velopes during CHeB, or even on the HG, exposing nuclear and the two grow together with the H-burning shell sup- processed material. If this occurs then a naked helium star plying most of the luminosity. During the following phase is produced and such stars, or stars about to become naked the helium shell is unstable, which can cause a helium shell helium stars, may be Wolf-Rayet stars. Wolf-Rayet stars are flash in which the helium shell will suddenly release a large massive objects which are found near the MS, are losing amount of luminosity. The energy released in the flash ex- mass at very high rates and show weak, or no, hydrogen pands the star resulting in the hydrogen shell cooling so lines in their spectra. Luminous Blue Variables (LBVs) are much that it is extinguished. Convection once again reaches extremely massive post-MS objects with enormous mass loss downward past the dead hydrogen shell. This mixes helium rates in a stage of evolution just prior to becoming a Wolf- to the surface, as well as carbon that was mixed out of the Rayet star. Naked helium stars can also be produced from helium shell by flash-driven convection. As the star subse- less massive stars in binaries as a consequence of mass trans- quently contracts the convection recedes and the hydrogen fer. shell re-ignites but has now moved inwards in mass due to Variations in composition can also affect the stellar evo- the envelope convection. This process is called third dredge- lution timescales as well as the appearance of the evolution up. The star continues its evolution up the AGB with the on the HRD, and even the ultimate fate of the star. A more hydrogen shell producing almost all of the luminosity. The detailed discussion of the various phases of evolution can be helium shell flash can repeat itself many times and the cycle found throughout this paper. is known as a thermal pulse. This is the thermally pulsing asymptotic giant branch (TPAGB). The stellar radius can grow to very large values on the 3 STELLAR MODELS AGB which lowers the surface gravity of the star, so that the The fitting formulae are based on the stellar models com- surface material is less tightly bound. Thus mass loss from puted by Pols et al. (1998). They computed a grid of evo- the stellar surface can become significant with the rate of lution tracks for masses M between 0.5 and 50 M⊙ and mass loss actually accelerating with time during continued for seven values of metallicity, Z = 0.0001, 0.0003, 0.001, evolution up the AGB. Unfortunately, our understanding of 0.004, 0.01, 0.02 and 0.03. They also considered the prob- the mechanisms that cause this mass loss is poor with pos- lem of enhanced mixing such as overshooting beyond the sible suggestions linking it to the helium shell flashes or to classical boundary of convective instability. Its effect was periodic envelope pulsations. Whatever the cause, the influ- modelled with a prescription based on a modification of the ence on the evolution of AGB stars is significant. Mass loss Schwarzschild stability criterion, introducing a free param- will eventually remove all of the stars envelope so that the eter δ (which differs from the more commonly used pa- hydrogen burning shell shines through. The star then leaves ov rameter relating the overshooting distance to the pressure the AGB and evolves to hotter Teff at nearly constant lumi- scale height; see Pols et al. 1998 for details). The tracks com- nosity. As the photosphere gets hotter the energetic photons puted with a moderate amount of enhanced mixing (given become absorbed by the material which was thrown off while by δ = 0.12 and labeled the OVS tracks by Pols et al. on the AGB. This causes the material to radiate and the star ov 1998) were found to best reproduce observations in a series may be seen as a planetary nebula. The core of the star then of sensitive tests involving open clusters and ecliping bina- begins to fade as the nuclear burning ceases. The star is now ries (see Schr¨oder, Pols & Eggleton 1997; Pols et al. 1997, a white dwarf (WD) and cools slowly at high temperature 1998). We consequently use these OVS tracks as the data to but low luminosity. which we fit our formulae. If the mass of the star is large enough, M > 7 M⊙, the For each Z, 25 tracks were computed spaced by approx- carbon-oxygen core is not degenerate and will∼ ignite car- imately 0.1 in log M, except between 0.8 and 2.0 M⊙where bon as it contracts, followed by a succession of nuclear re- four extra models were added to resolve the shape of the action sequences which very quickly produce an inner iron main-sequence which changes rapidly in this mass range. core. Any further reactions are endothermic and cannot con- Hence we dispose of a database of 175 evolution tracks, each tribute to the luminosity of the star. Photodisintegration of containing several thousand individual models. iron, combined with electron capture by protons and heavy A subset of the resulting OVS tracks in the HRD are nuclei, then removes most of the electron degeneracy pres- shown in Fig. 1 for Z = 0.02 and Fig 2 for Z = 0.001. sure that was supporting the core and it begins to collapse The considerable variation of model behaviour introduced rapidly. When the density becomes large enough the inner by changes in metallicity is illustrated by Fig. 3. Detailed core rebounds sending a shockwave outwards through the models of the same mass, M = 6.35 M⊙, are shown on the outer layers of the star that have remained suspended above HRD for three different metallicities, Z = 0.0001, 0.001 and the collapsing core. As a result the envelope of the star is 0.02. Not only does a change in composition move the track ejected in a supernova (SN) explosion so that the AGB is to a different position in the HRD but it also changes the effectively truncated at the start of carbon burning and the appearance of each track, as can be seen by considering the star has no TPAGB phase. The remnant in the inner core extent of the hook feature towards the end of the main se- will stablise to form a neutron star (NS) supported by neu- quence and the blue loops during core helium burning. Fur- tron degeneracy pressure unless the initial stellar mass is thermore, the Z = 0.0001 model ignites helium in its core
c 1999 RAS, MNRAS 000, 1–29 4 J. R. Hurley, O. R. Pols and C. A. Tout
Figure 1. Selected OVS evolution tracks for Z = 0.02, for masses Figure 2. Same as Fig. 1 for Z = 0.001. The 1.0 M⊙ post He 0.64, 1.0, 1.6, 2.5, 4.0, 6.35, 10, 16, 25 and 40 M⊙. flash track has been omitted for clarity. while on the Hertzsprung gap as opposed to the other mod- 6 = Thermally Pulsing Asymptotic Giant Branch (TPAGB) els which evolve up the giant branch before reaching a high 7 = Naked Helium Star MS (HeMS) enough core temperature to start helium burning. In addi- tion the nuclear burning lifetime of a star can change by as 8 = Naked Helium Star Hertzsprung Gap (HeHG) much as a factor of 2 owing to differences in composition, 9 = Naked Helium Star Giant Branch (HeGB) as shown in Fig. 4 for a set of 2.5 M⊙ models. This em- 10 = Helium White Dwarf (He WD) phasizes the need to present the results of stellar evolution 11 = Carbon/Oxygen White Dwarf (CO WD) calculations for an extensive range of metallicity. Mass loss from stellar winds was neglected in the de- 12 = Oxygen/Neon White Dwarf (ONe WD) tailed stellar models, mainly because the mass loss rates are 13 = Neutron Star(NS) uncertain by at least a factor of three. We do include mass 14 = Black Hole (BH) loss in our analytic formulae in an elegant way, as will be de- scribed in Section 7.1, which allows us to experiment easily 15 = massless remnant, with different mass loss rates and prescriptions. and we divide the MS into two phases to distinguish be- tween deeply or fully convective low-mass stars and stars of higher mass with little or no convective envelope as these will respond differently to mass loss. 4 PROCEDURE To begin with we take different features of the evolution We assign each evolution phase an integer type, k, where: in turn, e.g. MS lifetime, ZAHB luminosity, and first try to fit them as f (M) for a particular Z in order to get an idea of < 0 = MS star M 0.7 deeply or fully convective the functional form. We then extend the function to g (M, Z) ∼ 1 = MS star M > 0.7 using f (M) as a starting point. In this way we fit formulae ∼ 2 = Hertzsprung Gap (HG) to the end-points of the various evolutionary phases as well as to the timescales. We then fit the behaviour within each 3 = First Giant Branch (GB) phase as h (t,M,Z), e.g. LMS (t,M,Z). 4 = Core Helium Burning (CHeB) As a starting point we take the work of Tout et 5 = Early Asymptotic Giant Branch (EAGB) al. (1996) who fitted the zero-age main-sequence luminos-
c 1999 RAS, MNRAS 000, 1–29 Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity 5
ity (LZAMS) and radius (RZAMS) as a function of M and Z. Their aim, as is ours, was to find simple computationally efficent functions which are accurate, continuous and differ- entiable in M and Z, such as rational polynomials. This is acheived using least-squares fitting to the data after choos- ing the initial functional form. In most cases we determine the type of function, the value of the powers, and the num- ber of coefficients to be used, simply by inspecting the shape of the data, however in some cases, such as the luminosity- core-mass relation on the giant branch, the choice will be dictated by an underlying physical process. For the ZAMS, accuracy is very important because it fixes the star’s posi- tion in the HRD. Tout et al. (1996) acheived LZAMS accurate to 3% and RZAMS accurate to 1.2% over the entire range. For the remainder of the functions we aim for RMS errors less than 5% and preferably a maximum individual error less than 5% although this has to be relaxed for some later stages of the evolution where the behaviour varies greatly with Z but also where the model points are more uncertain owing to shortcomings in stellar evolution theory.
5 FITTING FORMULAE In this section we present our formulae describing the evo- lution as a function of mass M and age t. The explicit Z- dependence is in most cases not given here because it would clutter up the presentation. This Z-dependence is implicit whenever a coefficient of the form an or bn appears in any of the formulae. The explicit dependence of these coefficients on Z is given in the Appendix. Coefficients of the form cn, whose numerical values are given in this section, do not de- Figure 3. Detailed OVS evolution tracks for M = 6.35 M⊙, for metallicities 0.0001, 0.001 and 0.02. pend on Z. We adopt the following unit conventions: numerical val- ues of mass, luminosity and radius are in solar units, and values of timescales and ages are in units of 106 yr, unless otherwise specified. We begin by giving formulae for the most important critical masses, Mhook (the initial mass above which a hook appears in the main-sequence), MHeF (the maximum initial mass for which He ignites degenerately in a helium flash) and MFGB (the maximum initial mass for which He ignites on the first giant branch). Values for these masses are given in Table 1 of Pols et al. (1998) estimated from the detailed models for 7 metallicities. These values can be accurately fitted as a function of Z by the following formulae, where ζ = log(Z/0.02): 2 Mhook = 1.0185 + 0.16015ζ + 0.0892ζ , (1) 2 MHeF = 1.995 + 0.25ζ + 0.087ζ , (2) 13.048 (Z/0.02)0.06 MFGB = . (3) 1 + 0.0012 (0.02/Z)1.27 Based on the last two critical masses, we make a distinction into three mass intervals, which will be useful in the later descriptions:
(i) low-mass (LM) stars, with M < MHeF, develope de- Figure 4. Radius evolution as a function of stellar age for generate He cores on the GB and ignite He in a degenerate M = 2.5M⊙, for metallicities 0.0001, 0.001 and 0.02. Tracks are flash at the top of the GB; from the detailed models and run from the ZAMS to the point of (ii) intermediate-mass (IM) stars, with MHeF M < termination on the AGB. ≤ MFGB, which evolve to the GB without developing degener- ate He cores, also igniting He at the top of the GB;
c 1999 RAS, MNRAS 000, 1–29 6 J. R. Hurley, O. R. Pols and C. A. Tout
(iii) high-mass (HM) stars, with M > MFGB, ignite He in the HG before the GB is reached, and consequently do not have a GB phase. Note that this definition of IM and HM stars is different from the more often used one, based on whether or not carbon ignites non-degenerately.
5.1 Main-sequence and Hertzsprung gap To determine the base of the giant branch (BGB) we find where the mass of the convective envelope MCE first exceeds a set fraction of the envelope mass ME as MCE increases on the HG. From inspection the following fractions 2 MCE = ME M < MHeF 5 ∼ Figure 5. Time taken to reach the base of the giant branch as 1 > MCE = ME M MHeF a function of stellar mass as given by eq. (4), shown against the 3 ∼ detailed model points, for Z = 0.0001 and 0.03 which give the generally give a BGB point corresponding to the local mini- worst fit of all the metallicities. The maximum error over the mum in luminosity at the start of the GB. We define helium entire metallicity range is 4.8% and the RMS error is 1.9%. ignition as the point where LHe = 0.01L for the first time. For HM stars this will occur before the BGB point is found, ie. no GB, and thus we set tBGB = tHeI for the sake of defin- RGB (LBGB) M < MFGB REHG = . ing an end-point to the HG, so that BGB is more correctly RHeI M MFGB the end of the HG (EHG) as this is true over the entire mass ≥ range. The luminosity at the end of the MS is approximated The resultant lifetimes to the BGB are fitted as a func- by 3 4 a16+1.8 tion of M and Z by a11M + a12M + a13M LTMS = (8) 4 5.5 7 5 a16 a1 + a2M + a3M + M a14 + a15M + M tBGB = 2 7 . (4) a4M + a5M with a16 7.2. This proved fairly straightforward to fit but ≈ Figure 5 shows how eq. (4) fits the detailed model points for the behaviour of RTMS is not so smooth and thus requires Z = 0.0001 and 0.03 which are the metallicities which lead a more complicated function in order to fit it continuously. to the largest errors. Over the entire metallicity range the The resulting fit is function gives a rms error of 1.9% and a maximum error of a21 a18 + a19M 4.8%. In order that the time spent on the HG will always be RTMS = M a17 (9) a20 + Ma22 ≤ a small fraction of the time taken to reach the BGB, even 3 a26 a26+1.5 c1M + a23M + a24M for low-mass stars which don’t have a well defined HG, the ∗ RTMS = 5 M M , (9a) MS lifetimes are taken to be a25 + M ≥ with straight-line interpolation to connect eqs. (9) and (9a) t = max(t ,xt ) , (5) MS hook BGB between the endpoints, where where thook = µtBGB and M∗ = a17 + 0.1 , a17 1.4 ≈ x = max (0.95, min (0.95 0.03 (ζ + 0.30103) , 0.99)) (6) −2 − and c1 = 8.672073 10 , a21 1.47, a22 3.07, a26 − × ≈ ≈ ≈ a6 a9 5.50. Note that for low masses, M < 0.5, where the function µ = max 0.5, 1.0 0.01 max ,a8 + . (7) − Ma7 Ma10 is being extrapolated we add the condition Note that µ is ineffective for M < Mhook, ie. stars without RTMS = max(RTMS, 1.5RZAMS) a hook feature, and in this case the functions ensure that x > µ. to avoid possible trouble in the distant future. So we now have defined the time at the end of the MS, The luminosity at the base of the GB is approximated tMS and the time taken to reach the start of the GB (or end by of the HG), tBGB such that a31 c2 a27M + a28M LBGB = 3 32 , (10) t : 0.0 tMS MS evolution a29 + a30Mc + Ma → t : tMS tBGB HG evolution. with c2 = 9.301992, c3 = 4.637345, a31 4.60 and a32 → ≈ ≈ 6.68. The description of L , R and R is given in The starting values for L and R are the ZAMS points HeI GB HeI later sections. fitted by Tout et al. (1996). We fit the values at the end of the MS, LTMS and RTMS, as well as at the end of the HG,
LBGB M < MFGB 5.1.1 Main-sequence evolution LEHG = LHeI M MFGB ≥ On the MS we define a fractional timescale
c 1999 RAS, MNRAS 000, 1–29 Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity 7
Lb(t) for a typical detailed model is illustrated by Fig. 6. The luminosity pertubation is approximated by
0.0 M Mhook ≤ 0.4 M Mhook ∆L = B − Mhook 0.0 M Mhook ≤ 0.5 M Mhook a43 − Mhook < M a42 a42 Mhook ≤ − ∆R = (17) a44 M a42 a43 +(B a43) − a42