We have conducted a data study of leadership and participation in NASA's Astrophysics Explorer-class missions for the nine solicitations issued during the period 2008-2016, using gender as a marker of diversity. During this time, 102 Principal Investigators (PIs) submitted Explorer-class proposals; only four of these PIs were female. Among the 102 PIs, there were 61 unique PIs overall; of these, just three were female. The percentage of females in science teams in these proposals ranges from a low of 10% to a high of 19% across the various solicitations. Combining data from all these Explorer-class proposals, we find that the overall participation by females in science teams is 14%. Eighteen of the Explorer-class proposals had zero females in science roles, and this includes science teams with as many as 28 members. These results demonstrate that participation by women in the leadership of and, in many cases, on the science teams of proposals for Explorer-class missions is well below the representation of women in astronomy and astrophysics as a whole. In this white paper, we present our data and a discussion of our results, their context, and the ramifications for consideration by Astro2020 in its study of the state of the profession.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

The gravitational wave window onto the universe will open in roughly five years, when Advanced LIGO and Virgo achieve the first detections of high frequency gravitational waves, most likely coming from compact binary mergers. Electromagnetic follow-up of these triggers, using radio, optical, and high energy telescopes, promises exciting opportunities in multi-messenger time domain astronomy. In the next decade, space-based observations of low frequency gravitational waves from massive black hole mergers, and their electromagnetic counterparts, will open up further vistas for discovery. This two-part workshop at featured brief presentations and stimulating discussions on the challenges and opportunities presented by gravitational wave astronomy. Highlights from the workshop, with the emphasis on strategies for electromagnetic follow-up, are presented in this report.

We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms from equal-mass aligned-spin mergers, following an approach first presented in the complementary context of nonspinning black holes of varying mass ratio [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve together in lock-step through inspiral and merger, supporting the previous waveform description in terms of an adiabatically rigid rotator driving gravitational-wave emission-an implicit rotating source. We further apply the late-time merger-ringdown model for the rotational frequency introduced in [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)], along with an improved amplitude model appropriate for the dominant (2, \pm2) modes. This provides a quantitative description of the merger-ringdown waveforms, and suggests that the major features of these waveforms can be described with reference only to the intrinsic parameters associated with the state of the final black hole formed in the merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this model agrees to fitting factors better than 95% with the original numerical waveforms for system masses above \sim150M\odot. This model may be directly applicable to gravitational-wave detection of intermediate-mass black-hole mergers.

The gravitational wave window onto the universe is expected to open in ~ 5 years, when ground-based detectors make the first detections in the high-frequency regime. Gravitational waves are ripples in spacetime produced by the motions of massive objects such as black holes and neutron stars. Since the universe is nearly transparent to gravitational waves, these signals carry direct information about their sources - such as masses, spins, luminosity distances, and orbital parameters - through dense, obscured regions across cosmic time. This article explores gravitational waves as cosmic messengers, highlighting key sources, detection methods, and the astrophysical payoffs across the gravitational wave spectrum.

Recent breakthroughs in the field of numerical relativity have led to dramatic progress in understanding the predictions of General Relativity for the dynamical interactions of two black holes in the regime of very strong gravitational fields. Such black-hole binaries are important astrophysical systems and are a key target of current and developing gravitational-wave detectors. The waveform signature of strong gravitational radiation emitted as the black holes fall together and merge provides a clear observable record of the process. After decades of slow progress, these mergers and the gravitational-wave signals they generate can now be routinely calculated using the methods of numerical relativity. We review recent advances in understanding the predicted physics of events and the consequent radiation, and discuss some of the impacts this new knowledge is having in various areas of astrophysics.

Understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long-standing unsolved problem in theoretical physics. Black-hole mergers are monumental astrophysical events, releasing tremendous amounts of energy in the form of gravitational radiation, and are key sources for both ground- and space-based gravitational-wave detectors. The black-hole merger dynamics and the resulting gravitational waveforms can only be calculated through numerical simulations of Einstein's equations of general relativity. For many years, numerical relativists attempting to model these mergers encountered a host of problems, causing their codes to crash after just a fraction of a binary orbit could be simulated. Recently, however, a series of dramatic advances in numerical relativity has allowed stable, robust black-hole merger simulations. This remarkable progress in the rapidly maturing field of numerical relativity, and the new understanding of black-hole binary dynamics that is emerging is chronicled. Important applications of these fundamental physics results to astrophysics, to gravitational-wave astronomy, and in other areas are also discussed.

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum ``recoil kick'' of up to $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\sim\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters.

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%.

Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for ground-based interferometers as well as the space-based interferometer LISA. Recent progress in numerical relativity now makes it possible to calculate the waveforms from the strong-field dynamical merger and is revolutionizing our understanding of these systems. We review these dramatic developments, emphasizing applications to issues in gravitational wave observations. These new capabilities also make possible accurate calculations of the recoil or kick imparted to the final remnant black hole when the merging components have unequal masses, or unequal or unaligned spins. We highlight recent work in this area, focusing on results of interest to astrophysics.

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation, and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. We have performed multiple runs with different initial separations, orbital angular momenta, resolutions, extraction radii, and gauges. The full range of our kick speeds is 86--116 km s$^{-1}$, and the most reliable runs give kicks between 86 and 97 km s$^{-1}$. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts $z\gtrsim 10$, halos with masses $\lesssim 10^9 M_\odot$ will have difficulty retaining coalesced black holes after major mergers.

A remarkable series of breakthroughs in numerical relativity modeling of black hole binary mergers has occurred over the past few years. This paper provides a general overview of these exciting developments, focusing on recent progress in merger simulations and calculations of the resulting gravitational waveforms.

We study dynamics and radiation generation in the last few orbits and merger of a binary black hole system, applying recently developed techniques for simulations of moving black holes. Our analysis of the gravitational radiation waveforms and dynamical black hole trajectories produces a consistent picture for a set of simulations with black holes beginning on circular-orbit trajectories at a variety of initial separations. We find profound agreement at the level of one percent among the simulations for the last orbit, merger and ringdown. We are confident that this part of our waveform result accurately represents the predictions from Einstein's General Relativity for the final burst of gravitational radiation resulting from the merger of an astrophysical system of equal-mass non-spinning black holes. The simulations result in a final black hole with spin parameter a/m=0.69. We also find good agreement at a level of roughly 10 percent for the radiation generated in the preceding few orbits.

We present new techniqes for evolving binary black hole systems which allow the accurate determination of gravitational waveforms directly from the wave zone region of the numerical simulations. Rather than excising the black hole interiors, our approach follows the "puncture" treatment of black holes, but utilzing a new gauge condition which allows the black holes to move successfully through the computational domain. We apply these techniques to an inspiraling binary, modeling the radiation generated during the final plunge and ringdown. We demonstrate convergence of the waveforms and good conservation of mass-energy, with just over 3% of the system's mass converted to gravitional radiation.

LISA will be able to detect gravitational waves from inspiralling massive black hole (MBH) binaries out to redshifts z > 10. If the binary masses and luminosity distances can be extracted from the LISA data stream, this information can be used to reveal the merger history of MBH binaries and their host galaxies in the evolving universe. Since this parameter extraction generally requires that LISA observe the inspiral for a significant fraction of its yearly orbit, carrying out this program requires adequate sensitivity at low frequencies, f < 10^{-4} Hz. Using several candidate low frequency sensitivities, we examine LISA's potential for characterizing MBH binary coalescences at redshifts z > 1.

Calculations of the final merger stage of binary black hole evolution can only be carried out using full scale numerical relativity simulations. This article provides a general overview of these calculations, highlighting recent progress and current challenges.

Observations of low frequency gravitational waves by the space-based LISA mission will open a new observational window on the early universe and the emergence of structure. LISA will observe the dynamical coalescence of massive black hole binaries at high redshifts, giving an unprecedented look at the merger history of galaxies and the reionization epoch. LISA will also observe gravitational waves from the collapse of supermassive stars to form black holes, and will map the spacetime in the central regions of galaxy cusps at high precision.

