File:Hr pre main sequence 1msun 1.svg
Original file (SVG file, nominally 1,138 × 731 pixels, file size: 60 KB)
Captions
Summary
[edit]DescriptionHr pre main sequence 1msun 1.svg |
English: Pre-main sequence Sun-like star evolution in HR diagram, 0.1-40 million years.
Teff/K versus L/LSun. |
Date | |
Source | Own work |
Author | Merikanto |
Source of data is Universita di Pisa database, pre-main sequence star
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
Python source code
- visualize star/stellar evolution data
- single star pre-main sequence hr diagram
- python3 code
- 15.2.2023 0000.0001
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="TRK_M1.00_Z0.01377_Y0.2533_XD2E5_ML1.68_AS05.DAT"
- 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=';')
- 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(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)
- plt.rcParams.update({'font.size': 17})
font = {'family' : 'DejaVu Sans',
'weight' : 'bold',
'size' : 16}
- 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()
Licensing
[edit]- 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.
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 11:19, 15 February 2023 | 1,138 × 731 (60 KB) | Merikanto (talk | contribs) | Update of time labels, simpler layout | |
10:41, 15 February 2023 | 982 × 676 (49 KB) | Merikanto (talk | contribs) | Uploaded own work with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
Metadata
This file contains additional information such as Exif metadata which may have been added by the digital camera, scanner, or software program used to create or digitize it. If the file has been modified from its original state, some details such as the timestamp may not fully reflect those of the original file. The timestamp is only as accurate as the clock in the camera, and it may be completely wrong.
Width | 910.08pt |
---|---|
Height | 584.64pt |