X-ray reflection from ionised accretion discs a new XSPEC model - - PowerPoint PPT Presentation

x ray reflection from ionised accretion discs a new xspec
SMART_READER_LITE
LIVE PREVIEW

X-ray reflection from ionised accretion discs a new XSPEC model - - PowerPoint PPT Presentation

X-ray reflection from ionised accretion discs a new XSPEC model Michal Dov ciak Astronomical Institute Academy of Sciences of the Czech Republic Prague Ji r Svoboda Matteo Guainazzi European Space Astronomy Centre,


slide-1
SLIDE 1

X-ray reflection from ionised accretion discs – a new XSPEC model

Michal Dovˇ ciak

Astronomical Institute Academy of Sciences of the Czech Republic Prague

Jiˇ r´ ı Svoboda Matteo Guainazzi

European Space Astronomy Centre, Villafranca del Castillo

From the Dolomites to the event horizon: Sledging down the Black Hole potential well Sesto Val Pusteria, 15th – 19th July 2013

slide-2
SLIDE 2

StrongGravity logo explanation

slide-3
SLIDE 3

Scheme of the lamp-post geometry

r

in

r

  • ut

δi δe

a

corona accretion disc black hole

  • bserver

∆Φ

h M

◮ spin a ◮ inclination θo ◮ height h ◮ photon index Γ ◮ luminosity L/Ledd ◮ mass M/M8

(M8 = 108 M⊙)

◮ density nH

slide-4
SLIDE 4

Illumination geometry

10-12 10-10 10-8 10-6 10-4 10-2 100 1 10 100 1000 Ninc [arbitrary units] r [GM/c2] a = 1.0, Γ = 2.0

Heights 100 15 6 3 2 1.5 1.1

1 2 3 4 5 6 7 8 1 10 100 1000 q r [GM/c2] a = 1.0, Γ = 2.0 Heights 1.1 1.5 2 3 6 15 100

◮ Wilkins DR & Fabian AC (2011) MNRAS, 414, 1269 ◮ Svoboda J, Dovˇ ciak M, Goosmann RW, Jethwa P , Karas V, Miniutti G & Guainazzi M (2012) A&A, 545, A106 ◮ Wilkins DR & Fabian AC (2012) MNRAS, 424, 1284

slide-5
SLIDE 5

Emission directionality

a = 1, θ = 30◦, Γ = 2 incident emission angle µi angle µe

4 8 12 16 20 1 10 100 h [GM/c2] r [GM/c2] 1 0.2 0.4 0.6 0.8 1 −10 −5 5 10 x −10 −5 5 10 y 0.00 0.25 0.50 0.75 1.00

◮ Svoboda J, Dovˇ ciak M, Goosmann RW & Karas V (2009) A&A, 507, 1 emission directionality M (µi,µe)

slide-6
SLIDE 6

Emission directionality

a = 1, θo = 30◦, h = 3, Γ = 2 G — transfer function M — angular directionality

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

0.00 0.25 0.50 0.75 1.00

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

3 4 5 6 7

relativistic effects local re-processing

slide-7
SLIDE 7

Emission directionality

a = 1, θo = 30◦, h = 3, Γ = 2 G — transfer function M — angular directionality

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

0.00 0.25 0.50 0.75 1.00

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

3 4 5 6 7

relativistic effects local re-processing

−10 −5 5 10 x −10 −5 5 10 y −8.00 −6.75 −5.50 −4.25 −3.00

G×R ×M

slide-8
SLIDE 8

Lamp-post geometry versus broken power law

0.1 0.2 0.3 1 2 3 4 5 6 7 8 Photon flux E [keV] h = 1.5 GM/c2, θo=30°, q = 6.2, rb = 5 GM/c2 LPI BPI LPN 0.1 0.2 0.3 0.4 1 2 3 4 5 6 7 8 Photon flux E [keV] h = 2 GM/c2, θo=30°, q = 4.3, rb = 10 GM/c2 LPI BPI LPN

For low heights: → broken power-law is not a good approximation of lamp-post geometry → line shape is greatly influenced by the emission directionality → this is mainly due to its dependence on the incident angle

slide-9
SLIDE 9

