by Vallisneri, Michele and Galley, Chad R.
24 pages, 7 PDF figures. Mathematica code at http://www.vallis.org/publications/sensitivity

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, and we use it to generate a sampling distribution of sensitivities for a given source population. In effect, we derive error bars for the sky-averaged sensitivity curve. 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 distribution of the sensitivity for a moving detector turn out to be similar to the isotropic case.

by O’Shaughnessy, R. and Healy, J. and London, L. and Meeks, Z. and Shoemaker, D.
Submitted to PRD

The gravitational wave signature emitted from a merging binary depends on the orientation of an observer relative to the binary. Previous studies suggest that emission along the total initial or total final angular momenta leads to both the strongest and simplest signal from a precessing compact binary. In this paper we describe a concrete counterexample: a binary with $latex m_1/m_2=4$, $latex a_1=0.6 \hat{x} = -a_2$, placed in orbit in the x,y plane. We extract the gravitational wave emission along several proposed emission directions, including the initial (Newtonian) orbital angular momentum; the final (~ initial) total angular momentum; and the dominant principal axis of $latex _M$. Using several diagnostics, we show that the suggested preferred directions are not representative. For example, only for a handful of other directions (0.95). We conclude that the information available in just one direction (or mode) does not adequately encode the complexity of orientation-dependent emission for even short signals from merging black hole binaries. Future investigations of precessing, unequal-mass binaries should carefully explore and model their orientation-dependent emission.

by Ohme, Frank
12 pages, 2 figures, 1 table, NRDA2011/Amaldi 9 proceedings

Models of gravitational waveforms from coalescing black-hole binaries play a crucial role in the efforts to detect and interpret those signatures in the data of large-scale interferometers. Here we summarize recent models that combine information both from analytical approximations and numerical relativity. We briefly lay out and compare the strategies employed to build such complete models and we recapitulate the errors associated with various aspects of the modelling process.

by McWilliams, Sean T. and Lang, Ryan N. and Baker, John G. and Thorpe, James Ira
9 pages, 5 figures, submitted to Phys. Rev. D

We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble of systems near the peak of LISA’s sensitivity band, with total rest mass of 2\times10^6 M\odot, a redshift of z = 1, and randomly chosen orientations and sky positions. We find median sky localization errors of approximately \sim3 arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic signals associated with merging massive black holes. We investigate the way in which parameter errors decrease with measurement time, focusing specifically on the additional information provided during the merger-ringdown segment of the signal. We find that this information improves all parameter estimates directly, rather than through diminishing correlations with any subset of well- determined parameters. Although we have employed the baseline LISA design for this study, many of our conclusions regarding the information provided by mergers will be applicable to alternative mission designs as well.

by Petiteau, Antoine and Babak, Stanislav and Sesana, Alberto
12 pages, 8 figures, revised version to address referee’s comments, submitted to ApJ

Gravitational wave signals from coalescing Massive Black Hole (MBH) binaries could be used as standard sirens to measure cosmological parameters. The future space based gravitational wave observatory Laser Interferometer Space Antenna (LISA) will detect up to a hundred of those events, providing very accurate measurements of their luminosity distances. To constrain the cosmological parameters we also need to measure the redshift of the galaxy (or cluster of galaxies) hosting the merger. This requires the identification of a distinctive electromagnetic event associated to the binary coalescence. However, putative electromagnetic signatures may be too weak to be observed. Instead, we study here the possibility of constraining the cosmological parameters by enforcing statistical consistency between all the possible hosts detected within the measurement error box of a few dozen of low redshift (z<3) events. We construct MBH populations using merger tree realizations of the dark matter hierarchy in a LambdaCDM Universe, and we use data from the Millennium simulation to model the galaxy distribution in the LISA error box. We show that, assuming that all the other cosmological parameters are known, the parameter w describing the dark energy equation of state can be constrained to a 4-8% level (2sigma error), competitive with current uncertainties obtained by type Ia supernovae measurements, providing an independent test of our cosmological model.

by Marsh, T. R.
11 pages, 2 figures, in press Classical and Quantum Gravity. Summary of talk given at 8th LISA Symposium, Stanford, July 2010

