Source author record

Salvatore Vitale

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

40 published item(s)

preprint2022arXiv

A Standard Siren Cosmological Measurement from the Potential GW190521 Electromagnetic Counterpart ZTF19abanrhr

The identification of the electromagnetic counterpart candidate ZTF19abanrhr to the binary black hole merger GW190521 opens the possibility to infer cosmological parameters from this standard siren with a uniquely identified host galaxy. The distant merger allows for cosmological inference beyond the Hubble constant. Here we show that the three-dimensional spatial location of ZTF19abanrhr calculated from the electromagnetic data remains consistent with the updated sky localization of GW190521 provided by the LIGO-Virgo Collaboration. If ZTF19abanrhr is associated with the GW190521 merger and assuming a flat wCDM model we find that $H_0 =48^{+24}_{-10}$ km/s/Mpc, $Ω_m =0.39^{+0.38}_{-0.29}$, and $w_0 = -1.29^{+0.63}_{-0.50}$ (median and 68% credible interval). If we use the Hubble constant value inferred from another gravitational-wave event, GW170817, as a prior for our analysis, together with assumption of a flat $Λ$CDM and the model-independent constraint on the physical matter density $ω_m$ from Planck, we find $H_0 = 69.1^{8.7}_{-6.0}$ km/s/Mpc.

preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

preprint2022arXiv

An Infrared Search for Kilonovae with the WINTER Telescope. I. Binary Neutron Star Mergers

The Wide-Field Infrared Transient Explorer (WINTER) is a new 1 $\text{deg}^2$ seeing-limited time-domain survey instrument designed for dedicated near-infrared follow-up of kilonovae from binary neutron star (BNS) and neutron star-black hole mergers. WINTER will observe in the near-infrared Y, J, and short-H bands (0.9-1.7 microns, to $\text{J}_{AB}=21$ magnitudes) on a dedicated 1-meter telescope at Palomar Observatory. To date, most prompt kilonova follow-up has been in optical wavelengths; however, near-infrared emission fades more slowly and depends less on geometry and viewing angle than optical emission. We present an end-to-end simulation of a follow-up campaign during the fourth observing run (O4) of the LIGO, Virgo, and KAGRA interferometers, including simulating 625 BNS mergers, their detection in gravitational waves, low-latency and full parameter estimation skymaps, and a suite of kilonova lightcurves from two different model grids. We predict up to five new kilonovae independently discovered by WINTER during O4, given a realistic BNS merger rate. Using a larger grid of kilonova parameters, we find that kilonova emission is $\approx$2 times longer-lived and red kilonovae are detected $\approx$1.5 times further in the infrared than in the optical. For 90% localization areas smaller than 150 (450) $\rm{deg}^{2}$, WINTER will be sensitive to more than 10% of the kilonova model grid out to 350 (200) Mpc. We develop a generalized toolkit to create an optimal BNS follow-up strategy with any electromagnetic telescope and present WINTER's observing strategy with this framework. This toolkit, all simulated gravitational-wave events, and skymaps are made available for use by the community.

preprint2022arXiv

Constraining high-redshift stellar-mass primordial black holes with next-generation ground-based gravitational-wave detectors

The possible existence of primordial black holes in the stellar mass window has received considerable attention because their mergers may contribute to current and future gravitational-wave detections. Primordial black hole mergers, together with mergers of black holes originating from Population~III stars, are expected to dominate at high redshifts ($z\gtrsim 10$). However the primordial black hole merger rate density is expected to rise monotonically with redshift, while Population~III mergers can only occur after the birth of the first stars. Next-generation gravitational-wave detectors such as Cosmic Explorer~(CE) and Einstein Telescope~(ET) can access this distinctive feature in the merger rates as functions of redshift, allowing for a direct measurement of the abundance of the two populations, and hence for robust constraints on the abundance of primordial black holes. We simulate four-months worth of data observed by a CE-ET detector network and perform hierarchical Bayesian analysis to recover the merger rate densities. We find that if the Universe has no primordial black holes with masses of $\mathcal{O}(10M_{\odot})$, the projected upper limit on their abundance $f_{\rm PBH}$ as a fraction of dark matter energy density may be as low as $f_{\rm PBH}\sim \mathcal{O}({10^{-5}})$, about two orders of magnitude lower than current upper limits in this mass range. If instead $f_{\rm PBH}\gtrsim 10^{-4}$, future gravitational wave observations would exclude $f_{\rm PBH}=0$ at the 95\% credible interval.

preprint2022arXiv

Evidence of large recoil velocity from a black hole merger signal

The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to $v_f \sim 1542^{+747}_{-1098}$ km/s or $v_f \gtrsim 698$ km/s (one-sided limit), at 90\% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than $0.48\%$ ($7.7\%$) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.

