Source author record

Walter Del Pozzo

Walter Del Pozzo 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

24works
8topics
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

24 published item(s)

preprint2026arXiv

The Sequential Monte Carlo goes NUTS: Boosting Gravitational-Wave Inference

Sequential Monte Carlo (SMC) methods have recently been applied to gravitational-wave inference as a powerful alternative to standard sampling techniques, such as Nested Sampling. At the same time, gradient-based Markov Chain Monte Carlo algorithms, most notably the No-U-Turn Sampler (NUTS), provide an efficient way to explore high-dimensional parameter spaces. In this work we present SHARPy, a Bayesian inference framework that combines the parallelism and evidence-estimation capabilities of SMC with the state-of-the-art sampling performance of NUTS. Moreover, SHARPy exploits the local geometric structure of the posterior to further improve efficiency. Built on JAX and accelerated on GPUs, SHARPy performs gravitational-wave inference on binary black-hole events in around ten minutes, yielding posterior samples and Bayesian evidence estimates that are consistent with those obtained through Nested Sampling. This work sets a new milestone in GW inference with likelihood-based methods and paves the way for model comparison tasks to be accomplished in minutes.

preprint2022arXiv

Constraints on Kerr-Newman black holes from merger-ringdown gravitational-wave observations

We construct a template to model the post-merger phase of a binary black hole coalescence in the presence of a remnant $U(1)$ charge. We include the quasi-normal modes typically dominant during a binary black hole coalescence, $(\ell,m,n) = \{(2,2,0), (2,2,1)\}$ and also present analytical fits for the quasinormal mode frequencies of a Kerr-Newman black hole in terms of its spin and charge, here also including the $(3,3,0)$ mode. Aside from astrophysical electric charge, our template can accommodate extensions of the Standard Model, such as a dark photon. Applying the model to LIGO-Virgo detections, we find that we are unable to distinguish between the charged and uncharged hypotheses from a purely post-merger analysis of the current events. However, restricting the mass and spin to values compatible with the analysis of the full signal, we obtain a 90th percentile bound $\bar{q} < 0.33$ on the black hole charge-to-mass ratio, for the most favorable case of GW150914. Under similar assumptions, by simulating a typical loud signal observed by the LIGO-Virgo network at its design sensitivity, we assess that this model can provide a robust measurement of the charge-to-mass ratio only for values $\bar{q} \gtrsim 0.5$; here we also assume that the mode amplitudes are similar to the uncharged case in creating our simulated signal. Lower values, down to $\bar{q} \sim 0.3$, could instead be detected when evaluating the consistency of the pre-merger and post-merger emission.

preprint2022arXiv

Maximally informed Bayesian modelling of disc galaxies

Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement.

preprint2022arXiv

Multiband gravitational wave cosmology with stellar origin black hole binaries

Massive stellar origin black hole binaries (SBHBs), originating from stars above the pair-instability mass gap, are primary candidates for multiband gravitational wave (GW) observations. Here we study the possibility to use them as effective dark standard sirens to constrain cosmological parameters. The long lasting inspiral signal emitted by these systems is accessible by the future $Laser \; Interferometer \; Space \; Antenna$ (LISA), while the late inspiral and merger are eventually detected by third generation ground-based telescopes such as the $Einstein \; Telescope$ (ET). The direct measurement of the luminosity distance and the sky position to the source, together with the inhomogeneous redshift distribution of possible host galaxies, allow us to infer cosmological parameters by probabilistic means. The efficiency of this statistical method relies in high parameter estimation performances. We show that this multiband approach allows a precise determination of the Hubble constant H$_0$ with just ${\cal O}(10)$ detected sources. For selected SBHB population models, assuming $4$ ($10$) years of LISA observations, we find that H$_0$ is typically determined at $\sim 2\%$ ($\sim 1.5\%$), whereas $Ω_m$ is only mildly constrained with a typical precision of $30\%$ ($20\%$). We discuss the origin of some outliers in our final estimates and we comment on ways to reduce their presence.

preprint2022arXiv

Rapid localization of gravitational wave hosts with FIGARO

The copious scientific literature produced after the detection of GW170817 electromagnetic counterpart demonstrated the importance of a prompt and accurate localization of the gravitational wave within the co-moving volume. In this letter, we present FIGARO, a ready to use and publicly available software that relies on Bayesian non-parametrics. FIGARO is designed to run in parallel with parameter estimation algorithms to provide updated three-dimensional volume localization information. Differently from any existing algorithms, the analytical nature of the FIGARO reconstruction allows a ranking of the entries of galaxy catalogues by their probability of being the host of a gravitational wave event, hence providing an additional tool for a prompt electromagnetic follow up of gravitational waves. We illustrate the features of FIGARO on binary black holes as well as on GW170817. Finally, we demonstrate the robustness of FIGARO by producing so-called pp-plots and we present a method based on information entropy to assess when, during the parameter estimation run, it is reasonable to begin releasing skymaps.

