Source author record

Charles R. Evans

Charles R. Evans 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

11works
6topics
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

11 published item(s)

preprint2022arXiv

High-order post-Newtonian expansion of the redshift invariant for eccentric-orbit non-spinning extreme-mass-ratio inspirals

We calculate the eccentricity dependence of the high-order post-Newtonian (PN) series for the generalized redshift invariant $\langle u^t \rangle_τ$ for eccentric-orbit extreme-mass-ratio inspirals on a Schwarzschild background. These results are calculated within first-order black hole perturbation theory (BHPT) using Regge-Wheeler-Zerilli (RWZ) gauge. Our \textsc{Mathematica} code is based on a familiar procedure, using PN expansion of the Mano-Suzuki-Takasugi (MST) analytic function formalism for $l$ modes up to a certain maximum and then using a direct general-$l$ PN expansion of the RWZ equation for arbitrarily high $l$. We calculate dual expansions in PN order and in powers of eccentricity, reaching 10PN relative order and $e^{20}$. Detailed knowledge of the eccentricity expansion at each PN order allows us to find within the eccentricity dependence numerous closed-form expressions and multiple infinite series with known coefficients. We find leading logarithm sequences in the PN expansion of the redshift invariant that reflect a similar behavior in the PN expansion of the energy flux to infinity. A set of flux terms and special functions that appear in the energy flux, like the Peters-Mathews flux itself, are shown to reappear in the redshift PN expansion.

preprint2022arXiv

Post-Newtonian expansion of the spin-precession invariant for eccentric-orbit non-spinning extreme-mass-ratio inspirals to 9PN and $e^{16}$

We calculate the eccentricity dependence of the high-order post-Newtonian (PN) expansion of the spin-precession invariant $ψ$ 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 $ψ$ 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 $ψ$ 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$.

preprint2020arXiv

Determination of new coefficients in the angular momentum and energy fluxes at infinity to 9PN for eccentric Schwarzschild extreme-mass-ratio inspirals using mode-by-mode fitting

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.

preprint2016arXiv

Eccentric-orbit EMRI gravitational wave energy fluxes to 7PN order

We present new results through 7PN order on the energy flux from eccentric extreme-mass-ratio binaries. The black hole perturbation calculations are made at very high accuracy (200 decimal places) using a Mathematica code based on the Mano-Suzuki-Takasugi (MST) analytic function expansion formalism. All published coefficients in the expansion through 3PN order are confirmed and new analytic and numeric terms are found to high order in $e^2$ at orders between 3.5PN and 7PN. We also show original work in finding (nearly) arbitrarily accurate expansions for hereditary terms at 1.5PN, 2.5PN, and 3PN orders. We fit to a model where at each PN order an eccentricity singular function is factored out, improving substantially the fit even as $e \to 1$.

preprint2016arXiv

Highly eccentric inspirals into a black hole

We model the inspiral of a compact stellar-mass object into a massive nonrotating black hole including all dissipative and conservative first-order-in-the-mass-ratio effects on the orbital motion. The techniques we develop allow inspirals with initial eccentricities as high as $e\sim0.8$ and initial separations as large as $p\sim 50$ to be evolved through many thousands of orbits up to the onset of the plunge into the black hole. The inspiral is computed using an osculating elements scheme driven by a hybridized self-force model, which combines Lorenz-gauge self-force results with highly accurate flux data from a Regge-Wheeler-Zerilli code. The high accuracy of our hybrid self-force model allows the orbital phase of the inspirals to be tracked to within $\sim0.1$ radians or better. The difference between self-force models and inspirals computed in the radiative approximation is quantified.

preprint2016arXiv

Mirror Reflections of a Black Hole

An exact correspondence between a black hole and an accelerating mirror is demonstrated. It is shown that for a massless minimally coupled scalar field the same Bogolubov coefficients connecting the "in" and "out" states occur for a (1+1)D flat spacetime with a particular perfectly reflecting accelerating boundary trajectory and a (1+1)D curved spacetime in which a null shell collapses to form a black hole. Generalization of the latter to the (3+1)D case is discussed. The spectral dynamics is computed in both (1+1)-dimensional spacetimes along with the energy flux in the spacetime with a mirror. It is shown that the approach to equilibrium is monotonic, asymmetric in terms of the rate, and there is a specific time which characterizes the system when it is the most out-of-equilibrium.

preprint2015arXiv

Fast spectral source integration in black hole perturbation calculations

