## When will we observe binary black holes precessing?

Stephen Fairhurst, 1 Rhys Green, 1 Mark Hannam, 1 and Charlie Hoy 1

1 School of Physics and Astronomy, Cardiff University, Cardiff, CF24 3AA, United Kingdom (Dated: July 1, 2020)

After eleven gravitational-wave detections from compact-binary mergers, we are yet to observe the striking general-relativistic phenomenon of orbital precession. Measurements of precession would provide valuable insights into the distribution of black-hole spins, and therefore into astrophysical binary formation mechanisms. Using our recent two-harmonic approximation of precessing-binary signals [1], we introduce the 'precession signal-to-noise ratio', ρ p . We demonstrate that this can be used to clearly identify whether precession was measured in an observation (by comparison with both current detections and simulated signals), and can immediately quantify the measurability of precession in a given signal, which currently requires computationally expensive parameter-estimation studies. ρ p has numerous potential applications to signal searches, source-property measurements, and population studies. We give one example: assuming one possible astrophysical spin distribution, we predict that precession has a one in ∼ 25 chance of being observed in any detection.

The Advanced Laser Interferometer GravitationalWave Observatory (aLIGO) [2] and Advanced Virgo (AdV) [3], provide a unique method of observing mergers of black holes and/or neutron stars. Observations to date already provide insights into the mass and spin distributions of black holes [4, 5].

One important general relativistic effect that has not yet been observed is orbital precession. This arises when the black-hole spins are not aligned with the binary's orbital angular momentum. In contrast to Newtonian mechanics, where all angular momenta are individually conserved, in general relativity the binary's total angular momentum is (approximately) conserved, and the orbital angular momentum (and hence the orbital plane) and spins precess around it [6, 7]. This leads to modulations in the amplitude and phase of the gravitational wave (GW) signal. These are in general small effects and, in addition, whether they can be measured depends not only on the black-hole masses and spin magnitudes and directions, but also on the binary's orientation relative to the detector, and the observed GW polarization. For this reason, until now there was no straightforward way to determine how significantly precession would be imprinted onto a given waveform. The usual approach is to perform computationally expensive Bayesian analyses (see e.g. [4, 8]), but even then, the misaligned spin components (which signify whether the binary is precessing) are degenerate with other parameters, and do not provide a direct measure of precession features in the signal. This makes it difficult to infer the impact of precession measurements on the properties of astrophysical binary populations and their formation mechanisms.

The two-harmonic approximation: We take advantage of a new characterization of precession, where the gravitational waveform for a precessing bi- nary can be expressed as the sum of five harmonics, which form a power series in

$$b = \tan ( \beta / 2 )$$

where β is the opening angle between total and orbital angular momentum [1]. Large values of β are possible only when the in-plane component of spin, χ p [9] is large; the mass ratio of the system is large, or there is a significant spin component anti-aligned with the orbital angular momentum - both of which serve to reduce the relative magnitude of the orbit-aligned angular momentum. For most configurations b glyph[lessmuch] 1, and we may restrict attention to the first two harmonics and ignore terms of order b 2 or higher. It is only for binaries with β ≥ 45 , ° or equivalently b ≥ 0 4, . where the remaining harmonics become significant, and the two-harmonic approximation breaks down.

We may now write the observed waveform as a function of GW frequency as,

$$h _ { 2 \text{harm} } ( f ) \approx \mathcal { A } _ { 0 } h ^ { 0 } ( f ) + \mathcal { A } _ { 1 } h ^ { 1 } ( f ) \,, \quad \ ( 2 )$$

where A 0 and A 1 are complex, orientation dependent amplitudes and h 0 ( f ) and h 1 ( f ) are the waveforms of the two harmonics. Full details of the two-harmonic approximation are given in Ref. [1]. What is relevant to the current work is that the individual harmonics do not exhibit precession modulations; these show up due to beating between the harmonics, whose respective dephasing can be related to the frequency with which the orbital angular momentum precesses around the total angular momentum. The amplitudes A 0 and A 1 depend on the binary's orientation and polarisation, but are approximately constant throughout a given signal.

Observability of precession: Since the individual harmonics are indistinguishable from nonprecessing waveforms, it is only when two precession

harmonics can be independently observed that precession can be unambiguously identified. For precession to be observable, we therefore require that the expected signal-to-noise ratio (SNR) in both of the harmonics is above some threshold.

The expected SNR for a signal h embedded in data from a detector with a noise power spectral density S f ( ) is given as ˆ ρ 2 = ( h h | ), where the inner product is,