This Resource Letter provides a guide to the literature on the physics and astrophysics of gravitational waves. Journals, books, reports, archives, and websites are provided as basic resources and for current research frontiers in detectors, data analysis, and astrophysical source modeling.

One interesting class of gravitational radiation sources includes rapidly rotating astrophysical objects that encounter dynamical instabilities. We have carried out a set of simulations of rotationally induced instabilities in differentially rotating polytropes. An $n$=1.5 polytrope with the Maclaurin rotation law will encounter the $m$=2 bar instability at $T/|W| \gtrsim 0.27$. Our results indicate that the remnant of this instability is a persistent bar-like structure that emits a long-lived gravitational radiation signal. Furthermore, dynamical instability is shown to occur in $n$=3.33 polytropes with the $j$-constant rotation law at $T/|W| \gtrsim 0.14$. In this case, the dominant mode of instability is $m$=1. Such instability may allow a centrifugally-hung core to begin collapsing to neutron star densities on a dynamical timescale. If it occurs in a supermassive star, it may produce gravitational radiation detectable by LISA.

The development of both ground- and space-based gravitational wave detectors provides new opportunities to observe the radiation from binaries containing neutron stars and black holes. Numerical simulations in 3-D are essential for calculating the coalescence waveforms, and comprise some of the most challenging problems in astrophysics today. This article briefly reviews the current status of efforts to calculate black hole and neutron star coalescences, and highlights challenges for the future.

Dynamical instability is shown to occur in differentially rotating polytropes with N = 3.33 and $T/|W| \gtrsim 0.14$. This instability has a strong m=1 mode, although the m=2, 3, and 4 modes also appear. Such instability may allow a centrifugally-hung core to begin collapsing to neutron star densities on a dynamical timescale. The gravitational radiation emitted by such unstable cores may be detectable with advanced ground-based detectors, such as LIGO II. If the instability occurs in a supermassive star, it may produce gravitational radiation detectable by the space-based detector LISA.

Adaptive techniques are crucial for successful numerical modeling of gravitational waves from astrophysical sources such as coalescing compact binaries, since the radiation typically has wavelengths much larger than the scale of the sources. We have carried out an important step toward this goal, the evolution of weak gravitational waves using adaptive mesh refinement in the Einstein equations. The 2-level adaptive simulation is compared with unigrid runs at coarse and fine resolution, and is shown to track closely the features of the fine grid run.

Compact astrophysical objects that rotate rapidly may encounter the dynamical ``bar instability.'' The bar-like deformation induced by this rotational instability causes the object to become a potentially strong source of gravitational radiation. We have carried out a set of long-duration simulations of the bar instability with two Eulerian hydrodynamics codes. Our results indicate that the remnant of this instability is a persistent bar-like structure that emits a long-lived gravitational radiation signal.

The 3-level leapfrog time integration algorithm is an attractive choice for numerical relativity simulations since it is time-symmetric and avoids non-physical damping. In Newtonian problems without velocity dependent forces, this method enjoys the advantage of long term stability. However, for more general differential equations, whether ordinary or partial, delayed onset numerical instabilities can arise and destroy the solution. A known cure for such instabilities appears to have been overlooked in many application areas. We give an improved cure ("deloused leapfrog") that both reduces memory demands (important for 3+1 dimensional wave equations) and allows for the use of adaptive timesteps without a loss in accuracy. We show both that the instability arises and that the cure we propose works in highly relativistic problems such as tightly bound geodesics, spatially homogeneous spacetimes, and strong gravitational waves. In the gravitational wave test case (polarized waves in a Gowdy spacetime) the deloused leapfrog method was five to eight times less CPU costly at various accuracies than the implicit Crank-Nicholson method, which is not subject to this instability.