Close pairs of white dwarfs are potential progenitors of Type~Ia supernovae and they are common, with of order 100 — 300 million in the Galaxy. As such they will be significant, probably dominant, sources of the gravitational waves detectable by LISA. In the context of LISA’s goals for fundamental physics, double white dwarfs are a source of noise, but from an astrophysical perspective, they are of considerable interest in their own right. In this paper I discuss our current knowledge of double white dwarfs and their close relatives (and possible descendants) the AM~CVn stars. LISA will add to our knowledge of these systems by providing the following unique constraints: (i) an almost direct measurement of the Galactic merger rate of DWDs from the detection of short period systems and their period evolution, (ii) an accurate and precise normalisation of binary evolution models at the shortest periods, (iii) a determination of the evolutionary pathways to the formation of AM~CVn stars, (iv) measurements of the influence of tidal coupling in white dwarfs and its significance for stabilising mass transfer, and (v) discovery of numerous examples of eclipsing white dwarfs with the potential for optical follow-up to test models of white dwarfs.

by Lang, Ryan N. and Hughes, Scott A. and Cornish, Neil J.
19 pages, 16 figures, submitted to PRD

The future space-based gravitational wave detector LISA will be able to measure parameters of coalescing massive black hole binaries, often to extremely high accuracy. Previous work has demonstrated that the black hole spins can have a strong impact on the accuracy of parameter measurement. Relativistic spin-induced precession modulates the waveform in a manner which can break degeneracies between parameters, in principle significantly improving how well they are measured. Recent studies have indicated, however, that spin precession may be weak for an important subset of astrophysical binary black holes: those in which the spins are aligned due to interactions with gas. In this paper, we examine how well a binary’s parameters can be measured when its spins are partially aligned and compare results using waveforms that include higher post-Newtonian harmonics to those that are truncated at leading quadrupole order. We find that the weakened precession can substantially degrade parameter estimation, particularly for the “extrinsic” parameters sky position and distance. Absent higher harmonics, LISA typically localizes the sky position of a nearly aligned binary about an order of magnitude less accurately than one for which the spin orientations are random. Our knowledge of a source’s sky position will thus be worst for the gas-rich systems which are most likely to produce electromagnetic counterparts. Fortunately, higher harmonics of the waveform can make up for this degradation. By including harmonics beyond the quadrupole in our waveform model, we find that the accuracy with which most of the binary’s parameters are measured can be substantially improved. In some cases, the improvement is such that they are measured almost as well as when the binary spins are randomly aligned.

by Buonanno, Alessandra and Kidder, Lawrence E. and Mroué, Abdul H. and Pfeiffer, Harald P. and Taracchini, Andrea
17 pages, 11 figures

Building initial conditions for generic binary black-hole evolutions without initial spurious eccentricity remains a challenge for numerical-relativity simulations. This problem can be overcome by applying an eccentricity-removal procedure which consists in evolving the binary for a couple of orbits, estimating the eccentricity, and then correcting the initial conditions. The presence of spins can complicate this procedure. As predicted by post-Newtonian theory, spin-spin interactions and precession prevent the binary from moving along an adiabatic sequence of spherical orbits, inducing oscillations in the radial separation and in the orbital frequency. However, spin-induced oscillations occur at approximately twice the orbital frequency, therefore they can be distinguished from the initial spurious eccentricity, which occurs at approximately the orbital frequency. We develop a new removal procedure based on the derivative of the orbital frequency and find that it is successful in reducing the eccentricity measured in the orbital frequency to less than 0.0001 when moderate spins are present. We test this new procedure using numerical-relativity simulations of binary black holes with mass ratios 1.5 and 3, spin magnitude 0.5 and various spin orientations. The numerical simulations exhibit spin-induced oscillations in the dynamics at approximately twice the orbital frequency. Oscillations of similar frequency are also visible in the gravitational-wave phase and frequency of the dominant mode.

by Amaro-Seoane, Pau and Freitag, Marc Dewi
A small note of 5 pages, submitted to MNRAS letts

Two coalescing black holes (BHs) represent a conspicuous source of gravitational waves (GWs). The merger involves 17 parameters in the general case of Kerr BHs, so that a successful identification and parameter extraction of the information encoded in the waves will provide us with a detailed description of the physics of BHs. A search based on matched-filtering for characterization and parameter extraction requires the development of some $latex 10^{15}$ waveforms. If a third additional BH perturbed the system, the waveforms would not be applicable, and we would need to increase the number of templates required for a valid detection. In this letter, we calculate the probability that more than two BHs interact in the regime of strong relativity in a dense stellar cluster. We determine the physical properties necessary in a stellar system for three black holes to have a close encounter in this regime and also for an existing binary of two BHs to have a strong interaction with a third hole. In both cases the event rate is negligible. While dense stellar systems such as galactic nuclei, globular clusters and nuclear stellar clusters are the breeding grounds for the sources of gravitational waves that ground-based and space-borne detectors like Advanced LIGO and LISA will be exploring, the analysis of the waveforms in full general relativity needs only to evaluate the two-body problem. This reduces the number of templates of waveforms to create by orders of magnitude.

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.