AGN Feedback. Andrew King. Dept of Physics and Astronomy, University of Leicester, UK Anton Pannekoek Institute, Amsterdam, NL Sterrewacht Leiden, NL

Transcription

1 AGN Feedback Andrew King Dept of Physics and Astronomy, University of Leicester, UK Anton Pannekoek Institute, Amsterdam, NL Sterrewacht Leiden, NL Potsdam, 2018

2 SMBH affects galaxy bulge more than enough energy to unbind bulge only a fraction used galaxy must notice presence of hole

3 1. how do SMBH grow? Soltan => gas accretion (low z) disc formation is unavoidable all accreting gas has enough angular momentum to orbit the hole, so a disc always forms large disc mass => fragmentation, star formation, mass loss... disc probably never in steady state: Bondi is not a good estimate

4 2. disc accretion is slow spreads on viscous timescale t visc = R2 = 1 R H 2 t dyn where t dyn is the dynamical timescale R/v K =(R 3 /GM) 1/2 this is long: t visc ' yr for R 1pc can we get gas closer in - cancel angular momentum? either borrow some a.m. from SMBH (via Lense-Thirring) to cancel gas a.m. (`disc tearing ), or use radial SMBH feedback to give energy but not a.m. => eccentric => collisions.

5 SMBH affects galaxy bulge more than enough energy to unbind bulge only a fraction used galaxy must notice presence of hole

6 SMBH affects galaxy bulge more than enough energy to unbind bulge only a fraction used galaxy must notice presence of hole

7 SMBH affects galaxy bulge more than enough energy to unbind bulge only a fraction used galaxy must notice presence of hole how?

8 PG (Pounds & Reeves, 2009) P Cygni profile of iron K- alpha: outflow with v ' 0.1c `ultrafast outflow -- `UFO : measured ionization parameter => Ṁ out =4 bm p NR 2 v 1M yr 1 ṀEdd

9 R in intermittent outflow observed X ray column fixed by inner boundary of flow ṁ 3 N H ' b (R in/100r s ) cm 2 so if outflow stopped a time t o ago, we have R in t o ' 0.2 ṁ3 M 8 b N 23 yr recent! 9

10 a continuous Eddington wind would be Compton thick outflows are variable on timescales of weeks or less!

11 UFOs establish M relation by momentum driving at M = M outflows become energy-driven huge increase in lengthscales: pc to kpc most powerful feedback is Compton-thick : invisible? AGN driving is highly variable - correlations tricky! do shocks cool? - rarity of NGC 4051? UFO - molecular outflow correlations.? what determines the SMBH feeding rate?

12 Fig. 12. Correlation between the kinetic power of the outflow and the AGN bolometric luminosity. Symbols and colour-coding as in Fig. 8. The grey line represents the theoretical expectation of models of AGN feedback, for which P K,OF = 5%L AGN.Thereddashedlinerepresents the linear fit to our data, excluding the upper limits. The error bar shown at the bottom-right of the plot corresponds to an average error of ±0.5 dex. Cicone

13

14 outflow shock outflow must collide with bulge gas, and shock what happens? either (a) shocked gas cools: or (b) shocked gas does not cool: `momentum driven flow negligible thermal pressure `energy driven flow thermal pressure > ram pressure Compton cooling by quasar radiation field very effective out to large bulge radii (cf Ciotti & Ostriker, 1997, 2001) expansion into bulge gas is driven by momentum

15 swept-up ambient gas, mildly shocked wind shock outer shock driven into ambient gas SMBH ambient gas Eddington wind, v 0.1c

