Contribution of Radio-mode AGN Heating in Clusters. Cheng-Jiun Ma Brian McNamara, Paul Nulsen (SAO and University of Waterloo) arxiv:1206.

Transcription

1 Contribution of Radio-mode AGN Heating in Clusters Cheng-Jiun Ma Brian McNamara, Paul Nulsen (SAO and University of Waterloo) arxiv:

2 Heating in clusters o o Break of self-similar scaling relation. Cooling flow problem Makevitch+98

3 Continual AGN heating o Preheating (Kaiser 1991) About 1keV/particle can explain deviations from self-similar scaling (Wu et al. 2000). o Continual heating -- Radio mode " Linked to X-ray cavities: It has been used primarily used to explain quenching of cooling flows.

4 Sample: 8 cluster surveys + NVSS 0.1 < z < 0.6 and covered by NVSS ~400 NVSS sources (>2mJy) matched within 250 kpc Ø 8 X-ray cluster surveys ~700 clusters Require an efficient method to estimate AGN jet powers erg/s LX, LX,bol 1044 erg/s Ø 400SD 160SD WARPS MACS BCS NORAS REFLEX CIZA z z 1.0

5 X-ray cavity powers are costly o pdv work to inflate cavity Pcav = ΔH/τ = (E+pV)/τ = 4pV/τ Requires deep, high resolution X-ray data. About 70 cavity systems known Mostly at low z Hlavacek-Larrondo et al. 2012"

6 Jet powers determined from radio Large scatter: Pjet = 0.78dex < P 1.4 < erg/s erg/s < L X < erg/s Mostly FR-I sources Or references: Rafferty et al. 2006, Best et al P cav (10 42 erg/s) Hlavacek-Larrondo et al Cavagnolo et al Birzan et al. 2006, P 1.4 (10 42 erg/s) Log (P jet /10 42 erg/s) = 0.75(±0.14) Log (P 1.4 /10 40 erg/s) (±0.18) Cavagnolo et al. 2010

7 Average jet powers of central AGN Ø <Pjet> doesn t change significantly over L X <- large uncertainty Ø Sensitive to most powerful, rare, radio sources. Ø A net energy gain in clusters with L X < 3x10 44 erg/s <P jet > (10 44 erg/s) 10 1 w 43/85 89/ / / Lx (10 44 erg/s) including powerful radio sources saturated scaling excluding powerful radio sources R500

8 AGN energy deposited per particle preheating rate AGN heating (kev/par) w including powerful radio sources saturated scaling excluding powerful radio sources 0.1 < Z < Lx (10 44 erg/s) R<250 kpc