preprint2022arXiv

Hints of spin-orbit resonances in the binary black hole population

Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component black holes within the orbital plane, $ϕ_1$ and $ϕ_2$, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that $ϕ_1$ and $ϕ_2$ are best measured near the merger of the two black holes. In this work, we use these spin measurements to provide the first constraints on the full six-dimensional spin distribution of merging binary black holes. In particular, we find that there is a preference for $Δϕ= ϕ_1 - ϕ_2 \sim \pm π$ in the population, which can be a signature of spin-orbit resonances. We also find a preference for $ϕ_1 \sim -π/4$ with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, the strength of these preferences depends on our prior choices, and we are unable to constrain the widths of the $ϕ_1$ and $Δϕ$ distributions. Therefore, more observations are necessary to confirm the features we find. Finally, we derive constraints on the distribution of recoil kicks in the population, and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters. We make our spin and kick population constraints publicly available.

preprint2022arXiv

Measuring binary black hole orbital-plane spin orientations

Binary black hole spins are among the key observables for gravitational wave astronomy. Among the spin parameters, their orientations within the orbital plane, $ϕ_1$, $ϕ_2$ and $Δϕ=ϕ_1-ϕ_2$, are critical for understanding the prevalence of the spin-orbit resonances and merger recoils in binary black holes. Unfortunately, these angles are particularly hard to measure using current detectors, LIGO and Virgo. Because the spin directions are not constant for precessing binaries, the traditional approach is to measure the spin components at some reference stage in the waveform evolution, typically the point at which the frequency of the detected signal reaches 20 Hz. However, we find that this is a poor choice for the orbital-plane spin angle measurements. Instead, we propose measuring the spins at a fixed dimensionless time or frequency near the merger. This leads to significantly improved measurements for $ϕ_1$ and $ϕ_2$ for several gravitational wave events. Furthermore, using numerical relativity injections, we demonstrate that $Δϕ$ will also be better measured near the merger for louder signals expected in the future. Finally, we show that numerical relativity surrogate models are key for reliably measuring the orbital-plane spin orientations, even at moderate signal-to-noise ratios like $\sim 30-45$.

preprint2022arXiv

Searching for a subpopulation of primordial black holes in LIGO/Virgo gravitational-wave data

With several dozen binary black hole events detected by LIGO/Virgo to date and many more expected in the next few years, gravitational-wave astronomy is shifting from individual-event analyses to population studies. Using the GWTC-2 catalog, we perform a hierarchical Bayesian analysis that for the first time combines several state-of-the-art astrophysical formation models with a population of primordial black holes (PBHs) and constrains the fraction of a putative subpopulation of PBHs in the data. We find that this fraction depends significantly on the set of assumed astrophysical models. While a primordial population is statistically favored against certain competitive astrophysical channels, such as globular clusters and nuclear stellar clusters, a dominant contribution from the stable-mass-transfer isolated formation channel drastically reduces the need for PBHs, except for explaining the rate of mass-gap events like GW190521. The tantalizing possibility that black holes formed after inflation are contributing to LIGO/Virgo observations could only be verified by further reducing uncertainties in astrophysical and primordial formation models, and it may ultimately be confirmed by third-generation interferometers.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Future Gravitational-Wave Detector Facilities

The next generation of gravitational-wave observatories can explore a wide range of fundamental physics phenomena throughout the history of the universe. These phenomena include access to the universe's binary black hole population throughout cosmic time, to the universe's expansion history independent of the cosmic distance ladders, to stochastic gravitational-waves from early-universe phase transitions, to warped space-time in the strong-field and high-velocity limit, to the equation of state of nuclear matter at neutron star and post-merger densities, and to dark matter candidates through their interaction in extreme astrophysical environments or their interaction with the detector itself. We present the gravitational-wave detector concepts than can drive the future of gravitational-wave astrophysics. We summarize the status of the necessary technology, and the research needed to be able to build these observatories in the 2030s.

preprint2022arXiv

The binary black hole spin distribution likely broadens with redshift

The population-level distributions of the masses, spins, and redshifts of binary black holes (BBHs) observed using gravitational waves can shed light on how these systems form and evolve. Because of the complex astrophysical processes shaping the inferred BBH population, models allowing for correlations among these parameters will be necessary to fully characterize these sources. We hierarchically analyze the BBH population detected by LIGO and Virgo with a model allowing for correlations between the effective aligned spin and the primary mass and redshift. We find that the width of the effective spin distribution grows with redshift at 98.6% credibility. We determine this trend to be robust under the application of several alternative models and additionally verify that such a correlation is unlikely to be spuriously introduced using a simulated population. We discuss the possibility that this correlation could be due to a change in the natal black hole spin distribution with redshift.

preprint2022arXiv

