Tracing the growth history of the active BH population with the black hole mass function Andreas Schulze Kavli Institute for Astronomy and Astrophysics at Peking University (KIAA) Galaxies vs Black Holes, Castellammare del Golfo 24.06.2013 A. Bongiorno, I. Gavignaud, J. Silverman, M. Schramm, A. Merloni, G. Zamorani
Black hole - galaxy coevolution M σ relation cosmic BH accretion history parallel to SF history McConnell & Ma (2013) Zheng et al. (2009) link between black hole growth and galaxy evolution how are black holes growing? Andreas Schulze (KIAA) Galaxies vs Black Holes 2 / 16
How can we trace black hole growth? AGN Luminosity function is main demographic quantity Physical quantities of black holes: black hole mass M accretion rate / Eddington ratio λ = L bol /L Edd Andreas Schulze (KIAA) Galaxies vs Black Holes 3 / 16
How can we trace black hole growth? AGN Luminosity function is main demographic quantity Physical quantities of black holes: black hole mass M accretion rate / Eddington ratio λ = L bol /L Edd L λm implies degeneracy between M and λ additional M information able to break this degeneracy Andreas Schulze (KIAA) Galaxies vs Black Holes 3 / 16
How can we trace black hole growth? AGN Luminosity function is main demographic quantity Physical quantities of black holes: black hole mass M accretion rate / Eddington ratio λ = L bol /L Edd L λm implies degeneracy between M and λ additional M information able to break this degeneracy Active black hole mass function - Φ (M ) Eddington ratio distribution function - Φ λ (λ) well-defined AGN sample "virial" black hole mass estimates from single-epoch spectra: M L γ contfwhm 2 Andreas Schulze (KIAA) Galaxies vs Black Holes 3 / 16
Application of the BHMF Constraints on theoretical models comparison with galaxy evolution models Evolution of BH-bulge relations requires knowledge of distribution functions Hirschmann et al. (2012) Schulze & Wisotzki (2011) Andreas Schulze (KIAA) Galaxies vs Black Holes 4 / 16
Determination of the distribution functions binned DF via luminosity weighted V max + non-parametric + model independent - unaccounted sample incompleteness - BHMF not complete to lowest mass bin Andreas Schulze (KIAA) Galaxies vs Black Holes 5 / 16
Determination of the distribution functions binned DF via luminosity weighted V max + non-parametric + model independent - unaccounted sample incompleteness - BHMF not complete to lowest mass bin model DF via fitting of bivariate distribution function of M and λ Black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) determined jointly by fitting probability distribution in M λ-plane via Maximum likelihood method (Schulze & Wisotzki 2010) or via Bayesian framework (Kelly et al. 2009) Andreas Schulze (KIAA) Galaxies vs Black Holes 5 / 16
The local active black hole mass function and Eddington ratio distribution function Local (z < 0.3) BHMF and ERDF from the Hamburg/ESO Survey Active black hole mass function Eddington ratio distribution function Schulze & Wisotzki (2010) No evidence for downturn at low black hole mass or at low Eddington ratio Andreas Schulze (KIAA) Galaxies vs Black Holes 6 / 16
The local active black hole mass function and Eddington ratio distribution function Local (z < 0.3) BHMF and ERDF from the Hamburg/ESO Survey Active black hole mass function Eddington ratio distribution function Schulze & Wisotzki (2010) No evidence for downturn at low black hole mass or at low Eddington ratio Andreas Schulze (KIAA) Galaxies vs Black Holes 6 / 16
Active fraction of local black holes compare to quiescent BHMF of Marconi et al. 2004 log Φ [Mpc 3 logm 1 ] -2-8 total active active fraction -1-2 -3 6 7 8 9 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 7 / 16
Active fraction of local black holes compare to quiescent BHMF of Marconi et al. 2004 significant decrease of active fraction toward higher M indication for cosmic downsizing in black hole mass log Φ [Mpc 3 logm 1 ] active fraction -2-8 -1-2 -3 total active 6 7 8 9 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 7 / 16
Active BHMF and ERDF at higher redshifts at z > 0.4 BHMF and ERDF determined from SDSS QSO sample evidence for black hole mass downsizing only high mass end of BHMF, high λ end of ERDF Kelly & Shen (2013) Andreas Schulze (KIAA) Galaxies vs Black Holes 8 / 16
BH demographics from VVDS, zcosmos and SDSS combine bright, large area surveys (SDSS) with deep, small area AGN surveys (VVDS, zcosmos) SDSS: i < 19.1 Ω eff = 6248 deg 2 color selection VVDS: wide: I AB < 22.5 Ω eff = 4.5 deg 2 deep: I AB < 24.0 Ω eff = 0.6 deg 2 random selection zcosmos: I AB < 22.5 Ω eff = 1.6 deg 2 random + X-ray selection Andreas Schulze (KIAA) Galaxies vs Black Holes 9 / 16
BH demographics from VVDS, zcosmos and SDSS 1.0<z<2.15 use MgII BH masses SDSS: 31600 AGN log L bol [erg s 1 ] 47 46 45 44 log λ 0-1 -2 43 7 8 9 10 log M [M ] 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 10 / 16
BH demographics from VVDS, zcosmos and SDSS 1.0<z<2.15 use MgII BH masses SDSS: 31600 AGN VVDS: 91 + 62 AGN log L bol [erg s 1 ] 47 46 45 44 log λ 0-1 -2 43 7 8 9 10 log M [M ] 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 10 / 16
BH demographics from VVDS, zcosmos and SDSS 1.0<z<2.15 use MgII BH masses SDSS: 31600 AGN VVDS: 91 + 62 AGN zcosmos: 163 AGN log L bol [erg s 1 ] 47 46 45 44 log λ 0-1 -2 43 7 8 9 10 log M [M ] 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 10 / 16
The active black hole mass function at 1 < z < 2 Preliminary results! log Φ -5 SDSS DR7-7 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 11 / 16
The active black hole mass function at 1 < z < 2 Preliminary results! log Φ -5 SDSS DR7 VVDS -7 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 11 / 16
The active black hole mass function at 1 < z < 2 Preliminary results! log Φ -5 SDSS DR7 VVDS zcosmos (prelim.) -7 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 11 / 16
The active black hole mass function at 1 < z < 2 Preliminary results! log Φ -5-7 SDSS DR7 VVDS zcosmos (prelim.) SXDS (Nobuta+12) 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 11 / 16
The active black hole mass function at 1 < z < 2 Preliminary results! log Φ -5-7 7 8 9 10 log M [M ] SDSS DR7 VVDS zcosmos (prelim.) SXDS (Nobuta+12) combined Andreas Schulze (KIAA) Galaxies vs Black Holes 11 / 16
The Eddington ratio distribution function at 1 < z < 2 Preliminary results! log Φ λ -5 SDSS DR7-7 -2-1 0 log λ Andreas Schulze (KIAA) Galaxies vs Black Holes 12 / 16
The Eddington ratio distribution function at 1 < z < 2 Preliminary results! log Φ λ -5 SDSS DR7 VVDS -7-2 -1 0 log λ Andreas Schulze (KIAA) Galaxies vs Black Holes 12 / 16
The Eddington ratio distribution function at 1 < z < 2 Preliminary results! log Φ λ -5 SDSS DR7 VVDS zcosmos (prelim.) -7-2 -1 0 log λ Andreas Schulze (KIAA) Galaxies vs Black Holes 12 / 16
The Eddington ratio distribution function at 1 < z < 2 Preliminary results! log Φ λ -5-7 SDSS DR7 VVDS zcosmos (prelim.) SXDS (Nobuta+12) -2-1 0 log λ Andreas Schulze (KIAA) Galaxies vs Black Holes 12 / 16
The Eddington ratio distribution function at 1 < z < 2 Preliminary results! log Φ λ -5-7 -2-1 0 log λ SDSS DR7 VVDS zcosmos (prelim.) SXDS (Nobuta+12) combined Andreas Schulze (KIAA) Galaxies vs Black Holes 12 / 16
Evolution of the active black hole mass function and Eddington ratio distribution function comparison with local distribution functions AGN LF active BHMF ERDF log Φ [Mpc 3 mag 1 ] -8-10 -20-22 M B z =0-24 -26 [mag] -28 log Φ -5-7 z =0 7 8 9 10 log M [M ] log Φ λ -5-7 z =0-2 -1 0 log λ strong downsizing in the active BHMF only weak evolution in the ERDF Andreas Schulze (KIAA) Galaxies vs Black Holes 13 / 16
Conclusions active BHMF and ERDF provide additional observational constraints on BH growth and galaxy evolution active BHMF and ERDF have to be determined jointly as bivariate distribution function established local zero point from Hamburg/ESO Survey active BHMF and ERDF at z 1.5 determined from VVDS, zcosmos and SDSS downsizing in AGN LF driven by downsizing in the BHMF Andreas Schulze (KIAA) Galaxies vs Black Holes 14 / 16
Comparison with previous work on SDSS this work Max. likelihood Shen & Kelly (2012) Bayesian Kelly & Shen (2013) Bayesian log Φ -5 log Φ λ -5-7 -7 7 8 9 10 log M [M ] -2-1 0 log λ maximum likelihood estimate and Bayesian results consistent Andreas Schulze (KIAA) Galaxies vs Black Holes 15 / 16
Comparison with previous work on SDSS this work Max. likelihood Shen & Kelly (2012) Bayesian Kelly & Shen (2013) Bayesian log Φ -5 log Φ λ -5-7 -7 7 8 9 10 log M [M ] -2-1 0 log λ maximum likelihood estimate and Bayesian results consistent Andreas Schulze (KIAA) Galaxies vs Black Holes 15 / 16
Active black hole fraction at z 1.5 compare to quiescent BHMF derived from stellar mass function log Φ [Mpc 3 logm 1 ] -2-8 total active apparent increase of broad line AGN active fraction toward higher M active fraction 0-1 -2-3 7 8 9 10 log M [M ] Andreas Schulze (KIAA) Galaxies vs Black Holes 16 / 16