Researcher profile

Daniel E. Holz

Daniel E. Holz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

False positive control in time series coincidence detection

We study the problem of coincidence detection in time series data, where we aim to determine whether the appearance of simultaneous or near-simultaneous events in two time series is indicative of some shared underlying signal or synchronicity, or might simply be due to random chance. This problem arises across many applications, such as astrophysics (e.g., detecting astrophysical events such as gravitational waves, with two or more detectors) and neuroscience (e.g., detecting synchronous firing patterns between two or more neurons). In this work, we consider methods based on time-shifting, where the timeline of one data stream is randomly shifted relative to another, to mimic the types of coincidences that could occur by random chance. Our theoretical results establish rigorous finite-sample guarantees controlling the probability of false positives, under weak assumptions that allow for dependence within the time series data, providing reassurance that time-shifting methods are a reliable tool for inference in this setting. Empirical results with simulated and real data validate the strong performance of time-shifting methods in dependent-data settings.

preprint2022arXiv

Avoiding a Cluster Catastrophe: Retention Efficiency and the Binary Black Hole Mass Spectrum

The population of binary black hole mergers identified through gravitational waves has uncovered unexpected features in the intrinsic properties of black holes in the universe. One particularly surprising and exciting result is the possible existence of black holes in the pair-instability mass gap, $\sim50-120~M_\odot$. Dense stellar environments can populate this region of mass space through hierarchical mergers, with the retention efficiency of black hole merger products strongly dependent on the escape velocity of the host environment. We use simple toy models to represent hierarchical merger scenarios in various dynamical environments. We find that hierarchical mergers in environments with high escape velocities $(\gtrsim300~\mathrm{km\,s}^{-1})$ are efficiently retained. If such environments dominate the binary black hole merger rate, this would lead to an abundance of high-mass mergers that is potentially incompatible with the empirical mass spectrum from the current catalog of binary black hole mergers. Models that efficiently generate hierarchical mergers, and contribute significantly to the observed population, must therefore be tuned to avoid a "cluster catastrophe" of overproducing binary black hole mergers within and above the pair-instability mass gap.

preprint2022arXiv

Bridging the Gap: Categorizing Gravitational-Wave Events at the Transition Between Neutron Stars and Black Holes

We search for features in the mass distribution of detected compact binary coalescences which signify the transition between neutron stars and black holes. We analyze all gravitational wave detections by LIGO-Virgo-KAGRA made through the end of the first half of the third observing run, and find clear evidence for two different populations of compact objects based solely on gravitational wave data. We confidently (99.3%) find a steepening relative to a single power law describing neutron stars and low-mass black holes below $2.4^{+0.5}_{-0.5}\text{ M}_{\odot}$, which is consistent with many predictions for the maximum neutron star mass. We find suggestions of the purported lower mass gap between the most massive neutron stars and the least massive black holes, but are unable to conclusively resolve it with current data. If it exists, we find the lower mass gap's edges to lie at $2.2^{+0.7}_{-0.5}\text{ M}_{\odot}$ and $6.0^{+2.4}_{-1.4}\text{ M}_{\odot}$. We re-examine events that have been deemed "exceptional" by the LIGO-Virgo-KAGRA collaborations in the context of these features. We analyze GW190814 self-consistently in the context of the full population of compact binaries, finding support for its secondary to be either a neutron star or a lower mass gap object, consistent with previous claims. Our models are the first to accommodate this event, which is an outlier with respect to the binary black hole population. We find that GW200105 and GW200115 probe the edges of, and may have components within, the lower mass gap. As future data improve global population models, the classification of these events will also improve.

preprint2022arXiv

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Fundamental Physics and Beyond the Standard Model

Gravitational wave detectors are formidable tools to explore strong-field gravity, especially black holes and neutron stars. These compact objects are extraordinarily efficient at producing electromagnetic and gravitational radiation. As such, they are ideal laboratories for fundamental physics and have an immense discovery potential. The detection of black hole binaries by third-generation Earth-based detectors, space-based detectors and pulsar timing arrays will provide exquisite tests of general relativity. Loud "golden" events and extreme mass-ratio inspirals can strengthen the observational evidence for horizons by mapping the exterior spacetime geometry, inform us on possible near-horizon modifications, and perhaps reveal a breakdown of Einstein's gravity. Measurements of the black-hole spin distribution and continuous gravitational-wave searches can turn black holes into efficient detectors of ultralight bosons across ten or more orders of magnitude in mass. A precise monitoring of the phase of inspiralling binaries can constrain the existence of additional propagating fields and characterize the environment in which the binaries live, bounding the local dark matter density and properties. Gravitational waves from compact binaries will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address fundamental issues in our current understanding of the cosmos.

