Optimal scheduling of exoplanet direct imaging single-visit observations of a blind search survey

Abstract. We present an algorithm, effective over a broad range of planet populations and instruments, for optimizing integration times of an exoplanet direct imaging observation schedule, to maximize the number of unique exoplanet detections under realistic mission constraints. Our planning process uses “completeness” as a reward metric and the nonlinear combination of optimal integration time per target and constant overhead time per target as a cost metric constrained by a total mission time. We validate our planned target list and integration times for a specific telescope by running a Monte Carlo of full mission simulations using EXOSIMS, a code base for simulating telescope survey missions. These simulations encapsulate dynamic details such as time-varying local zodiacal light for each star, planet keep-out regions, exoplanet positions, and strict enforcement of observatory use over time. We test our methods on the Wide-Field Infrared Survey Telescope (WFIRST) coronagraphic instrument (CGI). We find that planet, Sun, and solar panel keep-out regions limit some target per-annum visibility to <28  %   and that the mean local zodiacal light flux for optimally scheduled observations is 22.79 mag arcsec  −  2. Both these values are more pessimistic than previous approximations and impact the simulated mission yield. We find that the WFIRST CGI detects 5.48  ±  0.17 and 16.26  ±  0.51 exoplanets, on average, when observing two different planet populations based on Kepler Q1-Q6 data and the full Kepler data release, respectively. Optimizing our planned observations using completeness derived from the more pessimistic planet population (in terms of overall planet occurrence rates) results in a more robust yield than optimization based on the more optimistic planet population. We also find optimization based on the more pessimistic population results in more small planet detections than optimization with the more optimistic population.

[1]  Dmitry Savransky,et al.  Parameterizing the Search Space of Starshade Fuel Costs for Optimal Observation Schedules , 2019 .

[2]  Margaret C Turnbull.  ExoCat-1: The Nearby Stellar Systems Catalog for Exoplanet Imaging Missions , 2015 .

[3]  Dmitry Savransky,et al.  Blind Search Single-Visit Exoplanet Direct Imaging Yield for Space Based Telescopes , 2019 .

[4]  John E. Krist.  End-to-end numerical modeling of AFTA coronagraphs , 2014, Astronomical Telescopes and Instrumentation.

[5]  Dmitry Savransky,et al.  PARAMETER DISTRIBUTIONS OF KEPLERIAN ORBITS , 2010, 1012.1526.

[6]  Edward J. Wollack,et al.  Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report , 2015, 1503.03757.

[7]  Stuart Bowyer,et al.  The 1997 reference of diffuse night sky brightness , 1998 .

[8]  Paul D. Spudis,et al.  Evaluating science return in space exploration initiative architectures , 1993 .

[9]  B. V. Semenov,et al.  A look towards the future in the handling of space science mission geometry , 2018 .

[10]  Robert W. Farquhar,et al.  The Utilization of Halo Orbits in Advanced Lunar Operations , 1971 .

[11]  Olivier Guyon,et al.  Science yield estimate with the Wide-Field Infrared Survey Telescope coronagraph , 2016 .

[12]  Dmitry Savransky,et al.  SEQUENTIAL COVARIANCE CALCULATION FOR EXOPLANET IMAGE PROCESSING , 2015, 1501.00991.

[13]  Jie Li,et al.  THE DISTRIBUTION OF TRANSIT DURATIONS FOR KEPLER PLANET CANDIDATES AND IMPLICATIONS FOR THEIR ORBITAL ECCENTRICITIES , 2011, 1102.0547.

[14]  Dmitry Savransky,et al.  Scheduling and target selection optimization for exoplanet imaging spacecraft , 2018, Astronomical Telescopes + Instrumentation.

[15]  E. Mamajek,et al.  INTRINSIC COLORS, TEMPERATURES, AND BOLOMETRIC CORRECTIONS OF PRE-MAIN-SEQUENCE STARS , 2013, 1307.2657.

[16]  Bijan Nemati.  Detector selection for the WFIRST-AFTA coronagraph integral field spectrograph , 2014, Astronomical Telescopes and Instrumentation.

[17]  Dmitry Savransky,et al.  WFIRST-AFTA coronagraph science yield modeling with EXOSIMS , 2015, 1511.02869.

[18]  Brian D. Mason,et al.  DATA MINING FOR DOUBLE STARS IN ASTROMETRIC CATALOGS , 2006 .

