Source author record

Carlos F. Sopuerta

Carlos F. Sopuerta 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

37works
17topics
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

37 published item(s)

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2021arXiv

Darboux Covariance: A Hidden Symmetry of Perturbed Schwarzschild Black Holes

Starting from the infinite set of possible master equations for the perturbations of Schwarzschild black holes, with master functions linear in the metric perturbations and their first-order derivatives, we show that of all them are connected via Darboux transformations. These transformations preserve physical quantities like the quasinormal mode frequencies and the infinite hierarchy of Korteweg-de Vries conserved quantities, revealing a new hidden symmetry in the description of the perturbations of Schwarzschild black holes: Darboux covariance.

preprint2020arXiv

Motion of localized sources in general relativity: gravitational self-force from quasilocal conservation laws

An idealized "test" object in general relativity moves along a geodesic. However, if the object has a finite mass, this will create additional curvature in the spacetime, causing it to deviate from geodesic motion. If the mass is nonetheless sufficiently small, such an effect is usually treated perturbatively and is known as the gravitational self-force due to the object. This issue is still an open problem in gravitational physics today, motivated not only by basic foundational interest, but also by the need for its direct application in gravitational-wave astronomy. In particular, the observation of extreme-mass-ratio inspirals by the future space-based detector LISA will rely crucially on an accurate modeling of the self-force driving the orbital evolution and gravitational wave emission of such systems. In this paper, we present a novel derivation, based on conservation laws, of the basic equations of motion for this problem. They are formulated with the use of a quasilocal (rather than matter) stress-energy-momentum tensor---in particular, the Brown-York tensor---so as to capture gravitational effects in the momentum flux of the object, including the self-force. Our formulation and resulting equations of motion are independent of the choice of the perturbative gauge. We show that, in addition to the usual gravitational self-force term, they also lead to an additional "self-pressure" force not found in previous analyses, and also that our results correctly recover known formulas under appropriate conditions. Our approach thus offers a fresh geometrical picture from which to understand the self-force fundamentally, and potentially useful new avenues for computing it practically.

preprint2020arXiv

Numerical Relativity with Arbitrary Precision Arithmetic: Applications to Gravitational Collapse

Numerical Relativity is a mature field with many applications in Astrophysics, Cosmology and even in Fundamental Physics. As such, we are entering a stage in which new sophisticated methods adapted to open problems are being developed. In this paper, we advocate the use of Pseudo-Spectral Collocation (PSC) methods in combination with high-order precision arithmetic for Numerical Relativity problems with high accuracy and performance requirements. The PSC method provides exponential convergence (for smooth problems, as is the case in many problems in Numerical Relativity) and we can use different bit precision without the need of changing the structure of the numerical algorithms. Moreover, the PSC method provides high-compression storage of the information. We introduce a series of techniques for combining these tools and show their potential in two problems in relativistic gravitational collapse: (i) The classical Choptuik collapse, estimating with arbitrary precision the location of the apparent horizon. (ii) Collapse in asympotically anti-de Sitter spacetimes, showing that the total energy is preserved by the numerical evolution to a very high degree of precision.

preprint2016arXiv

Entropy theorems in classical mechanics, general relativity, and the gravitational two-body problem

In classical Hamiltonian theories, entropy may be understood either as a statistical property of canonical systems, or as a mechanical property, that is, as a monotonic function of the phase space along trajectories. In classical mechanics, there are theorems which have been proposed for proving the non-existence of entropy in the latter sense. We explicate, clarify and extend the proofs of these theorems to some standard matter (scalar and electromagnetic) field theories in curved spacetime, and then we show why these proofs fail in general relativity; due to properties of the gravitational Hamiltonian and phase space measures, the second law of thermodynamics holds. As a concrete application, we focus on the consequences of these results for the gravitational two-body problem, and in particular, we prove the non-compactness of the phase space of perturbed Schwarzschild-Droste spacetimes. We thus identify the lack of recurring orbits in phase space as a distinct sign of dissipation and hence entropy production.

preprint2016arXiv

Moving Closer to the Collapse of a Massless Scalar Field in Spherically Symmetric Anti-de Sitter Spacetimes