preprint2022arXiv

Spectral sirens: cosmology from the full mass distribution of compact binaries

We explore the use of the mass spectrum of neutron stars and black holes in gravitational-wave compact binary sources as a cosmological probe. These standard siren sources provide direct measurements of luminosity distance. In addition, features in the mass distribution, such as mass gaps or peaks, will redshift, and thus provide independent constraints on their redshift distribution. We argue that the entire mass spectrum should be utilized to provide cosmological constraints. For example, we find that the mass spectrum of LIGO--Virgo--KAGRA events introduces at least five independent mass ``features'': the upper and lower edges of the pair instability supernova (PISN) gap, the upper and lower edges of the neutron star--black hole gap, and the minimum neutron star mass. We find that although the PISN gap dominates the cosmological inference with current detectors (2G), as shown in previous work, it is the lower mass gap that will provide the most powerful constraints in the era of Cosmic Explorer and Einstein Telescope (3G). By using the full mass distribution, we demonstrate that degeneracies between mass evolution and cosmological evolution can be broken, unless an astrophysical conspiracy shifts all features of the full mass distribution simultaneously following the (non-trivial) Hubble diagram evolution. We find that this self-calibrating ``spectral siren'' method has the potential to provide precision constraints of both cosmology and the evolution of the mass distribution, with 2G achieving better than $10\%$ precision on $H(z)$ at $z\lesssim1$ within a year, and 3G reaching $\lesssim1\%$ at $z\gtrsim2$ within one month.

preprint2022arXiv

Target of Opportunity Observations of Gravitational Wave Events with Vera C. Rubin Observatory

The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events.

preprint2021arXiv

Distance measures in gravitational-wave astrophysics and cosmology

