File:Hr pre main sequence 1msun 1.svg

Pre-Main Sequence models

Tognelli, Prada Moroni, Degl'Innocenti 2011, Astronomy & Astrophysics, 533, A109

http://astro.df.unipi.it/stellar-models/index.php?m=3

References for database

References

   Asplund, M., Grevesse, N., Sauval, A.J., & Scott, P. 2009, ARA&A, 47, 481
   Hammer, J.W., Fey, M., Kunz, R. et al. 2005, Nuclear Physics A, 758, 363

Dataset

Z=0.01377 Y=0.2533

mixing-length= 1.68 Hp XD= 2 · 10-5 Tracks

http://astro.df.unipi.it/stellar-models/pre-MS/html/Z0.01377_Y0.2533_XD2E5_ML1.68_AS05.html

1. 2. visualize star/stellar evolution data
   3. single star pre-main sequence hr diagram
   4. python3 code
   6. 15.2.2023 0000.0001

import pandas as pd import numpy as np import matplotlib.pyplot as plt import matplotlib from scipy import interpolate

1. 1. modified pisa pre-sequence data file
3. http://astro.df.unipi.it/stellar-models/index.php?m=3

1. filename1="TRK_M1.00_Z0.01377_Y0.2533_XD2E5_ML1.68_AS05.DAT"
2. TRK_M1.00_Z0.01377_Y0.2533_XD2E5_ML1.68_AS05

filename1="mass1_solar1.csv"

df0 = pd.read_csv(filename1,header=1, dtype=np.float64,sep=';')

1. df0 = pd.read_csv(filename1,header=6, dtype=np.float64,sep=" ")

1. quit(-1)

logage1 = df0.iloc[:,1].to_numpy() logl1 = df0.iloc[:,3].to_numpy() logte1 = df0.iloc[:,4].to_numpy()

age1=np.power(10, logage1) lum1=np.power(10, logl1) teff1=np.power(10, logte1)

f1 = interpolate.interp1d(age1, lum1) f2 = interpolate.interp1d(age1, teff1)

1. 1. first 40 Ma

agenew2 = np.arange(100000, 40*1000*1000, 1000*100) lum2 = f1(agenew2) teff2 = f2(agenew2)

plt.gca().invert_xaxis()

plt.plot(teff1, lum1, ":")

len1=int(len(agenew2)/40)

for i in range(1, len1, 1): x=teff2[i] y=lum2[i] agema1=agenew2[i]/1000000 agema2=round(agema1,1) agema3=str(agema2) plt.plot(x,y,color="blue",marker='x') plt.text(x,y,agema3, color="green",size=15, alpha=0.7)

len2=int(len(agenew2)/10)

for i in range(9, len2, 5): x=teff2[i] y=lum2[i] agema1=agenew2[i]/1000000 agema2=round(agema1,1) agema3=str(agema2) plt.plot(x,y,color="blue",marker='x') plt.text(x,y,agema3, color="green",size=15, alpha=0.7)

len3=int(len(agenew2)/1)

for i in range(89-10, len3, 20): x=teff2[i] y=lum2[i] agema1=agenew2[i]/1000000 agema2=round(agema1,0) agema3=str(int(agema2)) plt.plot(x,y,color="blue",marker='x') plt.text(x,y,agema3, color="green",size=15, alpha=0.7)

1. plt.rcParams.update({'font.size': 17})

font = {'family' : 'DejaVu Sans',

       'weight' : 'bold',
       'size'   : 16}

1. matplotlib.rc('font', **font)

plt.suptitle("Sun-like pre main sequence star", size=17) plt.title("Evolution track 0.1 - 40 million years", color="green", alpha=0.7, size=15)

plt.xlabel("Temperature Teff K", size=17) plt.ylabel("Luminosity L/L_sun", size=17) plt.xticks(fontsize=15) plt.yticks(fontsize=15) plt.text(5800,11,"M=1.00 X=0.73293 Y=0.2533 Z=0.01377 ", color="blue", alpha=0.8, size=15) plt.text(5200,2, "Age Ma", color="green", alpha=0.9, size=18) plt.grid(color='r', linestyle='-', linewidth=1, alpha=0.2) plt.axhline(y = 1, color = 'y', linestyle = 'dashed', alpha=0.2, label = "Solar luminosity")

plt.show()