by Cornish, Neil J.
4 pages, no figures

Theoretical studies in gravitational wave astronomy often require the calculation of Fisher Information Matrices and Likelihood functions, which in a direct approach entail the costly step of computing gravitational waveforms. Here I describe an alternative technique that sidesteps the need to compute full waveforms, resulting in significant computational savings. I describe how related techniques can be used to speed up Bayesian inference applied to real gravitational wave data.

by Key, Joey Shapiro and Cornish, Neil J.
11 pages, 19 figures

The Laser Interferometer Space Antenna (LISA) is designed to detect gravitational wave signals from astrophysical sources, including those from coalescing binary systems of compact objects such as black holes. Colliding galaxies have central black holes that sink to the center of the merged galaxy and begin to orbit one another and emit gravitational waves. Some galaxy evolution models predict that the binary black hole system will enter the LISA band with significant orbital eccentricity, while other models suggest that the orbits will already have circularized. Using a full seventeen parameter waveform model that includes the effects of orbital eccentricity, spin precession and higher harmonics, we investigate how well the source parameters can be inferred from simulated LISA data. Defining the reference eccentricity as the value one year before merger, we find that for typical LISA sources, it will be possible to measure the eccentricity to an accuracy of parts in a thousand. The accuracy with which the eccentricity can be measured depends only very weakly on the eccentricity, making it possible to distinguish circular orbits from those with very small eccentricities. LISA measurements of the orbital eccentricity can provide strong constraints on theories of galaxy mergers in the early universe.

by Finn, Lee Samuel and Lommen, Andrea N.
43 pages, 13 figures, submitted to ApJ.

Efforts to detect gravitational waves by timing an array of pulsars have focused traditionally on stationary gravitational waves: e.g., stochastic or periodic signals. Gravitational wave bursts — signals whose duration is much shorter than the observation period — will also arise in the pulsar timing array waveband. Sources that give rise to detectable bursts include the formation or coalescence of supermassive black holes (SMBHs), the periapsis passage of compact objects in highly elliptic or unbound orbits about a SMBH, or cusps on cosmic strings. Here we describe how pulsar timing array data may be analyzed to detect and characterize these bursts. Our analysis addresses, in a mutually consistent manner, a hierarchy of three questions: \emph{i}) What are the odds that a dataset includes the signal from a gravitational wave burst? \emph{ii}) Assuming the presence of a burst, what is the direction to its source? and \emph{iii}) Assuming the burst propagation direction, what is the burst waveform’s time dependence in each of its polarization states? Applying our analysis to synthetic data sets we find that we can \emph{detect} gravitational waves even when the radiation is too weak to either localize the source of infer the waveform, and \emph{detect} and \emph{localize} sources even when the radiation amplitude is too weak to permit the waveform to be determined. While the context of our discussion is gravitational wave detection via pulsar timing arrays, the analysis itself is directly applicable to gravitational wave detection using either ground or space-based detector data.

by Cohen, Michael I. and Cutler, Curt and Vallisneri, Michele
Submitted to CQG; 28 pages, 10 figures; higher-resolution plots available at http://www.vallis.org/publications/cosmicstrings

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.

by Manca, Gian Mario and Vallisneri, Michele
RevTeX4, 21 pages, 9 PDF figures

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.

by Petiteau, Antoine and Shang, Yu and Babak, Stanislav and Feroz, Farhan
25 pages, 9 figures

Coalescing massive Black Hole binaries are the strongest and probably the most important gravitational wave sources in the LISA band. The spin and orbital precessions bring complexity in the waveform and make the likelihood surface richer in structure as compared to the non-spinning case. We introduce an extended multimodal genetic algorithm which utilizes the properties of the signal and the detector response function to analyze the data from the third round of mock LISA data challenge (MLDC 3.2). The performance of this method is comparable, if not better, to already existing algorithms. We have found all five sources present in MLDC 3.2 and recovered the coalescence time, chirp mass, mass ratio and sky location with reasonable accuracy. As for the orbital angular momentum and two spins of the Black Holes, we have found a large number of widely separated modes in the parameter space with similar maximum likelihood values.

by Zanolin, M. and Vitale, S. and Makris, N.

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.

by Zanolin, M. and Vitale, S. and Makris, N.

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.

by Feroz, Farhan and Gair, Jonathan R. and Graff, Philip and Hobson, Michael P and Lasenby, Anthony
21 pages, 11 figures, submitted to CQG