This paper presents a new technique for achieving spectral accuracy and fast computational performance in a class of black hole perturbation and gravitational self-force calculations involving extreme mass ratios and generic orbits. Called \emph{spectral source integration} (SSI), this method should see widespread future use in problems that entail (i) point-particle description of the small compact object, (ii) frequency domain decomposition, and (iii) use of the background eccentric geodesic motion. Frequency domain approaches are widely used in both perturbation theory flux-balance calculations and in local gravitational self-force calculations. Recent self-force calculations in Lorenz gauge, using the frequency domain and method of extended homogeneous solutions, have been able to accurately reach eccentricities as high as $e \simeq 0.7$. We show here SSI successfully applied to Lorenz gauge. In a double precision Lorenz gauge code, SSI enhances the accuracy of results and makes a factor of three improvement in the overall speed. The primary initial application of SSI--for us its \emph{raison d'être}--is in an arbitrary precision \emph{Mathematica} code that computes perturbations of eccentric orbits in the Regge-Wheeler gauge to extraordinarily high accuracy (e.g., 200 decimal places). These high accuracy eccentric orbit calculations would not be possible without the exponential convergence of SSI. We believe the method will extend to work for inspirals on Kerr, and will be the subject of a later publication. SSI borrows concepts from discrete-time signal processing and is used to calculate the mode normalization coefficients in perturbation theory via sums over modest numbers of points around an orbit. A variant of the idea is used to obtain spectral accuracy in solution of the geodesic orbital motion.

preprint2013arXiv

Metric perturbations from eccentric orbits on a Schwarzschild black hole: I. Odd-parity Regge-Wheeler to Lorenz gauge transformation and two new methods to circumvent the Gibbs phenomenon

We calculate the odd-parity, radiative ($\ell \ge 2$) parts of the metric perturbation in Lorenz gauge caused by a small compact object in eccentric orbit about a Schwarzschild black hole. The Lorenz gauge solution is found via gauge transformation from a corresponding one in Regge-Wheeler gauge. Like the Regge-Wheeler gauge solution itself, the gauge generator is computed in the frequency domain and transferred to the time domain. The wave equation for the gauge generator has a source with a compact, moving delta-function term and a discontinuous non-compact term. The former term allows the method of extended homogeneous solutions to be applied (which circumvents the Gibbs phenomenon). The latter has required the development of new means to use frequency domain methods and yet be able to transfer to the time domain while avoiding Gibbs problems. Two new methods are developed to achieve this: a partial annihilator method and a method of extended particular solutions. We detail these methods and show their application in calculating the odd-parity gauge generator and Lorenz gauge metric perturbations. A subsequent paper will apply these methods to the harder task of computing the even-parity parts of the gauge generator.

preprint2013arXiv

Relativistic effects in the tidal interaction between a white dwarf and a massive black hole in Fermi normal coordinates

We consider tidal encounters between a white dwarf and an intermediate mass black hole. Both weak encounters and those at the threshold of disruption are modeled. The numerical code combines mesh-based hydrodynamics, a spectral method solution of the self-gravity, and a general relativistic Fermi normal coordinate system that follows the star and debris. Fermi normal coordinates provide an expansion of the black hole tidal field that includes quadrupole and higher multipole moments and relativistic corrections. We compute the mass loss from the white dwarf that occurs in weak tidal encounters. Secondly, we compute carefully the energy deposition onto the star, examining the effects of nonradial and radial mode excitation, surface layer heating, mass loss, and relativistic orbital motion. We find evidence of a slight relativistic suppression in tidal energy transfer. Tidal energy deposition is compared to orbital energy loss due to gravitational bremsstrahlung and the combined losses are used to estimate tidal capture orbits. Heating and partial mass stripping will lead to an expansion of the white dwarf, making it easier for the star to be tidally disrupted on the next passage. Finally, we examine angular momentum deposition. By including the octupole tide, we are able for the first time to calculate deflection of the center of mass of the star and debris. With this observed deflection, and taking into account orbital relativistic effects, we compute directly the change in orbital angular momentum and show its balance with computed spin angular momentum deposition.

preprint2013arXiv

Time Dependence of Particle Creation from Accelerating Mirrors

Particle production due to a quantized, massless, minimally coupled scalar field in two-dimensional flat spacetime with an accelerating mirror is investigated, with a focus on the time dependence of the process. We analyze first the classes of trajectories previously investigated by Carlitz and Willey and by Walker and Davies. We then analyze four new classes of trajectories, all of which can be expressed analytically and for which several ancillary properties can be derived analytically. The time dependence is investigated through the use of wave packets for the modes of the quantized field that are in the out vacuum state. It is shown for most of the trajectories studied that good time resolution of the particle production process can be obtained.

preprint2010arXiv

Gravitational perturbations and metric reconstruction: Method of extended homogeneous solutions applied to eccentric orbits on a Schwarzschild black hole

We calculate the gravitational perturbations produced by a small mass in eccentric orbit about a much more massive Schwarzschild black hole and use the numerically computed perturbations to solve for the metric. The calculations are initially made in the frequency domain and provide Fourier-harmonic modes for the gauge-invariant master functions that satisfy inhomogeneous versions of the Regge-Wheeler and Zerilli equations. These gravitational master equations have specific singular sources containing both delta function and derivative-of-delta function terms. We demonstrate in this paper successful application of the method of extended homogeneous solutions, developed recently by Barack, Ori, and Sago, to handle source terms of this type. The method allows transformation back to the time domain, with exponential convergence of the partial mode sums that represent the field. This rapid convergence holds even in the region of $r$ traversed by the point mass and includes the time-dependent location of the point mass itself. We present numerical results of mode calculations for certain orbital parameters, including highly accurate energy and angular momentum fluxes at infinity and at the black hole event horizon. We then address the issue of reconstructing the metric perturbation amplitudes from the master functions, the latter being weak solutions of a particular form to the wave equations. The spherical harmonic amplitudes that represent the metric in Regge-Wheeler gauge can themselves be viewed as weak solutions. They are in general a combination of (1) two differentiable solutions that adjoin at the instantaneous location of the point mass (a result that has order of continuity $C^{-1}$ typically) and (2) (in some cases) a delta function distribution term with a computable time-dependent amplitude.