Search

In this work we investigate the presence of integrable hidden structures in the dynamics of perturbed non-rotating black holes (BHs). This can also be considered as a first step in a wider program of an effective identification of ``slow'' and ``fast'' degrees of freedom (DoFs) in the (binary) BH dynamics, following a wave-mean flow perspective. The slow DoFs would be associated with a nonlinear integrable dynamics, on which the fast ones propagate following an effective linear dynamics. BH perturbation theory offers a natural ground to test these properties. Indeed, the decoupling of Einstein equations into wave master equations with a potential provides an instance of such splitting into (frozen) slow DoFs (background potential) over which the linear dynamics of the fast ones (perturbation master functions) evolve. It has been recently shown that these wave equations possess an infinite number of symmetries that correspond to the flow of the infinite hierarchy of Korteweg-de Vries (KdV) equations. Starting from these results, we systematically investigate the presence of integrable structures in BH perturbation theory. We first study them in Cauchy slices and then extend the analysis to hyperboloidal foliations. This second step introduces a splitting of the master equation into bulk and boundary contributions, unveiling an underlying structural relation with the slow and fast DoFs. This insight represents a first step to establish the integrable structures associated to the slow DoFs as bulk symmetries of the dynamics of perturbed BHs., Comment: 23 pages, comments are welcome

LISA (Laser Interferometer Space Antenna) is a space mission that has recently entered the implementation phase. The main goal of LISA is to survey, for the first time, the observable universe by detecting low-frequency gravitational waves (GWs). One of the main GW sources that LISA will detect is the coalescence and merger of Massive Black Hole Binaries (MBHBs). In this work, we explore the application of Sequential Neural Likelihood, a simulation-based inference algorithm, to the detection and parameter estimation of MBHBs signals in the LISA data. Instead of sampling from the conventional likelihood function, which typically relies on the assumptions of Gaussianity and stationarity of the noise and requires a forward simulation for each evaluation, this method constructs a surrogate likelihood that is ultimately described by a neural network trained from a dataset of simulations (noise + MBHB). Since the likelihood is independent of the priors, we can iteratively train models that target specific observations in a fraction of the time and computational cost that other traditional and machine learning-based strategies would require. One key aspect of the method is the compression of the data. Because of the iterative nature of the method, we are able to train models to obtain qualitatively similar posteriors with less than 2% of the simulator calls that Markov Chain Monte Carlo methods would require. We compare these posteriors with those obtained from standard techniques and discuss the differences that appear, some of which can be identified as caused by the loss of information from the data compression we use. We also discuss future improvements to this method., Comment: RevTeX 4.2, 28 pages, 17 figures, 20 plots. Comments welcome

Perturbation theory of vacuum spherically-symmetric spacetimes (including the cosmological constant) has greatly contributed to the understanding of black holes, relativistic compact stars and even inhomogeneous cosmological models. The perturbative equations can be decoupled in terms of (gauge-invariant) master functions satisfying $1+1$ wave equations. In this work, building on previous work on the structure of the space of master functions and equations, we study the reconstruction of the metric perturbations in terms of the master functions. To that end, we consider the general situation in which the perturbations are driven by an arbitrary energy-momentum tensor. Then, we perform the metric reconstruction in a completely general perturbative gauge. In doing so, we investigate the role of Darboux transformations and Darboux covariance, responsible for the isospectrality between odd and even parity in the absence of matter sources and also of the physical equivalence between the descriptions based on all the possible master equations. We also show that the metric reconstruction can be carried out in terms of any of the possible master functions and that the expressions admit an explicitly covariant form., Comment: 27 pages, comments are welcome

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

