Source author record

Priscilla Canizares

Priscilla Canizares 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

14works
9topics
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

14 published item(s)

preprint2015arXiv

Accelerated gravitational-wave parameter estimation with reduced order modeling

Inferring the astrophysical parameters of coalescing compact binaries is a key science goal of the upcoming advanced LIGO-Virgo gravitational-wave detector network and, more generally, gravitational-wave astronomy. However, current parameter estimation approaches for such scenarios can lead to computationally intractable problems in practice. Therefore there is a pressing need for new, fast and accurate Bayesian inference techniques. In this letter we demonstrate that a reduced order modeling approach enables rapid parameter estimation studies. By implementing a reduced order quadrature scheme within the LIGO Algorithm Library, we show that Bayesian inference on the 9-dimensional parameter space of non-spinning binary neutron star inspirals can be sped up by a factor of 30 for the early advanced detectors' configurations. This speed-up will increase to about $150$ as the detectors improve their low-frequency limit to 10Hz, reducing to hours analyses which would otherwise take months to complete. Although these results focus on gravitational detectors, the techniques are broadly applicable to any experiment where fast Bayesian analysis is desirable.

preprint2014arXiv

Overcoming the Gauge Problem for the Gravitational Self-Force

The gravitational waves emitted by binary systems with extreme-mass ratios carry unique astrophysical information that can only be detected by space-based detectors like eLISA. To that end, a very accurate modelling of the system is required. The gravitational self-force program, which has been fully developed in the Lorenz gauge, is the best approach we have so far. However, the computations required would be done more efficiently if we could work in other gauges, like the Regge-Wheeler (RW) one in the case of Schwarzschild black holes. In this letter we present a new scheme, based on the Particle-without-Particle formulation of the field equations, where the gravitational self-force can be obtained from just solving individual wave-type equations like the master equations of the RW gauge. This approach can help to tackle the yet unsolved Kerr case.

preprint2013arXiv

Gravitational wave parameter estimation with compressed likelihood evaluations

One of the main bottlenecks in gravitational wave (GW) astronomy is the high cost of performing parameter estimation and GW searches on the fly. We propose a novel technique based on Reduced Order Quadratures (ROQs), an application and data-specific quadrature rule, to perform fast and accurate likelihood evaluations. These are the dominant cost in Markov chain Monte Carlo (MCMC) algorithms, which are widely employed in parameter estimation studies, and so ROQs offer a new way to accelerate GW parameter estimation. We illustrate our approach using a four dimensional GW burst model embedded in noise. We build an ROQ for this model, and perform four dimensional MCMC searches with both the standard and ROQs quadrature rules, showing that, for this model, the ROQ approach is around 25 times faster than the standard approach with essentially no loss of accuracy. The speed-up from using ROQs is expected to increase for more complex GW signal models and therefore has significant potential to accelerate parameter estimation of GW sources such as compact binary coalescences.

preprint2012arXiv

Constraining Gravity with LISA Detections of Binaries

General Relativity (GR) describes gravitation well at the energy scales which we have so far been able to achieve or detect. However, we do not know whether GR is behind the physics governing stronger gravitational field regimes, such as near neutron stars or massive black-holes (MBHs). Gravitational-wave (GW) astronomy is a promising tool to test and validate GR and/or potential alternative theories of gravity. The information that a GW waveform carries not only will allow us to map the strong gravitational field of its source, but also determine the theory of gravity ruling its dynamics. In this work, we explore the extent to which we could distinguish between GR and other theories of gravity through the detection of low-frequency GWs from extreme-mass-ratio inspirals (EMRIs) and, in particular, we focus on dynamical Chern-Simons modified gravity (DCSMG). To that end, we develop a framework that enables us, for the first time, to perform a parameter estimation analysis for EMRIs in DCSMG. Our model is described by a 15-dimensional parameter space, that includes the Chern-Simons (CS) parameter which characterises the deviation between the two theories, and our analysis is based on Fisher information matrix techniques together with a (maximum-mismatch) criterion to assess the validity of our results. In our analysis, we study a 5-dimensional parameter space, finding that a GW detector like the Laser Interferometer Space Antenna (LISA) or eLISA (evolved LISA) should be able to discriminate between GR and DCSMG with fractional errors below 5%, and hence place bounds four orders of magnitude better than current Solar System bounds.

preprint2012arXiv

Testing Chern-Simons Modified Gravity with Gravitational-Wave Detections of Extreme-Mass-Ratio Binaries

