BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Michele Vallisneri

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

ResearcherUnclaimed source record
gr-qcastro-ph.GAastro-ph.IMastro-ph.HEastro-ph.COhep-phphysics.data-anApplicationsastro-ph.SRMachine Learningphysics.comp-phphysics.geo-phphysics.hist-phphysics.ins-det

Catalog footprint

What is connected

32works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

32 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GipsyX/RTGx, A New Tool Set for Space Geodetic Operations and Research

GipsyX/RTGx is the Jet Propulsion Laboratory's (JPL) next generation software package for positioning, navigation, timing, and Earth science using measurements from three geodetic techniques: Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS); with Very Long Baseline Interferometry (VLBI) under development. The software facilitates combined estimation of geodetic and geophysical parameters using a Kalman filter approach on real or simulated data in both post-processing and in real-time. The estimated parameters include station coordinates and velocities, satellite orbits and clocks, Earth orientation, ionospheric and tropospheric delays. The software is also capable of full realization of a dynamic terrestrial reference through analysis and combination of time series of ground station coordinates. We present some key aspects of its new architecture, and describe some of its major applications, including Real-time orbit determination and ephemeris predictions in the U.S. Air Force Next Generation GPS Operational Control Segment (OCX), as well as in JPL's Global Differential GPS (GDGPS) System, supporting User Range Error (URE) of $<$ 5 cm RMS; precision post-processing GNSS orbit determination, including JPL's contributions to the International GNSS Service (IGS) with URE in the 2 cm RMS range; Precise point positioning (PPP) with ambiguity resolution, both statically and kinematically, for geodetic applications with 2 mm horizontal, and 6.5 mm vertical repeatability for static positioning; Operational orbit and clock determination for Low Earth Orbiting (LEO) satellites, such as NASA's Gravity Recovery and Climate Experiment (GRACE) mission with GRACE relative clock alignment at the 20 ps level.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Learning Bayesian posteriors with neural networks for gravitational-wave inference

We seek to achieve the Holy Grail of Bayesian inference for gravitational-wave astronomy: using deep-learning techniques to instantly produce the posterior $p(θ|D)$ for the source parameters $θ$, given the detector data $D$. To do so, we train a deep neural network to take as input a signal + noise data set (drawn from the astrophysical source-parameter prior and the sampling distribution of detector noise), and to output a parametrized approximation of the corresponding posterior. We rely on a compact representation of the data based on reduced-order modeling, which we generate efficiently using a separate neural-network waveform interpolant [A. J. K. Chua, C. R. Galley & M. Vallisneri, Phys. Rev. Lett. 122, 211101 (2019)]. Our scheme has broad relevance to gravitational-wave applications such as low-latency parameter estimation and characterizing the science returns of future experiments. Source code and trained networks are available online at https://github.com/vallis/truebayes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On parametric tests of relativity with false degrees of freedom

General relativity can be tested by comparing the binary-inspiral signals found in LIGO--Virgo data against waveform models that are augmented with artificial degrees of freedom. This approach suffers from a number of logical and practical pitfalls. 1) It is difficult to ascribe meaning to the stringency of the resultant constraints. 2) It is doubtful that the Bayesian model comparison of relativity against these artificial models can offer actual validation for the former. 3) It is unknown to what extent these tests might detect alternative theories of gravity for which there are no computed waveforms; conversely, when waveforms are available, tests that employ them will be superior.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Emergence of Gravitational Wave Science: 100 Years of Development of Mathematical Theory, Detectors, Numerical Algorithms, and Data Analysis Tools

On September 14, 2015, the newly upgraded Laser Interferometer Gravitational-wave Observatory (LIGO) recorded a loud gravitational-wave (GW) signal, emitted a billion light-years away by a coalescing binary of two stellar-mass black holes. The detection was announced in February 2016, in time for the hundredth anniversary of Einstein's prediction of GWs within the theory of general relativity (GR). The signal represents the first direct detection of GWs, the first observation of a black-hole binary, and the first test of GR in its strong-field, high-velocity, nonlinear regime. In the remainder of its first observing run, LIGO observed two more signals from black-hole binaries, one moderately loud, another at the boundary of statistical significance. The detections mark the end of a decades-long quest, and the beginning of GW astronomy: finally, we are able to probe the unseen, electromagnetically dark Universe by listening to it. In this article, we present a short historical overview of GW science: this young discipline combines GR, arguably the crowning achievement of classical physics, with record-setting, ultra-low-noise laser interferometry, and with some of the most powerful developments in the theory of differential geometry, partial differential equations, high-performance computation, numerical analysis, signal processing, statistical inference, and data science. Our emphasis is on the synergy between these disciplines, and how mathematics, broadly understood, has historically played, and continues to play, a crucial role in the development of GW science. We focus on black holes, which are very pure mathematical solutions of Einstein's gravitational-field equations that are nevertheless realized in Nature, and that provided the first observed signals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The LIGO Open Science Center

