Detecting extreme mass ratio inspirals with LISA using time-frequency methods II: search characterization

The inspirals of stellar-mass compact objects into supermassive black holes constitute some of the most important sources for LISA. Detection of these sources using fully coherent matched filtering is computationally intractable, so alternative approaches are required. In a previous paper (Wen and Gair 2005, gr-qc/0502100), we outlined a detection method based on looking for excess power in a time-frequency spectrogram of the LISA data. The performance of the algorithm was assessed using a single `typical' trial waveform and approximations to the noise statistics. In this paper we present results of Monte Carlo simulations of the search noise statistics and examine its performance in detecting a wider range of trial waveforms. We show that typical extreme mass ratio inspirals (EMRIs) can be detected at distances of up to 1--3 Gpc, depending on the source parameters. We also discuss some remaining issues with the technique and possible ways in which the algorithm can be improved.Comment: 15 pages, 9 figures, to appear in proceedings of GWDAW 9, Annecy, France, December 200

Detecting extreme mass ratio inspiral events in LISA data using the Hierarchical Algorithm for Clusters and Ridges (HACR)

One of the most exciting prospects for the Laser Interferometer Space Antenna (LISA) is the detection of gravitational waves from the inspirals of stellar-mass compact objects into supermassive black holes. Detection of these sources is an extremely challenging computational problem due to the large parameter space and low amplitude of the signals. However, recent work has suggested that the nearest extreme mass ratio inspiral (EMRI) events will be sufficiently loud that they might be detected using computationally cheap, template-free techniques, such as a time-frequency analysis. In this paper, we examine a particular time-frequency algorithm, the Hierarchical Algorithm for Clusters and Ridges (HACR). This algorithm searches for clusters in a power map and uses the properties of those clusters to identify signals in the data. We find that HACR applied to the raw spectrogram performs poorly, but when the data is binned during the construction of the spectrogram, the algorithm can detect typical EMRI events at distances of up to $\sim2.6$Gpc. This is a little further than the simple Excess Power method that has been considered previously. We discuss the HACR algorithm, including tuning for single and multiple sources, and illustrate its performance for detection of typical EMRI events, and other likely LISA sources, such as white dwarf binaries and supermassive black hole mergers. We also discuss how HACR cluster properties could be used for parameter extraction.Comment: 21 pages, 11 figures, submitted to Class. Quantum Gravity. Modified and shortened in light of referee's comments. Updated results consider tuning over all three HACR thresholds, and show 10-15% improvement in detection rat

Gravitational radiation timescales for extreme mass ratio inspirals

The capture and inspiral of compact stellar masses into massive black holes is an important source of low-frequency gravitational waves (with frequencies of ~1-100mHz), such as those that might be detected by the planned Laser Interferometer Space Antenna (LISA). Simulations of stellar clusters designed to study this problem typically rely on simple treatments of the black hole encounter which neglect some important features of orbits around black holes, such as the minimum radii of stable, non-plunging orbits. Incorporating an accurate representation of the orbital dynamics near a black hole has been avoided due to the large computational overhead. This paper provides new, more accurate, expressions for the energy and angular momentum lost by a compact object during a parabolic encounter with a non-spinning black hole, and the subsequent inspiral lifetime. These results improve on the Keplerian expressions which are now commonly used and will allow efficient computational simulations to be performed that account for the relativistic nature of the spacetime around the central black hole in the system.Comment: 19 pages, 4 figures. Changed in response to referee's report. Accepted for publication in Astrophysical Journa

Fuel cells for power generation and organic waste treatment on the island of Mull

In-situ use of biomass and organic waste streams have the potential to provide the key to energy self sustainability for islands and remote communities. Traditionally biogas fuels have been used in combustion engines for electric power generation. However, fuel cells offer the prospect of achieving higher generating efficiencies, and additionally, important environmental benefits can be achieved by way of mitigating greenhouse gas emissions, whilst providing a carbon sink. This paper presents the design details of a biogas gas plant and fuel cell installation that will provide a practical solution on an island (and be applicable in other remote and rural areas) where connection to the grid can be expensive, and where biofuels can be produced on site at no significant extra cost

Classifying LISA gravitational wave burst signals using Bayesian evidence

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.Comment: 21 pages, 11 figures, accepted by CQG; v2 has minor changes for consistency with accepted versio

Detection Strategies for Extreme Mass Ratio Inspirals

The capture of compact stellar remnants by galactic black holes provides a unique laboratory for exploring the near horizon geometry of the Kerr spacetime, or possible departures from general relativity if the central cores prove not to be black holes. The gravitational radiation produced by these Extreme Mass Ratio Inspirals (EMRIs) encodes a detailed map of the black hole geometry, and the detection and characterization of these signals is a major scientific goal for the LISA mission. The waveforms produced are very complex, and the signals need to be coherently tracked for hundreds to thousands of cycles to produce a detection, making EMRI signals one of the most challenging data analysis problems in all of gravitational wave astronomy. Estimates for the number of templates required to perform an exhaustive grid-based matched-filter search for these signals are astronomically large, and far out of reach of current computational resources. Here I describe an alternative approach that employs a hybrid between Genetic Algorithms and Markov Chain Monte Carlo techniques, along with several time saving techniques for computing the likelihood function. This approach has proven effective at the blind extraction of relatively weak EMRI signals from simulated LISA data sets.Comment: 10 pages, 4 figures, Updated for LISA 8 Symposium Proceeding

Cosmology with the lights off: Standard sirens in the Einstein Telescope era

We explore the prospects for constraining cosmology using gravitational-wave (GW) observations of neutron-star binaries by the proposed Einstein Telescope (ET), exploiting the narrowness of the neutron-star mass function. Double neutron-star (DNS) binaries are expected to be one of the first sources detected after "first-light" of Advanced LIGO and are expected to be detected at a rate of a few tens per year in the advanced era. However the proposed ET could catalog tens of thousands per year. Combining the measured source redshift distributions with GW-network distance determinations will permit not only the precision measurement of background cosmological parameters, but will provide an insight into the astrophysical properties of these DNS systems. Of particular interest will be to probe the distribution of delay times between DNS-binary creation and subsequent merger, as well as the evolution of the star-formation rate density within ET's detection horizon. Keeping H_0, \Omega_{m,0} and \Omega_{\Lambda,0} fixed and investigating the precision with which the dark-energy equation-of-state parameters could be recovered, we found that with 10^5 detected DNS binaries we could constrain these parameters to an accuracy similar to forecasted constraints from future CMB+BAO+SNIa measurements. Furthermore, modeling the merger delay-time distribution as a power-law, and the star-formation rate (SFR) density as a parametrized version of the Porciani and Madau SF2 model, we find that the associated astrophysical parameters are constrained to within ~ 10%. All parameter precisions scaled as 1/sqrt(N), where N is the number of cataloged detections. We also investigated how precisions varied with the intrinsic underlying properties of the Universe and with the distance reach of the network (which may be affected by the low-frequency cutoff of the detector).Comment: 24 pages, 11 figures, 6 tables. Minor changes to reflect published version. References updated and correcte

Cosmic Swarms: A search for Supermassive Black Holes in the LISA data stream with a Hybrid Evolutionary Algorithm

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 $5\sigma$ (as estimated from the Fisher Matrix) of the true answer, for both the actual and antipodal sky solutions.Comment: submitted to Classical & Quantum Gravity. 19 pages, 4 figure

Improved approximate inspirals of test-bodies into Kerr black holes

We present an improved version of the approximate scheme for generating inspirals of test-bodies into a Kerr black hole recently developed by Glampedakis, Hughes and Kennefick. Their original "hybrid" scheme was based on combining exact relativistic expressions for the evolution of the orbital elements (the semi-latus rectum p and eccentricity e) with approximate, weak-field, formula for the energy and angular momentum fluxes, amended by the assumption of constant inclination angle, iota, during the inspiral. Despite the fact that the resulting inspirals were overall well-behaved, certain pathologies remained for orbits in the strong field regime and for orbits which are nearly circular and/or nearly polar. In this paper we eliminate these problems by incorporating an array of improvements in the approximate fluxes. Firstly, we add certain corrections which ensure the correct behaviour of the fluxes in the limit of vanishing eccentricity and/or 90 degrees inclination. Secondly, we use higher order post-Newtonian formulae, adapted for generic orbits. Thirdly, we drop the assumption of constant inclination. Instead, we first evolve the Carter constant by means of an approximate post-Newtonian expression and subsequently extract the evolution of iota. Finally, we improve the evolution of circular orbits by using fits to the angular momentum and inclination evolution determined by Teukolsky based calculations. As an application of the improved scheme we provide a sample of generic Kerr inspirals and for the specific case of nearly circular orbits we locate the critical radius where orbits begin to decircularise under radiation reaction. These easy-to-generate inspirals should become a useful tool for exploring LISA data analysis issues and may ultimately play a role in source detection.Comment: 25 pages, 14 figures, some typos corrected, short section on conservative corrections added, minor changes for consistency with published versio

Observable Properties of Orbits in Exact Bumpy Spacetimes

We explore the properties of test-particle orbits in "bumpy" spacetimes - stationary, reflection-symmetric, asymptotically flat solutions of Einstein equations that have a non-Kerr (anomalous) higher-order multipole-moment structure but can be tuned arbitrarily close to the Kerr metric. Future detectors should observe gravitational waves generated during inspirals of compact objects into supermassive central bodies. If the central body deviates from the Kerr metric, this will manifest itself in the emitted waves. Here, we explore some of the features of orbits in non-Kerr spacetimes that might lead to observable signatures. As a basis for this analysis, we use a family of exact solutions proposed by Manko & Novikov which deviate from the Kerr metric in the quadrupole and higher moments, but we also compare our results to other work in the literature. We examine isolating integrals of the orbits and find that the majority of geodesic orbits have an approximate fourth constant of the motion (in addition to the energy, angular momentum and rest mass) and the resulting orbits are tri-periodic to high precision. We also find that this fourth integral can be lost for certain orbits in some oblately deformed Manko-Novikov spacetimes. However, compact objects will probably not end up on these chaotic orbits in nature. We compute the location of the innermost stable circular orbit (ISCO) and find that the behavior of orbtis near the ISCO can be qualitatively different depending on whether the ISCO is determined by the onset of an instability in the radial or vertical direction. Finally, we compute periapsis and orbital-plane precessions for nearly circular and nearly equatorial orbits in both the strong and weak field, and discuss weak-field precessions for eccentric equatorial orbits.Comment: 42 pages, 20 figures, accepted by Phys. Rev. D, v2 has minor changes to make it consistent with published versio