Time lags and reverberation in the lamp-post geometry of the compact - - PowerPoint PPT Presentation

▶
time lags and reverberation in the lamp post geometry of
SMART_READER_LITE
LIVE PREVIEW

Time lags and reverberation in the lamp-post geometry of the compact - - PowerPoint PPT Presentation

Nov 17, 2022 169 likes •425 views

Time lags and reverberation in the lamp-post geometry of the compact corona illuminating a black-hole accretion disc Michal Dov ciak Astronomical Institute Academy of Sciences of the Czech Republic, Prague Astrophysics group seminar,

slide-1
SLIDE 1

Time lags and reverberation in the lamp-post geometry of the compact corona illuminating a black-hole accretion disc

Michal Dovˇ ciak

Astronomical Institute Academy of Sciences of the Czech Republic, Prague Astrophysics group seminar, Department of Physics, University of Crete, Heraklion, Greece 30th May 2014

slide-2
SLIDE 2

Acknowledgement

StrongGravity (2013–1017) — EU research project funded under the Space theme of the 7th Framework Programme on Cooperation Title: Probing Strong Gravity by Black Holes Across the Range of Masses Institutes: AsU, CNRS, UNIROMA3, UCAM, CSIC, UCO, CAMK Webpages: http://stronggravity.eu/ http://cordis.europa.eu/projects/rcn/106556 en.html

slide-3
SLIDE 3

Active Galactic Nuclei – scheme

Urry C. M. & Padovani P . (1995) Unified Schemes for Radio-Loud Active Galactic Nuclei PASP , 107, 803

slide-4
SLIDE 4

Active Galactic Nuclei – X-ray spectrum

Fabian A.C. (2005) X-ray Reflections on AGN, in proceedings of “The X-ray Universe 2005”, El Escorial, Madrid, Spain, 26-30/9/2005

slide-5
SLIDE 5

Active Galactic Nuclei – lags

104 103 100 50 50 100 150 200 Ν Hz TLΝ s PLΝ0.62 Ν0.61 s M1.8106 M Α0.32 Θ27° h2.4 rg 1H 0707495 104 103 500 500 1000 Ν Hz TLΝ s PLΝ0.064Ν1.0 s M7.1106 M Α0.91 Θ51° h8.8 rg ESO 113G010

Emmanoulopoulos et al. (2014) General relativistic modelling of the negative reverberation X-ray time delays in AGN, MNRAS 439 3931

slide-6
SLIDE 6

Active Galactic Nuclei – lags

Lag Energy (keV) 1 10

Kara et al. (2014) The curious time lags of PG 1244+026: Discovery of the iron K reverberation lag, MNRAS 439 L26

slide-7
SLIDE 7

References

◮ Blandford & McKee (1982) ApJ 255 419 → reverberation of BLR

Lo(v,t) =

∞

−∞ dt′ Lp(t′)Ψ(v,t −t′) ◮ Stella (1990) Nature 344 747 → time dependent Fe Kα shape (a = 0) ◮ Matt & Perola (1992) MNRAS 259 433

→ Fe Kα response and black hole mass estimate → t ∼ GM/c3 → time dependent light curve, centroid energy and line equivalent width (h = 6, 10; a = 0; θo)

◮ Campana & Stella (1995) MNRAS 272 585

→ line reverberation for a compact and extended source (a = 0)

◮ Reynolds, Young, Begelman & Fabian (1999) ApJ 514 164

→ fully relativistic line reverberation (h = 10; a = 0, 1) → more detailed reprocessing, off-axis flares → ionized lines for Schwarzschild case, outward and inward echo, reappearance of the broad relativistic line

slide-8
SLIDE 8

References

More recent references:

◮ Chainakun & Young (2012) MNRAS 420 1145

→ fully relativistic, lamp-post geometry, ionized accretion disc

◮ Wilkins & Fabian (2013) MNRAS 430 247

