English

Implementation of New OPAL Tables in Eggleton's Stellar Evolution Code

Astrophysics 2009-11-11 v2

Abstract

Based on the work of OPAL (Iglesias & Rogers 1996) and Alexander & Ferguson (1994), we construct a series of opacity tables for various metallicities Z=0, 0.00001, 0.00003, 0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08 and 0.1. These tables can be easily used in Eggleton's stellar evolution code in place of the old tables without changing the code. The OPAL tables are used for log10(T/K)>3.95{\log}_{10}(T/{\rm K})> 3.95 and Alexander's for log10(T/K)<3.95{\log}_{10}(T/{\rm K})< 3.95. At log10(T/K)=3.95{\rm log}_{10}(T/{\rm K})=3.95, the two groups' data fit well for all hydrogen mass fractions. Conductive opacities are included by reciprocal addition according to the formulae of Yakovlev & Urpin (1980). For different metallicities 1M1M_\odot models are constructed with the new opacity tables. Comparison of 1 and 5M5M\odot models constructed with the older OPAL tables (Iglesias & Rogers 1992) shows that the new opacities have most effect in the late stages of evolution, the extension of the blue loop during helium burning for intermediate-mass and massive stars.

Cite

@article{arxiv.astro-ph/0610305,
  title  = {Implementation of New OPAL Tables in Eggleton's Stellar Evolution Code},
  author = {Xuefei Chen and Christopher A. Tout},
  journal= {arXiv preprint arXiv:astro-ph/0610305},
  year   = {2009}
}

Comments

8 pages, 5 figures