16 outflow dynamics fixed by cooling close to quasar shocked wind gas cooled by inverse Compton effect (`momentum driven flow ) strong evidence for cooling shock: ionization parameter decreases with outflow velocity, conserving mass flow rate NGC4051: 10x decrease in v, seen in 14 species (Pounds et al.), correlates with ionization

17 wind from SMBH shock structure cooling shocked wind snowplough interstellar gas temperature density velocity SMBH inner (wind) shock contact discontinuity outer (ISM) shock (King, 2010)

18 NGC 4051, Pounds & Vaughan, 2011

19 NGC 4051, Pounds & Vaughan,

20 Eddington outflows momentum outflow rate speed Ṁ out v = L Edd c v = L Edd Ṁ out c = c ṁ 0.1c energy outflow rate where ṁ = Ṁout/ṀEdd Ṁoutv = 2. c2 Ṁ out = 2 L Edd ' 0.05L Edd where ṁ = Ṁout/ṀEdd 1

21 force balance total mass (dark, stars, gas) inside radius R of unperturbed bulge is M tot (R) = 2 2 R G but swept-up gas mass M(R) = 2f g G forces on shell are wind ram pressure: 2 R F = L Edd c = 4 GM apple and weight of gas within R, ram pressure balances weight M = M W = GM(R)M tot(r) R 2 (F = W ) when = 4f g G = f gapple G 2 4 ' M M 8 = M/10 8 M, 200 = /200 km s 1 (K, 2003, 2005)

22 transition to energy-driven flow once M reached close to quasar shocked gas cooled by inverse Compton effect (momentum-driven flow) but once M>M, R can exceed R C : wind shock no longer cools wind shock is adiabatic: hot postshock gas does P dv work on surroundings : energy-driven outflow dramatic change of lengthscales: momentum-driven flows are confined to ~ few pc (micro) energy-driven flows can be ~ kpc (meso) cf K & Pounds, ARAA, 2015

23 density contrast => energy-driven outflow shock may be Rayleigh-Taylor unstable two phase medium: gamma rays and molecular emission mixed large--scale high speed molecular outflows, e.g. Mrk 231:

24 density contrast => energy-driven outflow shock may be Rayleigh-Taylor unstable two phase medium: gamma rays and molecular emission mixed large--scale high speed molecular outflows, e.g. Mrk 231: very cool gas at high speeds

25 outer shock runs ahead of contact discontinuity into ambient ISM: velocity jump across it is a factor ( + 1)/( 1): fixes velocity as v out = +1 2 Ṙ ' / /3 lfc km s 1 f g outflow rate of shocked interstellar gas is Ṁ out = dm(r out) dt = ( + 1)f g 2 G Ṙ Ṁ out ' /3 200 l1/3 M yr 1 (Zubovas & K, 2012

26 Eddington outflows momentum outflow rate speed Ṁ out v = L Edd c v = L Edd Ṁ out c = c ṁ 0.1c energy outflow rate where ṁ = Ṁout/ṀEdd Ṁoutv = 2. c2 Ṁ out = 2 L Edd ' 0.05L Edd where ṁ = Ṁout/ṀEdd 1

27 Velocity / kms energy--driven outflows rapidly converge to apple 1/3 2 fc v e ' 2 c ' 925 2/ f (f c/f g ) 1/3 km s 1 g Time / yr 2000 high velocity outflow at large radius also for other potentials: Zubovas & King, 2012b Velocity / kms Radius / kpc **

28 observed ratios > 1 are possible! Velocity / kms energy--driven outflows rapidly converge to apple 1/3 2 fc v e ' 2 c ' 925 2/ f (f c/f g ) 1/3 km s 1 g and persist even after central quasar turns off! high velocity outflow at large radius Velocity / kms Time / yr also for other potentials: Zubovas & King, 2012b Radius / kpc **

29 observed ratios > 1 are possible! Velocity / kms energy--driven outflows rapidly converge to apple 1/3 2 fc v e ' 2 c ' 925 2/ f (f c/f g ) 1/3 km s 1 g and persist even after central quasar turns off! high velocity outflow at large radius Velocity / kms Time / yr also for other potentials: Zubovas & King, 2012b 500 AGN variability changes scaling! Radius / kpc **

30 spirals: outflow pressure => star formation in disc expanding shocked bulge gas galaxy disc bulge outflow pressurizes central disc, and stimulates star formation

31 observational picture (Tachella )

32 outflows must be episodic, as AGN driving is variable K & Pringle 2007 `chaotic accretion : each accretion disc event limited by self-gravity to a mass M d. H R M BH ' 10 3 M BH so characteristic variation (`flicker ) timescale is t var M d Ṁ HM BH RṀ 105 yr duty cycle yr (most galaxies are not AGN, but all have SMBH) (K & Pringle, 2006; K & Nixon 2015; Schawinski + 15) progress of outflow may be slower than measured velocity

33 outflows must be episodic, as AGN driving is variable

34 swept-out cavity piled-up ISM intermittent shells

35 swept-out cavity piled-up ISM shock cooling events as shells arrive intermittent shells

36 outflows collide with swept-up ISM gas, and shock but shell time of flight R/v 5pc/0.1c 150 yr so incidence of shocks reflects activity of AGN in the past duration of shocks reflects duty cycle of AGN in the past similarly UFO molecular outflow connection: AGN can vary, but not outflow