Hallo MAC,
ich habe mal schnell ein Programm geschrieben was die verschiedenen IMF simuliert.
einschränkende Annahmen keinen SN Ia, aber auch kein neueinströmendes Gas
Naherung: Sterne in Intervallen vom Faktor ~2
Annahmen: Die SFR(star formation rate) nimmt von -7,5Gyr auf heute auf S exponentiell ab, vorher konstant. Quelle "Madau-Plot", (es sollte auch spezelle errechtete SFR für die Milchstraße geben)
Ann. Metallfreisetzung nur durch Supernovae (außer Ia), d.h durch Sterne >)8Msonne
Sterne ab ~1Msonne werden nach Alter A WD (N,BH)
Vereinfachende Annhamen über Massen der "Reststerne"
1 -> 0,55Ms WD
2 -> 0,65Ms WD
4 -> 0,85Ms WD
8 -> 1,27Ms WD/N grenze
16 -> 1,4Ms N
32 -> 1,5? Ms N/BH Grenze
64 -> 3??Ms BH
128 -> 6??Ms BH
150 ->10??Ms BH???
Maxmale Masse 150Ms
Weiterhin, 6/25 der Ursprungsmasse minus Masse des N/BH, wird bei Supernova als Metall ausgeworfen. (passt zumindest für ~25Msstern)
Parameter
A Alter der Scheibe in Gyr
M Masse der Sterne+Gas heute in Mrd Msonne
G Gas Masse heute in Mrd Msonne
S Abnahmefaktor der SFR von -7,5Gyr zu heute (vorher konstent)
IMF
Output
MG gesamte je gebildete Sternenmasse
Z Metallizität heute
SNR1 Supernovarate/1000Jahre bei gleichmäßiger Sternentstehung
SNR2 Supernovarate/1000Jahre bei SFR nach Madau-plot
M_BD Masse in BD 0,01-0,08Ms heute
M_M Masse in M-Zwergen 0,08-0,5Ms heute
M_s Masse in Sternen >0,5, +WD+N+BH (sonstige)
N ist jeweils Anzahl, NG gesamtanzahl der Sterne
A=9, M=60, G=10, s=0,1
Modell Kroupka 1
alpha=0,3 (0,01-0,08)
alpha=1,3 (0,08-0,5)
alpha=2,3 (>0,5)
MG=84,2 NG=218,5 Z=0,048 SNR1=94,9 SNR2=18,1
M_BD=3,2 M_M=22,1 M_st=24,7
N_BD=81,1 N_M=105,3 N_s=32,1
A=9, M=60, G=10, s=0,1
Modell Kroupka 2
alpha=0,3 (0,01-0,08)
alpha=1,3 (0,08-0,5)
alpha=2,3 (0,5-1)
alpha=2,7 (>1)
MG=67,2 NG=230,0 Z=0,013 SNR1=33,5 SNR2=6,4
M_BD=3,4 M_M=23,6 M_st=22,9
N_BD=86,7 N_M=112,5 N_s=30,8
A=9, M=60, G=10, s=0,1
Modell Miller-Scalo
alpha=0,3 (0,01-0,08)
alpha=1,3 (0,08-0,5)
alpha=2,3 (0,5-1)
alpha=2,7 (>1)
MG=78,2 NG=223,3 Z=0,018 SNR1=58,2 SNR2=11,1
M_BD=3,4 M_M=20,3 M_st=26,2
N_BD=91,8 N_M=97,1 N_s=34,5
A=9, M=60, G=10, s=0,1
Modell MAC 1
alpha=0,3 (0,01-0,08)
alpha=1,3 (0,08-0,5)
alpha=1,6 (>0,5)
MG=268,9 NG=137,9 Z=0,438 SNR1=671,9 SNR2=128,3
M_BD=1,7 M_M=11,9 M_st=36,3
N_BD=43,8 N_M=56,8 N_s=37,3
A=9, M=100, G=10, s=0,1
Modell MAC 2
alpha=0,3 (0,01-0,08)
alpha=1,3 (0,08-0,5)
alpha=1,6 (>0,5)
MG=484,1 NG=248,3 Z=0,463 SNR1=1209 SNR2=231
M_BD=3,1 M_M=21,5 M_st=65,4
N_BD=78,8 N_M=102,3 N_s=67,1
A=9, M=1000, G=10, s=0,1
Modell MAC 3
alpha=0,3 (0,01-0,08)
alpha=1,3 (0,08-0,5)
alpha=1,6 (>0,5)
MG=5324,9 NG=2730,8 Z=0,496 SNR1=13303 SNR2=2541
M_BD=34,5 M_M=236,1 M_st=719,4
N_BD=867,1 N_M=1125,3 N_s=738,4
Welche IMF hast du eigentlich bei deinen Berechnungen wirklich verwendet?
Soll ich noch andere Parameter probieren?
Grüße UMa