Search

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$., Comment: 28 pages, 8 figures

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., Comment: Updated to reflect published version

An exact correspondence is shown between a new moving mirror trajectory in (1+1)D and a spacetime in (1+1)D in which a black hole forms from the collapse of a null shell. It is shown that the Bogolubov coefficients between the "in" and "out" states are identical and the exact Bogolubov coefficients are displayed. Generalization to the (3+1)D black hole case is discussed., Comment: Contribution to MG14 Proceedings, 5 pages, 1 figure

There is an exact correspondence between the simplest solution to Einstein's equations describing the formation of a black hole and a particular moving mirror trajectory. In both cases the Bogolubov coefficients in 1+1 dimensions are identical and can be computed analytically. Particle creation is investigated by using wave packets. The entire particle creation history is computed, incorporating the early-time non-thermal emission due to the formation of the black hole (or the early-time acceleration of the moving mirror) and the evolution to a Planckian spectrum., Comment: Contribution to MG14 Proceedings, 5 pages, 4 figures

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'\^{e}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., Comment: 15 pages, 7 figures

We present an algorithm for calculating the metric perturbations and gravitational self-force for extreme-mass-ratio inspirals (EMRIs) with eccentric orbits. The massive black hole is taken to be Schwarzschild and metric perturbations are computed in Lorenz gauge. The perturbation equations are solved as coupled systems of ordinary differential equations in the frequency domain. Accurate local behavior of the metric is attained through use of the method of extended homogeneous solutions and mode-sum regularization is used to find the self-force. We focus on calculating the self-force with sufficient accuracy to ensure its error contributions to the phase in a long term orbital evolution will be $\delta\Phi \lesssim 10^{-2}$ radians. This requires the orbit-averaged force to have fractional errors $\lesssim 10^{-8}$ and the oscillatory part of the self-force to have errors $\lesssim 10^{-3}$ (a level frequently easily exceeded). Our code meets this error requirement in the oscillatory part, extending the reach to EMRIs with eccentricities of $e \lesssim 0.8$, if augmented by use of fluxes for the orbit-averaged force, or to eccentricities of $e \lesssim 0.5$ when used as a stand-alone code. Further, we demonstrate accurate calculations up to orbital separations of $a \simeq 100 M$, beyond that required for EMRI models and useful for comparison with post-Newtonian theory. Our principal developments include (1) use of fully constrained field equations, (2) discovery of analytic solutions for even-parity static modes, (3) finding a pre-conditioning technique for outer homogeneous solutions, (4) adaptive use of quad-precision and (5) jump conditions to handle near-static modes, and (6) a hybrid scheme for high eccentricities., Comment: Updated to more closely reflect published version

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., Comment: 21 pages, 5 figures

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., Comment: Updated to match published version, one reference added, one figure changed

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., Comment: 17 pages, 9 figures, Updated with one modified figure and minor changes to the text. Added DOI and Journal reference

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., Comment: 25 pages, 5 figures, Updated with minor changes

We consider combining two important methods for constructing quasi-equilibrium initial data for binary black holes: the conformal thin-sandwich formalism and the puncture method. The former seeks to enforce stationarity in the conformal three-metric and the latter attempts to avoid internal boundaries, like minimal surfaces or apparent horizons. We show that these two methods make partially conflicting requirements on the boundary conditions that determine the time slices. In particular, it does not seem possible to construct slices that are quasi-stationary and avoid physical singularities and simultaneously are connected by an everywhere positive lapse function, a condition which must obtain if internal boundaries are to be avoided. Some relaxation of these conflicting requirements may yield a soluble system, but some of the advantages that were sought in combining these approaches will be lost., Comment: 8 pages, LaTeX2e, 2 postscript figures

We have conducted a VLA search for radio pulsars at the positions of 44 nearby OB runaway stars. The observations involved both searching images for point sources of continuum emission and a time series analysis. Our mean flux sensitivity to pulsars slower than 50 ms was 0.2 mJy. No new pulsars were found in the survey. The size of the survey, combined with the high sensitivity of the observations, sets a significant constraint on the probability, $f_p$, of a runaway OB star having an observable pulsar companion. We find $f_p \le 6.5$\% with 95\% confidence, if the general pulsar luminosity function is applicable to OB star pulsar companions. If a pulsar beaming fraction of \onethird\ is assumed, then we estimate that fewer than 20\% of runaway OB stars have neutron star companions, unless pulsed radio emission is frequently obscured by the OB stellar wind. Our result is consistent with the dynamical (or cluster) ejection model for the formation of OB runaways. The supernova ejection model is not ruled out, but is constrained by these observations to allow only a small binary survival fraction, which may be accommodated if neutron stars acquire significant natal kicks. According to Leonard, Hills and Dewey (1994), a 20\% survival fraction corresponds to a 3-d kick velocity of 420 km s$^{-1}$. This value is in close agreement with recent revisions of the pulsar velocity distribution., Comment: Submitted to the Astronomical Journal. 16 pages. Latex uses aaspp4.sty. 3 postscript figures. Address correspondence to Colin Philp (tas@physics.unc.edu). Revision was to replace .ps file with latex file

We observe critical phenomena in spherical collapse of radiation fluid. A sequence of spacetimes $\cal{S}[\eta]$ is numerically computed, containing models ($\eta\ll 1$) that adiabatically disperse and models ($\eta\gg 1$) that form a black hole. Near the critical point ($\eta_c$), evolutions develop a self-similar region within which collapse is balanced by a strong, inward-moving rarefaction wave that holds $m(r)/r$ constant as a function of a self-similar coordinate $\xi$. The self-similar solution is known and we show near-critical evolutions asymptotically approaching it. A critical exponent $\beta \simeq 0.36$ is found for supercritical ($\eta>\eta_c$) models., Comment: 10 pages (LaTeX) (to appear in Phys. Rev. Lett.), TAR-039-UNC