This thesis describes and develops procedures for the generation of theoretical lightcurves that can be used to model gravitational microlensing events that involve multiple lenses. Of particular interest are the cases involving a single lens star with one or more orbiting planets, as this has proven to be an effective way of detecting extrasolar planets. Although there is an analytical expression for microlensing lightcurves produced by single lensing body, the generation of model lightcurves for more than one lensing body requires the use of numerical techniques. The method developed here, known as the semi-analytic method, involves the analytical rearrangement of the relatively simple ‘lens equation’ to produce a high-order complex lens polynomial. Root-finding algorithms are then used to obtain the roots of this ‘lens polynomial’ in order to locate the positions of the images and calculate their magnifications.By running example microlensing events through the root-finding algorithms, both the speed and accuracy of the Laguerre and Jenkins-Traub algorithms were investigated. It was discovered that, in order to correctly identify the image positions, a method involving solutions of several ‘lens polynomials’ corresponding to different coordinate origins needed to be invoked. Multipole and polygon approximations were also developed to include finite source and limb darkening effects. The semi-analytical method and the appropriate numerical techniques were incorporated into a C++ modelling code at VUW (Victoria University of Wellington) known as mlens2. The effectiveness of the semi-analytic method was demonstrated using mlens2 to generate theoretical lightcurves for the microlensing events MOA-2009-BLG-319 and OGLE-2006-BLG-109. By comparing these theoretical lightcurves with the observed photometric data and the published models, it was demonstrated that the semi-analytic method described in this thesis is a robust and efficient method for discovering extrasolar planets.
展开▼
