-
Analytic Light Curve for Mutual Transits of Two Bodies Across a Limb-darkened Star
Authors:
Tyler A. Gordon,
Eric Agol
Abstract:
We present a solution for the light curve of two bodies mutually transiting a star with polynomial limb darkening. The term "mutual transit" in this work refers to a transit of the star during which overlap occurs between the two transiting bodies. These could be an exoplanet with an exomoon companion, two exoplanets, an eclipsing binary and a planet, or two stars eclipsing a third in a triple sta…
▽ More
We present a solution for the light curve of two bodies mutually transiting a star with polynomial limb darkening. The term "mutual transit" in this work refers to a transit of the star during which overlap occurs between the two transiting bodies. These could be an exoplanet with an exomoon companion, two exoplanets, an eclipsing binary and a planet, or two stars eclipsing a third in a triple star system. We include analytic derivatives of the light curve with respect to the positions and radii of both bodies. We provide code that implements a photodynamical model for a mutual transit. We include two dynamical models, one for hierarchical systems in which a secondary body orbits a larger primary (e.g. an exomoon system) and a second for confocal systems in which two bodies independently orbit a central mass (e.g. two planets in widely separated orbits). Our code is fast enough to enable inference with MCMC algorithms, and the inclusion of derivatives allows for the use of gradient-based inference methods such as Hamiltonian Monte Carlo. While applicable to a variety of systems, this work was undertaken primarily with exomoons in mind. It is our hope that making this code publicly available will reduce barriers for the community to assess the detectability of exomoons, conduct searches for exomoons, and attempt to validate existing exomoon candidates. We also anticipate that our code will be useful for studies of planet-planet transits in exoplanetary systems, transits of circumbinary planets, and eclipses in triple-star systems.
△ Less
Submitted 13 July, 2022;
originally announced July 2022.
-
exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
Authors:
Daniel Foreman-Mackey,
Rodrigo Luger,
Eric Agol,
Thomas Barclay,
Luke G. Bouma,
Timothy D. Brandt,
Ian Czekala,
Trevor J. David,
Jiayin Dong,
Emily A. Gilbert,
Tyler A. Gordon,
Christina Hedges,
Daniel R. Hey,
Brett M. Morris,
Adrian M. Price-Whelan,
Arjun B. Savel
Abstract:
"exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functi…
▽ More
"exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functions and probability distributions required when fitting exoplanet datasets or other astronomical time series. While it has been used for other applications, such as the study of stellar variability, the primary purpose of "exoplanet" is the characterization of exoplanets or multiple star systems using time-series photometry, astrometry, and/or radial velocity. In particular, the typical use case would be to use one or more of these datasets to place constraints on the physical and orbital parameters of the system, such as planet mass or orbital period, while simultaneously taking into account the effects of stellar variability.
△ Less
Submitted 24 June, 2021; v1 submitted 5 May, 2021;
originally announced May 2021.
-
Stellar Rotation in the K2 Sample: Evidence for Modified Spindown
Authors:
Tyler A. Gordon,
James R. A. Davenport,
Ruth Angus,
Daniel Foreman-Mackey,
Eric Agol,
Kevin R. Covey,
Marcel Agüeros,
David Kipping
Abstract:
We analyze light curves of 284,834 unique K2 targets using a Gaussian process model with a quasi-periodic kernel function. By crossmatching K2 stars to observations from Gaia Data Release 2, we have identified 69,627 likely main-sequence stars. From these we select a subsample of 8,977 stars on the main-sequence with highly precise rotation period measurements. With this sample we recover the gap…
▽ More
We analyze light curves of 284,834 unique K2 targets using a Gaussian process model with a quasi-periodic kernel function. By crossmatching K2 stars to observations from Gaia Data Release 2, we have identified 69,627 likely main-sequence stars. From these we select a subsample of 8,977 stars on the main-sequence with highly precise rotation period measurements. With this sample we recover the gap in the rotation period-color diagram first reported by \cite{McQuillan2013}. While the gap was tentatively detected in \cite{Reinhold2020}, this work represents the first robust detection of the gap in K2 data for field stars. This is significant because K2 observed along many lines of sight at wide angular separation, in contrast to Kepler's single line of sight. Together with recent results for rotation in open clusters, we interpret this gap as evidence for a departure from the $t^{-1/2}$ Skumanich spin down law, rather than an indication of a bimodal star formation history. We provide maximum likelihood estimates and uncertainties for all parameters of the quasi-periodic light curve model for each of the 284,834 stars in our sample.
△ Less
Submitted 1 April, 2021; v1 submitted 19 January, 2021;
originally announced January 2021.