We carry out 3-D numerical simulations of the dynamical instability in rapidly rotating stars initially modeled as polytropes with n = 1.5, 1.0, and 0.5. The calculations are done with a SPH code using Newtonian gravity, and the gravitational radiation is calculated in the quadrupole limit. All models develop the global m=2 bar mode, with mass and angular momentum being shed from the ends of the bar in two trailing spiral arms. The models then undergo successive episodes of core recontraction and spiral arm ejection, with the number of these episodes increasing as n decreases: this results in longer-lived gravitational wave signals for stiffer models. This instability may operate in a stellar core that has expended its nuclear fuel and is prevented from further collapse due to centrifugal forces. The actual values of the gravitational radiation amplitudes and frequencies depend sensitively on the radius of the star R_{eq} at which the instability develops.

We calculate the gravitational radiation produced by the coalescence of inspiraling binary neutron stars in the Newtonian regime using 3-dimensional numerical simulations. The stars are modeled as polytropes and start out in the point-mass regime at wide separation. The hydrodynamic integration is performed using smooth particle hydrodynamics (SPH) with Newtonian gravity, and the gravitational radiation is calculated using the quadrupole approximation. We have run a number of simulations varying the neutron star radii, equations of state, spins, and mass ratio. The resulting gravitational waveforms and spectra are rich in information about the hydrodynamics of coalescence, and show characteristic dependence on GM/Rc^2, the equation of state, and the mass ratio.

We have carried out 3-D numerical simulations of the dynamical bar instability in a rotating star and the resulting gravitational radiation using both an Eulerian code written in cylindrical coordinates and a smooth particle hydrodynamics (SPH) code. The star is modeled initially as a polytrope with index $n = 3/2$ and $T_{\rm rot}/|W| \approx 0.30$, where $T_{\rm rot}$ is the rotational kinetic energy and $|W|$ is the gravitational potential energy. In both codes the gravitational field is purely Newtonian, and the gravitational radiation is calculated in the quadrupole approximation. We have run 3 simulations with the Eulerian code, varying the number of angular zones and the treatment of the boundary between the star and the vacuum. Using the SPH code we did 7 runs, varying the number of particles, the artificial viscosity, and the type of initial model. We compare the growth rate and rotation speed of the bar, the mass and angular momentum distributions, and the gravitational radiation quantities. We highlight the successes and difficulties of both methods, and make suggestions for future improvements.

We calculate the gravitational radiation produced by the merger and coalescence of inspiraling binary neutron stars using 3-dimensional numerical simulations. The stars are modeled as polytropes and start out in the point-mass limit at wide separation. The hydrodynamic integration is performed using smooth particle hydrodynamics (SPH) with Newtonian gravity, and the gravitational radiation is calculated using the quadrupole approximation. We have run several simulations, varying both the neutron star radius and the equation of state. The resulting gravitational wave energy spectra $dE/df$ are rich in information about the hydrodynamics of merger and coalescence. In particular, our results demonstrate that detailed information on both $GM/Rc^2$ and the equation of state can in principle be extracted from the spectrum.

We present the first calculations of the gravitational radiation produced by nonaxisymmetric dynamical instability in a rapidly rotating compact star. The star deforms into a bar shape, shedding $\sim 4\%$ of its mass and $\sim 17\%$ of its angular momentum. The gravitational radiation is calculated in the quadrupole approximation. For a mass $M \sim 1.4$ M$_{\odot}$ and radius $R \sim 10$ km, the gravitational waves have frequency $\sim 4$ kHz and amplitude $h \sim 2 \times 10^{-22}$ at the distance of the Virgo Cluster. They carry off energy $\Delta E/M \sim 0.1\%$ and radiate angular momentum $\Delta J/J \sim 0.7\%$.

We have developed a 3-D Eulerian hydrodynamics code to model sources of gravitational radiation. The code is written in cylindrical coordinates $(\varpi,z,\varphi)$ and has moving grids in the $\varpi$ and $z$-directions. We use Newtonian gravity and calculate the gravitational radiation in the quadrupole approximation. This code has been tested on a variety of problems to verify its accuracy and stability, and the results of these tests are reported here.