The LIGO Open Science Center (LOSC) fulfills LIGO's commitment to release, archive, and serve LIGO data in a broadly accessible way to the scientific community and to the public, and to provide the information and tools necessary to understand and use the data. In August 2014, the LOSC published the full dataset from Initial LIGO's "S5" run at design sensitivity, the first such large-scale release and a valuable testbed to explore the use of LIGO data by non-LIGO researchers and by the public, and to help teach gravitational-wave data analysis to students across the world. In addition to serving the S5 data, the LOSC web portal (losc.ligo.org) now offers documentation, data-location and data-quality queries, tutorials and example code, and more. We review the mission and plans of the LOSC, focusing on the S5 data release.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $Ω_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $Ω_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $Gμ<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bayesian inference for pulsar timing models

The extremely regular, periodic radio emission from millisecond pulsars makes them useful tools for studying neutron star astrophysics, general relativity, and low-frequency gravitational waves. These studies require that the observed pulse times of arrival be fit to complex timing models that describe numerous effects such as the astrometry of the source, the evolution of the pulsar's spin, the presence of a binary companion, and the propagation of the pulses through the interstellar medium. In this paper, we discuss the benefits of using Bayesian inference to obtain pulsar timing solutions. These benefits include the validation of linearized least-squares model fits when they are correct, and the proper characterization of parameter uncertainties when they are not; the incorporation of prior parameter information and of models of correlated noise; and the Bayesian comparison of alternative timing models. We describe our computational setup, which combines the timing models of Tempo2 with the nested-sampling integrator MultiNest. We compare the timing solutions generated using Bayesian inference and linearized least-squares for three pulsars: B1953+29, J2317+1439, and J1640+2224, which demonstrate a variety of the benefits that we posit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Low-rank approximations for large stationary covariance matrices, as used in the Bayesian and generalized-least-squares analysis of pulsar-timing data

Many data-analysis problems involve large dense matrices that describe the covariance of stationary noise processes; the computational cost of inverting these matrices, or equivalently of solving linear systems that contain them, is often a practical limit for the analysis. We describe two general, practical, and accurate methods to approximate stationary covariance matrices as low-rank matrix products featuring carefully chosen spectral components. These methods can be used to greatly accelerate data-analysis methods in many contexts, such as the Bayesian and generalized-least-squares analysis of pulsar-timing residuals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

New advances in the Gaussian-process approach to pulsar-timing data analysis

In this work we review the application of the theory of Gaussian processes to the modeling of noise in pulsar-timing data analysis, and we derive various useful and optimized representations for the likelihood expressions that are needed in Bayesian inference on pulsar-timing-array datasets. The resulting viewpoint and formalism lead us to two improved parameter-sampling schemes inspired by Gibbs sampling. The new schemes have vastly lower chain autocorrelation lengths than the Markov Chain Monte Carlo methods currently used in pulsar-timing data analysis, potentially speeding up Bayesian inference by orders of magnitude. The new schemes can be used for a full-noise-model analysis of the large datasets assembled by the International Pulsar Timing Array collaboration, which present a serious computational challenge to existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stealth Bias in Gravitational-Wave Parameter Estimation