We present a new hybrid Cauchy-characteristic evolution method that is particularly suited for the study of gravitational collapse in spherically-symmetric asymptotically (global) Anti-de Sitter (AdS) spacetimes. The Cauchy evolution allows us to track the scalar field through the different bounces off the AdS boundary while the characteristic method can bring us very close to the point of formation of an apparent horizon. Here, we describe all the details of the method, including the transition between the two evolution schemes and the details of the numerical implementation for the case of massless scalar fields. We use this scheme to provide more numerical evidence for a recent conjecture on the power-law scaling of the apparent horizon mass resulting from the collapse of subcritical configurations. We also compute the critical exponents and echoing periods for a number of critical points and confirm the expectation that their values should be the same as in the asymptotically-flat case.

preprint2016arXiv

New Features of Gravitational Collapse in Anti-de Sitter Spacetimes

Gravitational collapse of a massless scalar field in spherically symmetric anti-de Sitter (AdS) spacetimes presents a new phenomenology with a series of critical points whose dynamics is discretely self-similar as in the asymptotically flat case. Each critical point is the limit of a branch of scalar field configurations that have bounced off the AdS boundary a fixed number of times before forming an apparent horizon. We present results from a numerical study that focus on the interfaces between branches. We find that there is a mass gap between branches and that subcritical configurations near the critical point form black holes with an apparent horizon mass that follows a power law of the form $M_{AH}-M_{g} \propto (p_{c}-p)^ξ$, where $M_g$ is the mass gap and the exponent $ξ\simeq 0.7$ appears to be universal.

preprint2014arXiv

Bayesian Model Selection for LISA Pathfinder

The main goal of the LISA Pathfinder (LPF) mission is to fully characterize the acceleration noise models and to test key technologies for future space-based gravitational-wave observatories similar to the eLISA concept. The data analysis team has developed complex three-dimensional models of the LISA Technology Package (LTP) experiment on-board LPF. These models are used for simulations, but more importantly, they will be used for parameter estimation purposes during flight operations. One of the tasks of the data analysis team is to identify the physical effects that contribute significantly to the properties of the instrument noise. A way of approaching this problem is to recover the essential parameters of a LTP model fitting the data. Thus, we want to define the simplest model that efficiently explains the observations. To do so, adopting a Bayesian framework, one has to estimate the so-called Bayes Factor between two competing models. In our analysis, we use three main different methods to estimate it: The Reversible Jump Markov Chain Monte Carlo method, the Schwarz criterion, and the Laplace approximation. They are applied to simulated LPF experiments where the most probable LTP model that explains the observations is recovered. The same type of analysis presented in this paper is expected to be followed during flight operations. Moreover, the correlation of the output of the aforementioned methods with the design of the experiment is explored.

preprint2014arXiv

Design of a dual species atom interferometer for space

Atom interferometers have a multitude of proposed applications in space including precise measurements of the Earth's gravitational field, in navigation & ranging, and in fundamental physics such as tests of the weak equivalence principle (WEP) and gravitational wave detection. While atom interferometers are realized routinely in ground-based laboratories, current efforts aim at the development of a space compatible design optimized with respect to dimensions, weight, power consumption, mechanical robustness and radiation hardness. In this paper, we present a design of a high-sensitivity differential dual species $^{85}$Rb/$^{87}$Rb atom interferometer for space, including physics package, laser system, electronics and software. The physics package comprises the atom source consisting of dispensers and a 2D magneto-optical trap (MOT), the science chamber with a 3D-MOT, a magnetic trap based on an atom chip and an optical dipole trap (ODT) used for Bose-Einstein condensate (BEC) creation and interferometry, the detection unit, the vacuum system for $10^{-11}$ mbar ultra-high vacuum generation, and the high-suppression factor magnetic shielding as well as the thermal control system. The laser system is based on a hybrid approach using fiber-based telecom components and high-power laser diode technology and includes all laser sources for 2D-MOT, 3D-MOT, ODT, interferometry and detection. Manipulation and switching of the laser beams is carried out on an optical bench using Zerodur bonding technology. The instrument consists of 9 units with an overall mass of 221 kg, an average power consumption of 608 W (819 W peak), and a volume of 470 liters which would well fit on a satellite to be launched with a Soyuz rocket, as system studies have shown.

preprint2014arXiv

Electromagnetic waves from neutron stars and black holes driven by polar gravitational perturbations