Extreme-mass-ratio inspirals consist of binary systems of compact objects, with orders of magnitude differences in their masses, in the regime where the dynamics are driven by gravitational wave emission. The unique nature of extreme-mass-ratio inspirals facilitates the exploration of diverse and stimulating tests related to various aspects of black holes and General Relativity. These tests encompass investigations into the spacetime geometry, the dissipation of the binary system energy and angular momentum, and the impact of the intrinsic non-linear effects of the gravitational interaction. This review accounts for the manifold opportunities for testing General Relativity, ranging from examinations of black hole properties to cosmological and multimessenger tests., Comment: V2. Updated to match published version. V1: 67+19 pages, 2 figures. Review Chapter for the book "Recent Progress on Gravity Tests" (Eds. C. Bambi and A. C\'ardenas-Avenda\~no, Springer Singapore)

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome., Comment: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

This document presents a summary of the 2023 Terrestrial Very-Long-Baseline Atom Interferometry Workshop hosted by CERN. The workshop brought together experts from around the world to discuss the exciting developments in large-scale atom interferometer (AI) prototypes and their potential for detecting ultralight dark matter and gravitational waves. The primary objective of the workshop was to lay the groundwork for an international TVLBAI proto-collaboration. This collaboration aims to unite researchers from different institutions to strategize and secure funding for terrestrial large-scale AI projects. The ultimate goal is to create a roadmap detailing the design and technology choices for one or more km-scale detectors, which will be operational in the mid-2030s. The key sections of this report present the physics case and technical challenges, together with a comprehensive overview of the discussions at the workshop together with the main conclusions., Comment: Summary of the Terrestrial Very-Long-Baseline Atom Interferometry Workshop held at CERN: https://indico.cern.ch/event/1208783/

The waveform of a binary black hole coalescence appears to be both simple and universal. In this essay we argue that the dynamics should admit a separation into 'fast and slow' degrees of freedom, such that the latter are described by an integrable system of equations, accounting for the simplicity and universality of the waveform. Given that Painlev\'e transcendents are a smoking gun of integrable structures, we propose the Painlev\'e-II transcendent as the key structural element threading a hierarchy of asymptotic models aiming at capturing different (effective) layers in the dynamics. Ward's conjecture relating integrable and (anti)self-dual solutions can provide the avenue to encode background binary black hole data in (non-local) twistor structures., Comment: 17 pages, no figures. Essay awarded a honorable mention in the 2023 Gravity Research Foundation essay competition

It has recently been shown that the dynamics of perturbed non-rotating black holes (BHs) admits an infinite number of symmetries that are generated by the flow of the Korteweg-de Vries (KdV) equation. These symmetries lead to an infinite number of conserved quantities that can be obtained as integrals of differential polynomials in the potential appearing in the gauge-invariant master equations describing the BH perturbations, the KdV integrals. These conserved quantities are the same for all the possible potentials, which means that they are invariant under Darboux transformations, and they fully determine the BHs transmission amplitudes, or greybody factors, via a moment problem. In this paper we introduce a new semi-analytical method to obtain the greybody factors associated with BH scattering processes by solving the moment problem using only the KdV integrals. The method is based on the use of Pad\'e approximants and we check it first by comparing with results from the case of a P\"oschl-Teller potential, for which we have analytical expressions for the greybody factors. Then, we apply it to the case of a Schwarzschild BH and compare with results from computations based on the Wentzel-Kramers-Brillouin (WKB) approximation. It turns out that the new method provides accurate results for the BH greybody factors for all frequencies. The method is also computationally very efficient., Comment: 21 pages, 13 figures, 26 plots, picture in Fig.(12) (bottom panel) replaced as it was duplicated

The dynamics of perturbed non-rotating black holes (BHs) can be described in terms of master equations of the wave type with a potential. In the frequency domain, the master equations become time-independent Schr\"odinger equations with no discrete spectrum. It has been recently shown that these wave equations possess an infinite number of symmetries that correspond to the flow of the infinite hierarchy of Korteweg-de Vries (KdV) equations. As a consequence, the infinite set of associated conserved quantities, the KdV integrals, are the same for all the different master equations that we can consider. In this paper we show that the BH scattering reflection and transmission coefficients characterizing the continuous spectrum can be fully determined via a moment problem, in such a way that the KdV integrals provide the momenta of a distribution function depending only on the reflection coefficient. We also discuss the existence and uniqueness of solutions, strategies to solve the moment problem, and finally show the case of the P\"oschl-Teller potential where all the steps can be carried out analytically., Comment: 22 pages, 5 figures

An M-class mission proposal in response to the 2021 call in ESA's science programme with a broad range of objectives in fundamental physics, which include testing the Equivalence Principle and Lorentz Invariance, searching for Ultralight Dark Matter and probing Quantum Mechanics., Comment: As submitted to the M7 call in July 2022, except updated for the recent (Sept. 2022) MICROSCOPE results, and new section 2.5 summarizing the information provided to ESA during the September 2022 audition

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., Comment: Accepted in: Living Reviews in Relativity

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

We summarize the discussions at a virtual Community Workshop on Cold Atoms in Space concerning the status of cold atom technologies, the prospective scientific and societal opportunities offered by their deployment in space, and the developments needed before cold atoms could be operated in space. The cold atom technologies discussed include atomic clocks, quantum gravimeters and accelerometers, and atom interferometers. Prospective applications include metrology, geodesy and measurement of terrestrial mass change due to, e.g., climate change, and fundamental science experiments such as tests of the equivalence principle, searches for dark matter, measurements of gravitational waves and tests of quantum mechanics. We review the current status of cold atom technologies and outline the requirements for their space qualification, including the development paths and the corresponding technical milestones, and identifying possible pathfinder missions to pave the way for missions to exploit the full potential of cold atoms in space. Finally, we present a first draft of a possible road-map for achieving these goals, that we propose for discussion by the interested cold atom, Earth Observation, fundamental physics and other prospective scientific user communities, together with ESA and national space and research funding agencies., Comment: Summary of the Community Workshop on Cold Atoms in Space and corresponding Road-map: https://indico.cern.ch/event/1064855/

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., Comment: 6 pages, revtex 4-2

Perturbation theory of vacuum spherically-symmetric spacetimes is a crucial tool to understand the dynamics of black hole perturbations. Spherical symmetry allows for an expansion of the perturbations in scalar, vector, and tensor harmonics. The resulting perturbative equations are decoupled for modes with different parity and different harmonic numbers. Moreover, for each harmonic and parity, the equations for the perturbations can be decoupled in terms of (gauge-invariant) master functions that satisfy 1+1 wave equations. By working in a completely general perturbative gauge, in this paper we study what is the most general master function that is linear in the metric perturbations and their first-order derivatives and satisfies a wave equation with a potential. The outcome of the study is that for each parity we have two branches of solutions with similar features. One of the branches includes the known results: In the odd-parity case, the most general master function is an arbitrary linear combination of the Regge-Wheeler and the Cunningham-Price-Moncrief master functions whereas in the even-parity case it is an arbitrary linear combination of the Zerilli master function and another master function that is new to our knowledge. The other branch is very different since it includes an infinite collection of potentials which in turn lead to an independent collection master of functions which depend on the potential. The allowed potentials satisfy a non-linear ordinary differential equation. Finally, all the allowed master functions are gauge invariant and can be written in a fully covariant form., Comment: 22 pages, revtex 4-2. (v2: title update) (v3: new conventions for a master function; update of a reference)

We present the scientific motivation for future space tests of the equivalence principle, and in particular the universality of free fall, at the $10^{-17}$ level or better. Two possible mission scenarios, one based on quantum technologies, the other on electrostatic accelerometers, that could reach that goal are briefly discussed., Comment: White paper submitted to the ESA Voyage 2050 long term plan. Partly derived from Altschul et al., Advances in Space Research, 2015, 55, 501 - 524, arXiv:1404.4307

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., Comment: v1: 41 pages, 12 figures. To be submitted for publication, comments welcome; v2: 42 pages. Comments and references added; v3: Matches published version

Since the 2017 Nobel Prize in Physics was awarded for the observation of gravitational waves, it is fair to say that the epoch of gravitational wave astronomy (GWs) has begun. However, a number of interesting sources of GWs can only be observed from space. To demonstrate the feasibility of the Laser Interferometer Space Antenna (LISA), a future gravitational wave observatory in space, the LISA Pathfinder satellite was launched on December, 3rd 2015. Measurements of the spurious forces accelerating an otherwise free-falling test mass, and detailed investigations of the individual subsystems needed to achieve the free-fall, have been conducted throughout the mission. This overview article starts with the purpose and aim of the mission, explains satellite hardware and mission operations and ends with a summary of selected important results and an outlook towards LISA. From the LISA Pathfinder experience, we can conclude that the proposed LISA mission is feasible., Comment: Preprint of an article submitted for consideration in Lepton-Photon 2017 conference proceedings \copyright 2019 copyright World Scientific Publishing Company

The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources for the space-based observatory LISA (Laser Interferometer Space Antenna). LISA observations will enable high-precision measurements of the physical characteristics of these extreme-mass-ratio inspirals (EMRIs): redshifted masses, massive black hole spin and orbital eccentricity can be determined with fractional errors $\sim 10^{-4}$-$10^{-6}$, the luminosity distance with better than $\sim 10\%$ precision, and the sky localization to within a few square degrees. EMRIs will provide valuable information about stellar dynamics in galactic nuclei, as well as precise data about massive black hole populations, including the distribution of masses and spins. They will enable percent-level measurements of the multipolar structure of massive black holes, precisely testing the strong-gravity properties of their spacetimes. EMRIs may also provide cosmographical data regarding the expansion of the Universe if inferred source locations can be correlated with galaxy catalogs., Comment: 13 pages, 1 figure; white paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics); Section III is arXiv-only

