BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Reed Essick

Reed Essick contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEastro-ph.IMastro-ph.SRnucl-thphysics.ins-det

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
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

12 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Calibration Uncertainty's Impact on Gravitational-Wave Observations

Our ability to calibrate current kilometer-scale interferometers can potentially confound the inference of astrophysical signals. Current calibration uncertainties are well described by a Gaussian process. I exploit this description to analytically examine the impact of calibration uncertainty. I derive closed-form expressions for the conditioned likelihood of the calibration error given the observed data and an astrophysical signal (astrophysical calibration) as well as for the marginal likelihood for the data given a signal (integrated over the calibration uncertainty). I show that calibration uncertainty always reduces search sensitivity and the amount of information available about astrophysical signals. Additionally, calibration uncertainty will fundamentally limit the precision to which loud signals can be constrained, a crucial factor when considering the scientific potential of proposed third-generation interferometers. For example, I estimate that with $1\%$ uncertainty in the detector response's amplitude and phase, one will only be able to measure the leading-order tidal parameter ($\tildeΛ$) for a 1.4+1.4$\,M_\odot$ system to better than $\pm 1$ ($\sim 0.2\%$ relative uncertainty) for signals with signal-to-noise ratios $\gtrsim 10^4$. At this signal-to-noise ratio, calibration uncertainty increases $σ_{\tildeΛ}$ by a factor of $2$ compared to stationary Gaussian noise alone. Furthermore, 1\% calibration uncertainty limits the precision to always be $σ_{\tildeΛ} \gtrsim 0.5$. I also show how to best select the frequencies at which calibration should be precisely constrained in order to minimize the information lost about astrophysical parameters. It is not necessary to constrain the calibration errors to be small at all frequencies to perform precise astrophysical inference for individual signals.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Implicit correlations within phenomenological parametric models of the neutron star equation of state

The rapid increase in the number and precision of astrophysical probes of neutron stars in recent years allows for the inference of their equation of state. Observations target different macroscopic properties of neutron stars which vary from star to star, such as mass and radius, but the equation of state allows for a common description of all neutron stars. To connect these observations and infer the properties of dense matter and neutron stars simultaneously, models for the equation of state are introduced. Parametric models rely on carefully engineered functional forms that reproduce a large array of realistic equations of state. Such models benefit from their simplicity but are limited because any finite-parameter model cannot accurately approximate all possible equations of state. Nonparametric models overcome this by increasing model freedom at the cost of increased complexity. In this study, we compare common parametric and nonparametric models, quantify the limitations of the former, and study the impact of modeling on our current understanding of high-density physics. We show that parametric models impose strongly model-dependent, and sometimes opaque, correlations between density scales. Such interdensity correlations result in tighter constraints that are unsupported by data and can lead to biased inference of the equation of state and of individual neutron star properties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Precision Requirements for Monte Carlo Sums within Hierarchical Bayesian Inference

Hierarchical Bayesian inference is often conducted with estimates of the target distribution derived from Monte Carlo sums over samples from separate analyses of parts of the hierarchy or from mock observations used to estimate sensitivity to a target population. We investigate requirements on the number of Monte Carlo samples needed to guarantee the estimator of the target distribution is precise enough that it does not affect the inference. We consider probabilistic models of how Monte Carlo samples are generated, showing that the finite number of samples introduces additional uncertainty as they act as an imperfect encoding of the components of the hierarchical likelihood. Additionally, we investigate the behavior of estimators marginalized over approximate measures of the uncertainty, comparing their performance to the Monte Carlo point estimate. We find that correlations between the estimators at nearby points in parameter space are crucial to the precision of the estimate. Approximate marginalization that neglects these correlations will either introduce a bias within the inference or be more expensive (require more Monte Carlo samples) than an inference constructed with point estimates. We therefore recommend that hierarchical inferences with empirically estimated target distributions use point estimates.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Selection Effects in Periodic X-ray Data from Maximizing Detection Statistics

The Neutron Star Interior Composition Explorer (NICER) records data of exceptional quality on the energy-dependent X-ray pulse profile of pulsars. However, in searching for evidence of pulsations, Guillot et al. (2019) introduce a procedure to select an ordered subset of the data that maximizes a detection statistic (the H-test). I show that this procedure can degrade subsequent analyses using an idealized model with a stationary expected count rates from both noise and signal. Specifically, the data-selection procedure biases the inferred mean count rate to be too low, biases the inferred pulsation amplitude to be too high, and that the size of these biases scales strongly with the amount of data that is rejected and the true signal amplitude. The procedure also alters the null-distribution of the H-test rendering nominal detection significance estimates overly optimistic. While the idealized model does not capture all the complexities of real NICER data, it suggests that these biases could be important for NICER's observations of J0740+6620 and other faint pulsars (observations of J0030+0451 are likely less affected). I estimate that these selection effects may introduce a bias of $\mathcal{O}(10\%)$ on average in the inferred modulation depth of lightcurves like J0740+6620's, and may be as large as $\mathcal{O}(50\%)$ for fainter pulsars. However, the change for a single dataset like J0740+6620 is expected to be a shift between $-5\%$ and $+20\%$. This could imply that the lower limit on J0740+6620's radius is slightly larger than it should be, although preliminary investigations suggest the change in the radius constraints are $\mathcal{O}(1\%)$ with real J0740+6620 data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Coincidence Null Test for Poisson-Distributed Events

When transient events are observed with multiple sensors, it is often necessary to establish the significance of coincident events. We derive a universal null test for an arbitrary number of sensors motivated by the archetypal detection problem for independent Poisson-distributed events in gravitational-wave detectors such as LIGO and Virgo. In these detectors, transient events may be witnessed by myriad channels that record interferometric signals and the surrounding physical environment. We apply our null test to a broad set of simulated gravitational-wave events as well as to a real gravitational-wave detection to determine which auxiliary channels do and do not witness real gravitational waves, and therefore which are safe to use when constructing vetoes. We also describe how our approach can be used to study detector artifacts and their origin, as well as to quantify the statistical independence of candidate GW signals from noise artifacts observed in auxiliary channels.

0 citations0 reviewsNo code linkedOpen access
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'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 "mass-gap" 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enhancing Gravitational-Wave Science with Machine Learning

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

iDQ: Statistical Inference of Non-Gaussian Noise with Auxiliary Degrees of Freedom in Gravitational-Wave Detectors

Gravitational-wave detectors are exquisitely sensitive instruments and routinely enable ground-breaking observations of novel astronomical phenomena. However, they also witness non-stationary, non-Gaussian noise that can be mistaken for astrophysical sources, lower detection confidence, or simply complicate the extraction of signal parameters from noisy data. To address this, we present iDQ, a supervised learning framework to autonomously detect noise artifacts in gravitational-wave detectors based only on auxiliary degrees of freedom insensitive to gravitational waves. iDQ has operated in low latency throughout the advanced detector era at each of the two LIGO interferometers, providing invaluable data quality information about each detection to date in real-time. We document the algorithm, describing the statistical framework and possible applications within gravitational-wave searches. In particular, we construct a likelihood-ratio test that simultaneously accounts for the presence of non-Gaussian noise artifacts and utilizes information from both the observed gravitational-wave strain signal and thousands of auxiliary degrees of freedom. We also present several examples of iDQ's performance with modern interferometers, showing iDQ's ability to autonomously reproduce known data quality monitors and identify noise artifacts not flagged by other analyses.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access