The dynamics of precessing binary black holes (BBHs) in the post-Newtonian regime has a strong timescale hierarchy: the orbital timescale is very short compared to the spin-precession timescale which, in turn, is much shorter than the radiation-reaction timescale on which the orbit is shrinking due to gravitational-wave emission. We exploit this timescale hierarchy to develop a multiscale analysis of BBH dynamics elaborating on the analysis of Kesden et al. [Phys. Rev. Lett. 114, 081103 (2015)]. We solve the spin-precession equations analytically on the precession time and then implement a quasiadiabatic approach to evolve these solutions on the longer radiation-reaction time. This procedure leads to an innovative “precession-averaged” post-Newtonian approach to studying precessing BBHs. We use our new solutions to classify BBH spin precession into three distinct morphologies, then investigate phase transitions between these morphologies as BBHs inspiral. These precession-averaged post-Newtonian inspirals can be efficiently calculated from arbitrarily large separations, thus making progress towards bridging the gap between astrophysics and numerical relativity.

We derive an effective potential for binary black hole (BBH) spin precession at second post-Newtonian order. This effective potential allows us to solve the orbit-averaged spin-precession equations analytically for arbitrary mass ratios and spins. These solutions are quasiperiodic functions of time: after a fixed period, the BBH spins return to their initial relative orientations and jointly precess about the total angular momentum by a fixed angle. Using these solutions, we classify BBH spin precession into three distinct morphologies between which BBHs can transition during their inspiral. We also derive a precession-averaged evolution equation for the total angular momentum that can be integrated on the radiation-reaction time and identify a new class of spin-orbit resonances that can tilt the direction of the total angular momentum during the inspiral. Our new results will help efforts to model and interpret gravitational waves from generic BBH mergers and predict the distributions of final spins and gravitational recoils.

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the early inspirals of stellar-mass black holes that will ultimately venture into the ground-based detectors’ view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA’s discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This White Paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.

We demonstrate the flexibility and utility of the Berger–Rigoutsos adaptive mesh refinement (AMR) algorithm used in the open-source numerical relativity (NR) code GRChombo for generating gravitational waveforms from binary black-hole (BH) inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary BH waveforms through a code comparison with the established NR code Lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several ‘rules of thumb’ for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.

If binary black holes form following the successive core collapses of sufficiently massive binary stars, precessional dynamics may align their spins, ${\mathbf{S}}_{\mathbf{1}}$ and ${\mathbf{S}}_{\mathbf{2}}$, and the orbital angular momentum $\mathbf{L}$ into a plane in which they jointly precess about the total angular momentum $\mathbf{J}$. These spin orientations are known as spin-orbit resonances since ${\mathbf{S}}_{\mathbf{1}}$, ${\mathbf{S}}_{\mathbf{2}}$, and $\mathbf{L}$ all precess at the same frequency to maintain their planar configuration. Two families of such spin-orbit resonances exist, differentiated by whether the components of the two spins in the orbital plane are either aligned or antialigned. The fraction of binary black holes in each family is determined by the stellar evolution of their progenitors, so if gravitational-wave detectors could measure this fraction they could provide important insights into astrophysical formation scenarios for binary black holes. In this paper, we show that even under the conservative assumption that binary black holes are observed along the direction of $\mathbf{J}$ (where precession-induced modulations to the gravitational waveforms are minimized), the waveforms of many members of each resonant family can be distinguished from all members of the other family in events with signal-to-noise ratios $\rho \simeq 10$, typical of those expected for the first detections with Advanced LIGO and Virgo. We hope that our preliminary findings inspire a greater appreciation of the capability of gravitational-wave detectors to constrain stellar astrophysics and stimulate further studies of the distinguishability of spin-orbit resonant families in more expanded regions of binary black-hole parameter space.