The effect of spin mismodeling on gravitational-wave measurements of the binary neutron star mass distribution

The binary neutron star (BNS) mass distribution measured with gravitational-wave observations has the potential to reveal information about the dense matter equation of state, supernova physics, the expansion rate of the universe, and tests of General Relativity. As most current gravitational-wave analyses measuring the BNS mass distribution do not simultaneously fit the spin distribution, the implied population-level spin distribution is the same as the spin prior applied when analyzing individual sources. In this work, we demonstrate that introducing a mismatch between the implied and true BNS spin distributions can lead to biases in the inferred mass distribution. This is due to the strong correlations between the measurements of the mass ratio and spin components aligned with the orbital angular momentum for individual sources. We find that applying a low-spin prior which excludes the true spin magnitudes of some sources in the population leads to significantly overestimating the maximum neutron star mass and underestimating the minimum neutron star mass at the population level with as few as six BNS detections. The safest choice of spin prior that does not lead to biases in the inferred mass distribution is one which allows for high spin magnitudes and tilts misaligned with the orbital angular momentum.

preprint2022arXiv

The orientations of the binary black holes in GWTC-3

It is expected that the orbital planes of gravitational-wave (GW) sources are isotropically distributed. However, both physical and technical factors, such as alternate theories of gravity with birefringence, catalog contamination, and search algorithm limitations, could result in inferring a non-isotropic distribution. Showing that the inferred astrophysical distribution of the orbital orientations is indeed isotropic can thus be used to rule out some violations of general relativity, as a null test about the purity of the GW catalog sample, and as a check that selection effects are being properly accounted for. We augment the default mass/spins/redshift model used by the LIGO-Virgo-KAGRA Collaboration in their most recent analysis to also measure the astrophysical distribution of orbital orientations. We show that the 69 binary black holes in GWTC-3 are consistent with having random orbital orientations. The inferred distribution is highly symmetric around $π/2$, with skewness $\mathcal{S}_{\rm{post}}=0.01^{+0.17}_{-0.17}$. Meanwhile, the median of the inferred distribution has a Jensen-Shannon divergence of $1.4\times 10^{-4}$ bits when compared to the expected isotropic distribution.

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

Measuring the star formation rate with gravitational waves from binary black holes

A measurement of the history of cosmic star formation is central to understand the origin and evolution of galaxies. The measurement is extremely challenging using electromagnetic radiation: significant modeling is required to convert luminosity to mass, and to properly account for dust attenuation, for example. Here we show how detections of gravitational waves from inspiraling binary black holes made by proposed third-generation detectors can be used to measure the star formation rate (SFR) of massive stars with high precision up to redshifts of ~10. Depending on the time-delay model, the predicted detection rates ranges from ~2310 to ~56,740 per month with the current measurement of local merger rate density. With 30,000 detections, parameters describing the volumetric SFR can be constrained at the few percent level, and the volumetric merger rate can be directly measured to 3% at z ~ 2. Given a parameterized SFR, the characteristic delay time between binary formation and merger can be measured to ~60%.

preprint2021arXiv

On the single-event-based identification of primordial black hole mergers at cosmological distances

The existence of primordial black holes (PBHs), which may form from the collapse of matter overdensities shortly after the Big Bang, is still under debate. Among the potential signatures of PBHs are gravitational waves (GWs) emitted from binary black hole (BBH) mergers at redshifts $z\gtrsim 30$, where the formation of astrophysical black holes is unlikely. Future ground-based GW detectors, Cosmic Explorer and Einstein Telescope, will be able to observe equal-mass BBH mergers with total mass of $\mathcal{O}(10-100)~M_{\odot}$ at such distances. In this work, we investigate whether the redshift measurement of a single BBH source can be precise enough to establish its primordial origin. We simulate BBHs of different masses, mass ratios and orbital orientations. We show that for BBHs with total masses between $20~M_{\odot}$ and $40~M_{\odot}$ merging at $z \geq 40$ one can infer $z>30$ at up to 97\% credibility, with a network of one Einstein Telescope, one 40-km Cosmic Explorer in the US and one 20-km Cosmic Explorer in Australia. A smaller network made of one Einstein Telescope and one 40-km Cosmic Explorer in the US measures $z>30$ at larger than 90\% credibility for roughly half of the sources than the larger network. We then assess the dependence of this result on the Bayesian redshift priors used for the analysis, specifically on the relative abundance of the BBH mergers originated from the first stars, and the primordial BBH mergers.

preprint2020arXiv

Astrophysical implications of GW190412 as a remnant of a previous black-hole merger

