SLIDE 1 Massive Black Hole Growth and Formation: Implications for LISA
- 1. Supermassive Black Holes:
When, Where, and How?
Theoretical Issues Observational Constraints & Clues
- 2. Does a Galaxy Merger Imply a
Black Hole Merger?
Where Are the Binaries? Gas-Rich vs. Gas-Poor Mergers
Do They Exist? Pop. III Seeds
Fan et al. 2003
P.Coppi, Yale
SLIDE 2
The existence of massive black holes is not necessarily so surprising. Many roads lead to a massive black hole? (Gravity is one way.) Rees 1984
SLIDE 3 2 6 Salpeter 0.1 9 Hubble
vs. 4 / exponential growth on timescale t 45 10 yr 4 t ( 6) 10 yr, i.e., marginally sufficient number of
acc BH Edd BH p T T p
L M c L GM m c Gm z ε π σ εσ ε π = = → = ≈ ≤ &
6 8 acc
growth e-foldings possible (even for 100 M seed) Problem worse if t / ( / ) 10 10 yr
AGN Gal
N c H −
Timescale Problem: Need to pack a lot of gas into small region FAST!
SLIDE 4 Formation of the First Quasars
- Seed BH by direct collapse of primordial gas cloud
Mass ~ 109 Mo, R ~ 1 kpc zvir = 5, no DM
Stars Gas
(Loeb & Rasio 1994, ApJ, 432, 52)
- Problem:
- Gas cooling
- Fragmentation
No compact central object! Neglected
- Star Formation
- Negative Feedback (SNe)
SLIDE 5
How to make life easier?
Give up Eddington limit? Not so hard… accreting compact objects in our galaxy seem to do it! Remember Eddington limit only applies to isotropic configurations, while gas flows in strong radiation fields are not. If really throw a lot of gas onto hole, can trap radiation and advect it into hole (e.g., Begelman 1978) Pre-existing Massive Seeds! Are they big (106 solar masses) or small (10-100 solar masses)? When do they appear and how rare are they? Big impact on LISA event rate!
SLIDE 6 Bahcall et al. 2000
HST QSO hosts
SLIDE 7
A “boring” object in the sky: the nearby elliptical galaxy M87
Optical Radio
SLIDE 8
SLIDE 9 Soltan 1982-type argument/problem:
5 1 3 6 1 3 lg lg /
1.4 10 ( / 0.01) . 1.1 10 ( / 0.002)( / 0.002 )
acc B BH relic BH Bu e Bu e
f M Mpc vs M h h M Mpc M
σ
ρ ε ρ
− − − −
= × = × < > Ω
accretion irrelevant, mergers key? most of accretion activity missed, i.e, "dark" (dust obscuration, low radiative eff / 1? iciency
acc BH
ρ ρ ⇒ accretion solns: super-Eddington photon-trapping, ADAF, etc.)
(e.g., Natarajan 1999 review)
SLIDE 10
Gilli 2003 review – astro-ph/0303115 The X-Ray Background (mostly AGN, hard X-rays clean BH signal)
SLIDE 11
Urry & Padovani 1995
The “Unified” AGN Model: Type I Type 2 Orientation Effect?
SLIDE 12
Type I “Type II”
+
The standard ingredients for an XRB model (e.g., Comastri et al. 1995, Gilli et al. 2002)
SLIDE 13
Deep ROSAT (one week exposure) of Lockman Hole Region
1000-2000+ sources per square degree!
[Chandra spacecraft can do this in half a day!] An image in only the 0.3-2.0 keV energy band
SLIDE 14
CDF-N/GOODS, 2Msec(!) Chandra CDF-S
The acid test: what happens when you start adding hard X-rays?
SLIDE 15 Aside: Chandra < 1” angular resolution absolutely critical! At R=27+ (>40% faint Chandra sources),
- ptical source density is huge …
[counterpart confusion serious problem for ROSAT, and even XMM]
Real HST GEMS data, w/real Chandra (1.5”) + simulated XMM (8”) error circles superimposed.
SLIDE 16 Success?! NO!
Hasinger 2003
SLIDE 17
Hasinger et al. 2003 version
N.B. cosmic variance!
SLIDE 18 Broad Line Narrow Line ?Unknown
More narrow line low-power
at lower z?
Berger et al. 2003
SLIDE 19 Wilkes et al. 2002
If select quasars by non-standard technique, indeed find “weird” objects!
2MASS Red/NIR Quasar Survey
(very bright, nearby objects; analog of Hellas2XMM)
Some show broad optical emission lines but absorbed X-rays??
Standard, old result
SLIDE 20
Heckman (+ SDSS) 2003 Sample of 56,000 emission-line galaxies!
starbursts ?composite
Further complications… confusion w/starburst
SLIDE 21 Cutri et al. 2001; Smith et al. 2001
? IR Detection of AGN? Ready for SIRTF!
SLIDE 22
In general, nice ROSAT era correlations kaput …
SLIDE 23 (mass-velocity dispersion) (mass-bulge luminosity dispersion)
2x disagreement?
vs.
[Yu & Tremaine 2002]
SLIDE 24
Observational Debates & Clues
Rare long-lived AGN vs. many short-lived AGN? Seems to be tilting decisively towards relation, ( many relic SMBH) X-ray/2MASS counts ( many active AGN missed optically) No more Soltan/ problem? M M σ σ − ⇒ ⇒ −
Also, relation BH and galaxy know about each other!? Galaxy & BH formation same process? (Once correct for obscuration, redshift evolution s M σ − ⇒ imilar?) Mergers/gas are clearly important in at least AGN phase.
SLIDE 25 Owen, VLA
3C 75: Merger Starting? Where are the SMBH binaries?
SLIDE 26
“Smoking Gun?”
Ekers & Merrit, 2002 NGC 326
SLIDE 27
NGC 6240: current best case for an eventual merger?
SLIDE 28 Simulation of idealized gas-rich merger… Dynamical friction phase
SLIDE 29
What happens when a binary forms? Drag continues! (If there’s enough gas…)
SLIDE 30
Merger happens very fast!
SLIDE 31
Bender and Pollack 2003
SLIDE 32 Black holes in globular clusters?
One of best studied cases : M15 a 2000 solar mass black hole?
Guhathakurta et al. 1996 Gebhardt et al. 2000
?
SLIDE 33 ULXs and IMBHs M82
Fabbiano et al. 2001, CXO
SLIDE 34
SLIDE 35
How massive were the First Stars?
Previous estimates: 1 Mo < MPopIII < 106 Mo M ~ 106 Mo
Massive Black Hole Cluster of Stars
normal IMF Top-heavy IMF
SLIDE 36 Atomic cooling H_2 cooling
The Physics of Population III
No magnetic fields yet (?) No metals no dust Initial conditions given by CDM
Well-posed problem
How to cool primordial gas?
No metals different cooling Below 104 K, main coolant is H2
Cooling sensitive to H2 abundance H2 formed in non-equilibrium
Have to solve coupled set of rate equations Metals
Tvir for Pop III
SLIDE 37 Cosmological Initial Conditions
- Consider situation at z = 20
Gas density
~ 7 kpc
Primordial Object
SLIDE 38
The First Star-Forming Region
~ 7 kpc 1 kpc M ~ 106 Mo
SLIDE 39
SLIDE 40 A Physical Explanation:
- Gravitational instability
(Jeans 1902)
MJ~T1.5 n-0.5 T vs. n MJ vs. n
- Thermodynamics of primordial gas
- Two characteristic numbers in
microphysics of H2 cooling:
- Tmin ~ 200 K
- ncrit ~ 103 - 104 cm-3 (NLTE LTE)
- Corresponding Jeans mass: MJ ~ 103 Mo
SLIDE 41 The Crucial Role of Accretion
- Final mass depends on accretion from dust-free
Envelope
- Development of core-envelope structure
- Omukai & Nishi 1998 , Ripamonti et al. 2002
Mcore ~ 10-3 Mo very similar to Pop. I
- Accretion onto core very different!
- dM/dtacc ~ MJ / tff ~ T3/2 (Pop I: T ~ 10 K, Pop III: T ~ 300 K)
- Can the accretion be shut off in the absence of dust?
SLIDE 42
SLIDE 43
SLIDE 44
The Death of the First Stars:
(Heger et al. 2002)
Initial Stellar Mass Z Pop III Pop I PISN
SLIDE 45
What happens to pop III remnant BH? Madau et al. 2003?
SLIDE 46 First Dwarf Galaxies as Sites of BH Formation
- 2 sigma peak
- M ~ 108 M0, zvir ~ 10
- Tvir ~ 104 K
Cooling possible due to atomic H
- Photo-dissociation of H2:
H2 + h nu 2 H
h nu = 11.2 – 13.6 eV
T vs. log n Tvir ~ 104 K
(Bromm & Loeb 2003)
SLIDE 47 En Route to a Supermassive Black Hole?
- Consider gas distribution in central 100 pc
Single object: M ~ 106 Mo Low-spin High-spin Binary: M1,2 ~ 106 M0
SLIDE 48
Summary:
SMBH growth must be a rapid and relatively robust process. Can happen very early on. Probably intimately tied to galaxy merger induced activity, especially nuclear star formation (M-σ relation!). Observations of SMBH improving rapidly. Field in state of flux. Chandra + SIRTF especially powerful, overcome obscuration problem to uncover true AGN and star formation activity. SBMH growth by merger vs. accretion? Both? ☺ Today, looks like accretion may be dominant mode. SMBH growth greatly facilitated by pre-existing massive “seeds.” Nature of number of seeds is major uncertainty in expected LISA event rate. Primordial (Pop. III) seeds appear plausible => very high z mini-AGN, WMAP reionization? Mergers and GRBs? High overall LISA rate? Especially if M-σ relation holds for early AGN, LISA powerful probe of early structure formation, at z > 10!?
SLIDE 49 Would be reassuring to actually find some MBH binaries before LISA. Where are they? Maybe binaries don’t accrete efficiently? NGC 6240 currently best and perhaps cautionary example. Binary BH merge quickly and are
- bscured? If so, EM counterpart to LISA signal may be difficult to find?
Angular resolution key to finding correct counterpart. Would also be nice to find IMBH/LMBH. Evidence scant right now. IMHO, most ULXs are NOT IMBH but beamed/super-Eddington stellar mass objects (e.g., GRS 1915 in our galaxy). However, a few ULXs are best explained as massive objects and there are a couple known AGN with Large characteristic Jeans mass for fragmentation might actually occur today, e.g., in galactic nuclei. Strong ionizing radiation field (e.g., from AGN) wipes
- ut metal coolants, heats gas? Top heavy IMF? Most massive stars in our
galaxy near galactic center. More massive remnants + dissipative central gas => good for LISA! Effects of BH spin seems to be major LISA signal uncertainty. If MBH grow by accretion, easy to get maximally rotating hole. Don’t ignore!?
5
10 . M M ≤
SLIDE 50
Metal abundance higher than solar, everything Happens fast Abandon eddington limit Seed Small or not? Blob vs. hierarchical formation? Tightness of M-sigma? Easy to understand scaling Why the constant, ie., why scatter so small?
SLIDE 51
Seeds!!! Lesson from present day star formation Bender “LMBH”
SLIDE 52
Henry 1999
The Diffuse Extragalactic Background
Energetic Particles!
SLIDE 53 From the Dark Ages to the Cosmic Renaissance
- First Stars Transition from Simplicity to Complexity
SLIDE 54
Berger et al. 2003
SLIDE 55 Cooling Rate vs. T
Paradise Lost: The Transition to Population II
(Bromm, Ferrara, Coppi, & Larson 2001, MNRAS, 328, 969)
- Add trace amount of metals
- Limiting case of no H2
- Heating by photoelectric
effect on dust grains
Consider two identical (other than Z) simulations !
SLIDE 56 Effect of Metallicity:
Z = 10-4 Zo Z = 10-3 Zo
- Insufficient cooling
- Vigorous fragmentation
Critical metallicity: Zcrit ~ 5 x 10-4 Zo
SLIDE 57 The First Supernova-Explosion
Gas density
~ 1 kpc
Disruption (PISN)
SLIDE 58 Barger et al. 2003 (CDF-N) Extra objects at wrong redshift(s)!
SLIDE 59 Region of Primordial Star Formation
- Gravitational Evolution of DM
- Gas Microphysic:
- Can gas sufficiently cool?
- tcool < tff (Rees-Ostriker)
- Collapse of First Luminous Objects expected:
- at: zcoll = 20 – 30
- with total mass: M ~ 106 Mo
SLIDE 60 A Tale of Two Timescales
- Gas particles loiter at: n ~ 103 – 104 cm-3
- tcool ~ tff
- Quasi-hydrostatic phase
- Runaway collapse occurs
- s.t. tcool ~ tff
- Consider the cooling and freefall times:
Timescale vs. n
tff tcool
SLIDE 61
The First Supernova Explosions
(with N. Yoshida & L. Hernquist)
~ 7 kpc 1 kpc M ~ 106 Mo
SLIDE 62
HII Regions around the First Stars
1 kpc
SLIDE 63 Simulating the Formation of the First Stars:
(Bromm, Coppi, & Larson and Bromm & Hernquist)
(both DM and gas)
primordial gas
- Non-equilibrium chemistry
- Initial conditions: ΛCDM
- Modifications to SPH:
- sink particles
- particle splitting
SLIDE 64
The First Supernova-Explosion
Metal Distribution
~ 1 kpc
SLIDE 65
Thermodynamics and Structure
T vs. log n Phase Distribution
SLIDE 66
Dense-shell Formation
Timescale vs Radius
tcool tff
Inverse Compton cooling
tshock
SLIDE 67 The First Supernova-Explosion
Gas density
~ 1 kpc
Disruption (PISN)
SLIDE 68 Nucleosynthetic Evidence:
(Qian & Wasserburg 2002)
enrichment at [Fe/H] < -3
higher [Fe/H]
Heavy r-process abund. vs. [Fe/H]
Zcrit
SLIDE 69 Cosmic Star Formation History
(Mackey, Bromm & Hernquist 2003)
- 2 modes of SF:
- Pop III VMS
- Pop I / II normal stars
- Pop III SF possible
in halos with:
- Tvir < 104 K molecular cooling
- Tvir > 104 K
atomic H cooling
Comoving SFR vs. redshift
Pop I / II Pop III
(Springel & Hernquist 2003)
SLIDE 70 Cosmic Star Formation History
(Mackey, Bromm & Hernquist 2003)
expected in halos with:
Tvir > 104K atomic H cooling
suppresses SF in mini-halos
(radiative and mechanical)
Comoving SFR vs. redshift
Pop III
SLIDE 71
The Pop III Pop II Transition
(Mackey, Bromm & Hernquist 2003)
Metallicity SFR vs. redshift
ztran~ 15 - 20 Zcrit
50% 5%
SLIDE 72 Relic of the Dawn of Time:
- HE0107-5240: [Fe/H] = - 5.3 (Christlieb et al. 2002)
- What does this star tell us about Population III ?
SLIDE 73 Metal Poor Halo Stars and the First Stars:
(with Schneider, Ferrara, Salvaterra, & Omukai 2003, Nature in press)
- Abundance pattern:
- core-collapse SN
- PISN
- Break degeneracy:
- r-process elements
- Z < Zcrit ?
- role of dust
- shock-compression
- statistics
SLIDE 74 First Dwarf Galaxies as Sites of BH Formation
- 2 sigma peak
- M ~ 108 M0, zvir ~ 10
- Tvir ~ 104 K
Cooling possible due to atomic H
- Photo-dissociation of H2:
H2 + h nu 2 H
h nu = 11.2 – 13.6 eV
T vs. log n Tvir ~ 104 K
SLIDE 75 Gamma-Ray Bursts as Probes of the First Stars:
massive stars
trace cosmic SFH
- Swift mission:
- Launch in 2003
- Sensitivity
GRBs from z > 15
SLIDE 76 Expected Redshift Distribution of GRBs:
( Bromm & Loeb 2002, ApJ, 575, 111 )
(Cf. Barkana & Loeb 2000, ApJ, 539, 20)
SF History GRB Redshift Distribution
- Fraction of all burst from z > 5: ~ 50%
- Fraction of GRBs detected by Swift from z > 5: ~ 25%
SLIDE 77 Summary
- Primordial gas typically attains:
- T ~ 200 – 300 K
- n ~ 103 – 104 cm-3
- Corresponding Jeans mass: MJ ~ 10 3 Mo
- Pop III SF might have favored very massive stars
- Transition to Pop II driven by presence of metals
(ztrans ~ 15 – 20)
- PISNe completely disrupt mini-halos and enriches
surroundings
- Metal-poor halo stars as probes of the first stars
SLIDE 78 Perspectives:
- Further fate of clumps
- Feedback of protostar on its envelope
- Inclusion of opacity effects (radiative transfer)
- The ``Second Generation of Stars’’
- SN feedback and metal enrichment from the first stars
- How does a VMO evolve and die?
Observability (lensing?) and statistics of high-z SNe
SLIDE 79
132 node Beowulf cluster (AMD Athlon)
SLIDE 80 The Mass of a Population III Star
fall: M ~ 100 Mo
with isothermal density profile First stars were predominantly very massive
SLIDE 81 Implications of a Heavy IMF For the First Stars
(Bromm, Kudritzki, Loeb 2001, ApJ, 552, 464)
- Consider: 100 Mo < M < 1000 Mo (VMO)
- Structure determined by:
- Radiation pressure, Luminosity close to EDDINGTON limit
log L vs. log Teff
Teff ~ 110,000 K lambda peak ~ 250 A (close to He II ionization edge)
SLIDE 82 How Do VMOs Evolve ?
log L vs. log Teff
He ignition
3 x 106 yr
Mass loss ???
SLIDE 83 Spectral Signature
Strong NLTE effects
- Close to black-body form
- Lines of H I and He II
Flux vs. Wavelength
SLIDE 84 A Generic Spectrum
L nu / M vs. lambda
M > 300 Mo
- Predict composite spectrum
almost independent of IMF
production
suffice to reionize the Universe
SLIDE 85 Probing the Primordial IMF with NGST
- Observed spectrum: Heavy IMF vs. Salpeter IMF
Observed flux vs. Wavelength
- Observed spectrum from cluster with heavy
IMF is significantly bluer
Tumlinson & Shull 2000
SLIDE 86 Why Study Population III?
- The Quest for our Origins
- Importance for Cosmological Structure Formation
Reheat / Reionize the Universe Feedback effects on IGM Initial enrichment with metals Pure H/He out of BBNS Need stars to synthesize heavy elements Pop III remnants Baryonic DM (?)
CMB anisotropy probes (WMAP / Planck) Study imprint of first stars IR missions (SIRTF/ JWST) Direct imaging
SLIDE 87 The Crucial Role of Accretion
dM/dt vs. time M vs. time
55 .
d d
−
∝ t t M
45 .
t M ∝
∗
SLIDE 88
Gilli 2003 review (astro-ph/0303115)