-
Analytic Light Curves in Reflected Light: Phase Curves, Occultations, and Non-Lambertian Scattering for Spherical Planets and Moons
Authors:
Rodrigo Luger,
Eric Agol,
Fran Bartolić,
Daniel Foreman-Mackey
Abstract:
We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations,…
▽ More
We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations, as well as to photometry of solar system bodies. We derive our solutions for Lambertian bodies illuminated by a point source, but extend them to model illumination sources of finite angular size and rough surfaces with phase-dependent scattering. Our algorithm is implemented in Python within the open-source starry mapping framework and is designed with efficient gradient-based inference in mind. The algorithm is 4-5 orders of magnitude faster than direct numerical evaluation methods and about 10 orders of magnitude more precise. We show how the techniques developed here may one day lead to the construction of two-dimensional maps of terrestrial planet surfaces, potentially enabling the detection of continents and oceans on exoplanets in the habitable zone.
△ Less
Submitted 10 March, 2021;
originally announced March 2021.
-
Occultation mapping of Io's surface in the near-infrared I: Inferring static maps
Authors:
Fran Bartolić,
Rodrigo Luger,
Daniel Foreman-Mackey,
Robert R. Howell,
Julie A. Rathbun
Abstract:
Jupiter's moon Io is the most volcanically active body in the Solar System with hundreds of active volcanoes varying in intensity on different timescales. Io has been observed during occultations by other Galilean moons and Jupiter since the 1980s, using high-cadence near infrared photometry. These observations encode a wealth of information about the volcanic features on its surface. We built a g…
▽ More
Jupiter's moon Io is the most volcanically active body in the Solar System with hundreds of active volcanoes varying in intensity on different timescales. Io has been observed during occultations by other Galilean moons and Jupiter since the 1980s, using high-cadence near infrared photometry. These observations encode a wealth of information about the volcanic features on its surface. We built a generative model for the observed occultations using the code starry which enables fast, analytic, and differentiable computation of occultation light curves in emitted and reflected light. Our probabilistic Bayesian model is able to recover known hotspots on the surface of Io using only two light curves and without any assumptions on the locations, shapes or the number of spots. The methods we have developed are also directly applicable to the problem of mapping the surfaces of stars and exoplanets.
△ Less
Submitted 5 March, 2021;
originally announced March 2021.
-
A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and IV
Authors:
Nathan Golovich,
William A. Dawson,
Fran Bartolić,
Casey Y. Lam,
Jessica R. Lu,
Michael S. Medford,
Michael D. Schneider,
George Chapline,
Edward F. Schlafly,
Alex Drlica-Wagner,
Kerianne Pruett
Abstract:
Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light…
▽ More
Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light curves as well as over-simplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian Process model. We used light curves for nearly 10,000 OGLE-III and IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the timescale distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and previous analyses over-estimated the number of long timescale ($t_E>100$ days) events due to their over simplistic models ignoring parallax effects and stellar variability. We use our fits to identify hundreds of events that are likely black holes.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
VBBinaryLensing: a public package for microlensing light curve computation
Authors:
V. Bozza,
E. Bachelet,
F. Bartolić,
T. Heintz,
A. Hoag,
M. Hundertmark
Abstract:
The computation of microlensing light curves represents a bottleneck for the modeling of planetary events, making broad searches in the vast parameter space of microlensing extremely time-consuming. The release of the first version of VBBinaryLensing (based on the advanced contour integration method presented in Bozza(2010)) has represented a considerable advance in the field, with the birth of se…
▽ More
The computation of microlensing light curves represents a bottleneck for the modeling of planetary events, making broad searches in the vast parameter space of microlensing extremely time-consuming. The release of the first version of VBBinaryLensing (based on the advanced contour integration method presented in Bozza(2010)) has represented a considerable advance in the field, with the birth of several analysis platforms running on this code. Here we present the version 2.0 of VBBinaryLensing, which contains several upgrades with respect to the first version, including new decision trees that introduce important optimizations in the calculations.
△ Less
Submitted 15 May, 2018;
originally announced May 2018.