Two of the dominant channels to produce merging stellar-mass black-hole binaries are believed to be the isolated evolution of binary stars in the field and dynamical formation in star clusters. The first reported black-hole binary event from the third LIGO/Virgo observing run (GW190412) is unusual in that it has unequal masses, nonzero effective spin, and nonzero primary spin at 90\% confidence interval. We show that this event should be exceedingly rare in the context of both the field and cluster formation scenarios. Interpreting GW190412 as a remnant of a previous black-hole merger provides a promising route to explain its features. If GW190412 indeed formed hierarchically, we show that the region of the parameter space that is best motivated from an astrophysical standpoint (low natal spins and light clusters) cannot accommodate the observation. We analyze public GW190412 LIGO/Virgo data with a Bayesian prior where the more massive black hole resulted from a previous merger, and find that this interpretation is equally supported by the data. If the heavier component of GW190412 is indeed a merger remnant, then its spin magnitude is $χ_1=0.56_{-0.21}^{+0.19}$, which is higher than the value previously reported by the LIGO/Virgo collaboration.

preprint2020arXiv

Constraining black hole mimickers with gravitational wave observations

LIGO and Virgo have recently observed a number of gravitational wave (GW) signals that are fully consistent with being emitted by binary black holes described by general relativity. However, there are theoretical proposals of exotic objects that can be massive and compact enough to be easily confused with black holes. Nevertheless, these objects differ from black holes in having nonzero tidal deformabilities, which can allow one to distinguish binaries containing such objects from binary black holes using GW observations. Using full Bayesian parameter estimation, we investigate the possibility of constraining the parameter space of such "black hole mimickers" with upcoming GW observations. Employing perfect fluid stars with a polytropic equation of state as a simple model that can encompass a variety of possible black hole mimickers, we show how the observed masses and tidal deformabilities of a binary constrain the equation of state. We also show how such constraints can be used to rule out some simple models of boson stars.

preprint2020arXiv

Constraining short gamma-ray burst jet properties with gravitational waves and gamma rays

Gamma-ray burst (GRB) prompt emission is highly beamed, and understanding the jet geometry and beaming configuration can provide information on the poorly understood central engine and circum-burst environment. Prior to the advent of gravitational-wave astronomy, astronomers relied on observations of jet breaks in the multi-wavelength afterglow to determine the GRB opening angle, since the observer's viewing angle relative to the system cannot be determined from the electromagnetic data alone. Gravitational-wave observations, however, provide an independent measurement of the viewing angle. We describe a Bayesian method for determining the geometry of short GRBs using coincident electromagnetic and gravitational-wave observations. We demonstrate how an ensemble of multi-messenger detections can be used to measure the distributions of the jet energy, opening angle, Lorentz factor, and angular profile of short GRBs; we find that for a population of 100 such observations, we can constrain the mean of the opening angle distribution to within 10 degrees regardless of the angular emission profile. Conversely, the constraint on the energy distribution depends on the shape of the profile, which can be distinguished.

preprint2020arXiv

Machine-learning non-stationary noise out of gravitational wave detectors

Signal extraction out of background noise is a common challenge in high precision physics experiments, where the measurement output is often a continuous data stream. To improve the signal to noise ratio of the detection, witness sensors are often used to independently measure background noises and subtract them from the main signal. If the noise coupling is linear and stationary, optimal techniques already exist and are routinely implemented in many experiments. However, when the noise coupling is non-stationary, linear techniques often fail or are sub-optimal. Inspired by the properties of the background noise in gravitational wave detectors, this work develops a novel algorithm to efficiently characterize and remove non-stationary noise couplings, provided there exist witnesses of the noise source and of the modulation. In this work, the algorithm is described in its most general formulation, and its efficiency is demonstrated with examples from the data of the Advanced LIGO gravitational wave observatory, where we could obtain an improvement of the detector gravitational wave reach without introducing any bias on the source parameter estimation.

preprint2020arXiv

Quantifying the Effect of Power Spectral Density Uncertainty on Gravitational-Wave Parameter Estimation for Compact Binary Sources

In order to perform Bayesian parameter estimation to infer the source properties of gravitational waves from compact binary coalescences (CBCs), the noise characteristics of the detector must be understood. It is typically assumed that the detector noise is stationary and Gaussian, characterized by a power spectral density (PSD) that is measured with infinite precision. We present a new method to incorporate the uncertainty in the power spectral density estimation into the Bayesian inference of the binary source parameters and apply it to the first 11 CBC detections reported by the LIGO- Virgo Collaboration. We find that incorporating the PSD uncertainty only leads to variations in the positions and widths of the binary parameter posteriors on the order of a few percent. Our results are publicly available for download on git [1].

preprint2020arXiv

You Can't Always Get What You Want: The Impact of Prior Assumptions on Interpreting GW190412

