Foraging with MUSHROOMS: A Mixed-Integer Linear Programming Scheduler for Multimessenger Target of Opportunity Searches with the Zwicky Transient Facility

Abstract

Electromagnetic follow-up of gravitational wave detections is very resource intensive, taking up hours of limited observation time on dozens of telescopes. Creating more efficient schedules for follow-up will lead to a commensurate increase in counterpart location efficiency without using more telescope time. Widely used in operations research and telescope scheduling, mixed integer linear programming (MILP) is a strong candidate to produce these higher-efficiency schedules, as it can make use of powerful commercial solvers that find globally optimal solutions to provided problems . We detail a new target of opportunity scheduling algorithm designed with Zwicky Transient Facility in mind that uses mixed integer linear programming. We compare its performance to \texttt{gwemopt}, the tuned heuristic scheduler used by the Zwicky Transient Facility and other facilities during the third LIGO-Virgo gravitational wave observing run. This new algorithm uses variable-length observing blocks to enforce cadence requirements and ensure field observability, along with having a secondary optimization step to minimize slew time. \blue{We show that by employing a hybrid method utilizing both this scheduler and \texttt{gwemopt}, the previous scheduler used, in concert, we can achieve an average improvement in detection efficiency of 3\%-11\% over \texttt{gwemopt} alone} for a simulated binary neutron star merger data set consistent with LIGO-Virgo's third observing run, highlighting the potential of mixed integer target of opportunity schedulers for future multimessenger follow-up surveys.

arXiv comment10 pages, 5 figures

References in corpus (38)

• Astropy: A Community Python Package for Astronomy
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• GW170817: Measurements of Neutron Star Radii and Equation of State
• GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
• The Zwicky Transient Facility: System Overview, Performance, and First Results
• An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi-GBM Detection of GRB 170817A
• The Zwicky Transient Facility: Data Processing, Products, and Archive
• Swope Supernova Survey 2017a (SSS17a), the Optical Counterpart to a Gravitational Wave Source
• The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. II. UV, Optical, and Near-IR Light Curves and Comparison to Kilonova Models
• Spectroscopic identification of r-process nucleosynthesis in a double neutron star merger
• The X-ray counterpart to the gravitational wave event GW 170817
• A kilonova as the electromagnetic counterpart to a gravitational-wave source
• The Zwicky Transient Facility: Science Objectives
• Illuminating Gravitational Waves: A Concordant Picture of Photons from a Neutron Star Merger
• Swift and NuSTAR observations of GW170817: detection of a blue kilonova
• The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. IV. Detection of Near-infrared Signatures of r-process Nucleosynthesis with Gemini-South
• GW170817: Joint Constraint on the Neutron Star Equation of State from Multimessenger Observations
• A Radio Counterpart to a Neutron Star Merger
• Rapid Bayesian position reconstruction for gravitational-wave transients
• The Zwicky Transient Facility: Surveys and Scheduler
• Early Spectra of the Gravitational Wave Source GW170817: Evolution of a Neutron Star Merger
• The Zwicky Transient Facility: Observing System
• Electromagnetic Evidence that SSS17a is the Result of a Binary Neutron Star Merger
• Multi-messenger Bayesian parameter inference of a binary neutron-star merger
• Constraints on the neutron star equation of state from AT2017gfo using radiative transfer simulations
• POSSIS: predicting spectra, light curves and polarization for multi-dimensional models of supernovae and kilonovae
• Discovery and confirmation of the shortest gamma ray burst from a collapsar
• GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR
• Data-driven expectations for electromagnetic counterpart searches based on LIGO/Virgo public alerts
• Optical follow-up of the neutron star-black hole mergers S200105ae and S200115j
• Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers
• Implications of the search for optical counterparts during the first six months of the Advanced LIGO's and Advanced Virgo's third observing run: possible limits on the ejecta mass and binary properties
• Optimizing searches for electromagnetic counterparts of gravitational wave triggers
• astroplan: An Open Source Observation Planning Package in Python
• Implications of the search for optical counterparts during the second part of the Advanced LIGO's and Advanced Virgo's third observing run: lessons learned for future follow-up observations
• Teamwork Makes the Dream Work: Optimizing Multi-Telescope Observations of Gravitational-Wave Counterparts
• 2900 square degree search for the optical counterpart of short gamma-ray burst GRB 180523B with the Zwicky Transient Facility
• Dynamic Scheduling: Target of Opportunity Observations of Gravitational Wave Events

Cited by in corpus (4)

• A data science platform to enable time-domain astronomy
• Optimal Follow-Up of Gravitational-Wave Events with the UltraViolet EXplorer (UVEX)
• GRRIS: a real-time intra-site observation scheduling scheme for distributed survey telescope arrays
• Automated Scheduling of Doppler Exoplanet Observations at Keck Observatory