BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Reed Essick

Reed Essick appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.IMastro-ph.EPastro-ph.SRnucl-thphysics.ins-det

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

16 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
preprint2016arXiv

Orbital decay of hot Jupiters due to nonlinear tidal dissipation within solar-type hosts

We study the orbital evolution of hot Jupiters due to the excitation and damping of tidally driven $g$-modes within solar-type host stars. Linearly resonant $g$-modes (the dynamical tide) are driven to such large amplitudes in the stellar core that they excite a sea of other $g$-modes through weakly nonlinear interactions. By solving the dynamics of large networks of nonlinearly coupled modes, we show that the nonlinear dissipation rate of the dynamical tide is several orders of magnitude larger than the linear dissipation rate. We find stellar tidal quality factors $Q \approx 10^5-10^6$ for systems with planet mass $M_p \geq 0.5 M_\mathrm{J}$ and orbital period $P \leq 2$ days, which implies that such systems decay on timescales that are small compared to the main-sequence lifetime of their solar-type hosts. According to our results, there are $\approx 10$ currently known exoplanetary systems, including WASP-19b and HAT-P-36-b, with orbital decay timescales shorter than a Gyr. Rapid, tidally induced orbital decay may explain the observed paucity of planets with $M_p \ge M_\mathrm{J}$ and $P < 2$ days around solar-type hosts and could generate detectable transit-timing variations in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Localization of short duration gravitational-wave transients with the early advanced LIGO and Virgo detectors

The Laser Interferometer Gravitational wave Observatory (LIGO) and Virgo, advanced ground-based gravitational-wave detectors, will begin collecting science data in 2015. With first detections expected to follow, it is important to quantify how well generic gravitational-wave transients can be localized on the sky. This is crucial for correctly identifying electromagnetic counterparts as well as understanding gravitational-wave physics and source populations. We present a study of sky localization capabilities for two search and parameter estimation algorithms: \emph{coherent WaveBurst}, a constrained likelihood algorithm operating in close to real-time, and \emph{LALInferenceBurst}, a Markov chain Monte Carlo parameter estimation algorithm developed to recover generic transient signals with latency of a few hours. Furthermore, we focus on the first few years of the advanced detector era, when we expect to only have two (2015) and later three (2016) operational detectors, all below design sensitivity. These detector configurations can produce significantly different sky localizations, which we quantify in detail. We observe a clear improvement in localization of the average detected signal when progressing from two-detector to three-detector networks, as expected. Although localization depends on the waveform morphology, approximately 50% of detected signals would be imaged after observing 100-200 deg$^2$ in 2015 and 60-110 deg$^2$ in 2016, although knowledge of the waveform can reduce this to as little as 22 deg$^2$. This is the first comprehensive study on sky localization capabilities for generic transients of the early network of advanced LIGO and Virgo detectors, including the early LIGO-only two-detector configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Application of machine learning algorithms to the study of noise artifacts in gravitational-wave data

The sensitivity of searches for astrophysical transients in data from the LIGO is generally limited by the presence of transient, non-Gaussian noise artifacts, which occur at a high-enough rate such that accidental coincidence across multiple detectors is non-negligible. Furthermore, non-Gaussian noise artifacts typically dominate over the background contributed from stationary noise. These "glitches" can easily be confused for transient gravitational-wave signals, and their robust identification and removal will help any search for astrophysical gravitational-waves. We apply Machine Learning Algorithms (MLAs) to the problem, using data from auxiliary channels within the LIGO detectors that monitor degrees of freedom unaffected by astrophysical signals. The number of auxiliary-channel parameters describing these disturbances may also be extremely large; an area where MLAs are particularly well-suited. We demonstrate the feasibility and applicability of three very different MLAs: Artificial Neural Networks, Support Vector Machines, and Random Forests. These classifiers identify and remove a substantial fraction of the glitches present in two very different data sets: four weeks of LIGO's fourth science run and one week of LIGO's sixth science run. We observe that all three algorithms agree on which events are glitches to within 10% for the sixth science run data, and support this by showing that the different optimization criteria used by each classifier generate the same decision surface, based on a likelihood-ratio statistic. Furthermore, we find that all classifiers obtain similar limiting performance, suggesting that most of the useful information currently contained in the auxiliary channel parameters we extract is already being used.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimizing Vetoes for Gravitational-wave Transient Searches

Interferometric gravitational-wave detectors like LIGO, GEO600 and Virgo record a surplus of information above and beyond possible gravitational-wave events. These auxiliary channels capture information about the state of the detector and its surroundings which can be used to infer potential terrestrial noise sources of some gravitational-wave-like events. We present an algorithm addressing the ordering (or equivalently optimizing) of such information from auxiliary systems in gravitational-wave detectors to establish veto conditions in searches for gravitational-wave transients. The procedure was used to identify vetoes for searches for unmodelled transients by the LIGO and Virgo collaborations during their science runs from 2005 through 2007. In this work we present the details of the algorithm; we also use a limited amount of data from LIGO's past runs in order to examine the method, compare it with other methods, and identify its potential to characterize the instruments themselves. We examine the dependence of Receiver Operating Characteristic curves on the various parameters of the veto method and the implementation on real data. We find that the method robustly determines important auxiliary channels, ordering them by the apparent strength of their correlations to the gravitational-wave channel. This list can substantially reduce the background of noise events in the gravitational-wave data. In this way it can identify the source of glitches in the detector as well as assist in establishing confidence in the detection of gravitaional-wave transients.

0 citations0 reviewsNo code linkedOpen access