$$( g | h ) = 4 \, \text{Re} \int _ { f _ { o } } ^ { f _ { \max } } \frac { g ( f ) h ^ { * } ( f ) } { S ( f ) } d f \,, \quad \ ( 3 ) \quad \sin _ { \substack { \text{sign} \\ \text{in } t. } }$$

and glyph[star] denotes complex conjugation [10]. In cases with more than one precession cycle between f o and f max , the two harmonics will be close to orthogonal, ( h 0 | h 1 ) ≈ 0, and the precession SNR is simply defined as the expected SNR in the weaker harmonic:

$$\rho _ { \mathrm p } = \text{Min} ( | \mathcal { A } _ { 0 } h ^ { 0 } |, | \mathcal { A } _ { 1 } h ^ { 1 } | ) \,, \quad \quad ( 4 ) \quad \mathrm p o s _ { \mathrm p }.$$

where | h | := √ ( h h | ). When there is less than one precession cycle, it is necessary to also subtract the power aligned with the dominant harmonic, as described in detail in Ref. [1], before determining whether the second harmonic is observable.

It remains to determine a threshold above which ρ p can be considered as evidence for precession. Consider the situation where an event has been observed, so there is significant SNR in at least one harmonic. In the absence of measurable precession, and assuming well-modelled Gaussian noise, the SNR in the second harmonic will be χ 2 distributed with two degrees of freedom, where the two degrees of freedom correspond to the real and imaginary parts of the complex amplitude. Therefore, in the absence of precession, ρ p &gt; 2 1 . is expected in less than 10% of cases, and we use this as a simple threshold for observable precession. 1

In Fig. 1 we show the recovered distribution of χ p and ρ p for a number of signals, both real and simulated. For each signal, we use a nested sampling routine within the LALInference code [8] to obtain posterior probability distributions for the parameters. This is the same infrastructure that was used to measure the properties of the LIGO-Virgo observations, and we

1 A more detailed analysis would consider the volume of the binary parameter space consistent with the non-precessing and precessing parameters and use these to appropriately weight the likelihoods of the precessing and non-precessing signals. This will have some impact on the required ρ p threshold, but is unlikely to change it significantly.

present our results in the same form as in, for example, the GWTC-1 catalogue [4], by using the PESummary library [11]. The new feature is our calculation of ρ p .

First, we show the recovered χ p and ρ p distributions for a set of simulated signals, generated using the IMRPhenomPv2 model [12], each with the same choices of masses and spins - total mass M = 40 M glyph[circledot] , mass ratio 2:1, and an in-plane spin of χ p = 0 4 on . the large black hole only - but varying orientation, encoded by the angle θ between the total angular momentum and the line of sight . The distance to each signal is chosen to ensure a fixed expected SNR of 20 in the aLIGO detectors at the sensitivity of the second observing run (O2) [4], resulting in a distance variation by a factor of ≈ 3 5 between the least and most . inclined systems.

For binaries with total angular momentum closely aligned with the line of sight, θ &lt; 45 , ° the precessing SNR is consistent with no power in the h 1 . The posterior on χ p is consistent with the prior at low χ p but excludes χ p &gt; ∼ 0 7. . When θ &gt; 45 , ° the angular momentum is significantly mis-aligned with the line of sight and there is significant power in both harmonics, leading to a value of ρ p inconsistent with noise alone and little support for values of χ p &lt; ∼ 0 1. . However, using χ p alone, even after performing the parameter recovery there are no simple criteria to determine when precession is observed. A natural choice might require that the 90% confidence interval for χ p exclude zero, but this will always be the case, primarily due to the shape of the prior. Furthermore, even though we know all of the parameters a priori , it is impossible to determine whether precession will be observable without generating the waveform and performing the parameter recovery.

The precession SNR solves these problems. A value of ρ p &gt; 2 1 tells us immediately that precession is ob-. servable. Most significantly, the expected precession SNR ρ p (red lines on the middle plot) can be calculated directly from the signal parameters; no detailed parameter estimation analysis is necessary. Thus, for the first time, we are able to identify immediately whether precession would be measurable in a given configuration. We see in Fig. 1 that for each inclination, the true value for ρ p lies within the recovered 90% credible interval however the posterior is not centred around the true value. This is due to selection and prior effects. In Ref. [13], we investigate these selection effects as well as providing a detailed exploration of the observability of precession over the parameter space of masses, spins and binary orientation.