Inspiraling binaries of compact objects are primary targets for current and future gravitational-wave observatories. Waveforms computed in General Relativity are used to search for these sources, and will probably be used to extract source parameters from detected signals. However, if a different theory of gravity happens to be correct in the strong-field regime, source-parameter estimation may be affected by a fundamental bias: that is, by systematic errors induced due to the use of waveforms derived in the incorrect theory. If the deviations from General Relativity are not large enough to be detectable on their own and yet these systematic errors remain significant (i.e., larger than the statistical uncertainties in parameter estimation), fundamental bias cannot be corrected in a single observation, and becomes stealth bias. In this article we develop a scheme to determine in which cases stealth bias could be present in gravitational-wave astronomy. For a given observation, the answer depends on the detection signal-to-noise ratio and on the strength of the modified-gravity correction. As an example, we study three representative stellar-mass binary systems that will be detectable with second-generation ground-based observatories. We find that significant systematic bias can occur whether or not modified gravity can be positively detected, for correction strengths that are not currently excluded by any other experiment. Thus, stealth bias may be a generic feature of gravitational-wave detections, and it should be considered and characterized, using expanded models such as the parametrized post-Einstein framework, when interpreting the results of parameter-estimation analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing General Relativity with Low-Frequency, Space-Based Gravitational-Wave Detectors

We review the tests of general relativity that will become possible with space-based gravitational-wave detectors operating in the ~0.01mHz - 1Hz low-frequency band. The fundamental aspects of gravitation that can be tested include the presence of additional gravitational fields other than the metric; the number and tensorial nature of gravitational-wave polarization states; the velocity of propagation of gravitational waves; the binding energy and gravitational-wave radiation of binaries, and therefore the time evolution of binary inspirals; the strength and shape of the waves emitted from binary mergers and ringdowns; the true nature of astrophysical black holes; and much more. The strength of this science alone calls for the swift implementation of a space-based detector; the remarkable richness of astrophysics, astronomy, and cosmology in the low-frequency gravitational-wave band make the case even stronger.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Gravitational-Wave Discovery Space of Pulsar Timing Arrays

Recent years have seen a burgeoning interest in using pulsar timing arrays (PTAs) as gravitational-wave (GW) detectors. To date, that interest has focused mainly on three particularly promising source types: supermassive--black-hole binaries, cosmic strings, and the stochastic background from early-Universe phase transitions. In this paper, by contrast, our aim is to investigate the PTA potential for discovering unanticipated sources. We derive significant constraints on the available discovery space based solely on energetic and statistical considerations: we show that a PTA detection of GWs at frequencies above ~3.e-5 Hz would either be an extraordinary coincidence or violate "cherished beliefs;" we show that for PTAs GW memory can be more detectable than direct GWs, and that, as we consider events at ever higher redshift, the memory effect increasingly dominates an event's total signal-to-noise ratio. The paper includes also a simple analysis of the effects of pulsar red noise in PTA searches, and a demonstration that the effects of periodic GWs in the 10^-8 -- 10^-4.5 Hz band would not be degenerate with small errors in standard pulsar parameters (except in a few narrow bands).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational-wave emission from compact Galactic binaries

Compact Galactic binaries where at least one member is a white dwarf or neutron star constitute the majority of individually detectable sources for future low-frequency space-based gravitational-wave (GW) observatories; they also form an unresolved continuum, the dominant Galactic foreground at frequencies below a few mHz. Due to the paucity of electromagnetic observations, the majority of studies of Galactic-binary populations so far have been based on population-synthesis simulations. However, recent surveys have reported several new detections of white-dwarf binaries, providing new constraints for population estimates. In this article, we evaluate the impact of revised local densities of interacting white-dwarf binaries on future GW observations. Specifically: we consider five scenarios that explain these densities with different assumptions on the formation of interacting systems; we simulate corresponding populations of detached and interacting white-dwarf binaries; we estimate the number of individually detectable GW sources and the magnitude of the confusion-noise foreground, as observed by space-based detectors with 5- and 1-Mkm arms. We confirm earlier estimates of thousands of detached-binary detections, but project only few ten to few hundred detections of interacting systems. This reduction is partly due to our assessment of detection prospects, based on the iterative identification and subtraction of bright sources with respect to both instrument and confusion noise. We also confirm earlier estimates for the confusion-noise foreground, except in one scenario that explains smaller local densities of interacting systems with smaller numbers of progenitor detached systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-sky-averaged sensitivity curves for space-based gravitational-wave observatories

