560 likes | 580 Views
AGB Star Yields for Everyone. John Lattanzio 1 and Amanda Karakas 2 1: Centre for Stellar and Planetary Astrophysics, Monash University,Australia 2: Origins Institute, Department of Physics & Astronomy, McMaster University , Canada. Structure of Talk.
E N D
AGB Star Yields for Everyone John Lattanzio1 and Amanda Karakas2 1: Centre for Stellar and Planetary Astrophysics, Monash University,Australia 2: Origins Institute, Department of Physics & Astronomy, McMaster University , Canada
Structure of Talk • Reminder of AGB stars and vital role • Outline of what we want • What we can do • What we cannot yet do • State of the art models • Synthetic AGB models • Strengths and weaknesses • Best yields available now (and soon) • Distribution
Public Speaking for Scientists Attention Time
Not this one… Trivia Question 1 • Who said: “To the extent that its possible, it’s the isotopes that keep the theorists honest.” • Dave Arnett (U Arizona)
Massive stars produce most of the galactic magnesium, which is primarily 24Mg at low Z But 3 - 6 Msun AGB stars can produce large amounts of the heavy magnesium isotopes (Y. Fenner, A. Karakas, B. Gibson, J. Lattanzio)
AGB stars are needed to recover the observed 25,26M/24Mg ratios at low metallicity Limongi et al. (2002) calculations generate more 25,26Mg thanWoosley & Weaver (1995) (Y. Fenner, A. Karakas, B. Gibson, J. Lattanzio, PASA, 2003)
GCE of 19F – Renda et al (2003) x = Milky Way = LMC SN, WR and AGB = w Cen SN & WR SN only
23Na 22Ne 23Na 14N 23Na 22Ne Sodium 12C Note: some 23Na is primary and some is secondary!
AGB star nucleosynthesis • Note the interplay of H and He burning… • Note the dependence on mixing regions… • Note that the return depends on mass-loss • Note that these are not the most precise inputs to the models…
Trivia Question 2 • What is the original definition of a SHIP? • A three masted sailing vessel with square-rigging on all three masts Barque Brig
What do we want? • Quantitative predictions of surface abundances over time for all M and compositions • Is this too much to ask???
Trivia Question 3 • Which famous (Italian) stellar astronomer said: “We provide evolutionary tracks, yields and physics from the pre-main sequence through to the thermally pulsing AGB. Its certainly wrong, but its freely available.” • Alessandro Chieffi
So why can’t we do it? • Well we can, but… • Each star takes a lot of cpu. • Typical AGB star (with my code) takes between 2 days and 4 weeks depending on the mass, mass-loss, and nucleosynthesis • There are many uncertainties…
Main Problems • Mass-loss • Rate determines when AGB ends • Hence huge effect on yields • Dredge-Up • Remains hard to do properly • Surely its hydrodynamic… • Convective borders and overshoot? Supernovae have the Mass-Cut...
How to Calculate Dredge-Up? • Depends on how you apply Schwarzschild …. some apply it naively and get NO dredge-up! • We search for a neutral point… a bit like overshooting… • Depends on how you implement it…
Calculating Dredge-Up No dredge-up Results depend on how you mix and iterate Deep dredge-up
Convection: Dredge-Up • Its not simple • Its not one-dimensional • Its not static • Its not mixing-length • But it is very very hard to calculate…
Convective Boundaries Eg Herwig and O… This will alter yields!
Trivia Question 4 • Who said: “The ultimate goal of nuclear physics is to produce dust.” • Erwin Sedlmayr (Has his own brewery in Berlin…)
The Best on Offer? • Best are full, detailed models • Star by star… • From ZAMS to end of AGB… • With the best you can do for all the physics… • One star per node on large cluster • Elements considered usually stop before the s-process…
Amanda Karakas’ Thesis • Three compositions • Solar, LMC and SMC (Z=0.02, 0.008, 0.004) • Since then has added Z=0.0001 • Many masses • M=1.5 to 7 • Various bits published but not all • Contact Amanda for yields…
Other Full Model Yields • Herwig (2004) Z=0.0001, M=2 to 6 • Simon Campbell (in Fenner et al 2004) • Globular Cluster chemical evolution • M=1.25 to 6.5 • [Fe/H]=-1.4 and mix from Big Bang and SN • Ie “Z” = 0.0017 (Z=1-X-Y is not proportional to Fe) • Ventura et al (2002) • Different theory of convection (FST, not MLT) • M=3 to 6 (HBB range) • Z=0.0002 to 0.01 (5 values)
Serious Advantage of Full Models • The interplay between HBB and TDU • Which ends first? • Unknown until we made some full detailed models go all the way to end of AGB • Has large effect on final envelope composition
Trivia Question 5 • Who were described by Ned Kelly as: “…big, ugly, fat-necked, wombat-headed, fat-bellied, magpie-legged, narrow-hipped, splay-footed sons of Irish bailifs or English landlords…” • The Victorian Police Force
The Other/Older Approach: Synthetic AGB Models • AGB stars are complicated • AGB stars take a long time to compute • But Paczynski noted that the surface L varied very nearly linearly with core mass • …at least for lower masses…
Core-mass v Luminosity L/L = 59250(Mc/M – 0.5)
Other relations can be fitted… • Mass in intershell convection zone • Duration of intershell convective zone • Interpulse period • Etc • Entire population can be synthesized very quickly!
Synthetic AGB Models • But only as good as the inputs • And things are not simple really • Eg dredge-up law does depend on the envelope mass as well as the core and composition • There is “memory” of the earlier life of the star, at some level
Not easy to know when to include Depends on M and Z Some turn it “on” at specific M and Z But what T? What rho? How much mass is burning? Very hard… The Biggest Problems? Hot Bottom Burning Dredge-Up • Can fit the dredge-up parameter l as a function of everything • But that has to come from detailed models and is not known… • Beware of this!
Advantages • Well suited to yields! • Quick! • Can change reaction rates quickly • Can investigate mass-loss dependence • In principle…very powerful • But limited by the inputs of course
Synthetic Yields Beware!! • Long history, starting with • Renzini and Voli 1981 • DO NOT USE RENZINI & VOLI (1981)! • Many better yields around today! • Who? • Van den Hoek & Groenewegen (1997) • Marigo (2001) • Forestini and Charbonnel (1997) • Izzard et al (2004)
Comparison of synthetic with detailed yields: 4He Black = Karakas (full) Blue = Izzard et al (syn) Red = Marigo (syn) Green = Forestini & Charbonnel (syn) Magenta = vdHG (syn) Cyan = Ventura (full) Big Variation...
12C 13C 14N Z=0.02 Z=0.02 Z=0.02 Z=0.008 Z=0.008 Z=0.008 Z=0.004 Z=0.004 Z=0.004
New Synthetic Yields from Izzard • Rob Izzard worked very hard on this • Aim is to include binaries • They go to places that single stars do not • So simple fits will not work all the time • Good fits to (our) detailed models • Real nucleosynthesis in most cases • Fit T and rho and then do real burning
Izzard Yields Synthetic = red Detailed = green
Izzard Yields Massive stars too Red and green = full models Cyan and magenta = synthetic
Huge number of extra parameters to do with binaries Includes mass accretion mass-loss mass transfer Add binary population as well
4.3 Population Synthesis H H 4He 4He 12C 12C 13C 13C
14N 14N 15N 15N 16O 16O 17O 17O
27Al 27Al 32S 32S 36Ar 36Ar 40Ca 40Ca
Trivia Question 6 • Which WWII General said: “I can take care of the Germans, but I am not sure I can beat Montgomery and Eisenhower.” • George Patton Bonus: "If Ike stops holding Monty's hand and gives me the supplies I'll go through the Siegfried line like shit through a goose."
What is the Best Now? • Full models are probably best • Monash group • Karakas, Campbell, Lattanzio etc • Ventura et al 2002 (limited mass range) • Herwig 2004 (Z=0.0001 only, and overshoot)
Synthetic if you must…. • Izzard et al 2004 • Beware: gap for M=7 to 12 = SUPER-AGB • Online yield calculator: • http://www.ciqua.org/binaryyields/index.php • For single stars set second at M=0.1 and very wide