We present EXOFAST, a fast, robust suite of routines written in IDL that is designed to fit exoplanetary transits and radial velocity variations simultaneously or separately and characterize the parameter uncertainties and covariances with a differential evolution Markov chain Monte Carlo method. We describe how our code incorporates both data sets to derive simultaneously stellar parameters along with the transit and RV parameters, resulting in more self-consistent results on an example fit of the discovery data of HAT-P-3b that is well-mixed in under 5 minutes on a standard desktop computer. We describe in detail how our code works and outline ways in which the code can be extended to include additional effects or generalized for the characterization of other data sets—including non-planetary data sets. We discuss the pros and cons of several common ways to parameterize eccentricity, highlight a subtle mistake in the implementation of MCMC that could bias the inferred eccentricity of intrinsically circular orbits to significantly non-zero results, discuss a problem with IDL's built-in random number generator in its application to large MCMC fits, and derive a method to analytically fit the linear and quadratic limb darkening coefficients of a planetary transit. Finally, we explain how we achieved improved accuracy and over a factor of 100 improvement in the execution time of the transit model calculation. Our entire source code, along with an easy-to-use online interface for several basic features of our transit and radial velocity fitting, are available online at http://astroutils.astronomy.ohio-state.edu/exofast.
[1]
William H. Press,et al.
Numerical Recipes 3rd Edition: The Art of Scientific Computing
,
2007
.
[2]
Elazar Uchupi,et al.
Annual Review of Earth and Planetary Sciences: Volume 19 (1991), 484p. US $60.00 and Volume 20 (1992), 631 p. US $64.00
,
1993
.
[3]
Daniel Durand,et al.
Astronomical Data Analysis Software and Systems XVIII
,
2009
.
[4]
Keivan G. Stassun,et al.
KELT-1b: A STRONGLY IRRADIATED, HIGHLY INFLATED, SHORT PERIOD, 27 JUPITER-MASS COMPANION TRANSITING A MID-F STAR
,
2012
.
[5]
L. Lucy.
Bayesian inference for orbital eccentricities
,
2012,
1206.3462.
[6]
Daniel Foreman-Mackey,et al.
emcee: The MCMC Hammer
,
2012,
1202.3665.
[7]
Keivan G. Stassun,et al.
KELT-3b: A HOT JUPITER TRANSITING A V = 9.8 LATE-F STAR
,
2012,
1211.1031.
[8]
L. Walkowicz,et al.
PHOTOMETRIC VARIABILITY IN KEPLER TARGET STARS: THE SUN AMONG STARS—A FIRST LOOK
,
2010,
1001.0414.
[9]
L. Hebb,et al.
KELT-2Ab: A HOT JUPITER TRANSITING THE BRIGHT (V = 8.77) PRIMARY STAR OF A BINARY SYSTEM
,
2012,
1206.1592.
[10]
Las Cumbres Observatory Global Telescope Network,et al.
PLANETARY CANDIDATES OBSERVED BY KEPLER. III. ANALYSIS OF THE FIRST 16 MONTHS OF DATA
,
2012,
1202.5852.
[11]
Eric Agol,et al.
Transit Analysis Package: An IDL Graphical User Interface for Exoplanet Transit Photometry
,
2011,
1102.1036.