Algorithms And Programs For Strong Gravitational Lensing In Kerr Space-time Including Polarization

Abstract

Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), a photon's polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing strong lensing effects of Kerr black holes, including effects on polarization. Our algorithms are realized in a program "KERTAP" in two versions: MATLAB and Python. The key ingredients of KERTAP are: a graphic user interface, a {\it backward} ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using 'mpi4py' package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity-related computations using interpreted languages such as MATLAB and Python.

11 pages, 7 figures

References in corpus (32)

• Magnetic Stress at the Marginally Stable Orbit: Altered Disk Structure, Radiation, and Black Hole Spin Evolution
• Quantitative Interpretation of Quasar Microlensing Light Curves
• The Sizes of the X-ray and Optical Emission Regions of RXJ1131-1231
• A Fast New Public Code for Computing Photon Orbits in a Kerr Spacetime
• X-ray Polarization from Accreting Black Holes: II. The Thermal State
• GYOTO: a new general relativistic ray-tracing code
• X-ray Microlensing in RXJ1131-1231 and HE1104-1805
• The Structure of the X-ray and Optical Emitting Regions of the Lensed Quasar Q 2237+0305
• grmonty: a Monte Carlo Code for Relativistic Radiative Transport
• X-Ray and Optical Microlensing in the Lensed Quasar PG 1115+080
• Light Curves for Rapidly-Rotating Neutron Stars
• Further Evidence that Quasar X-Ray Emitting Regions Are Compact: X-Ray and Optical Microlensing in the Lensed Quasar Q J0158-4325
• A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. I. Quasi-Kerr Black Holes
• A Model for High Frequency Quasi-periodic Oscillations from Accreting Black Holes
• Inferring the Inclination of a Black Hole Accretion Disk from Observations of its Polarized Continuum Radiation
• Line Emission from an Accretion Disk around a Rotating Black Hole: Toward a Measurement of Frame Dragging
• Iron Line Profiles in Strong Gravity
• Radiation Transfer of Emission Lines in Curved Space-Time
• A Monte Carlo Code for Relativistic Radiation Transport Around Kerr Black Holes
• Polarization signatures of strong gravity in AGN accretion discs
• X-ray Monitoring of Gravitational Lenses With Chandra
• YNOGK: A new public code for calculating null geodesics in the Kerr spacetime
• The Optical, Ultraviolet, and X-ray Structure of the Quasar HE 0435-1223
• The Structure of HE 1104-1805 from Infrared to X-Ray
• X-Ray and Optical Flux Anomalies in the Quadruply Lensed QSO 1RXS J1131-1231
• Scattering of circularly polarized light by a rotating black hole
• Discovery of Energy Dependent X-ray Microlensing in Q2237+0305
• On the influence of resonant absorption on the iron emission line profiles from accreting black holes
• Notes on Spinoptics in a Stationary Spacetime
• Radiative Transfer Along Rays in Curved Spacetimes
• Effects of Kerr Strong Gravity on Quasar X-ray Microlensing
• Inclination-Dependent AGN Flux Profiles From Strong Lensing of the Kerr Space-Time

Cited by in corpus (16)

• A public code for general relativistic, polarised radiative transfer around spinning black holes
• Measurement of the spin of the M87 black hole from its observed twisted light
• Thin Accretion Disk onto slowly rotating black holes in Einstein-Æther theory
• General purpose ray-tracing and polarized radiative transfer in General Relativity
• Thin accretion disk luminosity and its image around rotating black holes in perfect fluid dark matter
• Gravitational Faraday rotation of gravitational waves by a Kerr black hole
• Constraining the Generalized Uncertainty Principle with the light twisted by rotating black holes and M87*
• General relativistic polarized radiative transfer with inverse Compton scatterings
• Spin wave optics for gravitational waves lensed by a Kerr black hole
• Linear Spectropolarimetric Analysis of Fairall 9 with VLT/FORS2
• Constraining Quasar Relativistic Reflection Regions and Spins with Microlensing
• FANTASY: User-Friendly Symplectic Geodesic Integrator for Arbitrary Metrics with Automatic Differentiation
• A New Fast Monte Carlo Code for Solving Radiative Transfer Equations based on Neumann Solution
• Probing the gravitational Faraday rotation using quasar X-ray microlensing
• Epicyclic oscillations and accretion disk around a special Buchdahl-inspired spacetime
• Gravitational Faraday rotation of light propagation in the Kerr-Newman-Taub-NUT space-time