Accretion on Black Holes: Testing the Lamp Post Geometry Thomas Dauser Remeis Observatory Bamberg & ECAP in collaboration with J. García, J. Wilms, M. Fink, T. Beuchert, and many others
10 1 kev 2 (Photons cm -2 s -1 kev -1 ) 10 0 10-1 10-2 Total Model Coronal Emission Relativistic Reflection Distant Reflection 10-3 1 10 100 Energy (kev) X-ray Reflection from Accretion Disks
Accretion Geometry: The Primary Source of Radiation Lamp Post Corona steep emissivity emissivity r 3 Thomas Dauser Remeis Observatory Bamberg & ECAP 2
Accretion Geometry: The Primary Source of Radiation Lamp Post Corona steep emissivity emissivity r 3 Usually emissivities steeper than r 3 are observed (described by a broken power law emissivity) (see, e.g., Fabian et al., 2004; Miller et al., 2013) naturally explained in the lamp post geometry agrees with measurements of the emissivity (see, e.g., Wilkins & Fabian, 2012; Wilkins & Gallo, 2015) Thomas Dauser Remeis Observatory Bamberg & ECAP 2
How can we test the nature of the primary source? analyze the continuum: usually continuum produced by a corona or a jet fits equally well to data (see, e.g., Nowak et al., 2011) measure the emissivity profile (see, e.g., Wilkins & Fabian, 2012) timing / reverberation lags... (see, e.g., Kara et al., 2016) compare a full model of primary and reflected radiation to data (including normalization) Thomas Dauser Remeis Observatory Bamberg & ECAP 3
Reflection Fraction in the Lamp Post Geometry 1 0.9 Black Hole Fraction of the Total Flux 0.8 0.7 0.6 0.5 0.4 0.3 0.2 Accretion Disk Infinity 0.1 0 2 5 10 20 50 100 Height [r g ] lamp post geometry: fraction of photons / flux hitting the disk compared to infinity depends on the height source Thomas Dauser Remeis Observatory Bamberg & ECAP 4
Reflection Fraction in the Lamp Post Geometry 1 0.9 Black Hole 20 a = 0.998 Fraction of the Total Flux 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 Accretion Disk Infinity Reflection Fraction Rf 10 5 2 1 0 2 5 10 20 50 100 Height [r g ] 2 5 10 20 50 100 Height [r g ] reflection fraction (R f ) = Flux(Accretion Disk) Flux(Infinity) Thomas Dauser Remeis Observatory Bamberg & ECAP 4
Reflection Fraction in the Lamp Post Geometry Fraction of the Total Flux 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 Black Hole Accretion Disk Infinity 2 5 10 20 50 100 Height [r g ] Reflection Fraction Rf 20 10 5 2 1 a = 0.998 a = 0.990 a = 0.900 a = 0.500 a = 0.000 2 5 10 20 50 100 Height [r g ] reflection fraction (R f ) = Flux(Accretion Disk) Flux(Infinity) Thomas Dauser Remeis Observatory Bamberg & ECAP 4
Different ways to define the reflected normalization Reflection Fraction R f Reflection Strength R s Ratio between intensity illuminating the disk to intensity reaching the observer System Intrinsic Ratio between reflected and direct flux at the Compton hump (20 40 kev) Observed Strength Motivation: in order to investigate the accretion geometry, need to understand how much primary radiation is incident on the accretion disk Thomas Dauser Remeis Observatory Bamberg & ECAP 5
Relativistic Reflection: Predicted Reflection Fraction reflection spectra: source at height h seen under inclination θ 100 θ = 25 θ = 50 θ = 75 10 Flux 1 0.1 h = 2 r g h = 3 r g h = 6 r g h = 100 r g 0.1 1 10 Energy [kev] 100 0.1 1 10 100 0.1 1 10 100 Energy [kev] Energy [kev] strong dependency of reflection on inclination angle large boost in reflection fraction possible for low heights largest boost for large angle and low height source Thomas Dauser Remeis Observatory Bamberg & ECAP 6
Observed Spectrum: Reflection Strength θ = 25 θ = 50 θ = 75 R s (reflection strength) a = 0.99 a = 0.00 Reflection Fraction / Strength 10 1 0.1 R f (reflection fraction) 2 3 5 10 20 Height [r g ] 50 2 3 5 10 20 50 2 3 5 10 20 Height [r g ] Height [r g ] 50 strength (R s ) strongly depends on the inclination greatest strength for large inclination angles Thomas Dauser Remeis Observatory Bamberg & ECAP 7
Observed Spectrum: Reflection Strength θ = 25 θ = 50 θ = 75 R s (reflection strength) a = 0.99 a = 0.00 Reflection Fraction / Strength 10 1 0.1 R f (reflection fraction) 2 3 5 10 20 Height [r g ] 50 2 3 5 10 20 50 2 3 5 10 20 Height [r g ] Height [r g ] 50 strength (R s ) strongly depends on the inclination greatest strength for large inclination angles only high reflection fraction for low height and large spin prediction: some sources should show strong reflection Thomas Dauser Remeis Observatory Bamberg & ECAP 7
Sample Study: Testing the Lamp Post Geometry Fink at al., in prep (Poster J10) Aim: Constraining the height of the primary source for a sample of sources does the lamp post geometry describe the data? use 16 of the sources with best spectral coverage (XMM-Newton and NuSTAR) show that the physical model is able to describe the data constrain the height of the primary source does the predicted reflection fraction agree with measurements? compare with previous measurements Thomas Dauser Remeis Observatory Bamberg & ECAP 8
Sample Study: Height vs. Reflection Fraction Fink at al., in prep (Poster J10) Thomas Dauser Remeis Observatory Bamberg & ECAP 9
Sample Study: Height vs. Reflection Fraction Fink at al., in prep (Poster J10) 50 20 h (rg) 10 5 Free Lamppost Self-consistent LP self-consistent normalization results in equally good fits 2 f max refl (h, a = 0.998) 1 0.5 1 2 5 10 20 50 f refl more physically sound results and better constraints (reflection fraction strongly depends on inclination and height) Thomas Dauser Remeis Observatory Bamberg & ECAP 9
Sample Study: Spin Measurements (preliminary) 1 0.95 0.9 Ark 120 Fairall 9 PDS 456 NGC 1365 MCG-6-30-15 Mrk 841 3C 382 Fink at al., in prep (Poster J10) 3C 390.3 Spin a 0.85 0.8 0.75 0.7 preliminary! 0.65 0.6 Free Lamppost Self-consistent Lamppost 0.55 Walton (2013) 0.5 0 2 4 Source 6 8 spin measurements largely in good agreement Thomas Dauser Remeis Observatory Bamberg & ECAP 10
The Primary Source: Cutoff Energy determines energetics of corona compactness (Fabian et al., 2015,... ) García et al. (2015a): cutoff energy important for relativistic reflection energy shifted due to GR effects 600 kev Thomas Dauser Remeis Observatory Bamberg & ECAP 11
The Primary Source: Cutoff Energy determines energetics of corona compactness (Fabian et al., 2015,... ) García et al. (2015a): cutoff energy important for relativistic reflection energy shifted due to GR effects E cut = 300 kev 600 kev Thomas Dauser Remeis Observatory Bamberg & ECAP 11
The Primary Source: Cutoff Energy determines energetics of corona compactness (Fabian et al., 2015,... ) García et al. (2015a): cutoff energy important for relativistic reflection energy shifted due to GR effects E cut = 300 kev 600 kev 900 kev 400 kev 300 kev Thomas Dauser Remeis Observatory Bamberg & ECAP 11
The Primary Source: Cutoff Energy determines energetics of corona compactness (Fabian et al., 2015,... ) current limitations of reflection modeling: García et al. one (2015a): input cutoff spectral energy shape important for the complete for relativistic diskreflection energy shifted due to GR effects problem: convolution of multiple Ezones cut = 300 of the kevaccretion disk necessary, but computationally very expensive 600 kev 300 kev Thomas Dauser Remeis Observatory Bamberg & ECAP 11
relxill 2.0 : overcoming this limitation fully re-written relxill model with a FFT algorithm for the convolution and additional optimizations faster despite multiple zones reflection spectrum for each zone relativistically smeared separately E1 cut E2 cut E3 cut E4 cut Thomas Dauser Remeis Observatory Bamberg & ECAP 12
relxill 2.0 : overcoming this limitation fully re-written relxill model with a FFT algorithm for the convolution and additional optimizations faster despite multiple zones reflection spectrum for each zone relativistically smeared separately E1 cut E2 cut E3 cut E4 cut download: www.sternwarte.uni-erlangen.de/research/relxill/ (other models unchanged, but faster: including nthcomp, coronal geometry, line models...) Thomas Dauser Remeis Observatory Bamberg & ECAP 12
relxill 2.0 : overcoming this limitation fully re-written relxill model with a FFT algorithm for the convolution and additional optimizations faster despite multiple zones reflection spectrum for each zone relativistically smeared separately E1 cut E2 cut E3 cut E4 cut download: www.sternwarte.uni-erlangen.de/research/relxill/ (other models unchanged, but faster: including nthcomp, coronal geometry, line models...) Thomas Dauser Remeis Observatory Bamberg & ECAP 12
Multi-Zone Model: Change in Cutoff Energy Ecut (fitted) [kev] 1000 500 200 100 50 10 ksec NuSTAR simulations of GX-339 (García et al., 2015b) simulate with multi-zone model and fit with the previously used 1-zone model model is combination of xillver and relxill Ratio 0.5 0 30 50 80 150 E cut [kev] 200 300 400 Thomas Dauser Remeis Observatory Bamberg & ECAP 13
Multi-Zone Model: Change in Cutoff Energy Ecut (fitted) [kev] 1000 500 200 100 50 10 ksec NuSTAR simulations of GX-339 (García et al., 2015b) simulate with multi-zone model and fit with the previously used 1-zone model model is combination of xillver and relxill previously cutoff energy 10% over-predicted Ratio 0.5 0 30 50 80 150 E cut [kev] 200 300 400 Thomas Dauser Remeis Observatory Bamberg & ECAP 13
Multi-Zone Model: Are other parameters affected? Spin a 1 0.8 0.6 0.4 0.2 10 ksec NuSTAR simulations of GX-339 (García et al., 2015b) simulate with multi-zone model and fit with the previously used 1-zone model model is combination of xillver and relxill 0.2 previously cutoff energy 10% over-predicted Ratio 0-0.2 30 50 80 150 E cut [kev] 200 300 400 other parameters (e.g., spin) largely unaffected Thomas Dauser Remeis Observatory Bamberg & ECAP 13
Summary reflection fraction contains information about the primary source reflection strength (measured normalization) depends on height, inclination,... sample study: self-consistent reflection fraction leads to physically consistent results lamp post geometry describes data well new relxill model: multiply radial zones for a proper modeling of the reflection (cutoff energy) only changes around 10% in the cutoff energy Thomas Dauser Remeis Observatory Bamberg & ECAP 14
Summary reflection fraction contains information about the primary source reflection strength (measured normalization) depends on height, inclination,... sample study: self-consistent reflection fraction leads to physically consistent results lamp post geometry describes data well new relxill model: multiply radial zones for a proper modeling of the reflection (cutoff energy) only changes around 10% in the cutoff energy Spin measurements: only reliable with a sound understanding/constraint of the system geometry? Thomas Dauser Remeis Observatory Bamberg & ECAP 14
References Fabian A.C., Lohfink A., Kara E., et al., 2015, MNRAS 451, 4375 Fabian A.C., Miniutti G., Gallo L., et al., 2004, MNRAS 353, 1071 García J.A., Dauser T., Steiner J.F., et al., 2015a, Astrophys. J., Lett. 808, L37 García J.A., Steiner J.F., McClintock J.E., et al., 2015b, ApJ 813, 84 Kara E., Alston W.N., Fabian A.C., et al., 2016, MNRAS 462, 511 Miller J.M., Parker M.L., Fuerst F., et al., 2013, Astrophys. J., Lett. 775, L45 Nowak M.A., Hanke M., Trowbridge S.N., et al., 2011, ApJ 728, 13 Wilkins D.R., Fabian A.C., 2012, MNRAS 424, 1284 Wilkins D.R., Gallo L.C., 2015, MNRAS 449, 129 Thomas Dauser Remeis Observatory Bamberg & ECAP 15