(abridged) The signal-to-noise ratio (SNR) is used in gravitational-wave observations as the basic figure of merit for detection confidence and, together with the Fisher matrix, for the amount of physical information that can be extracted from a detected signal. SNRs are usually computed from a sensitivity curve, which describes the gravitational-wave amplitude needed by a monochromatic source of given frequency to achieve a threshold SNR. For interferometric space-based detectors similar to LISA, which are sensitive to long-lived signals and have constantly changing position and orientation, exact SNRs need to be computed on a source-by-source basis. For convenience, most authors prefer to work with sky-averaged sensitivities, accepting inaccurate SNRs for individual sources and giving up control over the statistical distribution of SNRs for source populations. In this paper, we describe a straightforward end-to-end recipe to compute the non-sky-averaged sensitivity of interferometric space-based detectors of any geometry: in effect, we derive error bars for the sky-averaged sensitivity curve, which provide a stringent statistical interpretation for previously unqualified statements about sky-averaged SNRs. As a worked-out example, we consider isotropic and Galactic-disk populations of monochromatic sources, as observed with the "classic LISA" configuration. We confirm that the (standard) inverse-rms average sensitivity for the isotropic population remains the same whether or not the LISA orbits are included in the computation. However, detector motion tightens the distribution of sensitivities, so for 50% of sources the sensitivity is within 30% of its average. For the Galactic-disk population, the average and the distribution of the sensitivity for a moving detector turn out to be similar to the isotropic case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Testing general relativity with gravitational waves: a reality check

The observations of gravitational-wave signals from astrophysical sources such as binary inspirals will be used to test General Relativity for self consistency and against alternative theories of gravity. I describe a simple formula that can be used to characterize the prospects of such tests, by estimating the matched-filtering signal-to-noise ratio required to detect non-General-Relativistic corrections of a given magnitude. The formula is valid for sufficiently strong signals; it requires the computation of a single number, the fitting factor between the General-Relativistic and corrected waveform families; and it can be applied to all tests that embed General Relativity in a larger theory, including tests of individual theories such as Brans-Dicke gravity, as well as the phenomenological schemes that introduce corrections and extra terms in the post-Newtonian phasing expressions of inspiral waveforms. The formula suggests that the volume-limited gravitational-wave searches performed with second-generation ground-based detectors would detect alternative-gravity corrections to General-Relativistic waveforms no smaller than 1-10% (corresponding to fitting factors of 0.9 to 0.99).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Beyond Fisher: exact sampling distributions of the maximum-likelihood estimator in gravitational-wave parameter estimation

Gravitational-wave astronomers often wish to characterize the expected parameter-estimation accuracy of future observations. The Fisher matrix provides a lower bound on the spread of the maximum-likelihood estimator across noise realizations, as well as the leading-order width of the posterior probability, but it is limited to high signal strengths often not realized in practice. By contrast, Monte Carlo Bayesian inference provides the full posterior for any signal strength, but it is too expensive to repeat for a representative set of noises. Here I describe an efficient semianalytical technique to map the exact sampling distribution of the maximum-likelihood estimator across noise realizations, for any signal strength. This technique can be applied to any estimation problem for signals in additive Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cover art: issues in the metric-guided and metric-less placement of random and stochastic template banks

The efficient placement of signal templates in source-parameter space is a crucial requisite for exhaustive matched-filtering searches of modeled gravitational-wave sources. Unfortunately, the current placement algorithms based on regular parameter-space meshes are difficult to generalize beyond simple signal models with few parameters. Various authors have suggested that a general, flexible, yet efficient alternative can be found in randomized placement strategies such as random placement and stochastic placement, which enhances random placement by selectively rejecting templates that are too close to others. In this article we explore several theoretical and practical issues in randomized placement: the size and performance of the resulting template banks; the effects of parameter-space boundaries; the use of quasi-random (self avoiding) number sequences; most important, the implementation of these algorithms in curved signal manifolds with and without the use of a Riemannian signal metric, which may be difficult to obtain. Specifically, we show how the metric can be replaced with a discrete triangulation-based representation of local geometry. We argue that the broad class of randomized placement algorithms offers a promising answer to many search problems, but that the specific choice of a scheme and its implementation details will still need to be fine-tuned separately for each problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Searches for Cosmic-String Gravitational-Wave Bursts in Mock LISA Data