[abridged] The detection of gravitational waves from extreme-mass-ratio (EMRI) binaries, comprising a stellar-mass compact object orbiting around a massive black hole, is one of the main targets for low-frequency gravitational-wave detectors in space, like the Laser Interferometer Space Antenna (LISA or eLISA/NGO). The long-duration gravitational-waveforms emitted by such systems encode the structure of the strong field region of the massive black hole, in which the inspiral occurs. The detection and analysis of EMRIs will therefore allow us to study the geometry of massive black holes and determine whether their nature is as predicted by General Relativity and even to test whether General Relativity is the correct theory to describe the dynamics of these systems. To achieve this, EMRI modeling in alternative theories of gravity is required to describe the generation of gravitational waves. In this paper, we explore to what extent EMRI observations with LISA or eLISA/NGO might be able to distinguish between General Relativity and a particular modification of it, known as Dynamical Chern-Simons Modified Gravity. Our analysis is based on a parameter estimation study that uses approximate gravitational waveforms obtained via a radiative-adiabatic method and is restricted to a five-dimensional subspace of the EMRI configuration space. This includes a Chern-Simons parameter that controls the strength of gravitational deviations from General Relativity. We find that, if Dynamical Chern-Simons Modified Gravity is the correct theory, an observatory like LISA or even eLISA/NGO should be able to measure the Chern-Simons parameter with fractional errors below 5%. If General Relativity is the true theory, these observatories should put bounds on this parameter at the level xi^(1/4) < 10^4 km, which is four orders of magnitude better than current Solar System bounds.

preprint2012arXiv

Testing Chern-Simons modified gravity with observations of extreme-mass-ratio binaries

Extreme-Mass-Ratio Inspirals (EMRIs) are one of the most promising sources of gravitational waves (GWs) for space-based detectors like the Laser Interferometer Space Antenna (LISA). EMRIs consist of a compact stellar object orbiting around a massive black hole (MBH). Since EMRI signals are expected to be long lasting (containing of the order of hundred thousand cycles), they will encode the structure of the MBH gravitational potential in a precise way such that features depending on the theory of gravity governing the system may be distinguished. That is, EMRI signals may be used to test gravity and the geometry of black holes. However, the development of a practical methodology for computing the generation and propagation of GWs from EMRIs in theories of gravity different than General Relativity (GR) has only recently begun. In this paper, we present a parameter estimation study of EMRIs in a particular modification of GR, which is described by a four-dimensional Chern-Simons (CS) gravitational term. We focus on determining to what extent a space-based GW observatory like LISA could distinguish between GR and CS gravity through the detection of GWs from EMRIs.

preprint2011arXiv

Are Time-Domain Self-Force Calculations Contaminated by Jost Solutions?

The calculation of the self force in the modeling of the gravitational-wave emission from extreme-mass-ratio binaries is a challenging task. Here we address the question of the possible emergence of a persistent spurious solution in time-domain schemes, referred to as a {\em Jost junk solution} in the literature, that may contaminate self force calculations. Previous studies suggested that Jost solutions are due to the use of zero initial data, which is inconsistent with the singular sources associated with the small object, described as a point mass. However, in this work we show that the specific origin is an inconsistency in the translation of the singular sources into jump conditions. More importantly, we identify the correct implementation of the sources at late times as the sufficient condition guaranteeing the absence of Jost junk solutions.

preprint2011arXiv

Time-domain modelling of Extreme-Mass-Ratio Inspirals for the Laser Interferometer Space Antenna

When a stellar-mass compact object is captured by a supermassive black hole located in a galactic centre, the system losses energy and angular momentum by the emission of gravitational waves. Subsequently, the stellar compact object evolves inspiraling until plunging onto the massive black hole. These EMRI systems are expected to be one of the main sources of gravitational waves for the future space-based Laser Interferometer Space Antenna (LISA). However, the detection of EMRI signals will require of very accurate theoretical templates taking into account the gravitational self-force, which is the responsible of the stellar-compact object inspiral. Due to its potential applicability on EMRIs, the obtention of an efficient method to compute the scalar self-force acting on a point-like particle orbiting around a massive black hole is being object of increasing interest. We present here a review of our time-domain numerical technique to compute the self-force acting on a point-like particle and we show its suitability to deal with both circular and eccentric orbits.

preprint2011arXiv

Tuning Time-Domain Pseudospectral Computations of the Self-Force on a Charged Scalar Particle

The computation of the self-force constitutes one of the main challenges for the construction of precise theoretical waveform templates in order to detect and analyze extreme-mass-ratio inspirals with the future space-based gravitational-wave observatory LISA. Since the number of templates required is quite high, it is important to develop fast algorithms both for the computation of the self-force and the production of waveforms. In this article we show how to tune a recent time-domain technique for the computation of the self-force, what we call the Particle without Particle scheme, in order to make it very precise and at the same time very efficient. We also extend this technique in order to allow for highly eccentric orbits.

