A hierarchical search for gravitational waves from supermassive black hole binary mergers

We present a method to search for gravitational waves from coalescing supermassive binary black holes in LISA data. The search utilizes the $\mathcal{F}$-statistic to maximize over, and determine the values of, the extrinsic parameters of the binary system. The intrinsic parameters are searched over hierarchically using stochastically generated multi-dimensional template banks to recover the masses and sky locations of the binary. We present the results of this method applied to the mock LISA data Challenge 1B data set.Comment: 11 pages, 2 figures, for GWDAW-12 proceedings edition of CQ

Evolutionary quantum cosmology in a gauge-fixed picture

We study the classical and quantum models of a flat Friedmann-Robertson-Walker (FRW) space-time, coupled to a perfect fluid, in the context of the consensus and a gauge-fixed Lagrangian frameworks. It is shown that, either in the usual or in the gauge-fixed actions, the evolution of the universe based on the classical cosmology represents a late time power law expansion, coming from a big-bang singularity in which the scale factor goes to zero for the standard matter, and tending towards a big-rip singularity in which the scale factor diverges for the phantom fluid. We then employ the familiar canonical quantization procedure in the given cosmological setting to find the cosmological wave functions in the corresponding minisuperspace. Using a gauge-fixed (reduced) Lagrangian, we show that, it may lead to a Schr\"{o}dinger equation for the quantum-mechanical description of the model under consideration, the eigenfunctions of which can be used to construct the time dependent wave function of the universe. We use the resulting wave function in order to investigate the possibility of the avoidance of classical singularities due to quantum effects by means of the many-worlds and ontological interpretation of quantum cosmology.Comment: 15 pages, 10 figures, typos corrected, Refs. adde

Forced motion near black holes

We present two methods for integrating forced geodesic equations in the Kerr spacetime, which can accommodate arbitrary forces. As a test case, we compute inspirals under a simple drag force, mimicking the presence of gas. We verify that both methods give the same results for this simple force. We find that drag generally causes eccentricity to increase throughout the inspiral. This is a relativistic effect qualitatively opposite to what is seen in gravitational-radiation-driven inspirals, and similar to what is observed in hydrodynamic simulations of gaseous binaries. We provide an analytic explanation by deriving the leading order relativistic correction to the Newtonian dynamics. If observed, an increasing eccentricity would provide clear evidence that the inspiral was occurring in a non-vacuum environment. Our two methods are especially useful for evolving orbits in the adiabatic regime. Both use the method of osculating orbits, in which each point on the orbit is characterized by the parameters of the geodesic with the same instantaneous position and velocity. Both methods describe the orbit in terms of the geodesic energy, axial angular momentum, Carter constant, azimuthal phase, and two angular variables that increase monotonically and are relativistic generalizations of the eccentric anomaly. The two methods differ in their treatment of the orbital phases and the representation of the force. In one method the geodesic phase and phase constant are evolved together as a single orbital phase parameter, and the force is expressed in terms of its components on the Kinnersley orthonormal tetrad. In the second method, the phase constants of the geodesic motion are evolved separately and the force is expressed in terms of its Boyer-Lindquist components. This second approach is a generalization of earlier work by Pound and Poisson for planar forces in a Schwarzschild background.Comment: 28 pages, 2 figures, submitted to Phys. Rev. D; v2 has minor changes for consistency with published version, plus a new section discussing the relative advantages of the two approache

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

Template-based searches for gravitational waves: efficient lattice covering of flat parameter spaces

The construction of optimal template banks for matched-filtering searches is an example of the sphere covering problem. For parameter spaces with constant-coefficient metrics a (near-) optimal template bank is achieved by the A_n* lattice, which is the best lattice-covering in dimensions n <= 5, and is close to the best covering known for dimensions n <= 16. Generally this provides a substantially more efficient covering than the simpler hyper-cubic lattice. We present an algorithm for generating lattice template banks for constant-coefficient metrics and we illustrate its implementation by generating A_n* template banks in n=2,3,4 dimensions.Comment: 10 pages, submitted to CQG for proceedings of GWDAW1

Random template placement and prior information

In signal detection problems, one is usually faced with the task of searching a parameter space for peaks in the likelihood function which indicate the presence of a signal. Random searches have proven to be very efficient as well as easy to implement, compared e.g. to searches along regular grids in parameter space. Knowledge of the parameterised shape of the signal searched for adds structure to the parameter space, i.e., there are usually regions requiring to be densely searched while in other regions a coarser search is sufficient. On the other hand, prior information identifies the regions in which a search will actually be promising or may likely be in vain. Defining specific figures of merit allows one to combine both template metric and prior distribution and devise optimal sampling schemes over the parameter space. We show an example related to the gravitational wave signal from a binary inspiral event. Here the template metric and prior information are particularly contradictory, since signals from low-mass systems tolerate the least mismatch in parameter space while high-mass systems are far more likely, as they imply a greater signal-to-noise ratio (SNR) and hence are detectable to greater distances. The derived sampling strategy is implemented in a Markov chain Monte Carlo (MCMC) algorithm where it improves convergence.Comment: Proceedings of the 8th Edoardo Amaldi Conference on Gravitational Waves. 7 pages, 4 figure

Heterogeneous Bond Percolation on Multitype Networks with an Application to Epidemic Dynamics

Considerable attention has been paid, in recent years, to the use of networks in modeling complex real-world systems. Among the many dynamical processes involving networks, propagation processes -- in which final state can be obtained by studying the underlying network percolation properties -- have raised formidable interest. In this paper, we present a bond percolation model of multitype networks with an arbitrary joint degree distribution that allows heterogeneity in the edge occupation probability. As previously demonstrated, the multitype approach allows many non-trivial mixing patterns such as assortativity and clustering between nodes. We derive a number of useful statistical properties of multitype networks as well as a general phase transition criterion. We also demonstrate that a number of previous models based on probability generating functions are special cases of the proposed formalism. We further show that the multitype approach, by naturally allowing heterogeneity in the bond occupation probability, overcomes some of the correlation issues encountered by previous models. We illustrate this point in the context of contact network epidemiology.Comment: 10 pages, 5 figures. Minor modifications were made in figures 3, 4 and 5 and in the text. Explanations and references were adde

Exit time as a measure of ecological resilience

Ecological resilience is the magnitude of the largest perturbation from which a system can still recover to its original state. However, a transition into another state may often be invoked by a series of minor synergistic perturbations rather than a single big one. We show how resilience can be estimated in terms of average life expectancy, accounting for this natural regime of variability. We use time series to fit a model that captures the stochastic as well as the deterministic components. The model is then used to estimate the mean exit time from the basin of attraction. This approach offers a fresh angle to anticipating the chance of a critical transition at a time when high-resolution time series are becoming increasingly available.</p

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