GW190412 is the first observation of a black hole binary with definitively unequal masses. GW190412's mass asymmetry, along with the measured positive effective inspiral spin, allowed for inference of a component black hole spin: the primary black hole in the system was found to have a dimensionless spin magnitude between 0.17 and 0.59 (90% credible range). We investigate how the choice of priors for the spin magnitudes and tilts of the component black holes affect the robustness of parameter estimates for GW190412, and report Bayes factors across a suite of prior assumptions. Depending on the waveform family used to describe the signal, we find either marginal to moderate (2:1-6:1) or strong ($\gtrsim$ 20:1) support for the primary black hole being spinning compared to cases where only the secondary is allowed to have spin. We show how these choices influence parameter estimates, and find the asymmetric masses and positive effective inspiral spin of GW190412 to be qualitatively, but not quantitatively, robust to prior assumptions. Our results highlight the importance of both considering astrophysically motivated or population-based priors in interpreting observations and considering their relative support from the data.

preprint2016arXiv

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 2015. We study an astrophysically motivated population of sources (binary components with masses $1.2~\mathrm{M}_\odot$--$1.6~\mathrm{M}_\odot$ and spins of less than $0.05$) using the full LIGO analysis pipeline. While this simulated population covers the observed range of potential binary neutron-star sources, we do not exclude the possibility of sources with parameters outside these ranges; given the existing uncertainty in distributions of mass and spin, it is critical that analyses account for the full range of possible mass and spin configurations. We find that conservative prior assumptions on neutron-star mass and spin lead to average fractional uncertainties in component masses of $\sim 16\%$, with little constraint on spins (the median $90\%$ upper limit on the spin of the more massive component is $\sim 0.7$). Stronger prior constraints on neutron-star spins can further constrain mass estimates, but only marginally. However, we find that the sky position and luminosity distance for these sources are not influenced by the inclusion of spin; therefore, if LIGO detects a low-spin population of BNS sources, less computationally expensive results calculated neglecting spin will be sufficient for guiding electromagnetic follow-up.

preprint2015arXiv

Constraining the neutron star equation of state with gravitational wave signals from coalescing binary neutron stars

Recently exploratory studies were performed on the possibility of constraining the neutron star equation of state (EOS) using signals from coalescing binary neutron stars, or neutron star-black hole systems, as they will be seen in upcoming advanced gravitational wave detectors such as Advanced LIGO and Advanced Virgo. In particular, it was estimated to what extent the combined information from multiple detections would enable one to distinguish between different equations of state through hypothesis ranking or parameter estimation. Under the assumption of zero neutron star spins both in signals and in template waveforms and considering tidal effects to 1 post-Newtonian (1PN) order, it was found that O(20) sources would suffice to distinguish between a hard, moderate, and soft equation of state. Here we revisit these results, this time including neutron star tidal effects to the highest order currently known, termination of gravitational waveforms at the contact frequency, neutron star spins, and the resulting quadrupole-monopole interaction. We also take the masses of neutron stars in simulated sources to be distributed according to a relatively strongly peaked Gaussian, as hinted at by observations, but without assuming that the data analyst will necessarily have accurate knowledge of this distribution for use as a mass prior. We find that especially the effect of the latter is dramatic, necessitating many more detections to distinguish between different EOS and causing systematic biases in parameter estimation, on top of biases due to imperfect understanding of the signal model pointed out in earlier work. This would get mitigated if reliable prior information about the mass distribution could be folded into the analyses.

preprint2015arXiv

Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era

Advanced ground-based gravitational-wave (GW) detectors begin operation imminently. Their intended goal is not only to make the first direct detection of GWs, but also to make inferences about the source systems. Binary neutron-star mergers are among the most promising sources. We investigate the performance of the parameter-estimation \edit{(PE)} pipeline that will be used during the first observing run of the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) in 2015: we concentrate on the ability to reconstruct the source location on the sky, but also consider the ability to measure masses and the distance. Accurate, rapid sky-localization is necessary to alert electromagnetic (EM) observatories so that they can perform follow-up searches for counterpart transient events. We consider PE accuracy in the presence of \edit{non-stationary}, non-Gaussian noise. We find that the character of the noise makes negligible difference to the PE performance \edit{at a given signal-to-noise ratio}. The source luminosity distance can only be poorly constrained, the median $90\%$ ($50\%$) credible interval scaled with respect to the true distance is $0.85$ ($0.38$). However, the chirp mass is well measured. Our chirp-mass estimates are subject to systematic error because we used gravitational-waveform templates without component spin to carry out inference on signals with moderate spins, but the total error is typically less than $10^{-3} M_\odot$. The median $90\%$ ($50\%$) credible region for sky localization is $\sim600~\mathrm{deg^{2}}$ ($\sim150~\mathrm{deg^{2}}$), with $3\%$ ($30\%$) of detected events localized within $100~\mathrm{deg^{2}}$. Early aLIGO, with only two detectors, will have a sky-localization accuracy for binary neutron stars of hundreds of square degrees; this makes EM follow-up challenging, but not impossible.

preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

preprint2014arXiv

Effect of squeezing on parameter estimation of gravitational waves emitted by compact binary systems

The LIGO gravitational wave (GW) detectors will begin collecting data in 2015, with Virgo following shortly after. The use of squeezing has been proposed as a way to reduce the quantum noise without increasing the laser power, and has been successfully tested at one of the LIGO sites and at GEO in Germany. When used in Advanced LIGO without a filter cavity, the squeezer improves the performances of detectors above about 100 Hz, at the cost of a higher noise floor in the low frequency regime. Frequency-dependent squeezing, on the other hand, will lower the noise floor throughout the entire band. Squeezing technology will have a twofold impact: it will change the number of expected detections and it will impact the quality of parameter estimation for the detected signals. In this work we consider three different GW detector networks, each utilizing a different type of squeezer, all corresponding to plausible implementations. Using LALInference, a powerful Monte Carlo parameter estimation algorithm, we study how each of these networks estimates the parameters of GW signals emitted by compact binary systems, and compare the results with a baseline advanced LIGO-Virgo network. We find that, even in its simplest implementation, squeezing has a large positive impact: the sky error area of detected signals will shrink by about 30% on average, increasing the chances of finding an electromagnetic counterpart to the GW detection. Similarly, we find that the measurability of tidal deformability parameters for neutron stars in binaries increases by about 30%, which could aid in determining the equation of state of neutron stars. The degradation in the measurement of the chirp mass, as a result of the higher low-frequency noise, is shown to be negligible when compared to systematic errors.

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.

preprint2014arXiv

Measuring the spin of black holes in binary systems using gravitational waves

Compact binary coalescences are the most promising sources of gravitational waves (GWs) for ground based detectors. Binary systems containing one or two spinning black holes are particularly interesting due to spin-orbit (and eventual spin-spin) interactions, and the opportunity of measuring spins directly through GW observations. In this letter we analyze simulated signals emitted by spinning binaries with several values of masses, spins, orientation, and signal-to-noise ratio. We find that spin magnitudes and tilt angles can be estimated to accuracy of a few percent for neutron star--black hole systems and $\sim$ 5-30% for black hole binaries. In contrast, the difference in the azimuth angles of the spins, which may be used to check if spins are locked into resonant configurations, cannot be constrained. We observe that the best performances are obtained when the line of sight is perpendicular to the system's total angular momentum, and that a sudden change of behavior occurs when a system is observed from angles such that the plane of the orbit can be seen both from above and below during the time the signal is in band. This study suggests that the measurement of black hole spin by means of GWs can be as precise as what can be obtained from X-ray binaries.

preprint2014arXiv

Prospects for doubling the range of Advanced LIGO

In the coming years, the gravitational wave community will be optimizing detector performance for a variety of astrophysical sources that make competing demands on the detector sensitivity in different frequency bands. In this paper we describe a number of technologies that are being developed as anticipated upgrades to the Advanced LIGO detector, and quantify the potential sensitivity improvement they offer. Specifically, we consider squeezed light injection for reduction of quantum noise, detector design and materials changes which reduce thermal noise, and mirrors with significantly increased mass. We explore how each of these technologies impacts the detection of the most promising gravitational wave sources, and suggest an effective progression of upgrades which culminate in a factor of two broadband sensitivity improvement.

preprint2014arXiv

Sensitivity Comparison of Searches for Binary Black Hole Coalescences with Ground-based Gravitational-Wave Detectors

Searches for gravitational-wave transients from binary black hole coalescences typically rely on one of two approaches: matched filtering with templates and morphology-independent excess power searches. Multiple algorithmic implementations in the analysis of data from the first generation of ground-based gravitational wave interferometers have used different strategies for the suppression of non-Gaussian noise transients, and targeted different regions of the binary black hole parameter space. In this paper we compare the sensitivity of three such algorithms: matched filtering with full coalescence templates, matched filtering with ringdown templates and a morphology-independent excess power search. The comparison is performed at a fixed false alarm rate and relies on Monte-carlo simulations of binary black hole coalescences for spinning, non-precessing systems with total mass 25-350 solar mass, which covers the parameter space of stellar mass and intermediate mass black hole binaries. We find that in the mass range of 25 -100 solar mass the sensitive distance of the search, marginalized over source parameters, is best with matched filtering to full waveform templates, to within 10 percent at a false alarm rate of 3 events per year. In the mass range of 100-350 solar mass, the same comparison favors the morphology-independent excess power search to within 20 percent. The dependence on mass and spin is also explored.

preprint2014arXiv

The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