A network of observable, macroscopic cosmic (super-)strings may have formed in the early universe. If so, the cusps that generically develop on cosmic-string loops emit bursts of gravitational radiation that could be detectable by both ground- and space-based gravitational-wave interferometers. Here we report on two versions of a LISA-oriented string-burst search pipeline that we have developed and tested within the context of the Mock LISA Data Challenges. The two versions rely on the publicly available MultiNest and PyMC software packages, respectively. To reduce the effective dimensionality of the search space, our implementations use the F-statistic to analytically maximize over the signal's amplitude and polarization, A and psi, and use the FFT to search quickly over burst arrival times t_C. The standard F-statistic is essentially a frequentist statistic that maximizes the likelihood; we also demonstrate an approximate, Bayesian version of the F-statistic that incorporates realistic priors on A and psi. We calculate how accurately LISA can expect to measure the physical parameters of string-burst sources. To understand LISA's angular resolution for string-burst sources, we draw maps of the waveform fitting factor [maximized over (A psi, t_C)] as a function of sky position; these maps dramatically illustrate why (for LISA) inferring the correct sky location of the emitting string loop will often be practically impossible. We also identify and elucidate several symmetries that are imbedded in this search problem, and we derive the distribution of cut-off frequencies f_max for observable bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Optimal filtering of the LISA data

The LISA time-delay-interferometry responses to a gravitational-wave signal are rewritten in a form that accounts for the motion of the LISA constellation around the Sun; the responses are given in closed analytic forms valid for any frequency in the band accessible to LISA. We then present a complete procedure, based on the principle of maximum likelihood, to search for stellar-mass binary systems in the LISA data. We define the required optimal filters, the amplitude-maximized detection statistic (analogous to the F statistic used in pulsar searches with ground-based interferometers), and discuss the false-alarm and detection probabilities. We test the procedure in numerical simulations of gravitational-wave detection.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A physical template family for gravitational waves from precessing binaries of spinning compact objects: Application to single-spin binaries

The detection of the gravitational waves (GWs) emitted by precessing binaries of spinning compact objects is complicated by the large number of parameters (such as the magnitudes and initial directions of the spins, and the position and orientation of the binary with respect to the detector) that are required to model accurately the precession-induced modulations of the GW signal. In this paper we describe a fast matched-filtering search scheme for precessing binaries, and we adopt the physical template family proposed by Buonanno, Chen, and Vallisneri [Phys.Rev.D 67, 104025 (2003)] for ground-based interferometers. This family provides essentially exact waveforms, written directly in terms of the physical parameters, for binaries with a single significant spin, and for which the observed GW signal is emitted during the phase of adiabatic inspiral (for LIGO-I and VIRGO, this corresponds to a total mass M < 15Msun). We show how the detection statistic can be maximized automatically over all the parameters (including the position and orientation of the binary with respect to the detector), except four (the two masses, the magnitude of the single spin, and the opening angle between the spin and the orbital angular momentum), so the template bank used in the search is only four-dimensional; this technique is relevant also to the searches for GW from extreme--mass-ratio inspirals and supermassive blackhole inspirals to be performed using the space-borne detector LISA. Using the LIGO-I design sensitivity, we compute the detection threshold (~10) required for a false-alarm probability of 10^(-3)/year, and the number of templates (~76,000) required for a minimum match of 0.97, for the mass range (m1,m2)=[7,12]Msun*[1,3]Msun.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A quasi-physical family of gravity-wave templates for precessing binaries of spinning compact objects: Application to double-spin precessing binaries

The gravitational waveforms emitted during the adiabatic inspiral of precessing binaries with two spinning compact bodies of comparable masses, evaluated within the post-Newtonian approximation, can be reproduced rather accurately by the waveforms obtained by setting one of the two spins to zero, at least for the purpose of detection by ground-based gravitational-wave interferometers. Here we propose to use this quasi-physical family of single-spin templates to search for the signals emitted by double-spin precessing binaries, and we find that its signal-matching performance is satisfactory for source masses (m1,m2) in [3,15]Msun x [3,15]Msun. For this mass range, using the LIGO-I design sensitivity, we estimate that the number of templates required to yield a minimum match of 0.97 is ~320,000. We discuss also the accuracy to which the single-spin template family can be used to estimate the parameters of the original double-spin precessing binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Detecting gravitational waves from precessing binaries of spinning compact objects: Adiabatic limit