preprint2010arXiv

An Efficient Time-Domain Method to Model Extreme-Mass-Ratio Inspirals

The gravitational-wave signals emitted by Extreme-Mass-Ratio Inspirals will be hidden in the instrumental LISA noise and the foreground noise produced by galactic binaries in the LISA band. Then, we need accurate gravitational-wave templates to extract these signals from the noise and obtain the relevant physical parameters. This means that in the modeling of these systems we have to take into account how the orbit of the stellar-mass compact object is modified by the action of its own gravitational field. This effect can be described as the action of a local force, the self-force. We present a time-domain technique to compute the self-force for geodesic eccentric orbits around a non-rotating massive black hole. To illustrate the method we have applied it to a testbed model consisting of scalar charged particle orbiting a non-dynamical black hole. A key feature of our method is that it does not introduce a small scale associated with the stellar-mass compact object. This is achieved by using a multidomain framework where the particle is located at the interface between two subdomains. In this way, we just have to evolve homogeneous wave-like equations with smooth solutions that have to be communicated across the subdomain boundaries using appropriate junction conditions. The numerical technique that we use to implement this scheme is the pseudospectral collocation method. We show the suitability of this technique for the modeling of Extreme-Mass-Ratio Inspirals and show that it can provide accurate results for the self-force.

preprint2010arXiv

Modelling Extreme-Mass-Ratio Inspirals using Pseudospectral Methods

We introduce a new time-domain method for computing the self-force acting on a scalar particle in a Schwarzschild geometry. The principal feature of our method consists in the division of the spatial domain into several subdomains and locating the particle at the interface betweem two them. In this way, we avoid the need of resolving a small length scale associated with the presence of a particle in the computational domain and, at the same time, we avoid numerical problems due to the low differentiability of solutions of equations with point-like singular behaviour.

preprint2010arXiv

The LISA PathFinder DMU and Radiation Monitor

The LISA PathFinder DMU (Data Management Unit) flight model was formally accepted by ESA and ASD on 11 February 2010, after all hardware and software tests had been successfully completed. The diagnostics items are scheduled to be delivered by the end of 2010. In this paper we review the requirements and performance of this instrumentation, specially focusing on the Radiation Monitor and the DMU, as well as the status of their programmed use during mission operations, on which work is ongoing at the time of writing.

preprint2009arXiv

An Efficient Pseudospectral Method for the Computation of the Self-force on a Charged Particle: Circular Geodesics around a Schwarzschild Black Hole

The description of the inspiral of a stellar-mass compact object into a massive black hole sitting at a galactic centre is a problem of major relevance for the future space-based gravitational-wave observatory LISA (Laser Interferometer Space Antenna), as the signals from these systems will be buried in the data stream and accurate gravitational-wave templates will be needed to extract them. The main difficulty in describing these systems lies in the estimation of the gravitational effects of the stellar-mass compact object on his own trajectory around the massive black hole, which can be modeled as the action of a local force, the self-force. In this paper, we present a new time-domain numerical method for the computation of the self-force in a simplified model consisting of a charged scalar particle orbiting a nonrotating black hole. We use a multi-domain framework in such a way that the particle is located at the interface between two domains so that the presence of the particle and its physical effects appear only through appropriate boundary conditions. In this way we eliminate completely the presence of a small length scale associated with the need of resolving the particle. This technique also avoids the problems associated with the impact of a low differentiability of the solution in the accuracy of the numerical computations. The spatial discretization of the field equations is done by using the pseudospectral collocation method and the time evolution, based on the method of lines, uses a Runge-Kutta solver. We show how this special framework can provide very efficient and accurate computations in the time domain, which makes the technique amenable for the intensive computations required in the astrophysically-relevant scenarios for LISA.

preprint2008arXiv

Simulations of Extreme-Mass-Ratio Inspirals Using Pseudospectral Methods

Extreme-mass-ratio inspirals (EMRIs), stellar-mass compact objects (SCOs) inspiralling into a massive black hole, are one of the main sources of gravitational waves expected for the Laser Interferometer Space Antenna (LISA). To extract the EMRI signals from the expected LISA data stream, which will also contain the instrumental noise as well as other signals, we need very accurate theoretical templates of the gravitational waves that they produce. In order to construct those templates we need to account for the gravitational backreaction, that is, how the gravitational field of the SCO affects its own trajectory. In general relativity, the backreaction can be described in terms of a local self-force, and the foundations to compute it have been laid recently. Due to its complexity, some parts of the calculation of the self-force have to be performed numerically. Here, we report on an ongoing effort towards the computation of the self-force based on time-domain multi-grid pseudospectral methods.