File:Hr pre main sequence 0.2msun 1.svg

Source of data is

Universita di Pisa database, pre-main sequence star

Pre-Main Sequence models

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

# M = 0.20 Mo, Y = 0.2533, Z = 1.37700E-02, alpha = 1.68, X_D = 2.0E-05
# Solar Mixture: Asplund et al. (2005)  

Python source code

##
## visualize star/stellar evolution data
## single star pre-main sequence hr diagram
## python3 code
##
## 15.2.2023 0000.0002
##

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

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

#filename1="mass1_solar1.csv"
filename1="mass0p2_solar1.csv"

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

#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)

## first 40 Ma

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

plt.gca().invert_xaxis()

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

len1=int(len(agenew2)/12)
	
for i in range(0, len1, 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)	
	
len2=int(len(agenew2)/1.2)
	
for i in range(95, len2, 100):
	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)	 

 
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)	
 
 
#plt.rcParams.update({'font.size': 17}) 

font = {'family' : 'DejaVu Sans',
        'weight' : 'bold',
        'size'   : 16}

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

plt.suptitle("Red dwarf (0.2 Msun) pre main sequence star", size=17)   
plt.title("Evolution track 0.5 - 80 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(3800,11,"M=1.00 X=0.73293 Y=0.2533 Z=0.01377 ", color="blue", alpha=0.8, size=15) 
plt.text(3200,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()

I, the copyright holder of this work, hereby publish it under the following license:

This file is licensed under the Creative Commons Attribution-Share Alike 4.0 International license.

You are free:

• to share – to copy, distribute and transmit the work
• to remix – to adapt the work

Under the following conditions:

• attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
• share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.

https://creativecommons.org/licenses/by-sa/4.0CC BY-SA 4.0 Creative Commons Attribution-Share Alike 4.0 truetrue