The science operations of the LISA Pathfinder mission has demonstrated the feasibility of sub-femto-g free-fall of macroscopic test masses necessary to build a LISA-like gravitational wave observatory in space. While the main focus of interest, i.e. the optical axis or the $x$-axis, has been extensively studied, it is also of interest to evaluate the stability of the spacecraft with respect to all the other degrees of freedom. The current paper is dedicated to such a study, with a focus set on an exhaustive and quantitative evaluation of the imperfections and dynamical effects that impact the stability with respect to its local geodesic. A model of the complete closed-loop system provides a comprehensive understanding of each part of the in-loop coordinates spectra. As will be presented, this model gives very good agreements with LISA Pathfinder flight data. It allows one to identify the physical noise source at the origin and the physical phenomena underlying the couplings. From this, the performances of the stability of the spacecraft, with respect to its geodesic, are extracted as a function of frequency. Close to $1 mHz$, the stability of the spacecraft on the $X_{SC}$, $Y_{SC}$ and $Z_{SC}$ degrees of freedom is shown to be of the order of $5.0\ 10^{-15} m\ s^{-2}/\sqrt{Hz}$ for X and $4.0 \ 10^{-14} m\ s^{-2}/\sqrt{Hz}$ for Y and Z. For the angular degrees of freedom, the values are of the order $3\ 10^{-12} rad\ s^{-2}/\sqrt{Hz}$ for $\Theta_{SC}$ and $3\ 10^{-13} rad\ s^{-2}/\sqrt{Hz}$ for $H_{SC}$ and $\Phi_{SC}$., Comment: 16 pages, 10 figures

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

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., Comment: Revised version. 30 pages, 12 figures