preprint2021arXiv

Eigenvalue repulsions in the quasinormal spectra of the Kerr-Newman black hole

We study the gravito-electromagnetic perturbations of the Kerr-Newman (KN) black hole metric and identify the two $-$ photon sphere and near-horizon $-$ families of quasinormal modes (QNMs) of the KN black hole, computing the frequency spectra (for all the KN parameter space) of the modes with the slowest decay rate. We uncover a novel phenomenon for QNMs that is unique to the KN system, namely eigenvalue repulsion between QNM families. Such a feature is common in solid state physics where \eg it is responsible for energy bands/gaps in the spectra of electrons moving in certain Schrödinger potentials. Exploiting the enhanced symmetries of the near-horizon limit of the near-extremal KN geometry we also develop a matching asymptotic expansion that allows us to solve the perturbation problem using separation of variables and provides an excellent approximation to the KN QNM spectra near extremality. The KN QNM spectra here derived are required not only to account for the gravitational emission in astrophysical environments, such as the ones probed by LIGO, Virgo and LISA, but also allow to extract observational implications on several new physics scenarios, such as mini-charged dark-matter or certain modified theories of gravity, degenerate with the KN solution at the scales of binary mergers.

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

Machine Learning Gravitational Waves from Binary Black Hole Mergers

We apply machine learning methods to build a time-domain model for gravitational waveforms from binary black hole mergers, called mlgw. The dimensionality of the problem is handled by representing the waveform's amplitude and phase using a principal component analysis. We train mlgw on about $\mathcal{O}(10^3)$ TEOBResumS and SEOBNRv4 effective-one-body waveforms with mass ratios $q\in[1,20]$ and aligned dimensionless spins $s\in[-0.80,0.95]$. The resulting models are faithful to the training sets at the ${\sim}10^{-3}$ level (averaged on the parameter space). The speed up for a single waveform generation is a factor 10 to 50 (depending on the binary mass and initial frequency) for TEOBResumS and approximately an order of magnitude more for SEOBNRv4. Furthermore, mlgw provides a closed form expression for the waveform and its gradient with respect to the orbital parameters; such an information might be useful for future improvements in GW data analysis. As demonstration of the capabilities of mlgw to perform a full parameter estimation, we re-analyze the public data from the first GW transient catalog (GWTC-1). We find broadly consistent results with previous analyses at a fraction of the cost, although the analysis with spin aligned waveforms gives systematic larger values of the effective spins with respect to previous analyses with precessing waveforms. Since the generation time does not depend on the length of the signal, our model is particularly suitable for the analysis of the long signals that are expected to be detected by third-generation detectors. Future applications include the analysis of waveform systematics and model selection in parameter estimation.

preprint2020arXiv

Massive black hole binary systems and the NANOGrav 12.5 year results

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported evidence for the presence of a common stochastic signal across their array of pulsars. The origin of this signal is still unclear. One of the possibilities is that it is due to a stochastic gravitational wave background (SGWB) in the $\sim 1-10\,{\rm nHz}$ frequency region. Taking the NANOGrav observational result at face value, we show that this signal would be fully consistent with a SGWB produced by an unresolved population of in-spiralling massive black hole binaries (MBHBs) predicted by current theoretical models. Considering an astrophysically agnostic model we find that the MBHB merger rate is loosely constrained to the range $10^{-11} - 2$ $\mathrm{Mpc}^{-3}\,\mathrm{Gyr}^{-1}$. Including additional constraints from galaxy pairing fractions and MBH-bulge scaling relations, we find that the MBHB merger rate is $10^{-5} - 5\times10^{-4}$ $\mathrm{Mpc}^{-3}\,\mathrm{Gyr}^{-1}$, the MBHB merger time-scale is $\le 3\,\mathrm{Gyr}$ and the norm of the $M_\mathrm{BH}-M_\mathrm{bulge}$ relation $\ge 1.2\times 10^{8}\,M_\odot$ (all intervals quoted at 90\% confidence). Regardless of the astrophysical details of MBHB assembly, this result would imply that a sufficiently large population of massive black holes pair up, form binaries and merge within a Hubble time.

preprint2016arXiv

On tests of general relativity with binary radio pulsars