→ fully relativistic, extended corona, propagation effects

◮ Cackett et al. (2014) MNRAS 438 2980

→ Fe Kα reverberation in lamp-post model (spectral line)

◮ Emmanoulopoulos et al. (2014) MNRAS 439 3931

→ Fe Kα reverberation in lamp-post model (soft excess)

slide-9
SLIDE 9

Why to study toy model of lamp-post geometry?

◮ physical motivation: base of a jet ◮ many effects should be qualitatively similar with this simple geometry ◮ it can give us certain limits on the model (e.g. spin versus height) ◮ we can easily explore the dependence on many parameters

(height of the corona, ionization of the disc, ...)

◮ if we want to study the dependence on geometry, we should know how

  • ther parameters influence the results

(i.e. Is the idea of measuring geometry of the corona via reverberation feasible?)

slide-10
SLIDE 10

Scheme of the lamp-post geometry

r

in

r

  • ut

δi δe

a

Ω

corona accretion disc black hole

  • bserver

∆Φ

h M

◮ central black hole – mass, spin ◮ compact corona with isotropic emission

→ height, photon index

◮ accretion disc

→ Keplerian, geometrically thin, optically thick → ionisation due to illumination (Lp, h, M, a, nH, qn)

◮ local re-processing in the disc

→ REFLIONX with different directional emissivity prescriptions

◮ relativistic effects:

→ Doppler and gravitational energy shift → light bending (lensing) → aberration (beaming) → light travel time

slide-11
SLIDE 11

Time delay

5 10 15 20 25 30 35 5 10 15 t [GM/c3] r [GM/c2] a = 1 height 3 6 15

−10 −5 5 10 x −10 −5 5 10 y

Time delay

(a = 1.0000 , θo = 60°)

−20 −10 10 20 −10 −5 5 10 x −10 −5 5 10 y

Total time delay

(a = 1.0000 , θo = 60°, h = 3)

8 16 24 32

◮ total light travel time includes the lamp-to-disc and disc-to-observer part ◮ first photons arrive from the region in front of the black hole

which is further out for higher source

◮ contours of the total time delay shows the ring of reflection that develops

into two rings when the echo reaches the vicinity of the black hole

slide-12
SLIDE 12

Light curve

0.05 0.1 0.15 0.2 0.25 10 20 30 40 Photon flux t [GM/c3] a = 0, h = 3, ∆t = 1 inclination 30° 60° 85° 0.5 1 1.5 2 2.5 10 20 30 40 Photon flux t [GM/c3] a = 1, h = 3, ∆t = 1 inclination 30° 60° 85° ◮ the flux for Schwarzschild BH is much smaller than for Kerr BH due to the hole

below ISCO (no inner ring in Schwarzschild case)

◮ the shape of the light curve differs substantially for different spins ◮ the “duration” of the echo is quite similar ◮ the higher the inclination the sooner first photons will be observed ◮ magnification due to lensing effect at high inclinations

slide-13
SLIDE 13

Dynamic spectrum – spectral line

◮ signature of outer and inner

echo in dynamic spectra

◮ large amplification when

the two echos separate

◮ intrinsically narrow Kα line

can acquire weird shapes

slide-14
SLIDE 14

Dynamic spectrum – ionised disc

E2 ×F(E)

slide-15
SLIDE 15

Definition of the phase lag

Frefl(E,t) = Np(t)∗ψ(E,t) ⇒ ˆ F refl(E,f) = ˆ Np(f). ˆ ψ(E,f) where ˆ ψ(E,f) = A(E,f)eiφ(E,f) if Np(t) = cos(2πft) and ˆ ψ(E) = A(E)eiφ(E) then Frefl(E,t) = A(E)cos{2πf[t +τ(E)]} where τ(E) ≡ φ(E) 2πf F(E,t) ∼ Np(t)∗(ψr(E,t)+δ(t)) ⇒ ˆ F(E,f) ∼ ˆ Np(f).( ˆ ψr(E,f)+1) and tanφtot(E,f) = Ar(E,f)sinφr(E,f) 1+Ar(E,f)cosφr(E,f)