Disc ionization

2 4 6 8 10 0.01 0.1 1 10 100 Photons cm−2 s−1 keV−1 Energy (keV) Current Theoretical Model

dovciak 23−Feb−2013 18:17

reflionx ξ = 1000 ξ = 10 ◮ Ross RR & Fabian AC (2005), MNRAS, 358, 211 ◮ Svoboda J, Dovˇ ciak M, Goosmann RW, Jethwa P , Karas V, Miniutti G & Guainazzi M (2012) A&A, 545, A106

slide-10
SLIDE 10

Disc ionization

Dependence on height:

  • 4
  • 3
  • 2
  • 1

1 2 3 4 1 10 100 log ξ r [GM/c2] h [GM/c2] 1.5 3 10 100

ξ ∼ L/Ledd M nH L = 0.001Ledd M = 108M⊙ nH = 1015cm−3 a = 1, Γ = 2.0

slide-11
SLIDE 11

Disc ionization

Dependence on photon index:

  • 4
  • 3
  • 2
  • 1

1 2 3 4 1 10 100 log ξ r [GM/c2] Γ 1.4 2.0 2.6

ξ ∼ L/Ledd M nH L = 0.001Ledd M = 108M⊙ nH = 1015cm−3 a = 1, h = 3

slide-12
SLIDE 12

Disc ionization

Dependence on photon index:

  • 4
  • 3
  • 2
  • 1

1 2 3 4 1 10 100 log ξ r [GM/c2] Γ 1.4 2.0 2.6

ξ ∼ L/Ledd M nH L = 0.001Ledd M = 108M⊙ nH = 1015cm−3 a = 1, h = 10

slide-13
SLIDE 13

Disc ionization

Dependence on density profile:

  • 4
  • 3
  • 2
  • 1

1 2 3 4 1 10 100 log ξ r [GM/c2] Γ = 2.0, qn = 0 Γ = 2.0, qn = -2 Γ = 2.6, qn = -2

ξ ∼ L/Ledd M nH L = 0.001Ledd M = 108M⊙ nH = 1015cm−3 a = 1, h = 10 nH ∼ r qn

slide-14
SLIDE 14

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

slide-15
SLIDE 15

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

◮ scales the primary

flux (given in Ledd)

◮ scales the incident

flux (as D−1)

◮ scales the ionization ◮ scales the reflected

flux

slide-16
SLIDE 16

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

◮ affects the primary

flux (light bending model)

◮ affects the incident

flux (radial structure)

◮ affects the ionization ◮ affects the reflected

flux

slide-17
SLIDE 17

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

slide-18
SLIDE 18

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

◮ scales the primary

flux

◮ scales the incident

flux

◮ scales the ionization ◮ scales the reflected

flux

slide-19
SLIDE 19

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2. may be used to estimate discrepancy between the primary and reflected flux (e.g. due to the anisotropy

  • r obscuration of the

primary radiation)

slide-20
SLIDE 20

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

◮ affect the ionization

slide-21
SLIDE 21

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

slide-22
SLIDE 22

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

slide-23
SLIDE 23

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

slide-24
SLIDE 24

KYREFLIONX parameters

a/M GM/c 0.9982 0. 1. theta o deg 30. 0. 89. rin GM/c2 1. 1. 1000. ms 1. 0. 1. rout GM/c2 400. 1. 1000. M/M8 1. 1e-8 1e+3 height GM/c2 3. 1.1 100. PhoIndex 2. 1.4 3.3 L/Ledd 0.001 1e-10 1e+10 Np:Nr 0. 0. 10. density 1. 1e-8 1e+8 den prof 0.

  • 5.

0. abun 1. 0.1 20. zshift 0.

  • 0.999

10. limb 0. 0. 2. tab 2. 1. 2. sw 2. 1. 2.

slide-25
SLIDE 25

KYREFLIONX example

0.1 1 10 100 0.01 0.1 1 10 keV2 (Photons cm−2 s−1 keV−1) Energy (keV) Current Theoretical Model

dovciak 19−Jul−2013 08:47

slide-26
SLIDE 26

Dynamic spectrum – ionized reflection

E ×FE