Partial differential equations with distributional sources---in particular, involving (derivatives of) delta distributions---have become increasingly ubiquitous in numerous areas of physics and applied mathematics. It is often of considerable interest to obtain numerical solutions for such equations, but any singular ("particle"-like) source modeling invariably introduces nontrivial computational obstacles. A common method to circumvent these is through some form of delta function approximation procedure on the computational grid; however, this often carries significant limitations on the efficiency of the numerical convergence rates, or sometimes even the resolvability of the problem at all. In this paper, we present an alternative technique for tackling such equations which avoids the singular behavior entirely: the "Particle-without-Particle" method. Previously introduced in the context of the self-force problem in gravitational physics, the idea is to discretize the computational domain into two (or more) disjoint pseudospectral (Chebyshev-Lobatto) grids such that the "particle" is always at the interface between them; thus, one only needs to solve homogeneous equations in each domain, with the source effectively replaced by jump (boundary) conditions thereon. We prove here that this method yields solutions to any linear PDE the source of which is any linear combination of delta distributions and derivatives thereof supported on a one-dimensional subspace of the problem domain. We then implement it to numerically solve a variety of relevant PDEs: hyperbolic (with applications to neuroscience and acoustics), parabolic (with applications to finance), and elliptic. We generically obtain improved convergence rates relative to typical past implementations relying on delta function approximations., Comment: 41 pages, 11 figures; v2: references and clarifications added (mostly in the introduction), matches the published version in Journal of Scientific Computing