slide-16
SLIDE 16

Parameter values and integrated spectrum

100 1000 0.1 1 10 E2 x F(E) E [keV] a = 1, h = 3, θo = 30° reflected primary

M = 108M⊙ a = 1 (0) θo = 30◦ (60◦) h = 3 (1.5, 6, 15, 30) Lp = 0.001LEdd Γ = 2 (1.5, 3) nH = 0.1 (0.01, 50, 5, 0.2)×1015cm−3 qn = −2 (0, −5, −3) Energy bands: soft excess: 0.3−0.8 keV primary: 1−3 keV iron line: 3−9 keV Compton hump: 15−40 keV

slide-17
SLIDE 17

Phase lag dependence on geometry

0.02 0.04 0.06 0.08 0.1 0.12 10 20 30 40 50 relative photon flux t [ks] a = 1, θo = 30° height [GM/c2] 1.5 3 6 15 30 0.05 0.1 0.15 0.2 0.25 10 20 30 40 50 relative photon flux t [ks] a = 1, θo = 60° height [GM/c2] 1.5 3 6 15 30

  • 300
  • 250
  • 200
  • 150
  • 100
  • 50

50 1 10 100 phase lag [s] f [µHz] a = 1, θo = 30° height [GM/c2] 1.5 3 6 15 30

  • 160
  • 140
  • 120
  • 100
  • 80
  • 60
  • 40
  • 20

20 1 10 100 phase lag [s] f [µHz] a = 1, θo = 60° height [GM/c2] 1.5 3 6 15 30

◮ reflected photon flux decreases with height ◮ primary flux increases with height ◮ the delay of response is increasing with height ◮ the “duration” of the response is longer ◮ the phase lag increases with height, it depends mainly on the “average” response time and magnitude of relative photon flux ◮ the phase lag null points are shifted to lower frequencies for higher heights due to longer timescales of response ◮ relative photon flux and the phase lag increase with inclination for low heights ◮ the delay and duration of response do not change much with the inclination and thus the phase lag null points frequencies change

  • nly slightly
slide-18
SLIDE 18

Phase lag dependence on spin and energy band

0.005 0.01 0.015 0.02 5 10 15 20 relative photon flux t [ks] a = 0, h = 3, θo = 60° energy band 0.3-0.8 keV 1-3 keV 3-9 keV 15-40 keV 0.05 0.1 0.15 0.2 5 10 15 20 relative photon flux t [ks] a = 1, h = 3, θo = 60° energy band 0.3-0.8 keV 1-3 keV 3-9 keV 15-40 keV

  • 50
  • 45
  • 40
  • 35
  • 30
  • 25
  • 20
  • 15
  • 10
  • 5

5 1 10 100 phase lag [s] f [µHz] a = 0, h = 3, θo = 60° energy band 0.3-0.8 keV 3-9 keV 15-40 keV

  • 100
  • 80
  • 60
  • 40
  • 20

20 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 60° energy band 0.3-0.8 keV 3-9 keV 15-40 keV

◮ the relative flux in the energy

band where primary dominates may in some cases be larger than that in Kα and Compton hump energy bands

◮ the magnitude of the phase

lag in different energy bands differs (in extreme Kerr case the larger lag in SE is due to larger ionisation near BH)

◮ the magnitude of the phase

lag is smaller in Schwarzschild case due to the hole in the disc under the ISCO

◮ the null points of the phase lag

change only slightly with energy and spin

slide-19
SLIDE 19

Ionisation

1 10 100 1000 10000 1 10 100 ξ r [GM/c2]

  • 15
  • 10
  • 5

5 10 15 20 25 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 30° ionisation 1 2 3 4 5

