No Arabic abstract
Gravitational wave emission from extreme mass ratio binaries (EMRBs) should be detectable by the joint NASA-ESA LISA project, spurring interest in analytical and numerical methods for investigating EMRBs. We describe a discontinuous Galerkin (dG) method for solving the distributionally forced 1+1 wave equations which arise when modeling EMRBs via the perturbation theory of Schwarzschild blackholes. Despite the presence of jump discontinuities in the relevant polar and axial gravitational master functions, our dG method achieves global spectral accuracy, provided that we know the instantaneous position, velocity, and acceleration of the small particle. Here these variables are known, since we assume that the particle follows a timelike geodesic of the Schwarzschild geometry. We document the results of several numerical experiments testing our method, and in our concluding section discuss the possible inclusion of gravitational self-force effects.
A powerful technique to calculate gravitational radiation from binary systems involves a perturbative expansion: if the masses of the two bodies are very different, the small body is treated as a point particle of mass $m_p$ moving in the gravitational field generated by the large mass $M$, and one keeps only linear terms in the small mass ratio $m_p/M$. This technique usually yields finite answers, which are often in good agreement with fully nonlinear numerical relativity results, even when extrapolated to nearly comparable mass ratios. Here we study two situations in which the point-particle approximation yields a divergent result: the instantaneous flux emitted by a small body as it orbits the light ring of a black hole, and the total energy absorbed by the horizon when a small body plunges into a black hole. By integrating the Teukolsky (or Zerilli/Regge-Wheeler) equations in the frequency and time domains we show that both of these quantities diverge. We find that these divergences are an artifact of the point-particle idealization, and are able to interpret and regularize this behavior by introducing a finite size for the point particle. These divergences do not play a role in black-hole imaging, e.g. by the Event Horizon Telescope.
We describe the hyperboloidal compactification for Teukolsky equations in Kerr spacetime. We include null infinity on the numerical grid by attaching a hyperboloidal layer to a compact domain surrounding the rotating black hole and the orbit of an inspiralling point particle. This technique allows us to study, for the first time, gravitational waveforms from large- and extreme-mass-ratio inspirals in Kerr spacetime extracted at null infinity. Tests and comparisons of our results with previous calculations establish the accuracy and efficiency of the hyperboloidal layer method.
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.
An extreme mass ratio inspiral takes place when a compact stellar object is inspiraling into a supermassive black hole due to gravitational radiation reaction. Gravitational waves (GWs) from this system can be calculated using the Teukolsky equation (TE). In our case, we compute the asymptotic GW fluxes of a spinning body orbiting a Kerr black hole by solving numerically the TE both in time and frequency domain. Our ultimate goal is to produce GW templates for space-based detectors such as LISA.
Binary-black-hole orbits precess when the black-hole spins are mis-aligned with the binarys orbital angular momentum. The apparently complicated dynamics can in most cases be described as simple precession of the orbital angular momentum about an approximately fixed total angular momentum. However, the imprint of the precession on the observed gravitational-wave signal is yet more complicated, with a non-trivial time-varying dependence on black-hole dynamics, the binarys orientation and the detector polarization. As a result, it is difficult to predict under which conditions precession effects are measurable in gravitational-wave observations, and their impact on both signal detection and source characterization. We show that the observed waveform can be simplified by decomposing it as a power series in a new precession parameter $b = tan(beta/2)$, where $beta$ is the opening angle between the orbital and total angular momenta. The power series is made up of five harmonics, with frequencies that differ by the binarys precession frequency, and individually do not exhibit amplitude and phase modulations. In many cases, the waveform can be well approximated by the two leading harmonics. In this approximation we are able to obtain a simple picture of precession as caused by the beating of two waveforms of similar frequency. This enables us to identify regions of the parameter space where precession is likely to have an observable effect on the waveform, and to propose a new approach to searching for signals from precessing binaries, based upon the two-harmonic approximation.