Source author record

Edward K. Porter

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

20 published item(s)

preprint2016arXiv

Detecting compact galactic binaries using a hybrid swarm-based algorithm

Compact binaries in our galaxy are expected to be one of the main sources of gravitational waves for the future eLISA mission. During the mission lifetime, many thousands of galactic binaries should be individually resolved. However, the identification of the sources, and the extraction of the signal parameters in a noisy environment are real challenges for data analysis. So far, stochastic searches have proven to be the most successful for this problem. In this work we present the first application of a swarm-based algorithm combining Particle Swarm Optimization and Differential Evolution. These algorithms have been shown to converge faster to global solutions on complicated likelihood surfaces than other stochastic methods. We first demonstrate the effectiveness of the algorithm for the case of a single binary in a 1 mHz search bandwidth. This interesting problem gave the algorithm plenty of opportunity to fail, as it can be easier to find a strong noise peak rather than the signal itself. After a successful detection of a fictitious low-frequency source, as well as the verification binary RXJ0806.3+1527, we then applied the algorithm to the detection of multiple binaries, over different search bandwidths, in the cases of low and mild source confusion. In all cases, we show that we can successfully identify the sources, and recover the true parameters within a 99\% credible interval.

preprint2015arXiv

Fisher vs. Bayes : A comparison of parameter estimation techniques for massive black hole binaries to high redshifts with eLISA

Massive black hole binaries are the primary source of gravitational waves (GW) for the future eLISA observatory. The detection and parameter estimation of these sources to high redshift would provide invaluable information on the formation mechanisms of seed black holes, and on the evolution of massive black holes and their host galaxies through cosmic time. The Fisher information matrix has been the standard tool for GW parameter estimation in the last two decades. However, recent studies have questioned the validity of using the Fisher matrix approach. For example, the Fisher matrix approach sometimes predicts errors of $\geq100\%$ in the estimation of parameters such as the luminosity distance and sky position. With advances in computing power, Bayesian inference is beginning to replace the Fisher matrix approximation in parameter estimation studies. In this work, we conduct a Bayesian inference analysis for 120 sources situated at redshifts of between $0.1\leq z\leq 13.2$, and compare the results with those from a Fisher matrix analysis. The Fisher matrix results suggest that for this particular selection of sources, eLISA would be unable to localize sources at redshifts of $z\lesssim6$. In contrast, Bayesian inference provides finite error estimations for all sources in the study, and shows that we can establish minimum closest distances for all sources. The study further predicts that we should be capable with eLISA, out to a redshift of at least $z\leq13$, of predicting a maximum error in the chirp mass of $\lesssim 1\%$, the reduced mass of $\lesssim20\%$, the time to coalescence of 2 hours, and to a redshift of $z\sim5$, the inclination of the source with a maximum error of $\sim60$ degrees.

preprint2015arXiv

Supermassive Black Hole Tests of General Relativity with eLISA

Motivated by the parameterized post-Einsteinian (ppE) scheme devised by Yunes and Pretorius, which introduces corrections to the post-Newtonian coefficients of the frequency domain gravitational waveform in order to emulate alternative theories of gravity, we compute analytical time domain waveforms that, after a numerical Fourier transform, aim to represent (phase corrected only) ppE waveforms. In this formalism, alternative theories manifest themselves via corrections to the phase and frequency, as predicted by General Relativity (GR), at different post-Newtonian (PN) orders. In order to present a generic test of alternative theories of gravity, we assume that the coupling constant of each alternative theory is manifestly positive, allowing corrections to the GR waveforms to be either positive or negative. By exploring the capabilities of massive black hole binary GR waveforms in the detection and parameter estimation of corrected time domain ppE signals, using the current eLISA configuration (as presented for the ESA Cosmic Vision L3 mission), we demonstrate that for corrections arising at higher than 1PN order in phase and frequency, GR waveforms are sufficient for both detecting and estimating the parameters of alternative theory signals. However, for theories introducing corrections at the 0 and 0.5 PN order, GR waveforms are not capable of covering the entire parameter space, requiring the use of non-GR waveforms for detection and parameter estimation.

preprint2015arXiv

The effect of different eLISA-like configurations on massive black hole parameter estimation