Fig. 1 also shows the distribution of ρ p for the BBH merger signals that were observed in the aLIGO and

FIG. 1. For a set of simulated signals with fixed masses and spins (see text), we show the posterior and prior (white) distributions for χ p (left), and posterior distributions for ρ p (middle) for a range of different binary orientations, θ . The grey lines show the 90% confidence regions, the solid red lines show the true values of χ p and ρ p respectively and the thick black line indicates the threshold for observable precession i.e., ρ p &gt; . 2 1. The right panel shows the ρ p distribution for the ten binary-black-hole observations in O1 and O2 [4].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00555v2-figure-8db9a105-fc72-4d34-8e9c-764e2bd3c609.jpg)

AdV O1 and O2 runs [4]. No evidence of precession was found in these signals [5], as is made clear from the recovery of ρ p . There are several cases where the distribution extends to higher values, but the median never exceeds the 2.1 threshold. These results demonstrate the efficacy of ρ p .

|      |   Aligned |   Aligned | Isotropic   |   Precessing |   Precessing |
|------|-----------|-----------|-------------|--------------|--------------|
| Low  |     0.043 |     0.644 | 0.151 0.194 |        0.173 |        0.15  |
| Flat |     0.077 |     0.448 | 0.276 0.040 |        0.327 |        0.019 |
| High |     0.105 |     0.331 | 0.354 0.013 |        0.412 |        0.005 |

When will we observe precession? We can use the observation of precession to distinguish different binary formation scenarios. The precession SNR makes it straightforward to perform an in-depth investigation of various models and identify the fraction of signals for which precession effects will be observable. Such a study was not previously possible, due to the difficulty in classifying observability of precession. Instead, limited investigations of the parameter space have been performed [14], or inferences of the distributions for the spin magnitudes and orientations obtained [15, 16], again with a limited sample size.

We investigate nine astrophysically-motivated populations of black hole binaries, comprised of three distributions of spin magnitude, and three distributions of spin orientation for the individual black holes in the binary. The spin-magnitude distributions are those used in Refs. [17-19]: low and high are triangular, peaked either at zero or extremal spin, and flat is a uniform distribution between zero and one. The spinorientation is characterized by the distribution for the angle σ between the black hole's spin and the orbital angular momentum: aligned is a triangular distribution in cos σ , which peaks at 1 and can take values 0 85 . &lt; cos σ &lt; 1 0, ( . σ &lt; ∼ 30 ); ° precessing is triangular in cos σ peaked at 0, with values -0 15 . &lt; cos σ &lt; 0 15, . (80 ° &lt; σ &lt; ∼ ∼ 100 ); ° isotropic is uniform in cos σ between -1 and 1. For each population, we generate 10 5 binaries with masses drawn from a power law mass distribution with p m ( 1 ) ∝ m -2 35 . , and p m ( 2 ) uniform in

TABLE I. The probability of observing precession, ρ p &gt; 2 1, for an observed binary (white) from each spin distribu-. tion and the probability of not observing precession in 10 random draws (grey) from each spin distribution.

m 2 between 5 M glyph[circledot] and m 1 (as in [19]), and distributed uniformly in volume and binary orientation.

Table I shows the probability of observing precession, ρ p &gt; . 2 1, in a single event drawn from each of the nine populations, observed with O2 senstivity. As expected, we are most likely to observe precession when the black holes have high spins that lie preferentially in the orbital plane ( high-precessing configurations) and least likely for black holes with low spins, or with spins preferentially aligned with the orbital angular momentum ( low-aligned configurations). Given that precession has not been observed in GW detections to date, we are able to restrict the spin distribution. Table I shows the probability of detecting ten signals with no observable precession from each of the nine spin distributions. Based on precession measurements alone, we strongly disfavour all precessing distributions. Although these are already considered astrophysically unlikely, there are models that predict preferentially inplane spins [20, 21]. We also disfavour isotropic spins with flat or high magnitudes. Thus, the lack of observed precession points towards low spins, or spins preferentially aligned with the orbital angular momentum.

Previous constraints on spins have primarily been

FIG. 2. The distribution of χ p , θ and q for observable binaries (grey), and those with measurable precession (blue), assuming a low isotropic spin distribution. θ is the inclination angle folded to [0 , π/ 2]. The y-axis labels the number of observed events in each bin, out of 10 5 simulated signals with low isotropic spins.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00555v2-figure-fed71140-1c47-4654-953c-e523182edaca.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00555v2-figure-e884add2-4daf-4049-b99f-7eac03b1700e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00555v2-figure-11d20d7d-761f-493e-a301-5dde6fade9d2.jpg)