We consider the problem of characterisation of burst sources detected with the Laser Interferometer Space Antenna (LISA) using the multi-modal nested sampling algorithm, MultiNest. We use MultiNest as a tool to search for modelled bursts from cosmic string cusps, and compute the Bayesian evidence associated with the cosmic string model. As an alternative burst model, we consider sine-Gaussian burst signals, and show how the evidence ratio can be used to choose between these two alternatives. We present results from an application of MultiNest to the last round of the Mock LISA Data Challenge, in which we were able to successfully detect and characterise all three of the cosmic string burst sources present in the release data set. We also present results of independent trials and show that MultiNest can detect cosmic string signals with signal-to-noise ratio (SNR) as low as ~7 and sine-Gaussian signals with SNR as low as ~8. In both cases, we show that the threshold at which the sources become detectable coincides with the SNR at which the evidence ratio begins to favour the correct model over the alternative.

by Amaro-Seoane, Pau and Santamaria, Lucia
Submitted to ApJ; abstract abridged, figure 1 has a lower resolution

Even though the existence of intermediate-mass black holes has not yet been corroborated observationally, these objects are of high interest for astrophysics. Our understanding of formation and evolution of supermassive black holes (SMBHs), as well as galaxy evolution modeling and cosmography would dramatically change if an IMBH was observed. The prospect of detection and, possibly, observation and characterization of an IMBH has good chances in lower-frequency gravitational-wave (GW) astrophysics with ground-based detectors such as LIGO, Virgo and the future Einstein Telescope (ET). We present an analysis of the signal of a system of a binary of IMBHs based on a waveform model obtained with numerical relativity simulations coupled with post-Newtonian calculations at the highest available order so as to extend the waveform to lower frequencies. We find that initial LIGO and Virgo are in the position of detecting IMBHs with a signal-to-noise ratio (SNR) of $latex \sim 10$ for systems with total mass between 100 and $latex 500 M_{\odot}$ situated at a distance of 100 Mpc. Nevertheless, the event rate is too low and the possibility that these signals are mistaken with a glitch is, unfortunately, non-negligible. When going to second- and third-generation detectors, such as Advanced LIGO or the proposed ET, the event rate becomes much more promising (tens per year for the first and thousands per year for the latter) and the SNR at 100 Mpc is as high as 100 — 1000 and 1000 — $latex 10^{5}$ respectively. The prospects for IMBH detection and characterization with ground-based GW observatories would not only provide us with a robust test of general relativity, but would also corroborate the existence of these systems. Such detections would be a probe to the stellar environments of IMBHs and their formation.

by Yunes, Nicolas and Pretorius, Frans
25 pages, submitted to Phys. Rev. D

We consider the concept of fundamental bias in gravitational wave astrophysics as the assumption that general relativity is the correct theory of gravity during the entire wave-generation and propagation regime. Such an assumption is valid in the weak-field, as verified by precision experiments and observations, but it need not hold in the dynamical strong-field regime where tests are lacking. Fundamental bias can cause systematic errors in the detection and parameter estimation of signals, which can lead to a mischaracterization of the universe through incorrect inferences about source event rates and populations. We propose a remedy through the introduction of the parameterized post-Einsteinian framework, which consists of the enhancement of waveform templates via the inclusion of post-Einsteinian parameters. These parameters would ostensibly be designed to interpolate between templates constructed in general relativity and well-motivated alternative theories of gravity, and also include extrapolations that follow sound theoretical principles, such as consistency with conservation laws and symmetries. As an example, we construct parameterized post-Einsteinian templates for the binary coalescence of equal-mass, non-spinning compact objects in a quasi-circular inspiral. The parametrized post-Einsteinian framework should allow matched filtered data to select a specific set of post-Einsteinian parameters without a priori assuming the validity of the former, thus either verifying general relativity or pointing to possible dynamical strong-field deviations.

by Whelan, John T. and Prix, Reinhard and Khurana, Deepak
26 pages, 11 figures

We describe the application of an F-statistic search for continuous gravitational waves to the search for galactic white-dwarf binaries in the Second Mock LISA Data Challenge. The search method employs a hierarchical template-grid based exploration of the parameter space, using a coincidence step to distinguish between primary (”true”) and secondary maxima, followed by a final (multi-TDI) “zoom” stage to provide an accurate parameter estimation of the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters comparable to statistical expectations, although there is still some systematic excess above expected statistical errors due to Gaussian noise. An experimental iterative pipeline with seven rounds of subtraction and re-analysis allows us to increase the number of signals recovered, up to a total of 3419 with 29 false alarms.

by Amaro-Seoane, Pau and Eichhorn, Christoph and Porter, Ed and Spurzem, Rainer
21 pages, 12 figs, to appear in MNRAS, abstract abridged for arxiv