Neutron stars and black holes are the most compact astrophysical objects we can think of and as a consequence they are the main sources of gravitational waves. There are many astrophysically relevant scenarios in which these objects are immersed in or endowed with strong magnetic fields, in such a way that gravitational perturbations can couple to electromagnetic ones and can potentially trigger synergistic electromagnetic signatures. In a recent paper we derived the main equations for gravito-electromagnetic perturbations and studied in detail the case of polar electromagnetic perturbations driven by axial gravitational perturbations. In this paper we deal with the case of axial electromagnetic perturbations driven by polar black-hole or neutron stars oscillations, in which the energy emitted in case is considerably larger than in the previous case. In the case of neutron stars the phenomenon lasts considerably longer since the fluid acts as an energy reservoir that shakes the magnetic field for a timescale of the order of secs.

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.

preprint2014arXiv

Research Update on Extreme-Mass-Ratio Inspirals

The inspirals of stellar-mass mass compact objects into massive black holes in the centres of galaxies are one of the most important sources of gravitational radiation for space-based detectors like LISA or eLISA. These extreme-mass-ratio inspirals (EMRIs) will enable an ambitious research program with implications for astrophysics, cosmology, and fundamental physics. This article is a summary of the talks delivered at the plenary session on EMRIs at the 10th International LISA Symposium. It contains research updates on the following topics: astrophysics of EMRIs; EMRI science potential; and EMRI modeling.

preprint2013arXiv

Gravitationally Driven Electromagnetic Perturbations of Neutron Stars and Black Holes

Gravitational perturbations of neutron stars and black holes are well known sources of gravitational radiation. If the compact object is immersed in or endowed with a magnetic field, the gravitational perturbations would couple to electromagnetic perturbations and potentially trigger synergistic electromagnetic signatures. We present a detailed analytic calculation of the dynamics of coupled gravitational and electromagnetic perturbations for both neutron stars and black holes. We discuss the prospects for detecting the electromagnetic waves in these scenarios and the potential that these waves have for providing information about their source.

preprint2012arXiv

Approximate Waveforms for Extreme-Mass-Ratio Inspirals: The Chimera Scheme

We describe a new kludge scheme to model the dynamics of generic extreme-mass-ratio inspirals (EMRIs; stellar compact objects spiraling into a spinning supermassive black hole) and their gravitational-wave emission. The Chimera scheme is a hybrid method that combines tools from different approximation techniques in General Relativity: (i) A multipolar, post-Minkowskian expansion for the far-zone metric perturbation (the gravitational waveforms) and for the local prescription of the self-force; (ii) a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and (iii) a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The EMRI trajectory is made out of Kerr geodesic fragments joined via the method of osculating elements as dictated by the multipolar post-Minkowskian radiation-reaction prescription. We implemented the proper coordinate mapping between Boyer-Lindquist coordinates, associated with the Kerr geodesics, and harmonic coordinates, associated with the multipolar post-Minkowskian decomposition. The Chimera scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme to intermediate mass ratios, and hence, it can provide valuable information for future space-based gravitational-wave observatories, like LISA, and even for advanced ground detectors. The local character in time of our multipolar post-Minkowskian self-force makes this scheme amenable to study the possible appearance of transient resonances in generic inspirals.

preprint2012arXiv

Blocking low-eccentricity EMRIs: A statistical direct-summation N-body study of the Schwarzschild barrier

The capture of a compact object in a galactic nucleus by a massive black hole (MBH), an extreme-mass ratio inspiral (EMRI), is the best way to map space and time around it. Recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on low-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statistical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian and also for the first time in a direct-summation code a geodesic approximation for the relativistic orbits. The existence of the barrier prevents low-eccentricity EMRIs from approaching the central MBH, but high-eccentricity EMRIs, which have been wrongly classified as "direct plunges" until recently, ignore the presence of the barrier, because they are driven by two-body relaxation. Hence, since the rates are significantly affected in the case of low-eccentricity EMRIs, we predict that a LISA-like observatory such as eLISA will predominantly detect high-eccentricity EMRIs.

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

NR/HEP: roadmap for the future

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

preprint2012arXiv

Probing the strong gravity regime with eLISA: Progress on EMRIs

