Guillermo Haro Workshop, 2015 Current Paradigm of Direct Collapse BH formation Bhaskar Agarwal Yale University (prev: PhD at Max Planck for Extraterrestrial Physics)
What is the Problem? SDSS Quasars at z > 6 highest redshift: z = 7.085 (770 Myr) (Momjian et al. 2013) most massive: z = 6.30 (900 Myr), M = 1.3x10 12 (Wu et al. 2015) First Stars form at z = 20 Latest Planck results push the z lower! How do we grow the Quasars to supermassive scales? time scale is t(z 6 ) - t(z first_star ) ~ 700 Myr* *Take this time with a pinch of salt
Solving the problem z 20 9 6 t 200 Myr 550 Myr 1 Gyr seed' a BH grow the seed make it a SMBH - For a given BH growth scenario: Bondi-Hoyle, Eddington etc. - sensitive to two parameters M seed t grow
The Possible Paths pristine gas in minihalo withstand SF Pop III/II cluster Pop III star: BH: 100 M-sun halo grows: T vir = 10 4 K critical LW radiation NO fragmentation: SMS/QS Merge before 100 M sun go SN BH: 1000 M sun super-edd. BH growth BH: 10 4-5 M sun sub-edd. BH growth Quasar BH: 10 9 M sun
PopIII stellar BH seeds Heger et al. 2003
PopIII stellar BH seeds Volonteri et al. Bromm et al. Yoshida et al. Assuming 10 10 z 19.00 11.50 8.50 7.00 6.00 f edd = 1 rad. eff = 0.1 M seed = 100 M sun M BH [M sun ] 10 8 10 6 10 4 t grow : t 1 {z(20-6)} = 700 Myr t 2 {z(9-6)} = 200 Myr z 20 10 2 200 400 600 800 1000 t age [Myr] 6 t 200 Myr 1 Gyr seed' a BH make it a SMBH
Hirano et al. 2014 PopIII stellar BH seeds: IMF PopIII : IMF Highest mass: 1000 M sun Mstar > 100 M sun desirable seed mass! Doesn t have to be too' common
PopIII stellar BH seeds: where? DM Mass 10 5-6 10 6-7 10 7-8 10 9-13? BH + II* SMBH + II* BH + II* III* PopIII : first sources of light! Minihaloes: M = 10 5-6 M sun (DM) Effect on neighbouring haloes: -Moderate LW feedback can kill SF -X-rays feedback -SN 'blow out'
PopIII stellar BH seeds: problems highly optimistic scenario! ignores feedback Assuming f edd = 1 rad. eff = 0.1 where is the gas? constant availability of gas form in mini-haloes! merger scenario unclear merger b/w 2 mini-haloes - time-scale - BH present in both M seed = 100 M sun t grow : t 1 {z(20-6)} = 700 Myr t 2 {z(9-6)} = 200 Myr need to be super-edd! all the time!
Direct Collapse BH seeds 10 10 Eisenstein & Loeb Oh & Haiman Bromm & Loeb Koushiappas et al. Bullock & Dekel Lodato & Natarajan z 19.00 11.50 8.50 7.00 6.00 Assuming f edd < 1 rad. eff = 0.1 M seed = 10 4-5 M sun M BH [M sun ] 10 8 10 6 10 4 t grow : t 1 {z(20-6)} = 700 Myr t 2 {z(9-6)} = 200 Myr 10 2 200 400 600 800 1000 t age [Myr] z 20 9 6 t 200 Myr 550 Myr 1 Gyr seed' a BH make it a SMBH
DCBH seeds: how LW photons: 11.2-13.6 ev Enough LW Radiation H2 + {LW} H2 Delay Pop III Star formation H2 H + H Tvir < 10,000 K
DCBH seeds: how LW photons: 11.2-13.6 ev Critical (high) LW Radiation H H H2 Stop Pop III Star formation H + H H Tvir > 10,000 K Pop III n = 10 5 cm 3 T = 200 K 8000 DCBH n = 10 5 cm 3 T = 8000 K M BE ~ 100 M sun M BE ~ 10 5 M sun
DCBH seeds: problems Assuming f edd < 1 rad. eff = 0.1 M seed = 10 4-5 M sun What is the critical LW flux? HOW does the gas collapse proceed? Feedback from the DCBH once it forms t grow : t 1 {z(20-6)} = 700 Myr t 2 {z(9-6)} = 200 Myr Subsequent merger: with the LW producing galaxy z 20 9 6 t 200 Myr 550 Myr 1 Gyr seed' a BH make it a SMBH
Do you like movies? Films, cinema. NOT simulation movies My guilty pleasure: BOLLYWOOD! 1. Introduce the audience to your plot 2. Everything is OK, songs/dance 3. Something bad happens! 4. Oh No! 5. There is hope
DCBH seeds: where Close to a LW source: Pop II galaxy M = 10 7 M sun (DM) DCBH Continued exposure to LW radiation LW LWcrit Mhot Mhot tdyn tdyn Mcold Mcold tdyn tdyn Metal Free Pop III Polluted Pop II Mhot Outflows Mhot Outflows
DCBH seeds: The Hunt SAM Agarwal et al. 2012 FiBY Dalla Vecchia et al. Agarwal et al. 2014 Cosmological Volume Yes Yes Minihaloes resolved Yes Yes Halo histories Yes Yes Outflows Yes Yes Pop III + Pop II Yes Yes LW radiation: Spatial + Global Yes Yes Gas No Yes IGM Metal Dispersion No Yes
DCBH seeds: The Hunt (SAM) SAM Halo history: Merger Trees Subfind output: LCDM simulation Min halo mass required : at least 10 6 M solar Box Size Cosmology N particles Min. Halo Mass z final 4.86 Mpc WMAP 7 768 3 10 5 M solar 6
DCBH seeds: The Hunt (SAM) DCBH LW LWcrit Mhot Mhot tdyn tdyn Mcold Mcold tdyn tdyn Mhot Mhot Outflows Outflows
DCBH seeds: The Hunt (SAM) Reionisation feedback: 13.6 ev : photoionisation / heating of gas in haloes gas needs a larger potential well to re-collapse Vc = 20 km/s (Dijkstra 08) LW Escape Fractions: more photons = more feedback (Ahn et al. 09, Kitayama 04) Stars: Pop II stars (Starburst 99) sfe: 1 % Pop III stars (Schaerer et al.) Range of SAM models Each explores a different parameter set Most vs. Least optimistic scenarios
DCBH seeds: The Hunt (SAM) Agarwal et al. 2012
DCBH seeds: The Hunt (SAM) 10 6 10 4 PopII J max PopII J bg PopII J crit PopIII J max PopIII J bg PopIII J crit J bg 10 2 J LW 10 0 10-2 10-4 10 15 20 25 z Agarwal et al. 2012
DCBH seeds: The Hunt (SAM) 2.0 esc0.1 2.5 log(dn/dz) [Mpc 3 ] 3.0 3.5 2.0 esc0.5 2.5 3.0 3.5 2.0 2.5 3.0 esc1.0 3.5 2.0 esc0.5hsfe 2.5 3.0 3.5 2.0 esc0.5reion 2.5 3.0 3.5 7 8 9 10 11 12 13 1 + z DC sites: Few (<10) at z>6 in a ~100 cmpc 3 box Quasars: few / Gpc 3 at z>6 Discrepancy: 10 7 objects too many! (?) Agarwal et al. 2012
DCBH seeds: The Hunt (SAM) 100.0 z=11.9 10.1 10.0 1.0 3.66, 0.801 3.34, 0.737 Steeper Two point correlation ξ Halo(DC) / ξ Halo(NoDC) 100.0 10.0 1.0 0.1 2.13, 0.451 0.591, 0.081 z=9.98 8.10 z=7.99 6.09 DC sites prefer a more clustered neighbourhood Need to be close to a larger galaxy giving out critical LW flux 10.0 1.0 0.1 2.29, 0.461 0.381, 0.041 2.5 3.0 3.5 4.0 4.5 5.0 log(d phy ) [pc] Agarwal et al. 2012
DCBH seeds: The Hunt (FiBY)
DCBH seeds: The Hunt (FiBY) 20 30 20 20 4.00 10 10 10 Y [phy. kpc] Y [phy. kpc] 0-10 -20 20 10 0-10 30-20 -10 0 10 20 X [phy. kpc] Z [phy. kpc] Z [phy. kpc] 0-10 -20 20 10 0-10 30-20 -10 0 10 20 Y [phy. kpc] X [phy. kpc] Y [phy. kpc] 0-10 -20 20 10 0-10 30-20 -10 0 10 20 Z [phy. kpc] 1.50-1.00-1.00-4.50 Log J LW,II Log Z sun Agarwal et al. 2014a -20-20 -10 0 10 20 X [phy. kpc] -20-20 -10 0 10 20 Y [phy. kpc] -20-20 -10 0 10 20 X [phy. kpc] -8.00
DCBH seeds: The Hunt (FiBY) J LW 10 2 10 1 10 0 10-3 DC3 ON08 J LW,II J LW,III J LW,Total 10-1 DC 10 (sites) seems to be more probable than previously thought -2 10 7 8 10 12 14 16 18 20 22 redshift T vir =10000 K redshift 30 25 20 15 first prog. next prog. DC0 at z=10.50 10.00 (Agarwal et al. 12 vs. Dijkstra et al. 08) Discrepancy broken if a higher J crit used 7.50 Log(M DM ) [M sun ] (Dijkstra et al 2015) Agarwal et al. 2014a M DM [M sun ] 10 6 T vir =2000 K 10 8 10 12 14 16 18 20 22 redshift 5 5.00
DCBH seeds: The Evolution DCBH Merger In Situ Pop III Star Pop II Cluster Empty halo Growth of BH ± Stars
DCBH seeds: The Evolution Metal Free Pop III Polluted Pop II DCBH Outflows Mhot LW Mhot Mhot LWcrit tdyn Mcold DCBH tdyn tdyn Mcold tdyn Mcold tdyn OBG Phase Mhot flim Mhot Outflows Mcold + Mhot Outflows
DCBH seeds: The Evolution 10 9 10 8 O1, z g =6.18 O2, z g =8.93 O3, z g =7.34 O4, z g =6.55 M BH [M O ] 10 7 10 6 10 5 10 4 f lim = 1.5 f lim = 0.75 f lim = 0.1 f edd =1 10 5 10 6 10 7 10 8 10 9 M * [M O ] Agarwal et al. 2013
DCBH seeds: Observations 20 22 Rest wavelength [Å] 10 3 10 4 O1 O2 O3 O4 f lim = 1.5 f lim = 0.75 f lim = 0.1 24 m AB 26 28 30 32 34 10 4 10 5 Observed wavelength [Å] Agarwal et al. 2014
DCBH seeds: Revised Theory? L hv - Rates depend on the shape of incident spectrum of external galaxy - In literature the J crit values have been derived by assuming Pop III source: black body at 10 5 K 1000 Pop II source : black body at 10 4 K 30-100 (Omukai et al. 01, Shang et al. 10, Wolcott-Green et al. 12, Sugimura et al. 14, Dijkstra et al. 14, Agarwal et al. '15)
DCBH seeds: Revised Theory? Rates depend on the shape of the incident spectrum from the external galaxy: relative ratios of 1eV to 12.4 ev photons Agarwal et al. 2014b
DCBH seeds: THE J crit Agarwal et al. 2015 Parameter space reaction rates: kh -, kh2 DC Input the rates directly into the Enzo framework with S.Glover 15 chem. PopIII Find the rate parameter space where you find DC Compare it to rate parameter space from realistic SEDs: Starburst 99
DCBH seeds: THE J crit from Starburst99 Agarwal et al. 2015 4 ST99: galaxy is placed 5 kpc from the atomic cooling halo de: photo-detachment of H - di: photo-dissociation of H 2
DCBH seeds: THE J crit ST99: ISF, Z = 0.004 IMF: Salpeter [1-100] No unique" value of J crit : depends on the age, SFH of the galaxy Agarwal et al. 2015
DCBH seed: Found? Observed by: Sobral et al. 2015 Claimed as DCBH: Pallotini et al. 2015
DCBH seed: found? asterisks: B+C diamonds: A A: BH B+C: stellar Agarwal et al. 2015b in prep
DCBH seeds: Recently Latif et al. NO fragmentation Aykutalp et al. X-ray feedback Glover et al. relevant chemistry Scahuer et al. fesc of LW Sugimura et al. Jcrit
DCBH: Take Away Message DCBH works maybe too well still << frequent as Pop III seeds Local LW flux essential for z>6 galaxy formation models DCBH form in satellites of LW producing galaxies DCBH prefer clustered environments DCBH can merge with the galaxies later on: OBG Need to account for realistic SEDs : BB not enough QUESTIONS QUESTI All DCBH are not 'future' SMBH
Stellar cluster BH seeds Zwart et al. Volonteri et al. Katz et al. Cluster M = 10 4-6 M sun Merge within few Myr M BH = 10 4-5 M sun Assuming f edd < 1 rad. eff = 0.1 10 10 10 8 z 19.00 11.50 8.50 7.00 6.00 M seed = 1000 M sun t grow : t 1 {z(20-6)} = 700 Myr t 2 {z(9-6)} = 200 Myr M BH [M sun ] 10 6 10 4 10 2 200 400 600 800 1000 t age [Myr] z 20 9 6 t 200 Myr 550 Myr 1 Gyr seed' a BH make it a SMBH
Stellar cluster BH seeds: where? DM Mass 10 5-6 10 6-7 10 7-8 10 9-13? MBH + II* SMBH + II* MBH + II* II*/III* cluster A proto-galaxy type stellar cluster Metals M = 10 6-7 M sun (DM) III* Pre-enrichment : VERY early on
Stellar cluster BH seeds: problems Cluster M = 10 4-6 M sun Merge within few Myr M BH = 10 4-5 M sun Assuming f edd < 1 rad. eff = 0.1 M seed = 1000 M sun t grow : t 1 {z(20-6)} = 700 Myr t 2 {z(9-6)} = 200 Myr formation epoch: Pop II cluster- VERY early on N-body problem mergers: all stars must merge timescale < 3 Myr z 20 9 6 t 200 Myr 550 Myr 1 Gyr seed' a BH make it a SMBH