The dynamical evolution of binaries of intermediate-massive black holes (IMBHs, massive black holes with a mass ranging between $latex 10^2$ and $latex 10^4 M_{\odot}$) in stellar clusters has recently received an increasing amount of attention. This is at least partially due to the fact that if the binary is hard enough to evolve to the phase at which it will start emitting gravitational waves (GWs) efficiently, there is a good probability that it will be detectable by future space-borne detectors like LISA. We study this evolution in the presence of rotation in the cluster. The eccentricity is strongly connected to the initial IMBHs velocities, and values of $latex \sim 0.7$ up to 0.9 are reached for low initial velocities, while almost circular orbits result if the initial velocities are increased. A Monte Carlo study indicates that these sources will be detectable by a detector such as LISA with median signal to noise ratios of between 10 and 20 over a three year period, although some events had signal to noise ratios of 300 or greater. Furthermore, one should also be able to estimate the chirp-mass with median fractional errors of $latex 10^{-4}$, reduced mass on the order of $latex 10^{-3}$ and luminosity distance on the order of $latex 10^{-1}$. Finally, these sources will have a median angular resolution in the LISA detector of about 3 square degrees, putting events firmly in the field of view of future electromagnetic detectors such as LSST.

by Asada, Hideki
21 pages, 3 figures, 1 table

Continuing work initiated in an earlier publication [Torigoe et al. Phys. Rev. Lett. {\bf 102}, 251101 (2009)], gravitational wave forms for a three-body system in Lagrange’s orbit are considered especially in an analytic method. First, we derive an expression of the three-body wave forms at the mass quadrupole, octupole and current quadrupole orders. By using the expressions, we solve a gravitational-wave {\it inverse} problem of determining the source parameters to this particular configuration (three masses, a distance of the source to an observer, and the orbital inclination angle to the line of sight) through observations of the gravitational wave forms alone. For this purpose, the chirp mass to a three-body system in the particular configuration is expressed in terms of only the mass ratios by deleting initial angle positions. We discuss also whether and how a binary source can be distinguished from a three-body system in Lagrange’s orbit or others.

by Klein, Antoine and Jetzer, Philippe and Sereno, Mauro
20 pages, 20 figures. Submitted to Phys. Rev. D

Gravitational waves emitted by binary systems in the inspiral phase carry a complicated structure, consisting in a superposition of different harmonics of the orbital frequency, the amplitude of each of them taking the form of a Post-Newtonian series. In addition to that, spinning binaries experience couplings which induce a precession of the orbital angular momentum and of the individual spins. So far, all studies of the measurement accuracy of gravitational wave experiments for comparable-mass binary systems have considered either spinless binaries, or spinning binaries without subdominant harmonics in the waveform as well as no amplitude modulations. In this paper, we consider supermassive black hole binaries as expected to be observed with the planned space-based interferometer LISA, and study the measurement accuracy for several astrophysically interesting parameters obtainable taking into account the full 2PN waveform for spinning bodies, as well as spin-precession effects. We find that for binaries with a total mass in the range 10^5 M_Sun < M 10 for M < 10^7 M_Sun, 1.5 – 5 times higher than with the restricted waveform. We computed that the full waveform allows to use supermassive black hole binaries as standard sirens up to a redshift of z = 1.6, about 0.4 larger than what previous studies allowed.

by Prix, Reinhard and Krishnan, Badri
12 pages, presented at GWDAW13, to appear in CQG

We investigate the Bayesian framework for detection of continuous gravitational waves (GWs) in the context of targeted searches, where the phase evolution of the GW signal is assumed to be known, while the four amplitude parameters are unknown. We show that the orthodox maximum-likelihood statistic (known as F-statistic) can be rediscovered as a Bayes factor with an unphysical prior in amplitude parameter space. We introduce an alternative detection statistic (”B-statistic”) using the Bayes factor with a more natural amplitude prior, namely an isotropic probability distribution for the orientation of GW sources. Monte-Carlo simulations of targeted searches show that the resulting Bayesian B-statistic is more powerful in the Neyman-Pearson sense (i.e. has a higher expected detection probability at equal false-alarm probability) than the frequentist F-statistic.

by De Laurentis, M. and Capozziello, S.
6 pages, 3 figures; Multifrequency Behaviour of High-Energy Cosmic Sources, Vulcano Workshop 2009

The emission of gravitational waves from a system of massive objects interacting on elliptical, hyperbolic and parabolic orbits is studied in the quadrupole approximation. Analytical expressions are then derived for the gravitational wave luminosity, the total energy output and gravitational radiation amplitude. A crude estimate of the expected number of events towards peculiar targets (i.e. globular clusters) is also given. In particular, the rate of events per year is obtained for the dense stellar cluster at the Galactic Center.

by Yagi, Kent and Tanaka, Takahiro
37 pages, 16 figures