The timing of radio pulsars in binary systems provides a superb testing ground of general relativity. Here we propose a Bayesian approach to carry out these tests, and a relevant efficient numerical implementation, that has several conceptual and practical advantages with respect to traditional methods based on least-square-fits that have been used so far: (i) it accounts for the actual structure of the likelihood function - and it is not predicated on the Laplace approximation which is implicitly built in least-square fits that can potentially bias the inference - (ii) it provides the ratio of the evidences of any two models under consideration as the statistical quantity to compare different theories, and (iii) it allows us to put joint constraints from the monitoring of multiple systems, that can be expressed in terms of ratio of evidences or probability intervals of global (thus not system-dependent) parameters of the theory, if any exists. Our proposed approach optimally exploits the progress in timing of radio pulsars and the increase in the number of observed systems. We demonstrate the power of this framework using simulated data sets that are representative of current observations.

preprint2016arXiv

Testing general relativity using golden black-hole binaries

The coalescences of stellar-mass black-hole binaries through their inspiral, merger, and ringdown are among the most promising sources for ground-based gravitational-wave (GW) detectors. If a GW signal is observed with sufficient signal-to-noise ratio, the masses and spins of the black holes can be estimated from just the inspiral part of the signal. Using these estimates of the initial parameters of the binary, the mass and spin of the final black hole can be uniquely predicted making use of general-relativistic numerical simulations. In addition, the mass and spin of the final black hole can be independently estimated from the merger--ringdown part of the signal. If the binary black hole dynamics is correctly described by general relativity (GR), these independent estimates have to be consistent with each other. We present a Bayesian implementation of such a test of general relativity, which allows us to combine the constraints from multiple observations. Using kludge modified GR waveforms, we demonstrate that this test can detect sufficiently large deviations from GR, and outline the expected constraints from upcoming GW observations using the second-generation of ground-based GW detectors.

preprint2015arXiv

Astrophysical constraints on massive black hole binary evolution from Pulsar Timing Arrays

We consider the information that can be derived about massive black-hole binary populations and their formation history solely from current and possible future pulsar timing array (PTA) results. We use models of the stochastic gravitational-wave background from circular massive black hole binaries with chirp mass in the range $10^6 - 10^{11} M_\odot$ evolving solely due to radiation reaction. Our parameterised models for the black hole merger history make only weak assumptions about the properties of the black holes merging over cosmic time. We show that current PTA results place an upper limit on the black hole merger density which does not depend on the choice of a particular merger history model, however they provide no information about the redshift or mass distribution. We show that even in the case of a detection resulting from a factor of 10 increase in amplitude sensitivity, PTAs will only put weak constraints on the source merger density as a function of mass, and will not provide any additional information on the redshift distribution. Without additional assumptions or information from other observations, a detection cannot meaningfully bound the massive black hole merger rate above zero or any particular mass.

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

Estimating parameters of binary black holes from gravitational-wave observations of their inspiral, merger and ringdown

We characterize the expected statistical errors with which the parameters of black-hole binaries can be measured from gravitational-wave (GW) observations of their inspiral, merger and ringdown by a network of second-generation ground-based GW observatories. We simulate a population of black-hole binaries with uniform distribution of component masses in the interval $(3,80)~M_\odot$, distributed uniformly in comoving volume, with isotropic orientations. From signals producing signal-to-noise ratio $\geq 5$ in at least two detectors, we estimate the posterior distributions of the binary parameters using the Bayesian parameter estimation code LALInference. The GW signals will be redshifted due to the cosmological expansion and we measure only the "redshifted" masses. By assuming a cosmology, it is possible to estimate the gravitational masses by inferring the redshift from the measured posterior of the luminosity distance. We find that the measurement of the gravitational masses will be in general dominated by the error in measuring the luminosity distance. In spite of this, the component masses of more than $50\%$ of the population can be measured with accuracy better than $\sim 25\%$ using the Advanced LIGO-Virgo network. Additionally, the mass of the final black hole can be measured with median accuracy $\sim 18\%$. Spin of the final black hole can be measured with median accuracy $\sim 5\% ~(17\%)$ for binaries with non-spinning (aligned-spin) black holes. Additional detectors in Japan and India significantly improve the accuracy of sky localization, and moderately improve the estimation of luminosity distance, and hence, that of all mass parameters. We discuss the implication of these results on the observational evidence of intermediate-mass black holes and the estimation of cosmological parameters using GW observations.

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

Testing general relativity with compact coalescing binaries: comparing exact and predictive methods to compute the Bayes factor