We anticipate the first direct detections of gravitational waves (GWs) with Advanced LIGO and Virgo later this decade. Though this groundbreaking technical achievement will be its own reward, a still greater prize could be observations of compact binary mergers in both gravitational and electromagnetic channels simultaneously. During Advanced LIGO and Virgo's first two years of operation, 2015 through 2016, we expect the global GW detector array to improve in sensitivity and livetime and expand from two to three detectors. We model the detection rate and the sky localization accuracy for binary neutron star (BNS) mergers across this transition. We have analyzed a large, astrophysically motivated source population using real-time detection and sky localization codes and higher-latency parameter estimation codes that have been expressly built for operation in the Advanced LIGO/Virgo era. We show that for most BNS events the rapid sky localization, available about a minute after a detection, is as accurate as the full parameter estimation. We demonstrate that Advanced Virgo will play an important role in sky localization, even though it is anticipated to come online with only one-third as much sensitivity as the Advanced LIGO detectors. We find that the median 90% confidence region shrinks from ~500 square degrees in 2015 to ~200 square degrees in 2016. A few distinct scenarios for the first LIGO/Virgo detections emerge from our simulations.

preprint2014arXiv

TIGER: A data analysis pipeline for testing the strong-field dynamics of general relativity with gravitational wave signals from coalescing compact binaries

The direct detection of gravitational waves with upcoming second-generation gravitational wave detectors such as Advanced LIGO and Virgo will allow us to probe the genuinely strong-field dynamics of general relativity (GR) for the first time. We present a data analysis pipeline called TIGER (Test Infrastructure for GEneral Relativity), which is designed to utilize detections of compact binary coalescences to test GR in this regime. TIGER is a model-independent test of GR itself, in that it is not necessary to compare with any specific alternative theory. It performs Bayesian inference on two hypotheses: the GR hypothesis $\mathcal{H}_{\rm GR}$, and $\mathcal{H}_{\rm modGR}$, which states that one or more of the post-Newtonian coefficients in the waveform are not as predicted by GR. By the use of multiple sub-hypotheses of $\mathcal{H}_{\rm modGR}$, in each of which a different number of parameterized deformations of the GR phase are allowed, an arbitrarily large number of 'testing parameters' can be used without having to worry about a model being insufficiently parsimonious if the true number of extra parameters is in fact small. TIGER is well-suited to the regime where most sources have low signal-to-noise ratios, again through the use of these sub-hypotheses. Information from multiple sources can trivially be combined, leading to a stronger test. We focus on binary neutron star coalescences, for which sufficiently accurate waveform models are available that can be generated fast enough on a computer to be fit for use in Bayesian inference. We show that the pipeline is robust against a number of fundamental, astrophysical, and instrumental effects, such as differences between waveform approximants, a limited number of post-Newtonian phase contributions being known, the effects of neutron star spins and tidal deformability on the orbital motion, and instrumental calibration errors.

preprint2013arXiv

Demonstrating the feasibility of probing the neutron star equation of state with second-generation gravitational wave detectors

Fisher matrix and related studies have suggested that with second-generation gravitational wave detectors, it may be possible to infer the equation of state of neutron stars using tidal effects in binary inspiral. Here we present the first fully Bayesian investigation of this problem. We simulate a realistic data analysis setting by performing a series of numerical experiments of binary neutron star signals hidden in detector noise, assuming the projected final design sensitivity of the Advanced LIGO- Virgo network. With an astrophysical distribution of events (in particular, uniform in co-moving volume), we find that only a few tens of detections will be required to arrive at strong constraints, even for some of the softest equations of state in the literature. Thus, direct gravitational wave detection will provide a unique probe of neutron star structure.

preprint2013arXiv

How serious can the stealth bias be in gravitational wave parameter estimation?

The upcoming direct detection of gravitational waves will open a window to probing the strong-field regime of general relativity (GR). As a consequence, waveforms that include the presence of deviations from GR have been developed (e.g. in the parametrized post-Einsteinian approach). TIGER, a data analysis pipeline which builds Bayesian evidence to support or question the validity of GR, has been written and tested. In particular, it was shown recently that data from the LIGO and Virgo detectors will allow to detect deviations from GR smaller than can be probed with Solar System tests and pulsar timing measurements or not accessible with conventional tests of GR. However, evidence from several detections is required before a deviation from GR can be confidently claimed. An interesting consequence is that, should GR not be the correct theory of gravity in its strong field regime, using standard GR templates for the matched filter analysis of interferometer data will introduce biases in the gravitational wave measured parameters with potentially disastrous consequences on the astrophysical inferences, such as the coalescence rate or the mass distribution. We consider three heuristic possible deviations from GR and show that the biases introduced by assuming GR's validity manifest in various ways. The mass parameters are usually the most affected, with biases that can be as large as $30$ standard deviations for the symmetric mass ratio, and nearly one percent for the chirp mass, which is usually estimated with sub-percent accuracy. We conclude that statements about the nature of the observed sources, e.g. if both objects are neutron stars, depend critically on the explicit assumption that GR it the right theory of gravity in the strong field regime.