The capture of a stellar-mass compact object by a supermassive black hole and the subsequent inspiral (driven by gravitational radiation emission) constitute one of the most important sources of gravitational waves for space-based observatories like eLISA/NGO. In this article we describe their potential as high-precision tools that can be used to perform tests of the geometry of black holes and also of the strong field regime of gravity.

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.

preprint2012arXiv

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

preprint2012arXiv

The role of the supermassive black hole spin in the estimation of the EMRI event rate

A common result to all EMRI investigations on rates is that the possibility that a compact object merges with the MBH after only one intense burst of GWs is much more likely than a slow adiabatic inspiral, an EMRI. The later is referred to as a "plunge" because the compact object dives into the MBH. The event rates for plunges are orders of magnitude larger than slow inspirals. On the other hand, nature MBH's are most likely Kerr and the magnitude of the spin has been sized up to be high. We calculate the number of periapsis passages that a compact object set on to an extremely radial orbit goes through before being actually swallowed by the Kerr MBH and we then translate it into an event rate for a LISA-like observatory, such as the proposed ESA mission eLISA/NGO. We prove that a "plunging" compact object is conceptually indistinguishable from an adiabatic, slow inspiral; plunges spend on average up to hundred of thousands of cycles in the bandwidth of the detector for a two years mission. This has an important impact on the event rate, enhancing in some cases significantly, depending on the spin of the MBH and the inclination. Moreover, it has been recently proved that the production of low-eccentricity EMRIs is severely blocked by the presence of a blockade in the rate at which orbital angular momenta change takes place. This is the result of relativistic precession on to the stellar potential torques and hence affects EMRIs originating via resonant relaxation at distances of about $\sim 10^{-2}$ pc from the MBH. Since high-eccentricity EMRIs are a result of two-body relaxation, they are not affected by this phenomenon. Therefore we predict that eLISA EMRI event rates will be dominated by high-eccentricity binaries, as we present here.

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

Black-hole binaries go to eleven orbits

We analyse an eleven-orbit inspiral of a non-spinning black-hole binary with mass ratio q=M1/M2=4. The numerically obtained gravitational waveforms are compared with post-Newtonian (PN) predictions including several sub-dominant multipoles up to multipolar indices (l=5,m=5). We find that (i) numerical and post-Newtonian predictions of the phase of the (2,2) mode accumulate a phase difference of about 0.35 rad at the PN cut off frequency 0.1 for the Taylor T1 approximant when numerical and PN waveforms are matched over a window in the early inspiral phase; (ii) in contrast to previous studies of equal-mass and specific spinning binaries, we find the Taylor T4 approximant to agree less well with numerical results, provided the latter are extrapolated to infinite extraction radius; (iii) extrapolation of gravitational waveforms to infinite extraction radius is particularly important for subdominant multipoles with l unequal m; (iv) 3PN terms in post-Newtonian multipole expansions significantly improve the agreement with numerical predictions for sub-dominant multipoles.

preprint2011arXiv

New Kludge Scheme for the Construction of Approximate Waveforms for Extreme-Mass-Ratio Inspirals

We introduce a new kludge scheme to model the dynamics of generic extreme mass-ratio inspirals (stellar compact objects spiraling into a spinning supermassive black hole) and to produce the gravitational waveforms that describe the gravitational-wave emission of these systems. This scheme combines tools from different techniques in General Relativity: It uses a multipolar, post-Minkowskian (MPM) expansion for the far-zone metric perturbation (which provides the gravitational waveforms, here taken up to mass hexadecapole and current octopole order) and for the local prescription of the self-force (since we are lacking a general prescription for it); a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The orbital evolution is thus equivalent to solving the geodesic equations with time-dependent orbital elements, as dictated by the MPM radiation-reaction prescription. To complete the scheme, both the orbital evolution and wave generation require to map the Boyer-Lindquist coordinates of the orbits to the harmonic coordinates in which the different MPM quantities have been derived, a mapping that we provide explicitly in this paper. This new kludge scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme mass ratios to systems with more moderate mass ratios, and hence can provide valuable information for future space-based gravitational-wave observatories like LISA and even for advanced ground detectors. Finally, due to the local character in time of our MPM self-force, this scheme can be used to perform studies of the possible appearance of transient resonances in generic inspirals.

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