As the theme for the future L3 Cosmic Vision mission, ESA has recently chosen the `Gravitational Wave Universe'. Within this call, a mission concept called eLISA has been proposed. This observatory has a current initial configuration consisting of 4 laser links between the three satellites, which are separated by a distance of one million kilometers, constructing a single channel Michelson interferometer. However, the final configuration for the observatory will not be fixed until the end of this decade. With this in mind, we investigate the effect of different eLISA-like configurations on massive black hole detections. This work compares the results of a Bayesian inference study of 120 massive black hole binaries out to a redshift of $z\sim13$ for a $10^6$m arm-length eLISA with four and six links, as well as a $2\times10^6$m arm-length observatory with four links. We demonstrate that the original eLISA configuration should allow us to recover the luminosity distance of the source with an error of less than 10% out to a redshift of $z\sim4$, and a sky error box of $\leq10^2\,deg^2$ out to $z\sim0.1$. In contrast, both alternative configurations suggest that we should be able to conduct the same parameter recovery with errors of less than 10% in luminosity distance out to $z\sim12$ and $\leq10^2\,deg^2$ out to $z\sim0.4$. Using the information from these studies, we also infer that if we were able to construct a 2Gm, 6-link detector, the above values would shift to $z\sim20$ for luminosity distance and $z\sim0.9$ for sky error. While the final configuration will also be dependent on both technological and financial considerations, our study suggests that increasing the size of a two arm detector is a viable alternative to the inclusion of a third arm in a smaller detector. More importantly, this work further suggests no clear scientific loss between either choice.

preprint2014arXiv

A New Method of Accelerated Bayesian Inference for Comparable Mass Binaries in both Ground and Space-Based Gravitational Wave Astronomy

With the advance in computational resources, Bayesian inference is increasingly becoming the standard tool of practise in GW astronomy. However, algorithms such as Markov Chain Monte Carlo (MCMC) require a large number of iterations to guarantee convergence to the target density. Each chain demands a large number of evaluations of the likelihood function, and in the case of a Hessian MCMC, calculations of the Fisher information matrix for use as a proposal distribution. As each iteration requires the generation of at least one gravitational waveform, we very quickly reach a point of exclusion for current Bayesian algorithms, especially for low mass systems where the length of the waveforms is large and the waveform generation time is on the order of seconds. This suddenly demands a timescale of many weeks for a single MCMC. As each likelihood and Fisher information matrix calculation requires the evaluation of noise-weighted scalar products, we demonstrate that by using the linearity of integration, and the fact that more than 90% of the generation time is spent at frequencies less that one third of the maximum, we can construct composite integrals that speed up the MCMCs for comparable mass binaries by a factor of between 3.5 and 5.5, depending on the waveform length. This method is both source and detector type independent, and can be applied to any waveform that displays significant frequency evolution, such as stellar mass binaries with Advanced LIGO/Virgo, as well as supermassive black holes with eLISA

preprint2014arXiv

A Time Domain Waveform for Testing General Relativity

Gravitational-wave parameter estimation is only as good as the theory the waveform generation models are based upon. It is therefore crucial to test General Relativity (GR) once data becomes available. Many previous works, such as studies connected with the ppE framework by Yunes and Pretorius, rely on the stationary phase approximation (SPA) to model deviations from GR in the frequency domain. As Fast Fourier Transform algorithms have become considerably faster and in order to circumvent possible problems with the SPA, we test GR with corrected time domain waveforms instead of SPA waveforms. Since a considerable amount of work has been done already in the field using SPA waveforms, we establish a connection between leading-order-corrected waveforms in time and frequency domain, concentrating on phase-only corrected terms. In a Markov Chain Monte Carlo study, whose results are preliminary and will only be available later, we will assess the ability of the eLISA detector to measure deviations from GR for signals coming from supermassive black hole inspirals using these corrected waveforms.

preprint2014arXiv

The Challenges in Gravitational Wave Astronomy for Space-Based Detectors

The Gravitational Wave (GW) universe contains a wealth of sources which, with the proper treatment, will open up the universe as never before. By observing massive black hole binaries to high redshifts, we should begin to explore the formation process of seed black holes and track galactic evolution to the present day. Observations of extreme mass ratio inspirals will allow us to explore galactic centers in the local universe, as well as providing tests of General Relativity and constraining the value of Hubble's constant. The detection of compact binaries in our own galaxy may allow us to model stellar evolution in the Milky Way. Finally, the detection of cosmic (super)strings and a stochastic background would help us to constrain cosmological models. However, all of this depends on our ability to not only resolve sources and carry out parameter estimation, but also on our ability to define an optimal data analysis strategy. In this presentation, I will examine the challenges that lie ahead in GW astronomy for the ESA L3 Cosmic Vision mission, eLISA.

preprint2012arXiv

Bayesian Inference for LISA Pathfinder using Markov Chain Monte Carlo Methods

We present a parameter estimation procedure based on a Bayesian framework by applying a Markov Chain Monte Carlo algorithm to the calibration of the dynamical parameters of a space based gravitational wave detector. The method is based on the Metropolis-Hastings algorithm and a two-stage annealing treatment in order to ensure an effective exploration of the parameter space at the beginning of the chain. We compare two versions of the algorithm with an application to a LISA Pathfinder data analysis problem. The two algorithms share the same heating strategy but with one moving in coordinate directions using proposals from a multivariate Gaussian distribution, while the other uses the natural logarithm of some parameters and proposes jumps in the eigen-space of the Fisher Information matrix. The algorithm proposing jumps in the eigen-space of the Fisher Information matrix demonstrates a higher acceptance rate and a slightly better convergence towards the equilibrium parameter distributions in the application to LISA Pathfinder data . For this experiment, we return parameter values that are all within $\sim1σ$ of the injected values. When we analyse the accuracy of our parameter estimation in terms of the effect they have on the force-per-unit test mass noise estimate, we find that the induced errors are three orders of magnitude less than the expected experimental uncertainty in the power spectral density.

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.

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.

preprint2011arXiv

Calculating the continued fraction coefficients of a sub-diagonal Padé approximant at arbitrary order

The inspiral of two compact objects in gravitational wave astronomy is described by a post-Newtonian expansion in powers of $(v/c)$. In most cases, it is believed that the post-Newtonian expansion is asymptotically divergent. A standard technique for accelerating the convergence of a power series is to re-sum the series by means of a rational polynomial called a Padé approximation. If we liken this approximation to a matrix, the best convergence is achieved by staying close to a diagonal Padé approximation. This broadly presents two subsets of the approximation : a super-diagonal approximation $P^M_N$ and a sub-diagonal approximation $P_M^N$, where $M = N+ε$, and $ε$ takes the values of 0 or 1. Left as rational polynomials, the coefficients in both the numerator and denominator need to be re-calculated as the order of the initial power series approximation is increased. However, the sub-diagonal Padé approximant is computationally advantageous as it can be expressed in terms of a Gauss-like continued fraction. Once in this form, each coefficient in the continued fraction is uniquely determined at each order. This means that as we increase the order of approximation of the original power series, we now have only one new additional coefficient to calculate in the continued fraction. While it is possible to provide explicit expressions for the continued fraction coefficients, they rapidly become unwieldy at high orders of approximation. It is also possible to numerically calculate the coefficients by means of ratios of Hankel determinants. However, these determinants can be ill-conditioned and lead to numerical instabilities. In this article, we present a method for calculating the continued fraction coefficients at arbitrary orders of approximation.

preprint2010arXiv

The Effect of Data Gaps on LISA Galactic Binary Parameter Estimation

In the last few years there has been an enormous effort in parameter estimation studies for different sources with the space based gravitational wave detector, LISA. While these studies have investigated sources of differing complexity, the one thing they all have in common is they assume continuous data streams. In reality, the LISA data stream will contain gaps from such possible events such as repointing of the satellite antennae, to discharging static charge build up on the satellites, to disruptions due to micro-meteor strikes. In this work we conduct a large scale Monte Carlo parameter estimation simulation for galactic binaries assuming data streams containing gaps. As the expected duration and frequency of the gaps are currently unknown, we have decided to focus on gaps of approximately one hour, occurring either once per day or once per week. We also study the case where, as well as the expected periodic gaps, we have a data drop-out of one continuous week. Our results show that for for galactic binaries, a gap of once per week introduces a bias of between 0.5% and 1% in the estimation of parameters, for the most important parameters such as the sky position, amplitude and frequency. This number rises to between 3% and 7% for the case of one gap a day, and to between 4% and 9% when we have one gap a day and a spurious gap of a week. A future study will investigate the effect of data gaps on supermassive black hole binaries and extreme mass ratio inspirals.

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.

preprint2009arXiv

Binaries of massive black holes in rotating clusters: Dynamics, gravitational waves, detection and the role of eccentricity

The dynamical evolution of binaries of intermediate-massive black holes (IMBHs, massive black holes with a mass ranging between $10^2$ and $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 $\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 $10^{-4}$, reduced mass on the order of $10^{-3}$ and luminosity distance on the order of $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.

preprint2009arXiv

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σ$ (as estimated from the Fisher Matrix) of the true answer, for both the actual and antipodal sky solutions.

preprint2009arXiv

Data Analysis Challenges for the Einstein Telescope

The Einstein Telescope is a proposed third generation gravitational wave detector that will operate in the region of 1 Hz to a few kHz. As well as the inspiral of compact binaries composed of neutron stars or black holes, the lower frequency cut-off of the detector will open the window to a number of new sources. These will include the end stage of inspirals, plus merger and ringdown of intermediate mass black holes, where the masses of the component bodies are on the order of a few hundred solar masses. There is also the possibility of observing intermediate mass ratio inspirals, where a stellar mass compact object inspirals into a black hole which is a few hundred to a few thousand times more massive. In this article, we investigate some of the data analysis challenges for the Einstein Telescope such as the effects of increased source number, the need for more accurate waveform models and the some of the computational issues that a data analysis strategy might face.

preprint2009arXiv

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

preprint2009arXiv

Use of the MultiNest algorithm for gravitational wave data analysis

We describe an application of the MultiNest algorithm to gravitational wave data analysis. MultiNest is a multimodal nested sampling algorithm designed to efficiently evaluate the Bayesian evidence and return posterior probability densities for likelihood surfaces containing multiple secondary modes. The algorithm employs a set of live points which are updated by partitioning the set into multiple overlapping ellipsoids and sampling uniformly from within them. This set of live points climbs up the likelihood surface through nested iso-likelihood contours and the evidence and posterior distributions can be recovered from the point set evolution. The algorithm is model-independent in the sense that the specific problem being tackled enters only through the likelihood computation, and does not change how the live point set is updated. In this paper, we consider the use of the algorithm for gravitational wave data analysis by searching a simulated LISA data set containing two non-spinning supermassive black hole binary signals. The algorithm is able to rapidly identify all the modes of the solution and recover the true parameters of the sources to high precision.

preprint2007arXiv

The Search for Massive Black Hole Binaries with LISA

In this work we focus on the search and detection of Massive black hole binary (MBHB) systems, including systems at high redshift. As well as expanding on previous works where we used a variant of Markov Chain Monte Carlo (MCMC), called Metropolis-Hastings Monte Carlo, with simulated annealing, we introduce a new search method based on frequency annealing which leads to a more rapid and robust detection. We compare the two search methods on systems where we do and do not see the merger of the black holes. In the non-merger case, we also examine the posterior distribution exploration using a 7-D MCMC algorithm. We demonstrate that this method is effective in dealing with the high correlations between parameters, has a higher acceptance rate than previously proposed methods and produces posterior distribution functions that are close to the prediction from the Fisher Information matrix. Finally, after carrying out searches where there is only one binary in the data stream, we examine the case where two black hole binaries are present in the same data stream. We demonstrate that our search algorithm can accurately recover both binaries, and more importantly showing that we can safely extract the MBHB sources without contaminating the rest of the data stream.

preprint2003arXiv

Comparison of filters for detecting gravitational wave bursts in interferometric detectors

Filters developed in order to detect short bursts of gravitational waves in interferometric detector outputs are compared according to three main points. Conventional Receiver Operating Characteristics (ROC) are first built for all the considered filters and for three typical burst signals. Optimized ROC are shown for a simple pulse signal in order to estimate the best detection efficiency of the filters in the ideal case, while realistic ones obtained with filters working with several ``templates'' show how detection efficiencies can be degraded in a practical implementation. Secondly, estimations of biases and statistical errors on the reconstruction of the time of arrival of pulse-like signals are then given for each filter. Such results are crucial for future coincidence studies between Gravitational Wave detectors but also with neutrino or optical detectors. As most of the filters require a pre-whitening of the detector noise, the sensitivity to a non perfect noise whitening procedure is finally analysed. For this purpose lines of various frequencies and amplitudes are added to a Gaussian white noise and the outputs of the filters are studied in order to monitor the excess of false alarms induced by the lines. The comparison of the performances of the different filters finally show that they are complementary rather than competitive.