Fast Multipole Method for Gravitational Lensing. Application to High Magnification Quasar Microlensing

Abstract

We introduce the use of the Fast Multipole Method (FMM) to speed up gravitational lensing ray tracing calculations. The method allows very fast calculation of ray deflections when a large number of deflectors, $N_{\ast }$, is involved, while keeping rigorous control on the errors. In particular, we apply this method, in combination with the Inverse Polygon Mapping technique (IPM), to quasar microlensing to generate microlensing magnification maps with very high workloads (high magnification, large size and/or high resolution) that require a very large number of deflectors. Using, FMM-IPM, the computation time can be reduced by a factor $\sim 10^5$ with respect to standard Inverse Ray Shooting, making the use of this algorithm on a personal computer comparable to the use of standard IRS on GPUs. We also provide a flexible web interface for easy calculation of microlensing magnification maps using FMM-IPM\footnote{http://gloton.ugr.es/microlensing/}. We exemplify the power of this new method by applying it to some challenging interesting astrophysical scenarios, including clustered primordial black holes, or extremely magnified stars close to the giant arcs of galaxy clusters. We also show the performance/use of FMM to calculate ray deflection for a halo resulting from cosmological simulations composed by a large number ($N\gtrsim 10^7$) of elements.

17 pages, 10 figures. Accepted for publication in The Astrophysical Journal. Accompanying web for microlensing map calculation at https://gloton.ugr.es/microlensing/

References in corpus (13)

• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• A Highly Magnified Star at Redshift 6.2
• Microlensing of Extremely Magnified Stars near Caustics of Galaxy Clusters
• Teraflop per second gravitational lensing ray-shooting using graphics processing units
• Microlensing of the broad emission lines in 27 gravitationally lensed quasars. Broad line region structure and kinematics
• An Efficient Method for Simulating Light Curves of Cosmological Microlensing and Caustic Crossing Events
• Gravitational Microlensing: A parallel, large-data implementation
• An Improved GPU-Based Ray-Shooting Code For Gravitational Microlensing
• Spectroscopy and polarimetry of the gravitationally lensed quasar SDSS J1004+4112 with the 6m SAO RAS telescope
• Fast simulations of extragalactic microlensing
• Abundance of LIGO/Virgo Black Holes from Microlensing Observations of Quasars with Reverberation Mapping Size Estimates
• Limiting the Abundance of LIGO/Virgo Black Holes with Microlensing Observations of Quasars of Finite Size
• A fast algorithm for simulating light curves of stars at extreme magnification affected by microlensing

Cited by in corpus (5)

• Revealing the structure of the lensed quasar Q 0957+561 III. Constraints on the size of the broad-line region
• Free-floating "planets'' in the macrolensed quasar Q2237+0305
• Amoeba: An AGN Model of Optical Emissions Beyond steady-state Accretion discs
• A GPU Code for Finding Microlensing Critical Curves and Caustics
• A Hamiltonian, post-Born, three-dimensional, on-the-fly ray tracing algorithm for gravitational lensing