2010ApJ...716.1408B

I present a previously unpublished method for calculating and modeling multiple lens microlensing events that is based on the image centered ray-shooting approach of Bennett & Rhie. It has been used to model a wide variety of binary and triple lens systems, but it is designed to efficiently model high-magnification planetary microlensing events, because these high-magnification events are, by far, the most challenging events to model. It is designed to be efficient enough to handle complicated microlensing events, which include more than two lens masses and lens orbital motion. This method uses a polar coordinate integration grid with a smaller grid spacing in the radial direction than in the angular direction, and it employs an integration scheme specifically designed to handle limb-darkened sources. I present tests that show that these features achieve second-order accuracy for the light curves of a number of high-magnification planetary events. They improve the precision of the calculations by a factor of >100 compared to first-order integration schemes with the same grid spacing in both directions (for a fixed number of grid points). This method also includes a χ² minimization method, based on the Metropolis algorithm, that allows the jump function to vary in a way that allows quick convergence to χ² minima. Finally, I introduce a global parameter space search strategy that allows a blind search of parameter space for light curve models without requiring χ² minimization over a large grid of fixed parameters. Instead, the parameter space is explored on a grid of initial conditions for a set of χ² minimizations using the full parameter space. While this method may be somewhat faster than methods that find the χ² minima over a large grid of parameters, I argue that the main strength of this method is for events with the signals of multiple planets, where a much higher dimensional parameter space must be explored to find the correct light curve model.