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 - - 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,
StrongGravity logo explanation
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
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
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)
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
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
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
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
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
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
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
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
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.
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
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
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.
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
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)
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
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.
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.
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.
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.
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