We present quantities which characterize the sensitivity of gravitational-wave observatories to sources at cosmological distances. In particular, we introduce and generalize the horizon, range, response, and reach distances. These quantities incorporate a number of important effects, including cosmologically well-defined distances and volumes, cosmological redshift, cosmological time dilation, and rate density evolution. In addition, these quantities incorporate unique aspects of gravitational wave detectors, such as the variable sky sensitivity of the detectors and the scaling of the sensitivity with inverse distance. An online calculator (https://users.rcc.uchicago.edu/~dholz/gwc/) and python notebook (https://github.com/hsinyuc/distancetool) to determine GW distances are available. We provide answers to the question: "How far can gravitational-wave detectors hear?"

preprint2021arXiv

Please repeat: Strong lensing of gravitational waves as a probe of compact binary and galaxy populations

Strong gravitational lensing of gravitational wave sources offers a novel probe of both the lens galaxy and the binary source population. In particular, the strong lensing event rate and the time delay distribution of multiply-imaged gravitational-wave binary coalescence events can be used to constrain the mass distribution of the lenses as well as the intrinsic properties of the source population. We calculate the strong lensing event rate for a range of second (2G) and third generation (3G) detectors, including Advanced LIGO/Virgo, A+, Einstein Telescope (ET), and Cosmic Explorer (CE). For 3G detectors, we find that {$\sim0.1\%$} of observed events are expected to be strongly lensed. We predict detections of {$\sim 1$} lensing pair per year with A+, and {$\sim 50$} pairs {per year} with ET/CE. These rates are highly sensitive to the characteristic galaxy velocity dispersion, $σ_*$, implying that observations of the rates will be a sensitive probe of lens properties. We explore using the time delay distribution between multiply-imaged gravitational-wave sources to constrain properties of the lenses. We find that 3G detectors would constrain $σ_*$ to {$\sim21\%$ after 5 years}. Finally, we show that the presence or absence of strong lensing {within the detected population} provides useful insights into the source redshift and mass distribution out to redshifts beyond the peak of the star formation rate, which can be used to constrain formation channels and their relation to the star formation rate and delay time distributions for these systems.

preprint2021arXiv

Probing Extremal Gravitational-Wave Events with Coarse-Grained Likelihoods

As catalogs of gravitational-wave transients grow, new records are set for the most extreme systems observed to date. The most massive observed black holes probe the physics of pair instability supernovae while providing clues about the environments in which binary black hole systems are assembled. The least massive black holes, meanwhile, allow us to investigate the purported neutron star-black hole mass gap, and binaries with unusually asymmetric mass ratios or large spins inform our understanding of binary and stellar evolution. Existing outlier tests generally implement leave-one-out analyses, but these do not account for the fact that the event being left out was by definition an extreme member of the population. This results in a bias in the evaluation of outliers. We correct for this bias by introducing a coarse-graining framework to investigate whether these extremal events are true outliers or whether they are consistent with the rest of the observed population. Our method enables us to study extremal events while testing for population model misspecification. We show that this ameliorates biases present in the leave-one-out analyses commonly used within the gravitational-wave community. Applying our method to results from the second LIGO--Virgo transient catalog, we find qualitative agreement with the conclusions of Abbott et al, ApJL 913 L7 (2021). GW190814 is an outlier because of its small secondary mass. We find that neither GW190412 nor GW190521 are outliers.

preprint2021arXiv

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for &#39;&#39;red&#39;&#39; kilonovae and rule out &#39;&#39;blue&#39;&#39; kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

preprint2020arXiv

Black Hole Coagulation: Modeling Hierarchical Mergers in Black Hole Populations

Data from the LIGO and Virgo detectors has confirmed that stellar-mass black holes can merge within a Hubble time, leaving behind massive remnant black holes. In some astrophysical environments such as globular clusters and AGN disks, it may be possible for these remnants to take part in further compact-object mergers, producing a population of hierarchically formed black holes. In this work, we present a parameterized framework for describing the population of binary black hole mergers, while self-consistently accounting for hierarchical mergers. The framework casts black holes as particles in a box which can collide based on an effective cross-section, but allows inputs from more detailed astrophysical simulations. Our approach is relevant to any population which is comprised of second or higher generation black holes, such as primordial black holes or dense cluster environments. We describe some possible inputs to this generic model and their effects on the black hole merger populations, and use the model to perform Bayesian inference on the catalog of black holes from LIGO and Virgo&#39;s first two observing runs. We find that models with a high rate of hierarchical mergers are disfavored, consistent with previous population analyses. Future gravitational-wave events will further constrain the inputs to this generic hierarchical merger model, enabling a deeper look into the formation environments of binary black holes.

preprint2020arXiv

Cosmological Inference using Gravitational Wave Standard Sirens: A Mock Data Challenge

The observation of binary neutron star merger GW170817, along with its optical counterpart, provided the first constraint on the Hubble constant $H_0$ using gravitational wave standard sirens. When no counterpart is identified, a galaxy catalog can be used to provide the necessary redshift information. However, the true host might not be contained in a catalog which is not complete out to the limit of gravitational-wave detectability. These electromagnetic and gravitational-wave selection effects must be accounted for. We describe and implement a method to estimate $H_0$ using both the counterpart and the galaxy catalog standard siren methods. We perform a series of mock data analyses using binary neutron star mergers to confirm our ability to recover an unbiased estimate of $H_0$. Our simulations used a simplified universe with no redshift uncertainties or galaxy clustering, but with different magnitude-limited catalogs and assumed host galaxy properties, to test our treatment of both selection effects. We explore how the incompleteness of catalogs affects the final measurement of $H_0$, as well as the effect of weighting each galaxy&#39;s likelihood of being a host by its luminosity. In our most realistic simulation, where the simulated catalog is about three times denser than the density of galaxies in the local universe, we find that a 4.4\% measurement precision can be reached using galaxy catalogs with 50\% completeness and $\sim 250$ binary neutron star detections with sensitivity similar to that of Advanced LIGO&#39;s second observing run.

preprint2020arXiv

Does Matter Matter? Using the mass distribution to distinguish neutron stars and black holes

Gravitational-wave detectors have opened a new window through which we can observe black holes (BHs) and neutron stars (NSs). Analyzing the 11 detections from LIGO/Virgo&#39;s first gravitational-wave catalog, GWTC-1, we investigate whether the power-law fit to the BH mass spectrum can also accommodate the binary neutron star (BNS) event GW170817, or whether we require an additional feature, such as a mass gap, in between the NS and BH populations. We find that with respect to the power-law fit to binary black hole (BBH) masses, GW170817 is an outlier at the 0.13\% level, suggesting a distinction between NS and BH masses. A single power-law fit across the entire mass range is in mild tension with: (a) the detection of one source in the BNS mass range ($\sim 1$--$2.5 \,M_\odot$), (b) the absence of detections in the &#34;mass-gap&#34; range ($\sim 2.5$--$5 \,M_\odot$), and (c) the detection of 10 sources in the BBH mass range ($\gtrsim 5 \,M_\odot$). Instead, the data favor models with a feature between NS and BH masses, including a mass gap (Bayes factor of 4.6) and a break in the power law, with a steeper slope at NS masses compared to BH masses (91\% credibility). We estimate the merger rates of compact binaries based on our fit to the global mass distribution, finding $\mathcal{R}_\mathrm{BNS} = 871^{+3015}_{-805} \ \mathrm{Gpc}^{-3} \ \mathrm{yr}^{-1}$ and $\mathcal{R}_\mathrm{BBH} = 47.5^{+57.9}_{-28.8} \ \mathrm{Gpc}^{-3} \ \mathrm{yr}^{-1}$. We conclude that, even in the absence of any prior knowledge of the difference between NSs and BHs, the gravitational-wave data alone already suggest two distinct populations of compact objects.

preprint2020arXiv

Picky Partners: The Pairing of Component Masses in Binary Black Hole Mergers

We examine the relationship between individual black hole (BH) masses in merging binary black hole (BBH) systems. Analyzing the ten BBH detections from LIGO/Virgo&#39;s first two observing runs, we find that the masses of the component BHs comprising each binary are unlikely to be randomly drawn from the same underlying distribution. Instead, the two BHs of a given binary prefer to be of comparable mass. We show that it is $\sim 5$ times more likely that the component BHs in a given binary are always equal (to within 5\%) than that they are randomly paired. If we assume that the probability of a merger between two BHs scales with the mass ratio $q$ as $q^β$, so that $β=0$ corresponds to random pairings, we find $β>0$ is favored at credibility 0.987. By modeling the mass distribution, we find that the median mass ratio is $q_{50\%} = 0.91^{+0.05}_{-0.17}$ at 90\% credibility. While the pairing between BHs depends on their mass ratio, we find no evidence that it depends on the total mass of the system: it is $\sim 6$ times more likely that the pairing depends purely on the mass ratio than on the total mass. We predict that 99\% of BBHs detected by LIGO/Virgo will have mass ratios $q > 0.5$. We conclude that merging black holes do not form random pairings; instead they are selective about their partners, preferring to mate with black holes of a similar mass. The details of these selective pairings provide insight into the underlying formation channels of merging binaries.

preprint2019arXiv

Counting on Short Gamma-Ray Bursts: Gravitational-Wave Constraints of Jet Geometry

The detection of GW170817 in gravitational waves and gamma rays revealed that short gamma-ray bursts are associated with the merger of neutron-stars. Gamma rays are thought to result from the formation of collimated jets, but the details of this process continue to elude us. One fundamental observable is the emission profile of the jet as a function of viewing angle. We present two methods to measure the effective angular width, $θ_B$, of short gamma-ray burst (sGRB) jets using gravitational wave and gamma-ray data, assuming all sGRBs have the same angular dependence for their luminosities. The first is a counting experiment, where we combine the known detection thresholds of the LIGO/Virgo and Fermi Gamma Ray Burst Monitor detectors to infer parameters of systems that are detected in gravitational waves. This method requires minimal knowledge about each event, beyond whether or not they were detected in gamma-rays. The second method uses additional information from the gravitational-wave and electromagnetic data to estimate parameters of the source, and thereby improve constraints on jet properties. Applying our methods to GW170817, we find only weak constraints on the sGRB luminosity profile, with statistical uncertainty dominating differences between models. We also analyze simulated events from future observing runs, and find that with 5 and 100 BNS detections, the counting method constrains the relative uncertainty in $θ_B$ to within 51% and 12%, respectively. Incorporating gravitational-wave parameter estimation would further tighten these constraints to 43% and 9.6%. In the limit of many detections, incorporating parameter estimation achieves only marginal improvements; we conclude that the majority of the information about jet structure comes from the relative sensitivities of gravitational-wave and gamma-ray detectors as encoded in simple counting experiments.

preprint2019arXiv

Nonparametric Inference of Neutron Star Composition, Equation of State, and Maximum Mass with GW170817

The detection of GW170817 in gravitational waves provides unprecedented constraints on the equation of state (EOS) of the ultra-dense matter within the cores of neutron stars (NSs). We extend the nonparametric analysis first introduced in Landry & Essick (2019), and confirm that GW170817 favors soft EOSs. We infer macroscopic observables for a canonical 1.4 $M_{\odot}$ NS, including the tidal deformability $Λ_{1.4} = 211^{+312}_{-137}$ ($491^{+216}_{-181}$) and radius $R_{1.4}= 10.86^{+2.04}_{-1.42}$ ($12.51^{+1.00}_{-0.88}$) km, as well as the maximum mass for nonrotating NSs, $M_{max} = 2.064^{+0.260}_{-1.34}$ ($2.017^{0.238}_{-0.087}$) $M_\odot$, with nonparametric priors loosely (tightly) constrained to resemble candidate EOSs from the literature. Furthermore, we find weak evidence that GW170817 involved at least one NS based on gravitational-wave data alone ($B^{NS}_{BBH}= 3.3 \pm 1.4$), consistent with the observation of electromagnetic counterparts. We also investigate GW170817&#39;s implications for the maximum spin frequency of millisecond pulsars, and find that the fastest known pulsar is spinning at more than 50% of its breakup frequency at 90% confidence. We additionally find modest evidence in favor of quark matter within NSs, and GW170817 favors the presence of at least one disconnected hybrid star branch in the mass--radius relation over a single stable branch by a factor of 2. Assuming there are multiple stable branches, we find a suggestive posterior preference for a sharp softening around nuclear density followed by stiffening around twice nuclear density, consistent with a strong first-order phase transition. While the statistical evidence in favor of new physics within NS cores remains tenuous with GW170817 alone, these tantalizing hints reemphasize the promise of gravitational waves for constraining the supranuclear EOS.

preprint2019arXiv

The Most Massive Binary Black Hole Detections and the Identification of Population Outliers

Advanced LIGO and Virgo detected ten binary black holes (BBHs) in their first two observing runs (O1 and O2). Analysis of these events found strong evidence for a dearth of BBHs with component masses greater than $\sim45 \ M_\odot$, as would be expected from a pair-instability mass gap. Meanwhile, a standalone analysis of the merger GW170729 found its primary mass $m_1 = {51.2^{+16.2}_{-11.0} \ M_\odot}$, with the majority of its posterior support at $m_1 > 45 \ M_\odot$. Although this appears to be in contradiction with the existence of a limit at $\sim45\ M_\odot$, we argue that individual events cannot be evaluated without reference to the entire population. When GW170729 is analyzed jointly with the rest of the detections, as part of a full hierarchical population analysis, its inferred primary mass tightens considerably, to $m_1 = {38.9^{+7.3}_{-4.5} \ M_\odot}$. For a large sample of events in the presence of noise, apparent outliers in the detected distribution are inevitable, even if the underlying population forbids outliers. We discuss methods of distinguishing between statistical fluctuations and population outliers using posterior predictive tests. Applying these tests to the primary mass distribution in O1 and O2, we find that the ten detections are consistent with even the simplest power-law plus maximum-mass model considered by the LVC. This supports the claim that GW170729 is not a population outlier. We also provide non-parametric constraints on the rate of high-mass mergers and conservatively bound the rate of mergers with $m_1 > 45 \ M_\odot$ at $2.8^{+5.4}_{-2.0}\%$ of the total merger rate. After 100 detections like those of O1 and O2 from a population with a maximum primary mass of $45 \, M_\odot$, it would be common for the most massive system to have an observed maximum-likelihood mass $m_1 \gtrsim 70 \, M_\odot$.