170 likes | 263 Views
New calculation method of multiple gravitational lensing system. F. Abe Nagoya University. 18 th International Conference on Gravitational Microlensing , Santa Barbara, 21 st Jan 2014. Contents. Introduction Lensing equation Magnification map on the lens plane Transfer matrix
E N D
New calculation method of multiple gravitational lensing system F. Abe Nagoya University 18th International Conference on Gravitational Microlensing, Santa Barbara, 21st Jan 2014
Contents • Introduction • Lensing equation • Magnification map on the lens plane • Transfer matrix • Approximation • Iteration • Problems • Summary
Quasar microlensing (Garsden, Bate, Lewis, 2011, NRAS 418, 1012) Multiple lenses cause complex magnification pattern!!
Calculation methods • Single lens • Simple quadratic equation (Lieb 1964) • Binary lens • Quintic equation (Witt & Mao 1995, Asada 2002) • Inverse ray shooting (Schneider & Weiss 1987) • Triple lens and more • 10th order polynomial equation (Rhie 2002) • Inverse ray shooting (Schneider & Weiss 1987) • Perturbation (Han 2005, Asada 2008)
Lensing configuration , j = 1, m m: number of images Lensplane Sourceplane and are normalized by Lensing equation θy βy Image Source Lensqi ? Single source makes multiple images Lensing equation is difficult to solve Observer θx βx DL DS
Lensing equation Lensing equation Scalar potential Lensing Straight projection
Jacobian matrix Jacobian matrix
Jacobian determinant and magnification Jacobian determinant Magnification
Magnification map on the lens plane θy , j = 1, m m: number of images = θx
Linear expression : infinitesimally small , Inverse matrix
Calculation of image position : reference point on the source plane exactly traced from a point on the lensing plane : a point on the source plane close to : first approximation of the image position corresponding to : second approximation of the image position corresponding to Iteration
Calculation of image position Lensplane Sourceplane Lensing equation θy βy Image Source 1 0 Lensqi 1 2 t Observer θx βx DL DS
Problems in and • This method only finds an image close to . • To find other images, we must try other . • One of the methods is to calculate and for mesh points to use as and . • The other method is to solve polynomial equation for similar configuration. • If steps over caustic, calculation become divergent. So we need to select other .
Summary • Analytic form of Jacobian matrix is derived for general multiple lens system • Using Jacobian determinant, magnification on the lens plane can be calculated • Approximate image position can be calculated from a close reference source point which is exactly traced from lens plane • Calculation to get image position converges in 3-5 times iteration • Although there are problems to get reference point, this method may be useful for future multiple lens analyses