Re: [mesa-users] ionization fraction or kap table ?

[RICHARD H D T. <tow...@as...>]

Hi Jing —

> On Jun 24, 2017, at 9:55 AM, Jing Luan <jin...@gm...> wrote:
> 
> Dear mesa users,
> 
> 1, Are there parameters in mesa that tells the ionization fraction of (HI to HII) and of (HeI to HeII) please? I search around mesa/star/public/star_data.inc, but did not find any. But maybe I miss something there. 
> 

I would suspect not. Nowadays, with tabulated EOS’s, it is unusual to actually solve the Saha equation during a stellar evolution calculation. However, MESA *can* tell you the free electron number density, and so you could calculate the ionization fractions yourself, using the Saha equation (for this, I suggest you ignore excited states of the H/He).

> 2, For calculating kappa mechanism due to H or He opacity bumps for driving stellar oscillation, which kappa table and eos table should be adopted please? The default in mesa-r8118 are 
> 
> kappa_file_prefix = 'gn93'
> 
> eos_file_prefix = 'mesa'
> eosDT_Z1_suffix = '_CO_1'
> eosPT_Z1_suffix = ‘_
> CO_1'
> 

For the blue edges of the classical (Cepheid) instability strip, plotted in Fig. 14 of Paxton et al (2015), I used:

initial_zfracs = 6 ! Asplund et al. (2009)
kappa_file_prefix = 'a09' ! Asplund et al. (2009)

…in the &star_job namelist, and

initial_Y = 0.2703 ! Asplund et al. (2009) / proto-solar
initial_Z = 0.0142 ! Asplund et al. (2009) / proto-solar

…in the &controls namelist. I expect the metallicity is less important to driving classical Cepheid pulsations than the relative H/He abundances.

cheers,

Rich

> Many thanks 😊😊😊😊
> 
> Sincerely,
> Jing
> ------------------------------------------------------------------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, Slashdot.org! http://sdm.link/slashdot_______________________________________________
> mesa-users mailing list
> mes...@li...
> https://lists.sourceforge.net/lists/listinfo/mesa-users