◮ the phase lag in Kα band is shown ◮ the reflection component of the spectra are steeper for higher ionisation ◮ the magnitude of the phase lag depend on ionisation ◮ the null points of the phase lag does not change with the ionisation

slide-20
SLIDE 20

Directionality and photon index

  • 80
  • 70
  • 60
  • 50
  • 40
  • 30
  • 20
  • 10

10 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 30° directionality darkening isotropic brightening

  • 120
  • 100
  • 80
  • 60
  • 40
  • 20

20 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 30° Np:Nr 0.5 1 2

  • 80
  • 70
  • 60
  • 50
  • 40
  • 30
  • 20
  • 10

10 20 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 30° photon index 1.5 2 3

◮ the phase lag in SE band is shown ◮ the magnitude of the phase lag changes in all three cases ◮ the null points of the phase lag does not change with different directionility

dependences or power-law photon index

slide-21
SLIDE 21

Phase lag energy dependence

for low f: Ar(E,f) ≃ AE(E)Af(f) φr(E,f) ≃ φr(f) and ∆τ(E,f) ≃ 1 2πf atan [Ar(E,f)−Ar(E0,f)] sinφr(f) 1+[Ar(E,f)+Ar(E0,f)] cosφr(f)+Ar(E,f)Ar(E0,f) and for f such that φr(f) = ±π 2 : ∆τ(E,f) ≃ 1 2πf [Ar(E,f)−Ar(E0,f)]

slide-22
SLIDE 22

Phase lag energy dependence

  • 100
  • 90
  • 80
  • 70
  • 60
  • 50
  • 40
  • 30
  • 20
  • 10

10 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 30° energy band 0.3-0.8 keV 3-9 keV 15-40 keV

  • 100
  • 80
  • 60
  • 40
  • 20

20 1 10 100 phase lag [s] f [µHz] a = 1, h = 3, θo = 60° energy band 0.3-0.8 keV 3-9 keV 15-40 keV 0.001 0.01 0.1 0.1 1 10 phase lag [s] E [keV] a = 1, h = 3, θo = 30° f = 40µHz f = 130µHz E2 x F(E) 0.001 0.01 0.1 0.1 1 10 phase lag [s] E [keV] a = 1, h = 3, θo = 60° f = 40µHz f = 120µHz E2 x F(E)

◮ the energy dependence of

the phase lag follows the spectral shape perfectly at particular frequencies

slide-23
SLIDE 23

Phase lag energy dependence

  • 70
  • 60
  • 50
  • 40
  • 30
  • 20
  • 10

10 1 10 100 phase lag [s] f [µHz] a = 0, h = 3, θo = 30° energy band 0.3-0.8 keV 3-9 keV 15-40 keV

  • 50
  • 45
  • 40
  • 35
  • 30
  • 25
  • 20
  • 15
  • 10
  • 5

5 1 10 100 phase lag [s] f [µHz] a = 0, h = 3, θo = 60° energy band 0.3-0.8 keV 3-9 keV 15-40 keV 0.0001 0.001 0.01 0.1 1 10 phase lag [s] E [keV] a = 0, h = 3, θo = 30° f = 30µHz f = 110µHz E2 x F(E)

  • 0.001

0.001 0.002 0.003 0.004 0.005 0.006 0.1 1 10 phase lag [s] E [keV] a = 0, h = 3, θo = 60° f = 30µHz f = 100µHz E2 x F(E)

◮ if the second phase lag

maximum is too small the phase lag energy dependence does not follow the spectral shape that well

slide-24
SLIDE 24

Summary

◮ the frequency dependence of the phase lag is mainly due to geometry

(height of the corona)

◮ the magnitude of the phase lag depends on many details of the model

(height, spin, ionisation, unisotropy, energy, ...)

◮ lag versus energy follows the spectral shape at the right frequencies ◮ extended corona

→ brings several new parameters (size, propagation speed, “ignition” position, inhomogeinities) → broadens the response of the disc