560 likes | 653 Views
Lessons on ISM modelling from kiloparsec scale simulations. Adrianne Slyz. University of Oxford. Zurich, September 18 th 2007. What physical processes regulate …. Kennicutt (1998). starbursts.
E N D
Lessons on ISM modelling from kiloparsec scale simulations Adrianne Slyz University of Oxford Zurich, September 18th 2007
What physical processes regulate … Kennicutt (1998) starbursts centers of normal disks logΣSFR (Msol yr-1 kpc-2) normal disks SFR Area ∑gas 1.4 ~ log Σgas(Msol pc-2) the rate at which gas turns into stars? Zurich, September 18th 2007
Link between density structure & star formation ? Zurich, September 18th 2007
Key insights from periodic box simulations in the 90’s: 1. Density structure of isothermal medium structured by supersonic, compressible turbulence is well described by a log-normal distribution whose dispersion reflects the Mach number of the medium (Vazquez-Semadeni 1994, Padoan, Nordlund & Jones 1997, Nordlund & Padoan 1999) 2. Isothermal and adiabatic turbulence decays quickly (within a sound crossing time) whether the medium is magnetized or not (Stone et al. 1998, Maclow et al. 1998, Padoan & Nordlund 1999) Zurich, September 18th 2007
Initial conditions Homogeneous gas ρ= 1 at/cm3 1.28 kpc T = 105 K Turbulent velocity field imposed on large scales P(k) ∝ k-4 1.28 kpc 1.28 kpc Periodic boundary conditions Zurich, September 18th 2007
Initial conditions Homogeneous gas ρ= 1 at/cm3 1.28 kpc T = 105 K Turbulent velocity field imposed on large scales P(k) ∝ k-4 1.28 kpc 1.28 kpc Periodic boundary conditions Zurich, September 18th 2007
Radiative cooling -21 -22 -23 -24 -25 -26 log10(T) ergs cm3/s 3 4 5 6 7 8 T (Kelvin) White and Sarazin (1987), Rosen and Bregman (1995) Zurich, September 18th 2007
PDF = 1___ (2 π)1/2σ - (lnρ- < lnρ>) 2 2σ2 exp with σ= ln 1 + (Mrms/2)2 log10(normalized PDF) log10gas - <log10gas> ( ) [ ] (Padoan & Nordlund 2002) Zurich, September 18th 2007
PDF = 1___ (2 π)1/2σ - (lnρ- < lnρ>) 2 2σ2 exp with σ= ln 1 + (Mrms/2)2 log10(normalized PDF) log10gas - <log10gas> ( ) [ ] (Padoan & Nordlund 2002) Zurich, September 18th 2007
PDF = 1___ (2 π)1/2σ - (lnρ- < lnρ>) 2 2σ2 exp with σ= ln 1 + (Mrms/2)2 log10(normalized PDF) log10gas - <log10gas> ( ) [ ] (Padoan & Nordlund 2002) Zurich, September 18th 2007
PDF = 1___ (2 π)1/2σ - (lnρ- < lnρ>) 2 2σ2 exp with σ= ln 1 + (Mrms/2)2 log10(normalized PDF) log10gas - <log10gas> ( ) [ ] (Padoan & Nordlund 2002) Zurich, September 18th 2007
PDF = 1___ (2 π)1/2σ - (lnρ- < lnρ>) 2 2σ2 exp with σ= ln 1 + (Mrms/2)2 log10(normalized PDF) log10gas - <log10gas> ( ) [ ] (Padoan & Nordlund 2002) Zurich, September 18th 2007
fc =∫ ρPDF dρ ρth ∞ ∫ ρPDF dρ 0 Open questions Universal PDF? (Elmegreen 2002 Krumholz & McKee 2005 Wada & Norman 2007) ∞ Is there a clear density threshold for star formation? fraction of mass above th Zurich, September 18th 2007
PDF = 1___ (2 π)1/2 σ - (ln ρ - < ln ρ >) 2 2σ2 exp Log-normal fit to high density end of run with stars, self-gravity, fbk ( ) log10(normalized PDF) < lnρ> = 3.9 ⇒ρpeak ≈ 50 at cm-3 σ= 1.22 ⇒Mrms= 3.1 log10gas (at/cm3) Zurich, September 18th 2007
PDF for SN driven stratified segment of a disk 1 X 1 X 20 kpc density Probability Distribution Function (PDF) 10 8 6 4 2 1 0.5 0 -0.5 -1 -2 -4 -6 -8 -10 x-y plane PDF Z (kpc) n (cm3) Avillez & Breitschwerdt 2005 Zurich, September 18th 2007
PDFs in different subbox sizes for SN driven stratified segment of a disk 0.5 X 0.5 X 10 kpc disk segment PDFs for gas near midplane 125 pc subbox 4 pc subbox Joung & MacLow 2006 Zurich, September 18th 2007
3D isolated disks, 25-50 pc resolution New generation of ISM simulations density temperature pressure stars face-on view 7.5 PDF log Temp 2.5 0.001 1.0 1000 (Msol pc3) edge-on view Tasker & Bryan 2006 Zurich, September 18th 2007
Different philosophies for adding supernovae explosions in ISM models A B Model observed supernovae rates & mimic their distribution (e.g. isolated, clustered) Model star formation m*=ρgasVcell ∆t/tdyn + Stellar Initial Mass Function e.g. SN frequency Milky Way Galaxy:1/330 yr-1 for Type I and 1/44 yr-1 for Type II (Tammann et al. 1994) Scale heights: Type I :325 pc (Heiles 1987) Type II :90 pc Power law distribution of superbubbles: dNB ~ n*-2dn* (Kennicutt et al. 1989, McKee & Williams 1997) Calculate energy and mass returned to interstellar medium via supernovae and stellar winds Zurich, September 18th 2007
B Supernovae feedback in Joung & Maclow 2006 1.) Identify supernovae site (stick to the observations) 2.) Grow a sphere at that site until it encloses 60 Msun. Radius of this sphere Rexp ~ 7 pc to 50 pc Rexp Mexp = 60 Msun 3.) Redistribute mass in that sphere so that it has uniform density = 3Mexp/(4 Rexp3) 4.) Inject thermal energy ESN = 1051 ergs evenly into the sphere NO mass ever removed to form a star Zurich, September 18th 2007
= 1 ~ 1 order of magnitude = 0.3 predictions Compare estimated SFR to input SN rates SFRs derived from input SN rates Identify Jeans unstable boxes Mbox/MJ > 1 where MJ = <J3 = avg density in box J = (/G<)1/2tot tot = (<cs>2 + 1/3<2)1/2 (Chandrasekhar 1951) (2 input supernovae rates: assuming 130 and 200 Msol required per SN) SFR = Mbox/tff where = 0.3 or 1 Joung & MacLow 2006 Zurich, September 18th 2007
Mrms=5.1 Mrms=11.8 Mrms=5.2 Mrms=5.8 no fbk, with s-g Is that because they ignore self-gravity in their model? Gotoh & Kraichnan (1993) found power law PDFs for 1D sims of Burgers flows⇒ infinitely compressible flows (Slyz et al. 2005) Zurich, September 18th 2007
r > rthresh ->dense T < T thresh->cold A First make stars . . . Heyer et al. 1998 (FCRAO CO survey) m*= ερgasVcell ∆t/tdyn if the gas satifies: Cen & Ostriker 1992 v < 0->contracting t cool< t dyn->cooling rapidly Zurich, September 18th 2007
1.28 kpc Then do stellar feedback . . . Cen & Ostriker 1992 Calculate a time dependent SFR: DmSF (t) =m*(t-t*)/т2 exp [-(t-t*)/τ] whereτ = max(tdyn, 10 Myr) Stellar winds: f ∆mSF returned to gas 100 pc @ 10 Myr if v=10 km/s Supernovae: ∆mSF c2 injected as thermal energy f, determined by IMF Zurich, September 18th 2007
Non-instantaneous feedback temp pressure density 4.5 Myr 22 Myr 41 Myr Slyz, Devriendt, Bryan, Silk (2005) Zurich, September 18th 2007
Instantaneous feedback temp pressure density 4.5 Myr 22 Myr 41 Myr Slyz, Devriendt, Bryan, Silk (2005) Zurich, September 18th 2007
Is this the same old story. . . ? When put supernova thermal energy ESN = 1051 ergs in dense regions most of the energy is quickly radiated away? Get neither thermal or dynamical heating (Katz 1992) Fixes:1)artificial time delay in cooling(Gerritsen 1997; Thacker & Couchman 2001; Governato et al. 2006) 2)assign explosion energy to fluid parcels as pure kinetic energy(Navarro & White 1993) 3)introduce a thermalization efficiency whereby assign some fraction of supernova energy as kinetic and some as thermal(Navarro & White 1993, Hernquist & Mihos 1995) 4)sub-grid models of multi-phase ISM (Yepes et al. 1999, Springel & Hernquist 2003) Zurich, September 18th 2007
Time evolution of density PDF green: inst fbk, black: non-inst fbk Zurich, September 18th 2007
Time evolution of energy spectra s-g no fbk no s-g no fbk no s-g fbk s-g fbk compressible ∇✘vcom= 0 solenoidal ∇⋅vsol = 0 ratio Zurich, September 18th 2007
Time evolution of energy spectra s-g no fbk Instantaneous feedback no s-g no fbk no s-g fbk s-g fbk compressible ∇✘vcom= 0 solenoidal ∇⋅vsol = 0 ratio Zurich, September 18th 2007
Comparison of inputs intoSilk prescription Q = SFR G-1/2 ρgas-3/2(σgas / σf ) -2.72 0.8 0.6 0.4 0.2 0.0 SFR (Msun/yr) no fbk, no gravity 5 4 3 2 1 Porosity no fbk, with gravity -1.5 -2.0 -2.5 -3.0 Different physics withnon-instantaneous fbk, no gravity log10<ρ> (Msun/pc3) withnon-instantaneous fbk, with gravity 40 30 20 10 0 <σ>MW (km/s) withinstantaneous fbk, with gravity 0 100 200 300 time (Myr) Zurich, September 18th 2007
Effect on star formation rate 4.5 Myr pressure temp density non-instantaneous feedback density 22 Myr pressure temp 0.8 0.6 41 Myr SFR (Msun/yr) 0.4 0.2 instantaneous feedback 0 100 200 300 time (Myr) Slyz, Devriendt, Bryan, Silk (2005) Zurich, September 18th 2007
Time evolution of density PDF green: inst fbk, black: non-inst fbk Zurich, September 18th 2007
How to erase a thermal instability… 1 kpc2 box P(k) ∝ k-4 Large scale forcing log (number of cells) vs Small scale forcing > thresh v < 0 heat for 6 X 106 yrs to mimic log « photo- ionization » Vazquez-Semadeni, Gazol, Scalo 2000 Zurich, September 18th 2007
Time evolution of density PDF non-instantaneous feedback run Zurich, September 18th 2007
Time evolution of density PDF non-instantaneous feedback run SFR Time→ ~85 Myr Zurich, September 18th 2007
Time evolution of density PDF non-instantaneous feedback run Zurich, September 18th 2007
Time evolution of density PDF green: inst fbk, black: non-inst fbk Zurich, September 18th 2007
Time evolution of thermal phase diagrams Initial conditions Zurich, September 18th 2007
Time evolution of thermal phase diagrams Initial conditions Zurich, September 18th 2007
Thermally unstable regime Zurich, September 18th 2007
Lines of constant pressure (kB-1 cm-3 K) Zurich, September 18th 2007
1D cut Accretion shock! 2D pressure map x-velocity Density, pressure X (kpc) Zurich, September 18th 2007
Different philosophies for adding supernovae explosions in ISM models A B Model observed supernovae rates & mimic their distribution (e.g. isolated, clustered) Model star formation m*=ρgasVcell ∆t/tdyn + Stellar Initial Mass Function e.g. SN frequency Milky Way Galaxy:1/330 yr-1 for Type I and 1/44 yr-1 for Type II (Tammann et al. 1994) Scale heights: Type I :325 pc (Heiles 1987) Type II :90 pc Power law distribution of superbubbles: dNB ~ n*-2dn* (Kennicutt et al. 1989, McKee & Williams 1997) Calculate energy and mass returned to interstellar medium via supernovae and stellar winds Zurich, September 18th 2007