preprint2013arXiv

Reconstructing the sky location of gravitational-wave detected compact binary systems: methodology for testing and comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

preprint2012arXiv

Effect of calibration errors on Bayesian parameter estimation for gravitational wave signals from inspiral binary systems in the advanced detectors era

By 2015 the advanced versions of the gravitational-wave detectors Virgo and LIGO will be online. They will collect data in coincidence with enough sensitivity to potentially deliver multiple detections of gravitation waves from inspirals of compact-object binaries. This work is focused on understanding the effects introduced by uncertainties in the calibration of the interferometers. We consider plausible calibration errors based on estimates obtained during LIGO's fifth and Virgo's third science runs, which include frequency-dependent amplitude errors of $\sim 10%$ and frequency-dependent phase errors of $\sim 3$ degrees in each instrument. We quantify the consequences of such errors estimating the parameters of inspiraling binaries. We find that the systematics introduced by calibration errors on the inferred values of the chirp mass and mass ratio are smaller than 20% of the statistical measurement uncertainties in parameter estimation for 90% of signals in our mock catalog. Meanwhile, the calibration-induced systematics in the inferred sky location of the signal are smaller than $\sim 50%$ of the statistical uncertainty. We thus conclude that calibration-induced errors at this level are not a significant detriment to accurate parameter estimation.

preprint2011arXiv

Application of asymptotic expansions for maximum likelihood estimators' errors to gravitational waves from binary mergers: the network case

This paper describes the most accurate analytical frequentist assessment to date of the uncertainties in the estimation of physical parameters from gravitational waves generated by non spinning binary systems and Earth-based networks of laser interferometers. The paper quantifies how the accuracy in estimating the intrinsic parameters mostly depends on the network signal to noise ratio (SNR), but the resolution in the direction of arrival also strongly depends on the network geometry. We compare results for 6 different existing and possible global networks and two different choices of the parameter space. We show how the fraction of the sky where the one sigma angular resolution is below 2 square degrees increases about 3 times when transitioning from the Hanford (USA), Livingston (USA) and Cascina (Italy) network to possible 5 sites ones (while keeping the network SNR fixed). The technique adopted here is an asymptotic expansion of the uncertainties in inverse powers of the signal to noise ratio where the first order is the inverse Fisher information matrix. We show that a common approach to use simplified parameter spaces and only the Fisher information matrix can largely underestimate the uncertainties (by a factor ~7 for the one sigma sky uncertainty in square degrees at a network SNR of ~15).

preprint2011arXiv

Application of asymptotic expansions of maximum likelihood estimators errors to gravitational waves from binary mergers: the single interferometer case

In this paper we describe a new methodology to calculate analytically the error for a maximum likelihood estimate (MLE) for physical parameters from Gravitational wave signals. All the existing litterature focuses on the usage of the Cramer Rao Lower bounds (CRLB) as a mean to approximate the errors for large signal to noise ratios. We show here how the variance and the bias of a MLE estimate can be expressed instead in inverse powers of the signal to noise ratios where the first order in the variance expansion is the CRLB. As an application we compute the second order of the variance and bias for MLE of physical parameters from the inspiral phase of binary mergers and for noises of gravitational wave interferometers . We also compare the improved error estimate with existing numerical estimates. The value of the second order of the variance expansions allows to get error predictions closer to what is observed in numerical simulations. It also predicts correctly the necessary SNR to approximate the error with the CRLB and provides new insight on the relationship between waveform properties SNR and estimation errors. For example the timing match filtering becomes optimal only if the SNR is larger than the kurtosis of the gravitational wave spectrum.

preprint2010arXiv

Parameter estimation from Gravitational waves generated by non-spinning binary black holes with laser interferometers: beyond the Fisher information

In this paper we apply to gravitational waves from non-spinning binary systems a recently intro- duced frequentist methodology to calculate analytically the error for a maximum likelihood estimate (MLE) of physical parameters. While existing literature focuses on using the Cramer Rao Lower bound (CRLB) and Monte Carlo simulations, we use a power expansion of the bias and covariance in inverse powers of the signal to noise ratio. The use of higher order derivatives of the likelihood function in the expansions makes the prediction also sensitive to the secondary lobes of the MLE probability distribution. We discuss conditions for validity of the CRLB and predict new features in regions of the parameter space currently not explored. For example, we see how the bias can become the most important contributor to the parameters' errors for high mass systems (200M and above).