9 lim 1 dnsrc (> P1.4 )Ncl is number of clusters 0.3) = 1.26, calculated for P (z) in population. The ex- = jet, (3) ϕ(p1.4 ) = cluster lim radio sources for Ncl d log P1.4 pected number of background P jet (z0 ) = each erg s 1. The is subtracted from Nsrc and re may be more than one radio is used calculation of ofavera source with in a cluster. that, in since normalization ϕ e Nsrc (> P1.4 ) is number of radio sources pow- Note gives meanand number of radio sources per 7. cluster it may be reater than P1.4 in cluster population of interest AVERAGE JET PO greater than is number of clusters in population. Theunity. ex- The distribution ϕ(p1.4 ) is plotted in The average Figure 7 for clusters in jet fourpowers, redshift$p ranges, jet %, sh ed number of background radio sourcesleft forpanel each of cluster Figure 8 were using of0.6. Following 3,estimated only clusters witha M 1. from Cumulative power formay be more btracted Nsrc and re than one radio luminosities in range 3 10 < L < erg s accounts for uncertainties in X 0.05 < z < 0.1 used to ce in a cluster. Note that, since normalization of ϕ are included. for ϕ(p1.4 )inare shown (1), only fordistribu powrameters Equation correct forof faint s mean number radio sources per cluster it may Values be ers above a in threshold corresponding to fluxintrinsic limit of 3scatter mjy dices, and large radio sources at Ma et al. ter than unity. The distribution ϕ(p1.4 ) is plotted 1.0 for redshifts z = 0.1, 0.2, 0.4,(1). andnote 0.6. that The distribution, Equation log of higher 7redshifts. panel of Figure for clusters in fourϕ(p redshift ranges, redshift, evolvingdistribution more significantly 1.4 ), increases with for a lognormal is greate Following 3, only clusters with at higher power end, withbyprevious findings 2. Fit a linear evolution 45 1 logconsistent that is given Equation (1). Fo nosities in range < LX < erg s (e.g., Galametz et al. 2009; Hart et al. 2011). using redshift bins of Figure 1 are chosen ncluded. model Values for44ϕ(p1.4 ) are shown Closely only forrelated pow- to ϕ(p define φ(p ), First, number 1.4 ), we jet Φz(Pjet=3x10 erg/s) evenly between bins. $Pjet 0.05<z<0.1 bove a threshold corresponding to flux limit of 3 mjy of radio jets per cluster per unit log P, by replacing radio jet at different redshifts. bin in Figure 1, n this is integrate 0.1<z<0.2 he redshifts z = 0.1, 0.2, 0.4, and 0.6. The distribution, power P in Equation (3) with P calculated from Equa1.4 jet 0.1 range of L to give time averaged X 0.2<z< significantly tion (1), 1.4 ), increases with redshift, evolving more " 0.4<z<0.6 e higher power end, consistent with previous findings 1 dnsrc (> Pjet ) (zi, i [$Pjet(4) ". φ(p ) = $Pjet %int =, Galametz et al. 2009; Hart et al. 2011). jet Ncl d log Pjet it osely related to ϕ(p1.4 ), we define φ(pjet ), number 0.1 radio dio jets per cluster per unit log Pjet, by replacing The cumulative jet power per jets more powerful where tzi is from time interval for reds 44 cluster P (10 erg/s) lim er P1.4 in Equation (3) with Pjet calculated Equa- jet correction factor calculated from Equ than Pfrom jet is n (1), range is 0.1 to 0.4 for lowe!shift to 0.6 for0.1 remainder. In bin 1 dnsrc (> Pjet ) lim Φz (P(4) φ(pjet )Pjet d log Pjet. (5). φ(pjet ) = jet ) = evolution of P seen in rig lim jet Pjet Ncl d log Pjet overwhelmed by large uncertainti scatter in relation of Equation (1) cumulative jet power per cluster from jets more powerful lim ble differences between serendipi Φz(Pjet) (1044 erg/s) z (Pjet) (1044 erg/s) Cumulative jet power per cluster

10 AGN energy deposited per particle preheating rate AGN heating (kev/par) w including powerful radio sources saturated scaling excluding powerful radio sources R<250 kpc Constant hea;ng from z=2 Evolu;on of radio LF from z= Lx (10 44 erg/s)

11 Note 1: Significant accreted mass of black hole from z=2 P jet = ηṁc2, =0.1 è M BH 10 9 M since z = 2 The black hole mass may increase by a factor of 2. Caveat: The accreted mass can be distributed to multiple black holes.

12 Note 2: Jet power is underestimated due to assembling of halos At higher z à less massive à more efficient jet heating Heating energy is conserved in groups, and will be merged into clusters. Hart at al. 2011: Jet powers increase by a factor of 10 from z=0.2 to z=1.2 considering growth of clusters.

13 Summary " A composite cluster sample with about 700 clusters at 0.1 < z < 0.6. " AGN jet power is calculated using correlation between radio power and cavity power in Cavagnolo et al " Average jet power doesn t increase significantly in luminous clusters. " The integrated AGN energy to z=2 may reach 1 kev per particle, energy required in preheating model.