provided by considering the measurable aligned-spin component [4, 17-19] and provide strong evidence against all but low aligned or iostropic distributions, with low isotropic spins preferred. Combining the aligned spin and precession results will further restrict the spin distribution consistent with GW observations, and will likely require spin magnitudes even smaller than our low distribution (see also Ref. [5]).

Where will we observe precession? We are able to identify, for the first time, the regions of parameter space that lead to signals with observable precession. In Fig. 2, we show the expected distribution of the precessing spin χ p , binary orientation θ and mass ratio q for observable binaries and binaries with observable precession, ρ p &gt; . 2 1, assuming a low isotropic distribution of spins. We identify clear regions of the parameter space where precession is more likely to be observed: large values of χ p , binaries that are close to edge-on, θ &gt; 45 , ◦ and systems with high mass ratio. Regions where the chance of observing precession is close to zero include binaries with χ p &lt; 0 2 or where . the total angular momentum is within 20 ° of the line of sight. These results are consistent with expectations based upon smaller studies using detailed parameter estimation techniques [14]. We also note that most observations of precession will be in comparable-mass binaries, i.e., q ≤ 2. This is a surprising, new result. It is well known that precession is more easily measured at higher mass ratios [6], which is confirmed by our study: precession is observed in &lt; 12% of detections with q &lt; 2, but &gt; 35% for q &gt; 2. However, with ∼ 90% of observations expected to have q &lt; 2, these vastly outnumber the higher-mass-ratio observations, and we find that ∼ 75% of precession observations will come from detections of binaries with q &lt; 2.

Discussion: In this letter we have used a simple method to identify when precession is measurable in a compact binary GW signal. The gravitational wave- form is well approximated by the first two harmonics in a power series expansion in the tangent of the halfopening angle [1], and the unambiguous observation of precession requires the identification of both of these harmonics in the data. The precession SNR ρ p is a simple measure of this observability. We have demonstrated the efficacy of ρ p through parameter estimation studies and also provided the distributions of ρ p for the aLIGO-AdV observations to date. Using our definition of precession SNR, we have identified how often precession will be observed for a variety of potential astrophysical spin distributions. For the most likely distribution, based on current observations ( low-aligned ) there is a 83% chance that precession will be measured after ∼ 40 observations, and is therefore likely to be observed during the current third aLIGO-AdV observing run (O3). The non-measurement of precession by the end of O3 would place much stronger constraints on spin orientations and magnitudes.

The precession SNR has many applications. Most immediately, it allows us to determine the measurability of precession in a system without performing computationally expensive parameter estimation . This allows us to, e.g., easily fold precession information into population analyses of black-hole binaries. In future work, we will explore whether the value of ρ p can be used to predict the measured χ p distribution. The precessing SNR also gives us a simple way to identify regions of the parameter space where precession is important, a necessary first step in extending existing GW searches to explicitly use precessing waveforms [22].

Note: Since the first submission of this Letter, ρ p has been used in the analysis of the aLIGO-AdV signals GW190412 [23] and GW190814 [24].

Acknowledgements We thank Eleanor Hamilton, Ian Harry, Alistair Muir, Francesco Pannarale, Vivien Raymond, Cory Thomas, and Vaibhav Tiwari for discussions. This work was supported by Sci-

ence and Technology Facilities Council (STFC) grant ST/L000962/1 and European Research Council Con-

- [1] Stephen Fairhurst, Rhys Green, Mark Hannam, Charlie Hoy, and Alistair Muir. The two-harmonic approximation for gravitational waveforms from precessing binaries. In preparation .
- [2] J. Aasi et al. Advanced LIGO. Class. Quant. Grav. , 32:074001, 2015.
- [3] F. Acernese et al. Advanced Virgo: a secondgeneration interferometric gravitational wave detector. Class. Quant. Grav. , 32(2):024001, 2015.
- [4] B. P. Abbott et al. GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs. Phys. Rev. , X9(3):031040, 2019.
- [5] B. P. Abbott et al. Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo. Astrophys. J. , 882(2):L24, 2019.
- [6] Theocharis A. Apostolatos, Curt Cutler, Gerald J. Sussman, and Kip S. Thorne. Spin induced orbital precession and its modulation of the gravitational wave forms from merging binaries. Phys. Rev. , D49:62746297, 1994.
- [7] T. A. Apostolatos. Search templates for gravitational waves from precessing, inspiraling binaries. Phys. Rev. , D52:605-620, 1995.
- [8] J. Veitch et al. Parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library. Phys. Rev. , D91(4):042003, 2015.
- [9] Patricia Schmidt, Frank Ohme, and Mark Hannam. Towards models of gravitational waveforms from generic binaries II: Modelling precession effects with a single effective precession parameter. Phys. Rev. , D91(2):024043, 2015.
- [10] Bruce Allen, Warren G. Anderson, Patrick R. Brady, Duncan A. Brown, and Jolien D. E. Creighton. FINDCHIRP: An Algorithm for detection of gravitational waves from inspiraling compact binaries. Phys. Rev. , D85:122006, 2012.
- [11] Charlie Hoy and Vivien Raymond. PESummary: the code agnostic Parameter Estimation Summary page builder. 6 2020.
- [12] Mark Hannam, et al. Simple Model of Complete Pre-