We calculate how strong one can put constraints on the alternative theories of gravities such as Brans-Dicke and massive graviton theories with LISA. We consider the inspiral gravitational waves from NS/IMBH binaries in Brans-Dicke theory and SMBH/BH binaries in massive graviton theories. We use the 2PN waveforms including spins. We also take both precession and small eccentricity of the orbit into account. We neglect the spin of one of the binary object so that we can apply the so-called \textit{simple precession}. We perform the Monte Carlo simulations of $latex 10^4$ binaries, whose parameters include the Brans-Dicke parameter $latex \omega_{\mathrm{BD}}$ and the graviton Compton length $latex \lambda_g$. We find that including both the spin-spin coupling $latex \sigma$ and the small eccentricity into the binary parameters reduces the determination accuracy by an order of magnitude for the Brans-Dicke case, whilst it has less influence on massive graviton theories. On the other hand, including precession enhances the constraint on $latex \omega_{\mathrm{BD}}$ only 20% but it increases the constraint on $latex \lambda_g$ by several factors. For $latex (1.4+1000)M_{\odot}$ NS/BH binaries of SNR=10, one can put $latex \omega_{\mathrm{BD}}>7040$, whilst for $latex (10^7+10^6)M_{\odot}$ BH/BH binaries at 3Gpc, one can put $latex \lambda_g>4.24\times10^{21}$cm, on average. This is four orders of magnitude stronger than the one obtained from the solar system experiment. From these results, it is understood that the effects of precession and eccentricity cannot be neglected in the parameter estimation analysis.

by Berti, Emanuele and Cardoso, Vitor and Starinets, Andrei O.
112 pages, 20 figures, 4 tables. Invited Topical Review for CQG. Comments and suggestions are welcome. QNM data and Mathematica notebooks are available at http://www.phy.olemiss.edu/~berti/qnms.html and http://gamow.ist.utl.pt/~vitor/ringdown

Quasinormal modes are eigenmodes of dissipative systems. Perturbations of classical gravitational backgrounds involving black holes or branes naturally lead to quasinormal modes. The analysis and classification of the quasinormal spectra requires solving non-Hermitian eigenvalue problems for the associated linear differential equations. Within the recently developed gauge-gravity duality, these modes serve as an important tool for determining the near-equilibrium properties of strongly coupled quantum field theories, in particular their transport coefficients, such as viscosity, conductivity and diffusion constants. In astrophysics, the detection of quasinormal modes in gravitational wave experiments would allow precise measurements of the mass and spin of black holes as well as new tests of general relativity. This review is meant as an introduction to the subject, with a focus on the recent developments in the field.

by Laguna, Pablo and Larson, Shane L. and Spergel, David and Yunes, Nicolas
4 pages, 1 figure

Gravitational waves are messengers carrying valuable information about their sources. For sources at cosmological distances, the waves will contain also the imprint left by the intervening matter. The situation is in close analogy with cosmic microwave photons, for which the large-scale structures the photons traverse contribute to the observed temperature anisotropies, in a process known as the integrated Sachs-Wolfe effect. We derive the gravitational wave counterpart of this effect for waves propagating on a Friedman-Robertson-Walker background with scalar perturbations. We find that the phase, frequency and amplitude of the gravitational waves experience Sachs-Wolfe type integrated effects, this in addition to the magnification effects on the amplitude from gravitational lensing. We show that for supermassive black hole binaries, the integrated effects could account for measurable changes on the frequency, chirp mass and luminosity distance of the binary, thus unveiling the presence of inhomogeneities, and potentially dark energy, in the Universe.

by Gair, Jonathan R. and Porter, Edward K.
submitted to Classical & Quantum Gravity. 19 pages, 4 figures

We describe a hybrid evolutionary algorithm that can simultaneously search for multiple supermassive black hole binary (SMBHB) inspirals in LISA data. The algorithm mixes evolutionary computation, Metropolis-Hastings methods and Nested Sampling. The inspiral of SMBHBs presents an interesting problem for gravitational wave data analysis since, due to the LISA response function, the sources have a bi-modal sky solution. We show here that it is possible not only to detect multiple SMBHBs in the data stream, but also to investigate simultaneously all the various modes of the global solution. In all cases, the algorithm returns parameter determinations within $latex 5\sigma$ (as estimated from the Fisher Matrix) of the true answer, for both the actual and antipodal sky solutions.

by Hopman, Clovis
Accepted to CQG, special LISA issue

We study the rate at which stars spiral into a massive black hole (MBH) due to the emission of gravitational waves (GWs), as a function of the mass M of the MBH. In the context of our model, it is shown analytically that the rate approximately depends on the MBH mass as $latex M^{-1/4}$. Numerical simulations confirm this result, and show that for all MBH masses, the event rate is highest for stellar black holes, followed by white dwarfs, and lowest for neutron stars. The Laser Interferometer Space Antenna (LISA) is expected to see hundreds of these extreme mass ratio inspirals per year. Since the event rate derived here formally diverges as M->0, the model presented here cannot hold for MBHs of masses that are too low, and we discuss what the limitations of the model are.