Search

We calculate the eccentricity dependence of the high-order post-Newtonian (PN) expansion of the spin-precession invariant $\psi$ for eccentric-orbit extreme-mass-ratio inspirals with a Schwarzschild primary. The series is calculated in first-order black hole perturbation theory through direct analytic expansion of solutions in the Regge-Wheeler-Zerilli formalism, using a code written in \textsc{Mathematica}. Modes with small values of $l$ are found via the Mano-Suzuki-Takasugi (MST) analytic function expansion formalism for solutions to the Regge-Wheeler equation. Large-$l$ solutions are found by applying a PN expansion ansatz to the Regge-Wheeler equation. Previous work has given $\psi$ to 9.5PN order and to order $e^2$ (i.e., the near circular orbit limit). We calculate the expansion to 9PN but to $e^{16}$ in eccentricity. It proves possible to find a few terms that have closed-form expressions, all of which are associated with logarithmic terms in the PN expansion. We also compare the numerical evaluation of our PN expansion to prior numerical calculations of $\psi$ in close orbits to assess its radius of convergence. We find that the series is not as rapidly convergent as the one for the redshift invariant at $r \simeq 10M$ but still yielding $\sim 1\%$ accuracy for eccentricities $e \lesssim 0.25$., Comment: 23 pages, 2 figures

In a recent paper we showed that for eccentric-orbit extreme-mass-ratio inspirals the analytic forms of the leading-logarithm energy and angular momentum post-Newtonian (PN) flux terms (radiated to infinity) can, to arbitrary PN order, be determined by sums over the Fourier spectrum of the Newtonian quadrupole moment. We further showed that an essential part of the eccentricity dependence of the related subleading-logarithm PN sequences, at lowest order in the symmetric mass ratio $\nu$, stems as well from the Newtonian quadrupole moment. Once that part is factored out, the remaining eccentricity dependence is more easily determined by black hole perturbation theory. In this paper we show how the sequences that are the 1PN corrections to the entire leading-logarithm series, namely terms that appear at PN orders $x^{3k+1} \log^k(x)$ and $x^{3k+5/2} \log^k(x)$ (for PN compactness parameter $x$ and integers $k\ge 0$), at lowest order in $\nu$, are determined by the Fourier spectra of the Newtonian mass octupole, Newtonian current quadrupole, and 1PN part of the mass quadrupole moments. We also develop a conjectured (but plausible) form for 1PN correction to the leading logs at second order in $\nu$. Further, in analogy to the first paper, we show that these same source multipole moments also yield nontrivial parts of the 1PN correction to the subleading-logarithm series, and that the remaining eccentricity dependence (at lowest order in $\nu$) can then more easily be determined using black hole perturbation theory. We use this method to determine the entire analytic eccentricity dependence of the perturbative (i.e., lowest order in $\nu$) 4PN non-log terms, $\mathcal{R}_4(e_t)$ and $\mathcal{Z}_4(e_t)$, for energy and angular momentum respectively., Comment: 36 pages, 1 figure

We present an extension of work in an earlier paper showing high precision comparisons between black hole perturbation theory and post-Newtonian (PN) theory in their region of overlapping validity for bound, eccentric-orbit, Schwarzschild extreme-mass-ratio inspirals. As before we apply a numerical fitting scheme to extract eccentricity coefficients in the PN expansion of the gravitational wave fluxes, which are then converted to exact analytic form using an integer-relation algorithm. In this work, however, we fit to individual $lmn$ modes to exploit simplifying factorizations that lie therein. Since the previous paper focused solely on the energy flux, here we concentrate initially on analyzing the angular momentum flux to infinity. A first step involves finding convenient forms for hereditary contributions to the flux at low-PN order, analogous to similar terms worked out previously for the energy flux. We then apply the upgraded techniques to find new PN terms through 9PN order and (at many PN orders) to $e^{30}$ in the power series in eccentricity. With the new approach applied to angular momentum fluxes, we return to the energy fluxes at infinity to extend those previous results. Like before, the underlying method uses a \textsc{Mathematica} code based on use of the Mano-Suzuki-Takasugi (MST) function expansion formalism to represent gravitational perturbations and spectral source integration (SSI) to find numerical results at arbitrarily high precision., 36 pages, 1 figure

We develop new high-order results for the post-Newtonian (PN) expansions of the energy and angular momentum fluxes at infinity for eccentric-orbit extreme-mass-ratio inspirals (EMRIs) on a Schwarzschild background. The series are derived through direct expansion of the MST solutions within the RWZ formalism for first-order black hole perturbation theory (BHPT). By utilizing factorization and a few computational simplifications, we are able to compute the fluxes to 19PN, with each PN term calculated as a power series in (Darwin) eccentricity to $e^{10}$. This compares favorably with the numeric fitting approach used in previous work. We also compute PN terms to $e^{20}$ through 10PN. Then, we analyze the convergence properties of the composite energy flux expansion by checking against numeric data for several orbits, both for the full flux and also for the individual 220 mode, with various resummation schemes tried for each. The match between the high-order series and numerical calculations is generally strong, maintaining relative error better than $10^{-5}$ except when $p$ (the semi-latus rectum) is small and $e$ is large. However, the full-flux expansion demonstrates superior fidelity (particularly at high $e$), as it is able to incorporate additional information from PN theory. For the orbit $(p=10, e=1/2)$, the full flux achieves a best error near $10^{-5}$, while the 220 mode exhibits error worse than $1\%$. Finally, we describe a procedure for transforming these expansions to the harmonic gauge of PN theory by analyzing Schwarzschild geodesic motion in harmonic coordinates. This will facilitate future comparisons between BHPT and PN theory., Comment: 29 pages, 3 figures

We present new results on the analytic eccentricity dependence of several sequences of gravitational wave flux terms at high post-Newtonian (PN) order for extreme-mass-ratio inspirals. These sequences are the leading logarithms, which appear at PN orders ${x}^{3k}{\mathrm{log}}^{k}(x)$ and ${x}^{3k+3/2}{\mathrm{log}}^{k}(x)$ for integers $k\ensuremath{\ge}0$ ($x$ is a PN compactness parameter), and the subleading logarithms, which appear at orders ${x}^{3k}{\mathrm{log}}^{k\ensuremath{-}1}(x)$ and ${x}^{3k+3/2}{\mathrm{log}}^{k\ensuremath{-}1}(x)$ ($k\ensuremath{\ge}1$), in both the energy and angular momentum radiated to infinity. For the energy flux leading logarithms, we show that to arbitrarily high PN order, their eccentricity dependence is determined by particular sums over the function $g(n,{e}_{t})$, derived from the Newtonian mass quadrupole moment, that normally gives the spectral content of the Peters-Mathews flux as a function of radial harmonic $n$. An analogous power spectrum $\stackrel{\texttildelow{}}{g}(n,{e}_{t})$ determines the leading logarithms of the angular momentum flux. For subleading logs, the quadrupole power spectra are again shown to play a role, providing a distinguishable part of the eccentricity dependence of these flux terms to high PN order. With the quadrupole contribution understood, the remaining analytic eccentricity dependence of the subleading logs can, in principle, be determined more easily using black hole perturbation theory. We show this procedure in action, deriving the complete analytic structure of the ${x}^{6}\mathrm{log}(x)$ subleading-log term and an analytic expansion of the ${x}^{9/2}$ subleading log to high order in a power series in eccentricity. We discuss how these methods might be extended to other sequences of terms in the PN expansion involving logarithms.