The Large Synoptic Survey Telescope is designed to provide an unprecedented
optical imaging dataset that will support investigations of our Solar System,
Galaxy and Universe, across half the sky and over ten years of repeated
observation. However, exactly how the LSST observations will be taken (the
observing strategy or "cadence") is not yet finalized. In this
dynamically-evolving community white paper, we explore how the detailed
performance of the anticipated science investigations is expected to depend on
small changes to the LSST observing strategy. Using realistic simulations of
the LSST schedule and observation properties, we design and compute diagnostic
metrics and Figures of Merit that provide quantitative evaluations of different
observing strategies, analyzing their impact on a wide range of proposed
science projects. This is work in progress: we are using this white paper to
communicate to each other the relative merits of the observing strategy choices
that could be made, in an effort to maximize the scientific value of the
survey. The investigation of some science cases leads to suggestions for new
strategies that could be simulated and potentially adopted. Notably, we find
motivation for exploring departures from a spatially uniform annual tiling of
the sky: focusing instead on different parts of the survey area in different
years in a "rolling cadence" is likely to have significant benefits for a
number of time domain and moving object astronomy projects. The communal
assembly of a suite of quantified and homogeneously coded metrics is the vital
first step towards an automated, systematic, science-based assessment of any
given cadence simulation, that will enable the scheduling of the LSST to be as
well-informed as possible.
Models of gravitational waveforms play a critical role in detecting and
characterizing the gravitational waves (GWs) from compact binary coalescences.
Waveforms from numerical relativity (NR), while highly accurate, are too
computationally expensive to produce to be directly used with Bayesian
parameter estimation tools like Markov-chain-Monte-Carlo and nested sampling.
We propose a Gaussian process regression (GPR) method to generate accurate
reduced-order-model waveforms based only on existing accurate (e.g. NR)
simulations. Using a training set of simulated waveforms, our GPR approach
produces interpolated waveforms along with uncertainties across the parameter
space. As a proof of concept, we use a training set of IMRPhenomD waveforms to
build a GPR model in the 2-d parameter space of mass ratio $q$ and
equal-and-aligned spin $\chi_1=\chi_2$. Using a regular, equally-spaced grid of
120 IMRPhenomD training waveforms in $q\in[1,3]$ and $\chi_1 \in [-0.5,0.5]$,
the GPR mean approximates IMRPhenomD in this space to mismatches below
$4.3\times 10^{-5}$. Our approach can alternatively use training waveforms
directly from numerical relativity. Beyond interpolation of waveforms, we also
present a greedy algorithm that utilizes the errors provided by our GPR model
to optimize the placement of future simulations. In a fiducial test case we
find that using the greedy algorithm to iteratively add simulations achieves
GPR errors that are $\sim 1$ order of magnitude lower than the errors from
using Latin-hypercube or square training grids.