Tables and their content

MetalFile

Mass of different elements, given in solar mass, produced by a star of mass M.

example, from Tsujimoto 95, Table 1 : tables/MetalEjection.dat:

# M     Fe56    Mg24    O16     Metals
10      0.0     0.0     0.0     0.0
13      1.50e-1 9.23e-3 1.51e-1 4.51e-1
15      1.44e-1 3.16e-2 3.55e-1 8.20e-1
18      7.57e-2 3.62e-2 7.92e-1 1.45
20      7.32e-2 1.47e-1 1.48    2.28
25      5.24e-2 1.59e-1 2.99    4.28
40      7.50e-2 3.54e-1 9.11    11.47
70      7.50e-2 7.87e-1 21.4    26.71

During the table construction, content of the table is processed by:

SNIIelts = libmetalsSNII.WriteParams(f,Mmin,Mmax,m_s,a_s,MetalFile,HeliumCoreFile,elts,nelts,Mmin1,Mmax1,Mmin2,Mmax2,dM,n)

Values initially in mass are converted in mass fraction. Values are then written in the block labeled #### Metal Parameters ####.

In the C/API, this table is read by:

Elts[it][i].label                     # label of elements
Elts[it][i].Mmin Elts[it][i].Step     # min mass, steps in mass
Elts[it][i].Metal[j]                  # mass fraction of element as a function of mass

it=i+2 for element i.

To retrieve the mass fraction of ejected elements for a 20 solar mass star, from python:

from PyChem import chemistry
chemistry.init_chimie("../chemistry.simple.dat")
chemistry.SNII_single_mass_ejection(20*chemistry.SolarMass_to_CodeUnits_Factor())[2:]

This recover the values of the original table:

20*chemistry.SNII_single_mass_ejection(20*chemistry.SolarMass_to_CodeUnits_Factor())[2:]

To display the tables (in mass fraction):

../scripts/pychem_plot_SNII_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 chemistry.simple.dat

HeliumCoreFile

Mass of the remnent object (Residu, \omega(m)) and helium core (\alpha(m)) mass, both in solar mass, a star of mass m.

example (from Kobayashi 2000, private comm.) : tables/HeliumCore.dat:

# Mass  Residu  Helium Core
0.80    0.459   0.459
1.00    0.473   0.473
1.50    0.480   0.480
1.75    0.487   0.487
2.00    0.404   0.404
2.50    0.501   0.501
3.00    0.508   0.508
3.50    0.515   0.515
4.00    0.550   0.550
4.30    0.620   0.620
4.60    0.675   0.675
4.80    0.730   0.730
5.00    0.870   0.870
6.00    1.010   1.010
8.00    1.120   1.120
10.0    1.150   1.150
13.0    1.500   3.300
15.0    1.500   4.000
18.0    1.500   5.000
20.0    1.600   6.000
25.0    2.000   8.000
40.0    2.300   16.00
70.0    2.500   32.00

During the table construction, content of the table is processed by:

SNIIelts = libmetalsSNII.WriteParams(f,Mmin,Mmax,m_s,a_s,MetalFile,HeliumCoreFile,elts,nelts,Mmin1,Mmax1,Mmin2,Mmax2,dM,n)

Values initially in mass are converted in mass fraction. In addition, we keep 1-\omega(m) and 1-\alpha(m), instead of \omega(m) \alpha(m). Values are then written in the block labeled #### Metal Parameters ####.

In the C/API, this table is read by:

Elts[it][i].label                     # label of elements
Elts[it][i].Mmin Elts[it][i].Step     # min mass, steps in mass
Elts[it][i].Metal[j]                  # mass fraction of element as a function of mass

it=0 for 1-\omega(m) (labeled Ej) and it=1 for 1-\alpha(m) (labeled Ejnp).

To retrieve the mass fraction of the ejected mass (1-\omega(m)) and the mass fraction of the non processed mass (1-\alpha(m)):

from PyChem import chemistry
chemistry.init_chimie("../chemistry.simple.dat")
chemistry.SNII_single_mass_ejection(20*chemistry.SolarMass_to_CodeUnits_Factor())[:2]

This recover the values of the original table:

20*(1-chemistry.SNII_single_mass_ejection(20*chemistry.SolarMass_to_CodeUnits_Factor())[:2])

To display the tables (in mass fraction):

../scripts/pychem_plot_SNII_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 chemistry.simple.dat

SNIaFile

Mass of different elements, in solar mass, ejected during a SNIa explosion involving a 1.38 white dwarf.

example (from Iwamoto et al 99) : tables/SNIa/SNIaEjection.dat:

# Fe56  Mg24    O16     Metals
6.13e-1 8.5e-3  1.43e-1 1.3743

During the table construction, content of the table is processed by:

SNIaelts = libmetalsSNIa.WriteParams(f,SNIaFile,MeanWDMass,elts,nelts)

Contrary to the SNII tables, here, the elements are not converted in mass fraction. Values are then still in solar mass. Their are then written in the block labeled # SNIa Metal Ejection.

In the C/API, this table is read by:

Elts[it][i].MSNIa                     # mass of element ejected by a SNIa

it=i+2 for element i. it=0 and it=1 are not used.

To retrieve the ejected masses in python:

from PyChem import chemistry
chemistry.init_chimie("../chemistry.simple.dat")
chemistry.SNIa_single_mass_ejection(20*chemistry.SolarMass_to_CodeUnits_Factor())[2:] * chemistry.CodeUnits_to_SolarMass_Factor()

This recover directly the values of the original table. Note that by default the output of SNIa_single_mass_ejection is in code units.

SolarAbundancesFile

example : tables/SolarAbundances/Grevesse98.txt:

# Z     Elt     X       A          Flag
01      H       12.00   1.00797    0
02      He      10.93   4.0026     0
03      Li       1.10   6.941      0
04      Be       1.40   9.0122     0
05      B        2.55   10.811     0
06      C        8.52   12.01115   0
07      N        7.92   14.0067    0
08      O        8.83   15.9994    0
09      F        4.56   18.9984    0
10      Ne       8.08   20.183     0
11      Na       6.33   22.9898    0
12      Mg       7.58   24.312     0
13      Al       6.47   26.9815    0
14      Si       7.55   28.086     0

During the table construction, content of the table is processed by:

libSolarAbundances.WriteParams(f,SolarAbundancesFile,elts)