One of the key astrophysical sources for the Laser Interferometer Space Antenna (LISA) are the inspirals of stellar-origin compact objects into massive black holes in the centres of galaxies. These extreme-mass-ratio inspirals (EMRIs) have great potential for astrophysics, cosmology and fundamental physics. In this paper we describe the likely numbers and properties of EMRI events that LISA will observe. We present the first results computed for the 2.5 Gm interferometer that was the new baseline mission submitted in January 2017 in response to the ESA L3 mission call. In addition, we attempt to quantify the astrophysical uncertainties in EMRI event rate estimates by considering a range of different models for the astrophysical population. We present both likely event rates and estimates for the precision with which the parameters of the observed sources could be measured. We finish by discussing the implications of these results for science using EMRIs., Comment: 12 pages, 3 figures, 5 tables. These proceedings are a summary of the results presented in arXiv:1703.09722

The space-based Laser Interferometer Space Antenna (LISA) will be able to observe the gravitational-wave signals from systems comprised of a massive black hole and a stellar-mass compact object. These systems are known as extreme-mass-ratio inspirals (EMRIs) and are expected to complete $\sim 10^4$-$10^5$ cycles in band, thus allowing exquisite measurements of their parameters. In this work, we attempt to quantify the astrophysical uncertainties affecting the predictions for the number of EMRIs detectable by LISA, and find that competing astrophysical assumptions produce a variance of about three orders of magnitude in the expected intrinsic EMRI rate. However, we find that irrespective of the astrophysical model, at least a few EMRIs per year should be detectable by the LISA mission, with up to a few thousands per year under the most optimistic astrophysical assumptions. We also investigate the precision with which LISA will be able to extract the parameters of these sources. We find that typical fractional statistical errors with which the intrinsic parameters (redshifted masses, massive black hole spin and orbital eccentricity) can be recovered are $\sim 10^{-6}$-$10^{-4}$. Luminosity distance (which is required to infer true masses) is inferred to about $10\%$ precision and sky position is localized to a few square degrees, while tests of the multipolar structure of the Kerr metric can be performed to percent-level precision or better., Comment: 13 figures, 22 pages; updated to match published version

We extend our previous analysis of holographic heavy ion collisions in non-conformal theories. We provide a detailed description of our numerical code. We study collisions at different energies in gauge theories with different degrees of non-conformality. We compare four relaxation times: the hydrodynamization time (when hydrodynamics becomes applicable), the EoSization time (when the average pressure approaches its equilibrium value), the isotropization time (when the longitudinal and transverse pressures approach each other) and the condensate relaxation time (when the expectation value of a scalar operator approaches its equilibrium value). We find that these processes can occur in several different orderings. In particular, the condensate can remain far from equilibrium even long after the plasma has hydrodynamized and EoSized. We also explore the rapidity distribution of the energy density at hydrodynamization. This is far from boost-invariant and its width decreases as the non-conformality increases. Nevertheless, the velocity field at hydrodynamization is almost exactly boost-invariant regardless of the non-conformality. This result may be used to constrain the initialization of hydrodynamic fields in heavy ion collisions., Comment: 45 pages, 17 figures. Expanded some sections. Matches published version in JHEP