The second generation of gravitational-wave detectors is scheduled to start operations in 2015. Gravitational-wave signatures of compact binary coalescences could be used to accurately test the strong-field dynamical predictions of general relativity. Computationally expensive data analysis pipelines, including TIGER, have been developed to carry out such tests. As a means to cheaply assess whether a particular deviation from general relativity can be detected, Cornish et al. and Vallisneri recently proposed an approximate scheme to compute the Bayes factor between a general-relativity gravitational-wave model and a model representing a class of alternative theories of gravity parametrised by one additional parameter. This approximate scheme is based on only two easy-to-compute quantities: the signal-to-noise ratio of the signal and the fitting factor between the signal and the manifold of possible waveforms within general relativity. In this work, we compare the prediction from the approximate formula against an exact numerical calculation of the Bayes factor using the lalinference library. We find that, using frequency-domain waveforms, the approximate scheme predicts exact results with good accuracy, providing the correct scaling with the signal-to-noise ratio at a fitting factor value of 0.992 and the correct scaling with the fitting factor at a signal-to-noise ratio of 20, down to a fitting factor of $\sim$ 0.9. We extend the framework for the approximate calculation of the Bayes factor which significantly increases its range of validity, at least to fitting factors of $\sim$ 0.7 or higher.

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.

preprint2012arXiv

A Mock Data Challenge for the Einstein Gravitational-Wave Telescope

Einstein Telescope (ET) is conceived to be a third generation gravitational-wave observatory. Its amplitude sensitivity would be a factor ten better than advanced LIGO and Virgo and it could also extend the low-frequency sensitivity down to 1--3 Hz, compared to the 10--20 Hz of advanced detectors. Such an observatory will have the potential to observe a variety of different GW sources, including compact binary systems at cosmological distances. ET's expected reach for binary neutron star (BNS) coalescences is out to redshift $z\simeq 2$ and the rate of detectable BNS coalescences could be as high as one every few tens or hundreds of seconds, each lasting up to several days. %in the sensitive frequency band of ET. With such a signal-rich environment, a key question in data analysis is whether overlapping signals can be discriminated. In this paper we simulate the GW signals from a cosmological population of BNS and ask the following questions: Does this population create a confusion background that limits ET's ability to detect foreground sources? How efficient are current algorithms in discriminating overlapping BNS signals? Is it possible to discern the presence of a population of signals in the data by cross-correlating data from different detectors in the ET observatory? We find that algorithms currently used to analyze LIGO and Virgo data are already powerful enough to detect the sources expected in ET, but new algorithms are required to fully exploit ET data.

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.

preprint2012arXiv

Inference of the cosmological parameters from gravitational waves: application to second generation interferometers

The advanced world-wide network of gravitational waves (GW) observatories is scheduled to begin operations within the current decade. Thanks to their improved sensitivity, they promise to yield a number of detections and thus to open a new observational windows for astronomy and astrophysics. Among the scientific goals that should be achieved, there is the independent measurement of the value of the cosmological parameters, hence an independent test of the current cosmological paradigm. Due to the importance of such task, a number of studies have evaluated the capabilities of GW telescopes in this respect. However, since GW do not yield information about the source redshift, different groups have made different assumptions regarding the means through which the GW redshift can be obtained. These different assumptions imply also different methodologies to solve this inference problem. This work presents a formalism based on Bayesian inference developed to facilitate the inclusion of all assumptions and prior information about a GW source within a single data analysis framework. This approach guarantees the minimisation of information loss and the possibility of including naturally event-specific knowledge (such as the sky position for a Gamma Ray Burst - GW coincident observation) in the analysis. The workings of the method are applied to a specific example, loosely designed along the lines of the method proposed by Schutz in 1986, in which one uses information from wide-field galaxy surveys as prior information for the location of a GW source. I show that combining the results from few tens of observations from a network of advanced interferometers will constrain the Hubble constant $H_0$ to an accuracy of $\sim 4 - 5$% at 95% confidence.

preprint2011arXiv

Testing General Relativity using Bayesian model selection: Applications to observations of gravitational waves from compact binary systems

Second generation interferometric gravitational wave detectors, such as Advanced LIGO and Advanced Virgo, are expected to begin operation by 2015. Such instruments plan to reach sensitivities that will offer the unique possibility to test General Relativity in the dynamical, strong field regime and investigate departures from its predictions, in particular using the signal from coalescing binary systems. We introduce a statistical framework based on Bayesian model selection in which the Bayes factor between two competing hypotheses measures which theory is favored by the data. Probability density functions of the model parameters are then used to quantify the inference on individual parameters. We also develop a method to combine the information coming from multiple independent observations of gravitational waves, and show how much stronger inference could be. As an introduction and illustration of this framework - and a practical numerical implementation through the Monte Carlo integration technique of nested sampling - we apply it to gravitational waves from the inspiral phase of coalescing binary systems as predicted by General Relativity and a very simple alternative theory in which the graviton has a non-zero mass. This method can trivially (and should) be extended to more realistic and physically motivated theories.