Gravitational Waves Notes, Issue #4 : "A Roadmap to Fundamental Physics from LISA EMRI Observations"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Carlos F. Sopuerta to write the highlight article of this issue on his thoughts on fundamental physics with LISA.

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

Testing Modified Gravity with Gravitational Wave Astronomy

The emergent area of gravitational wave astronomy promises to provide revolutionary discoveries in the areas of astrophysics, cosmology, and fundamental physics. One of the most exciting possibilities is to use gravitational-wave observations to test alternative theories of gravity. In this contribution we describe how to use observations of extreme-mass-ratio inspirals by the future Laser Interferometer Space Antenna to test a particular class of theories: Chern-Simons modified gravity.

preprint2010arXiv

Towards Tests of Alternative Theories of Gravity with LISA

The inspiral of stellar compact objects into massive black holes, usually known as extreme-mass-ratio inspirals (EMRIs), is one of the most important sources of gravitational-waves for the future Laser Interferometer Space Antenna (LISA). Intermediate-mass-ratio inspirals (IMRIs are also of interest to advance ground-based gravitational-wave observatories. We discuss here how modifications to the gravitational interaction can affect the signals emitted by these systems and their detectability by LISA. We concentrate in particular on Chern-Simons modified gravity, a theory that emerges in different quantum gravitational approaches.

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.

preprint2007arXiv

Gravitational Recoil from Binary Black Hole Mergers: the Close-Limit Approximation

[abridged] The coalescence of a binary black hole system is one of the main sources of gravitational waves that present and future detectors will study. Apart from the energy and angular momentum that these waves carry, for unequal-mass binaries there is also a net flux of linear momentum that implies a recoil velocity of the resulting final black hole in the opposite direction. We present a computation of the recoil velocity based on the close-limit approximation scheme, which gives excellent results for head-on and grazing collisions of black holes when compared to full numerical relativistic calculations. We obtain a maximum recoil velocity of ~ 64 km/s for a symmetric mass ratio eta = M_1 M_2/(M_1+M_2)^2 ~ 0.19 and an initial proper separation of 4 M, where M is the total ADM mass of the system. This separation is the maximum at which the close-limit approximation is expected to provide accurate results. If we supplement this estimate with PN calculations up to the innermost stable circular orbit, we obtain a lower bound for the recoil velocity, with a maximum around 84 km/s. This is a lower bound because it neglects the initial merger phase. We can however obtain a rough estimate by using PN methods or the close-limit approximation. Since both methods are known to overestimate the amount of radiation, we obtain in this way an upper bound for the recoil with maxima in the range of 220-265 km/s. We also provide non-linear fits to these estimated upper and lower bounds. These estimates are subject to uncertainties related to issues such as the choice of initial data and higher effects in perturbation theory. Nonetheless, our estimates are consistent with previous results in the literature and suggest a narrower range of possible recoil velocities.

preprint2003arXiv

Stacking a 4D geometry into an Einstein-Gauss-Bonnet bulk

In Einstein gravity there is a simple procedure to build D-dimensional spacetimes starting from (D-1)-dimensional ones, by stacking any (D-1)-dimensional Ricci-flat metric into the extra-dimension. We analyze this procedure in the context of Einstein-Gauss-Bonnet gravity, and find that it can only be applied to metrics with a constant Krestschmann scalar. For instance, we show that solutions of the black-string type are not allowed in this framework.

preprint2003arXiv

Two-parameter non-linear spacetime perturbations: gauge transformations and gauge invariance

An implicit fundamental assumption in relativistic perturbation theory is that there exists a parametric family of spacetimes that can be Taylor expanded around a background. The choice of the latter is crucial to obtain a manageable theory, so that it is sometime convenient to construct a perturbative formalism based on two (or more) parameters. The study of perturbations of rotating stars is a good example: in this case one can treat the stationary axisymmetric star using a slow rotation approximation (expansion in the angular velocity Omega), so that the background is spherical. Generic perturbations of the rotating star (say parametrized by lambda) are then built on top of the axisymmetric perturbations in Omega. Clearly, any interesting physics requires non-linear perturbations, as at least terms lambda Omega need to be considered. In this paper we analyse the gauge dependence of non-linear perturbations depending on two parameters, derive explicit higher order gauge transformation rules, and define gauge invariance. The formalism is completely general and can be used in different applications of general relativity or any other spacetime theory.