The gravitational waves emitted by binary systems with extreme mass ratios carry unique astrophysical information expected to be probed by the next generation of gravitational wave detectors such as LISA. The detection of these binaries rely on an accurate modeling of the gravitational self-force that drives their orbital evolution. Although the theoretical formalism to compute the self-force has been largely established, the mathematical tools needed to implement it are still under development, and the self-force computation remains an open problem. We present here a frequency-domain implementation of the particle-without-particle (PwP) technique previously developed for the computation of the scalar self-force -- a helpful testbed for the gravitational self-force., Comment: 2 pages, to appear in the proceedings of the 11th International LISA Symposium (Zurich, 5-9 September 2016), in Journal of Physics: Conference Series

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz., Comment: Submitted to ESA on January 13th in response to the call for missions for the L3 slot in the Cosmic Vision Programme

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., Comment: 39 pages, 3 figures; v2: version to appear in Phys. Rev. D, references added

We numerically simulate gravitational shock wave collisions in a holographic model dual to a non-conformal four-dimensional gauge theory. We find two novel effects associated to the non-zero bulk viscosity of the resulting plasma. First, the hydrodynamization time increases. Second, if the bulk viscosity is large enough then the plasma becomes well described by hydrodynamics before the energy density and the average pressure begin to obey the equilibrium equation of state. We discuss implications for the quark-gluon plasma created in heavy ion collision experiments., Comment: 10 pages, 3 figures; published version

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., Comment: 25 pages, 14 figures

We study the equilibrium and near-equilibrium properties of a holographic five-dimensional model consisting of Einstein gravity coupled to a scalar field with a non-trivial potential. The dual four-dimensional gauge theory is not conformal and, at zero temperature, exhibits a renormalisation group flow between two different fixed points. We quantify the deviations from conformality both in terms of thermodynamic observables and in terms of the bulk viscosity of the theory. The ratio of bulk over shear viscosity violates Buchel's bound. We study relaxation of small-amplitude, homogeneous perturbations by computing the quasi-normal modes of the system at zero spatial momentum. In this approximation we identify two different relaxation channels. At high temperatures, the different pressures first become approximately equal to one another, and subsequently this average pressure evolves towards the equilibrium value dictated by the equation of state. At low temperatures, the average pressure first evolves towards the equilibrium pressure, and only later the different pressures become approximately equal to one another., Comment: 22 pages, 29 figures; published version

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)^{\xi}$, where $M_g$ is the mass gap and the exponent $\xi\simeq 0.7$ appears to be universal., Comment: 5 pages, 5 figures. Corrects typos and adds clarifications to match the published version

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., Comment: 30 pages, 23 figures, accepted for publication in Experimental Astronomy

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., Comment: 17 pages, no figures. Proceedings of the LISA Symposium X, to be published at the Journal of Physics

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., Comment: 5 pages, 1 figure

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., Comment: accepted for publication in GRG

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.

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., Comment: accepted for publication in PRD

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., Comment: 10 pages, submitted

The capture of a compact object in a galactic nucleus by a massive black hole (MBH) is the best way to map space and time around it. Compact objects such as stellar black holes on a capture orbit with a very high eccentricity have been wrongly assumed to be lost for the system after an intense burst of radiation, which has been described as a "direct plunge". We prove that these very eccentric capture orbits spend actually a similar number of cycles in a LISA-like detector as those with lower eccentricities if the central MBH is spinning. Although the rates are higher for high-eccentricity EMRIs, the spin also enhances the rates of lower-eccentricity EMRIs. This last kind have received more attention because of the fact that high-eccentricity EMRIs were thought to be direct plunges and thus negligible. On the other hand, recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on these lower-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 statitical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian but also, and for the first time, a geodesic approximation for the relativistic orbits. The existence of the barrier prevents "traditional EMRIs" from approaching the central MBH, but if the central MBH is spinning the rate will be anyway dominated by highly-eccentric extreme-mass ratio inspirals, which insolently ignore the presence of the barrier, because they are driven by two-body relaxation., Comment: Some rephrasing of two sentences in this version. arXiv admin note: substantial text overlap with arXiv:1205.4713

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., Comment: 10 pages, 2 figures. LaTeX, ASP style. Plenary talk presented at the 9th International LISA Symposium, 21-25 May 2012, BnF-Paris, France