[19]  Dmitry Savransky,et al.  ANALYTICAL FORMULATION OF THE SINGLE-VISIT COMPLETENESS JOINT PROBABILITY DENSITY FUNCTION , 2016, 1607.01682.

[20]  F. Fressin,et al.  THE FALSE POSITIVE RATE OF KEPLER AND THE OCCURRENCE OF PLANETS , 2013, 1301.0842.

[21]  John E. Krist,et al.  Sensitivity of the WFIRST coronagraph performance to key instrument parameters , 2017, Optical Engineering + Applications.

[22]  N. J. Kasdin,et al.  Analyzing the Designs of Planet-Finding Missions , 2009, 0903.4915.

[23]  Howard Isaacson,et al.  The Occurrence and Mass Distribution of Close-in Super-Earths, Neptunes, and Jupiters , 2010, Science.

[24]  Dmitry Savransky,et al.  Space mission design for exoplanet imaging , 2013, Optics & Photonics - Optical Engineering + Applications.

[25]  Dmitry Savransky,et al.  Optimal starshade observation scheduling , 2018, Astronomical Telescopes + Instrumentation.

[26]  Andrew Cumming,et al.  The Keck Planet Search: Detectability and the Minimum Mass and Orbital Period Distribution of Extrasolar Planets , 2008, 0803.3357.

[27]  Dmitry Savransky,et al.  A Simple Depth-of-Search Metric for Exoplanet Imaging Surveys , 2017 .

[28]  Robert A. Brown.  Single-Visit Photometric and Obscurational Completeness , 2005, astro-ph/0503077.

[29]  Dmitry Savransky,et al.  Building Better Planet Populations for EXOSIMS , 2018 .

[30]  Christoph Mordasini,et al.  Hints for a Turnover at the Snow Line in the Giant Planet Occurrence Rate , 2018, The Astrophysical Journal.

[31]  Pini Gurfil,et al.  Quasi-periodic orbits of the restricted three-body problem made easy , 2007 .

[32]  K. Cahoy,et al.  EXOPLANET ALBEDO SPECTRA AND COLORS AS A FUNCTION OF PLANET PHASE, SEPARATION, AND METALLICITY , 2010, 1009.3071.

[33]  Dmitry Savransky,et al.  EXOSIMS: Exoplanet Open-Source Imaging Mission Simulator , 2017 .

[34]  Aki Roberge,et al.  MAXIMIZING THE ExoEarth CANDIDATE YIELD FROM A FUTURE DIRECT IMAGING MISSION , 2014, 1409.5128.

[35]  James P. Lloyd,et al.  Optimization of high-inclination orbits using planetary flybys for a zodiacal light-imaging mission , 2017, Optical Engineering + Applications.

[36]  Dmitry Savransky,et al.  Multi-mission modeling for space-based exoplanet imagers , 2017, Optical Engineering + Applications.

[37]  Shawn Domagal-Goldman,et al.  Exoplanet Classification and Yield Estimates for Direct Imaging Missions , 2018, 1802.09602.

[38]  Remi Soummer,et al.  NEW COMPLETENESS METHODS FOR ESTIMATING EXOPLANET DISCOVERIES BY DIRECT DETECTION , 2010 .

[39]  John E. Krist,et al.  Numerical modeling of the proposed WFIRST-AFTA coronagraphs and their predicted performances , 2015 .

[40]  Robin Lougee,et al.  The Common Optimization INterface for Operations Research: Promoting open-source software in the operations research community , 2003, IBM J. Res. Dev..

[41]  C. H. Acton,et al.  Ancillary data services of NASA's Navigation and Ancillary Information Facility , 1996 .

[42]  Dmitry Savransky,et al.  Planet Occurrence Rate Density Models Including Stellar Effective Temperature , 2018, Publications of the Astronomical Society of the Pacific.

[43]  Dmitry Savransky,et al.  WFIRST: Exoplanet Target Selection and Scheduling with Greedy Optimization , 2018 .

[44]  Stuart B. Shaklan,et al.  The lighter side of TPF-C: evaluating the scientific gain from a smaller mission concept , 2007, SPIE Optical Engineering + Applications.

[45]  Paul T. Boggs,et al.  Sequential Quadratic Programming , 1995, Acta Numerica.

[46]  Dmitry Savransky,et al.  The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics from 10 to 100 au , 2019, The Astronomical Journal.