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

Patrick R. Brady

Patrick R. Brady 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-phastro-ph.CO

Catalog footprint

What is connected

15works
5topics
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

15 published item(s)

preprint2022arXiv

Supplementing rapid Bayesian parameter estimation schemes with adaptive grids

In the age of multi-messenger astrophysics, low-latency parameter estimation of gravitational-wave signals is essential for electromagnetic follow-up observations. In this paper, we present a new edition of the Bayesian parameter estimation scheme for compact binaries known as Rapid PE. Rapid PE parallelizes parameter estimation by fixing the intrinsic parameters (such as the masses and spins of the binary) to a grid around a search-identified intrinsic point and by integrating over the extrinsic parameters via Monte Carlo sampling. The original version of Rapid PE calculated an effective Fisher matrix to fix the intrinsic parameters to an ellipsoidal grid. Here we use rectilinear gridding in intrinsic space. The use of rectilinear intrinsic grids provides the option to more easily include spin and allows for adaptive grid refinements to mitigate search biases and map out the posterior more completely without sacrificing sampling efficiency. We demonstrate that this parameter estimation method produces reliable results for simulated gravitational-wave signals from binary neutron star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inferring kilonova population properties with a hierarchical Bayesian framework I : Non-detection methodology and single-event analyses

We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of non-detections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility (ZTF) as a single-event test case for two different prior choices of model parameters -- (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with ${\tt simsurvey}$ -- an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Machine Learning Based Source Property Inference for Compact Binary Mergers

The detection of the binary neutron star (BNS) merger, GW170817, was the first success story of multi-messenger observations of compact binary mergers. The inferred merger rate along with the increased sensitivity of the ground-based gravitational-wave (GW) network in the present LIGO/Virgo, and future LIGO/Virgo/KAGRA observing runs, strongly hints at detection of binaries which could potentially have an electromagnetic (EM) counterpart. A rapid assessment of properties that could lead to a counterpart is essential to aid time-sensitive follow-up operations, especially robotic telescopes. At minimum, the possibility of counterparts require a neutron star (NS). Also, the tidal disruption physics is important to determine the remnant matter post merger, the dynamics of which could result in the counterparts. The main challenge, however, is that the binary system parameters such as masses and spins estimated from the real time, GW template-based searches are often dominated by statistical and systematic errors. Here, we present an approach that uses supervised machine-learning to mitigate such selection effects to report possibility of counterparts based on presence of a NS component, and presence of remnant matter post merger in real time.

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extracting equation of state parameters from black hole-neutron star mergers: aligned-spin black holes and a preliminary waveform model

Information about the neutron-star equation of state is encoded in the waveform of a black hole-neutron star system through tidal interactions and the possible tidal disruption of the neutron star. During the inspiral this information depends on the tidal deformability Lambda of the neutron star, and we find that Lambda is the best measured parameter during the merger and ringdown as well. We performed 134 simulations where we systematically varied the equation of state as well as the mass ratio, neutron star mass, and aligned spin of the black hole. Using these simulations we have developed an analytic representation of the full inspiral-merger-ringdown waveform calibrated to these numerical waveforms, and we use this analytic waveform to estimate the accuracy to which Lambda can be measured with gravitational-wave detectors. We find that although the inspiral tidal signal is small, coherently combining this signal with the merger-ringdown matter effect improves the measurability of Lambda by a factor of ~3 over using just the merger-ringdown matter effect alone. However, incorporating correlations between all the waveform parameters then decreases the measurability of Lambda by a factor of ~3. The uncertainty in Lambda increases with the mass ratio, but decreases as the black hole spin increases. Overall, a single Advanced LIGO detector can measure Lambda for mass ratios Q = 2--5, black hole spins J_BH/M_BH^2 = -0.5--0.75, neutron star masses M_NS = 1.2M_sun--1.45M_sun, and an optimally oriented distance of 100Mpc to an uncertainty of ~10%--100%. For the proposed Einstein Telescope, the uncertainty in Lambda is an order of magnitude smaller.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Detecting transient gravitational waves in non-Gaussian noise with partially redundant analysis methods

There is a broad class of astrophysical sources that produce detectable, transient, gravitational waves. Some searches for transient gravitational waves are tailored to known features of these sources. Other searches make few assumptions about the sources. Typically events are observable with multiple search techniques. This work describes how to combine the results of searches that are not independent, treating each search as a classifier for a given event. This will be shown to improve the overall sensitivity to gravitational-wave events while directly addressing the problem of consistent interpretation of multiple trials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Likelihood-ratio ranking of gravitational-wave candidates in a non-Gaussian background

We describe a general approach to detection of transient gravitational-wave signals in the presence of non-Gaussian background noise. We prove that under quite general conditions, the ratio of the likelihood of observed data to contain a signal to the likelihood of it being a noise fluctuation provides optimal ranking for the candidate events found in an experiment. The likelihood-ratio ranking allows us to combine different kinds of data into a single analysis. We apply the general framework to the problem of unifying the results of independent experiments and the problem of accounting for non-Gaussian artifacts in the searches for gravitational waves from compact binary coalescence in LIGO data. We show analytically and confirm through simulations that in both cases the likelihood ratio statistic results in an improved analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extracting equation of state parameters from black hole-neutron star mergers. I. Nonspinning black holes

The late inspiral, merger, and ringdown of a black hole-neutron star (BHNS) system can provide information about the neutron-star equation of state (EOS). Candidate EOSs can be approximated by a parametrized piecewise-polytropic EOS above nuclear density, matched to a fixed low-density EOS; and we report results from a large set of BHNS inspiral simulations that systematically vary two parameters. To within the accuracy of the simulations, we find that, apart from the neutron-star mass, a single physical parameter Lambda, describing its deformability, can be extracted from the late inspiral, merger, and ringdown waveform. This parameter is related to the radius, mass, and l=2 Love number, k_2, of the neutron star by Lambda = 2k_2 R^5/3M_{NS}^5, and it is the same parameter that determines the departure from point-particle dynamics during the early inspiral. Observations of gravitational waves from BHNS inspiral thus restrict the EOS to a surface of constant Lambda in the parameter space, thickened by the measurement error. Using various configurations of a single Advanced LIGO detector, we find that neutron stars are distinguishable from black holes of the same mass and that Lambda^{1/5} or equivalently R can be extracted to 10-40% accuracy from single events for mass ratios of Q=2 and 3 at a distance of 100 Mpc, while with the proposed Einstein Telescope, EOS parameters can be extracted to accuracy an order of magnitude better.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Black-hole threshold solutions in stiff fluid collapse

Numerical studies of the gravitational collapse of a stiff (P=rho) fluid have found the now familiar critical phenomena, namely scaling of the black hole mass with a critical exponent and continuous self-similarity at the threshold of black hole formation. Using the equivalence of an irrotational stiff fluid to a massless scalar field, we construct the critical solution as a scalar field solution by making a self-similarity ansatz. We find evidence that this solution has exactly one growing perturbation mode; both the mode and the critical exponent, gamma ~ 0.94, derived from its eigenvalue agree with those measured in perfect fluid collapse simulations. We explain why this solution is seen as a critical solution in stiff fluid collapse but not in scalar field collapse, and conversely why the scalar field critical solution is not seen in stiff fluid collapse, even though the two systems are locally equivalent.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

A power filter for the detection of burst sources of gravitational radiation in interferometric detectors

We present a filter for detecting gravitational wave signals from burst sources. This filter requires only minimal advance knowledge of the expected signal: i.e. the signal's frequency band and time duration. It consists of a threshold on the total power in the data stream in the specified signal band during the specified time. This filter is optimal (in the Neyman-Pearson sense) for signal searches where only this minimal information is available.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Phases of massive scalar field collapse

We study critical behavior in the collapse of massive spherically symmetric scalar fields. We observe two distinct types of phase transition at the threshold of black hole formation. Type II phase transitions occur when the radial extent $(λ)$ of the initial pulse is less than the Compton wavelength ($μ^{-1}$) of the scalar field. The critical solution is that found by Choptuik in the collapse of massless scalar fields. Type I phase transitions, where the black hole formation turns on at finite mass, occur when $λμ\gg 1$. The critical solutions are unstable soliton stars with masses $\alt 0.6 μ^{-1}$. Our results in combination with those obtained for the collapse of a Yang-Mills field~{[M.~W. Choptuik, T. Chmaj, and P. Bizon, Phys. Rev. Lett. 77, 424 (1996)]} suggest that unstable, confined solutions to the Einstein-matter equations may be relevant to the critical point of other matter models.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Searching for periodic sources with LIGO

We investigate the computational requirements for all-sky, all-frequency searches for gravitational waves from spinning neutron stars, using archived data from interferometric gravitational wave detectors such as LIGO. These sources are expected to be weak, so the optimal strategy involves coherent accumulaton of signal-to-noise using Fourier transforms of long stretches of data (months to years). Earth-motion-induced Doppler shifts, and intrinsic pulsar spindown, will reduce the narrow-band signal-to-noise by spreading power across many frequency bins; therefore, it is necessary to correct for these effects before performing the Fourier transform. The corrections can be implemented by a parametrized model, in which one does a search over a discrete set of parameter values. We define a metric on this parameter space, which can be used to determine the optimal spacing between points in a search; the metric is used to compute the number of independent parameter-space points Np that must be searched, as a function of observation time T. The number Np(T) depends on the maximum gravitational wave frequency and the minimum spindown age tau=f/(df/dt) that the search can detect. The signal-to-noise ratio required, in order to have 99% confidence of a detection, also depends on Np(T). We find that for an all-sky, all-frequency search lasting T=10^7 s, this detection threshhold is at a level of 4 to 5 times h(3/yr), where h(3/yr) is the corresponding 99% confidence threshhold if one knows in advance the pulsar position and spin period.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

The central density of a neutron star is unaffected by a binary companion at linear order in $μ/R$

Recent numerical work by Wilson, Mathews, and Marronetti [J. R. Wilson, G. J. Mathews and P. Marronetti, Phys. Rev. D 54, 1317 (1996)] on the coalescence of massive binary neutron stars shows a striking instability as the stars come close together: Each star's central density increases by an amount proportional to 1/(orbital radius). This overwhelms any stabilizing effects of tidal coupling [which are proportional to 1/(orbital radius)^6] and causes the stars to collapse before they merge. Since the claimed increase of density scales with the stars' mass, it should also show up in a perturbation limit where a point particle of mass $μ$ orbits a neutron star. We prove analytically that this does not happen; the neutron star's central density is unaffected by the companion's presence to linear order in $μ/R$. We show, further, that the density increase observed by Wilson et. al. could arise as a consequence of not faithfully maintaining boundary conditions.

0 citations0 reviewsNo code linkedOpen access