Black-hole (BH) binaries with single-BH masses m=5--20 Msun, moving on quasicircular orbits, are among the most promising sources for first-generation ground-based gravitational-wave (GW) detectors. Until now, the development of data-analysis techniques to detect GWs from these sources has been focused mostly on nonspinning BHs. The data-analysis problem for the spinning case is complicated by the necessity to model the precession-induced modulations of the GW signal, and by the large number of parameters needed to characterize the system, including the initial directions of the spins, and the position and orientation of the binary with respect to the GW detector. In this paper we consider binaries of maximally spinning BHs, and we work in the adiabatic-inspiral regime to build families of modulated detection templates that (i) are functions of very few physical and phenomenological parameters, (ii) model remarkably well the dynamical and precessional effects on the GW signal, with fitting factors on average >~ 0.97, but (iii) might require increasing the detection thresholds, offsetting at least partially the gains in the fitting factors. Our detection-template families are quite promising also for the case of neutron-star--black-hole binaries, with fitting factors on average ~ 0.93. For these binaries we also suggest (but do not test) a further template family, which would produce essentially exact waveforms written directly in terms of the physical spin parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Detection template families for gravitational waves from the final stages of binary--black-hole inspirals: Nonspinning case

We investigate the problem of detecting gravitational waves from binaries of nonspinning black holes with masses m = 5--20 Msun, moving on quasicircular orbits, which are arguably the most promising sources for first-generation ground-based detectors. We analyze and compare all the currently available post--Newtonian approximations for the relativistic two-body dynamics; for these binaries, different approximations predict different waveforms. We then construct examples of detection template families that embed all the approximate models, and that could be used to detect the true gravitational-wave signal (but not to characterize accurately its physical parameters). We estimate that the fitting factor for our detection families is >~0.95 (corresponding to an event-rate loss <~15%) and we estimate that the discretization of the template family, for ~10^4 templates, increases the loss to <~20%.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Detecting gravitational waves from precessing binaries of spinning compact objects. II. Search implementation for low-mass binaries

Detection template families (DTFs) are built to capture the essential features of true gravitational waveforms using a small set of phenomenological waveform parameters. Buonanno, Chen, and Vallisneri [Phys. Rev. D 67, 104025 (2003)] proposed the ``BCV2'' DTF to perform computationally efficient searches for signals from precessing binaries of compact stellar objects. Here we test the signal-matching performance of the BCV2 DTF for asymmetric--mass-ratio binaries, and specifically for double--black-hole binaries with component masses (m1,m2): (6~12Msun, 1~3Msun), and for black-hole--neutron-star binaries with component masses (m1,m2) = (10Msun, 1.4Msun); we take all black holes to be maximally spinning. We find a satisfactory signal-matching performance, with fitting factors averaging between 0.94 and 0.98. We also scope out the region of BCV2 parameters needed for a template-based search, we evaluate the template match metric, we discuss a template-placement strategy, and we estimate the number of templates needed for searches at the LIGO design sensitivity. In addition, after gaining more insight in the dynamics of spin--orbit precession, we propose a modification of the BCV2 DTF that is parametrized by physical (rather than phenomenological) parameters. We test this modified ``BCV2P'' DTF for the (10Msun, 1.4Msun) black-hole--neutron-star system, finding a signal-matching performance comparable to the BCV2 DTF, and a reliable parameter-estimation capability for target-binary quantities such as the chirp mass and the opening angle (the angle between the black-hole spin and the orbital angular momentum).

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Geometric Time Delay Interferometry

The space-based gravitational-wave observatory LISA, a NASA-ESA mission to be launched after 2012, will achieve its optimal sensitivity using Time Delay Interferometry (TDI), a LISA-specific technique needed to cancel the otherwise overwhelming laser noise in the inter-spacecraft phase measurements. The TDI observables of the Michelson and Sagnac types have been interpreted physically as the virtual measurements of a synthesized interferometer. In this paper, I present Geometric TDI, a new and intuitive approach to extend this interpretation to all TDI observables. Unlike the standard algebraic formalism, Geometric TDI provides a combinatorial algorithm to explore exhaustively the space of second-generation TDI observables (i.e., those that cancel laser noise in LISA-like interferometers with time-dependent armlengths). Using this algorithm, I survey the space of second-generation TDI observables of length (i.e., number of component phase measurements) up to 24, and I identify alternative, improved forms of the standard second-generation TDI observables. The alternative forms have improved high-frequency gravitational-wave sensitivity in realistic noise conditions (because they have fewer nulls in the gravitational-wave and noise response functions), and are less susceptible to instrumental gaps and glitches (because their component phase measurements span shorter time periods).

0 citations0 reviewsNo code linkedOpen access