solidator Grant 647839. Computational resources were provided by Cardiff University, funded by STFC support for UK involvement in aLIGO operations.

cessing Black-Hole-Binary Gravitational Waveforms. Phys. Rev. Lett. , 113(15):151101, 2014.

- [13] Rhys Green, et al. Identifying when precession can be measured in gravitational waveforms. In preparation .
- [14] Salvatore Vitale, Ryan Lynch, John Veitch, Vivien Raymond, and Riccardo Sturani. Measuring the spin of black holes in binary systems using gravitational waves. Phys. Rev. Lett. , 112(25):251101, 2014.
- [15] Daniel Wysocki, Jacob Lange, and Richard O'Shaughnessy. Reconstructing phenomenological distributions of compact binaries via gravitational wave observations. Phys. Rev. D , 100(4):043012, 2019.
- [16] Colm Talbot and Eric Thrane. Determining the population properties of spinning black holes. Phys. Rev. , D96(2):023012, 2017.
- [17] Will M. Farr, et al. Distinguishing Spin-Aligned and Isotropic Black Hole Populations With Gravitational Waves. Nature , 548:426, 2017.
- [18] Ben Farr, Daniel E. Holz, and Will M. Farr. Using Spin to Understand the Formation of LIGO and Virgos Black Holes. Astrophys. J. , 854(1):L9, 2018.
- [19] Vaibhav Tiwari, Stephen Fairhurst, and Mark Hannam. Constraining black-hole spins with gravitational wave observations. Astrophys. J. , 868(2):140, 2018.
- [20] Fabio Antonini, Carl L. Rodriguez, Cristobal Petrovich, and Caitlin L. Fischer. Precessional dynamics of black hole triples: binary mergers with near-zero effective spin. Mon. Not. Roy. Astron. Soc. , 480(1):L58L62, 2018.
- [21] Carl L. Rodriguez and Fabio Antonini. A Triple Origin for the Heavy and Low-Spin Binary Black Holes Detected by LIGO/Virgo. Astrophys. J. , 863(1):7, 2018.
- [22] Ian Harry, Stephen Privitera, Alejandro Boh, and Alessandra Buonanno. Searching for Gravitational Waves from Compact Binaries with Precessing Spins. Phys. Rev. , D94(2):024012, 2016.
- [23] GW190412: Observation of a Binary-Black-Hole Coalescence with Asymmetric Masses. 4 2020.
- [24] R. Abbott et al. GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object. Astrophys. J. , 896(2):L44, 2020.
Stephen Fairhurst Rhys Green Mark Hannam Charlie Hoy
When will we observe binary black holes precessing?
After eleven gravitational-wave detections from compact-binary mergers, we
are yet to observe the striking general-relativistic phenomenon of orbital
precession. Measurements of precession would provide valuable insights into the
distribution of black-hole spins, and therefore into astrophysical binary
formation mechanisms. Using our recent two-harmonic approximation of
precessing-binary signals~\cite{Fairhurst:2019_2harm}, we introduce the
``precession signal-to-noise ratio'', $\rho_p$. We demonstrate that this can be
used to clearly identify whether precession was measured in an observation (by
comparison with both current detections and simulated signals), and can
immediately quantify the measurability of precession in a given signal, which
currently requires computationally expensive parameter-estimation studies.
$\rho_p$ has numerous potential applications to signal searches,
source-property measurements, and population studies. We give one example:
assuming one possible astrophysical spin distribution, we predict that
precession has a one in $\sim 25$ chance of being observed in any detection.