Source author record

Samuel E. Gralla

Samuel E. Gralla 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

29works
9topics
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

29 published item(s)

preprint2022arXiv

How narrow is the M87* ring? I. The choice of closure likelihood function

Event Horizon Telescope (EHT) observations of the core of the galaxy M87 suggest an observational appearance dominated by a ring of approximately 40$μ$as in diameter. The thickness of the ring is less certain: imaging efforts constrained it to be less than half the diameter (consistent with an imaging resolution of 20$μ$as), while visibility-domain modeling suggested a variety of fractional widths, including as low as $10\%$ on some days. The fractional width is very interesting as it has the potential to discriminate between different astrophysical scenarios for the source; in fact, the $10$--$20\%$ range is so narrow as to be in tension with theoretical expectations. In the first of a series of papers on the width of the observed ring, we reproduce a subset of EHT visibility-domain modeling results and we explore whether alternative data analysis methods might favor thicker rings. We point out that the closure phase (and closure amplitude) likelihood function is not independent of residual station gain amplitudes, even at high signal-to-noise, and explore two approximations of practical interest: one standard in the field (and employed by the EHT collaboration), and a new one that we propose. Analyzing the public data, we find that the new likelihood approximation prefers somewhat thicker rings, more in line with theoretical expectations. Further analysis is needed, however, to determine which approximation is better for the EHT data.

preprint2020arXiv

Can the EHT M87 results be used to test general relativity?

No. All theoretical predictions for the observational appearance of an accreting supermassive black hole, as measured interferometrically by a sparse Earth-sized array at current observation frequencies, are sensitive to many untested assumptions about accretion flow and emission physics. There is no way to distinguish a violation of general relativity (GR) from the much more likely scenario that the relevant "gastrophysical" assumptions simply do not hold. Tests of GR will become possible with longer interferometric baselines (likely requiring a space mission) that reach the resolution where astrophysics-independent predictions of the theory become observable.

preprint2020arXiv

Consistent Blandford-Znajek Expansion

The Blandford-Znajek mechanism is the continuous extraction of energy from a rotating black hole via plasma currents flowing on magnetic field lines threading the horizon. In the discovery paper, Blandford and Znajek demonstrated the mechanism by solving the equations of force-free electrodynamics in a perturbative expansion valid at small black hole spin. Attempts to extend this perturbation analysis to higher order have encountered inconsistencies.We overcome this problem using the method of matched asymptotic expansions, taking care to resolve all of the singular surfaces (light surfaces) in the problem. Working with the monopole field configuration, we show explicitly how the inconsistencies are resolved in this framework and calculate the field configuration to one order higher than previously known. However, there is no correction to the energy extraction rate at this order. These results confirm the basic consistency of the split monopole at small spin and lay a foundation for further perturbative studies of the Blandford-Znajek mechanism.

preprint2020arXiv

Horizon Instability of the Extremal BTZ Black Hole

We study real-time propagation of a massive scalar field on the extremal BTZ black hole spacetime, focusing on the Aretakis instability of the event horizon. We obtain a simple time-domain expression for the $\textrm{AdS}_3$ retarded Green function with Dirichlet boundary conditions and construct the corresponding time-domain BTZ retarded Green function using the method of images. The field decays at different rates on and off the horizon, indicating that transverse derivatives grow with time on the horizon (Aretakis instability). We solve the null geodesic equation in full generality and show that the instability is associated with a class of null geodesics that orbit near the event horizon arbitrarily many times before falling in. In an appendix we also treat the problem in the frequency domain, finding consistency between the methods.

preprint2020arXiv

Lensing by Kerr Black Holes

Interpreting horizon-scale observations of astrophysical black holes demands a general understanding of null geodesics in the Kerr spacetime. These may be divided into two classes: "direct" rays that primarily determine the observational appearance of a given source, and highly bent rays that produce a nested sequence of exponentially demagnified images of the main emission: the so-called "photon ring". We develop heuristics that characterize the direct rays and study the highly bent geodesics analytically. We define three critical parameters $γ$, $δ$, and $τ$ that respectively control the demagnification, rotation, and time delay of successive images of the source, thereby providing an analytic theory of the photon ring. These observable parameters encode universal effects of general relativity, independent of the details of the emitting matter.

preprint2020arXiv

Measuring the shape of a black hole photon ring

General relativity predicts that gravitational lensing near black holes will produce narrow "photon rings" on images. Building on recent work of Johnson, Lupsasca et al. focusing on circular rings, I calculate the long-baseline interferometric signature of a narrow feature of arbitrary shape and intensity. The shape information is contained in the oscillation frequencies of the complex visibility as a function of baseline length, for each baseline angle. These results raise the possibility of measuring the detailed shape of a photon ring and comparing to the precise predictions of general relativity.

preprint2016arXiv

Electromagnetic Jets from Stars and Black Holes

We present analytic force-free solutions modeling rotating stars and black holes immersed in the magnetic field of a thin disk that terminates at an inner radius. The solutions are exact in flat spacetime and approximate in Kerr spacetime. The compact object produces a conical jet whose properties carry information about its nature. For example, the jet from a star is surrounded by a current sheet, while that of a black hole is smooth. We compute an effective resistance in each case and compare to the canonical values used in circuit models of energy extraction. These solutions illustrate all of the basic features of the Blandford-Znajek process for energy extraction and jet formation in a clean setting.

preprint2016arXiv

Horizon Instability of Extremal Kerr Black Holes: Nonaxisymmetric Modes and Enhanced Growth Rate

We show that the horizon instability of the extremal Kerr black hole is associated with a singular branch point in the Green function at the superradiant bound frequency. We study generic initial data supported away from the horizon and find an enhanced growth rate due to nonaxisymmetric modes. The growth is controlled by the conformal weight of each mode. We speculate on connections to near-extremal black holes and holographic duality.

preprint2016arXiv

Inspiral into Gargantua

We model the inspiral of a compact object into a more massive black hole rotating very near the theoretical maximum. We find that once the body enters the near-horizon regime the gravitational radiation is characterized by a constant frequency, equal to (twice) the horizon frequency, with an exponentially damped profile. This contrasts with the usual "chirping" behavior and, if detected, would constitute a "smoking gun" for a near-extremal black hole in nature.

preprint2016arXiv

Near-horizon Kerr Magnetosphere

We exploit the near-horizon conformal symmetry of rapidly spinning black holes to determine universal properties of their magnetospheres. Analytic expressions are derived for the limiting form of the magnetosphere in the near-horizon region. The symmetry is shown to imply that the black hole Meissner effect holds for free Maxwell fields but is generically violated for force-free fields. We further show that in the extremal limit, near-horizon plasma particles are infinitely boosted relative to accretion flow. Active galactic nuclei powered by rapidly spinning black holes are therefore natural sites for high-energy particle collisions.

preprint2016arXiv

Plasma Waves and Jets from Moving Conductors

We consider force-free plasma waves launched by the motion of conducting material through a magnetic field. We develop a spacetime-covariant formalism for perturbations of a uniform magnetic field and show how the transverse motion of a conducting fluid acts as a source. We show that fast-mode waves are sourced by the compressibility of the fluid, with incompressible fluids launching a pure-Alfvén outflow. Remarkably, this outflow can be written down in closed form for an arbitrary time-dependent, nonaxisymmetric incompressible flow. The instantaneous flow velocity is imprinted on the magnetic field and transmitted away at the speed of light, carrying detailed information about the conducting source at the time of emission. These results can be applied to transients in pulsar outflows and to jets from neutron stars orbiting in the magnetosphere of another compact object. We discuss jets from moving conductors in some detail.

preprint2016arXiv

Pulsar Magnetospheres: Beyond the Flat Spacetime Dipole

Most studies of the pulsar magnetosphere have assumed a pure magnetic dipole in flat spacetime. However, recent work suggests that the effects of general relativity are in fact of vital importance and that realistic pulsar magnetic fields will have a significant nondipolar component. We introduce a general analytical method for studying the axisymmetric force-free magnetosphere of a slowly-rotating star of arbitrary magnetic field, mass, radius and moment of inertia, including all the effects of general relativity. We confirm that spacelike current is generically present in the polar caps (suggesting a pair production region), irrespective of the stellar magnetic field. We show that general relativity introduces a ~60% correction to the formula for the dipolar component of the surface magnetic field inferred from spindown. Finally, we show that the location and shape of the polar caps can be modified dramatically by even modestly strong higher moments. This can affect emission processes occurring near the star and may help explain the modified beam characteristics of millisecond pulsars.

preprint2016arXiv

Transient Instability of Rapidly Rotating Black Holes

We analytically study the linear response of a near-extremal Kerr black hole to external scalar, electromagnetic, and gravitational field perturbations. We show that the energy density, electromagnetic field strength, and tidal force experienced by infalling observers exhibit transient growth near the horizon. The growth lasts arbitrarily long in the extremal limit, reproducing the horizon instability of extremal Kerr. We explain these results in terms of near-horizon geometry and discuss potential astrophysical implications.

preprint2015arXiv

Nonaxisymmetric Poynting Jets

The relativistic plasma jets from a misaligned black hole-accretion disk system will not be axially symmetric. Here we analyze nonaxisymmetric, stationary, translation invariant jets in the force-free approximation where the field energy dominates the particle energy. We derive a stream equation for these configurations involving the flux function $ψ$ for the transverse magnetic field, the linear velocity $v(ψ)$ of field lines along the jet, and the longitudinal magnetic field $B_z(ψ)$. The equations can be completely solved when $|v|=1$, and when $|v|<1$ the problem can be reduced to the pure magnetic case $v=0$ by a "field line dependent boost". We also find a large class of nonaxisymmetric solutions with arbitrary dependence on the longitudinal retarded time. A subclass of these constitutes a novel type of jet that has vanishing electromagnetic pressure $\tfrac{1}{2}(B^2-E^2)$ and requires no external pressure for confinement. We prove that such self-confinement is impossible when $B^2>E^2$. Finally, we write down specific solutions approximating numerical results for the nonaxisymmetric jet produced by a spinning black hole in an external, misaligned magnetic field.

preprint2015arXiv

Note on Bunching of Field Lines in Black Hole Magnetospheres

Numerical simulations of Blandford-Znajek energy extraction at high spin have revealed that field lines tend to bunch near the poles of the event horizon. We show that this behavior can be derived analytically from the assumption of fixed functional dependence of current and field line rotation on magnetic flux. The argument relies crucially on the existence of the Znajek condition, which offers non-trivial information about the fields on the horizon without requiring a full force-free solution. We also provide some new analytic expressions for the parabolic field configuration.

preprint2015arXiv

Particle on the Innermost Stable Circular Orbit of a Rapidly Spinning Black Hole

We compute the radiation emitted by a particle on the innermost stable circular orbit of a rapidly spinning black hole both (a) analytically, working to leading order in the deviation from extremality and (b) numerically, with a new high-precision Teukolsky code. We find excellent agreement between the two methods. We confirm previous estimates of the overall scaling of the power radiated, but show that there are also small oscillations all the way to extremality. Furthermore, we reveal an intricate mode-by-mode structure in the flux to infinity, with only certain modes having the dominant scaling. The scaling of each mode is controlled by its conformal weight, a quantity that arises naturally in the representation theory of the enhanced near-horizon symmetry group. We find relationships to previous work on particles orbiting in precisely extreme Kerr, including detailed agreement of quantities computed here with conformal field theory calculations performed in the context of the Kerr/CFT correspondence.

preprint2015arXiv

QED Plasma and Magnetars

Magnetars are surrounded by diffuse plasma in magnetic field strengths well above the quantum electrodynamic critical value. We derive equations of "quantum force-free electrodynamics" for this plasma using effective field theory arguments. We argue that quantum effects do not modify the large scale structure of the magnetosphere, and in particular that the spin-down rate does not deviate significantly from the classical result. We provide definite evolution equations that can be used to explore potentially important small-scale corrections, such as shock formation, which has been proposed as a mechanism for both burst and quiescent emission from magnetars.

preprint2014arXiv

On the Magnetosphere of an Accelerated Pulsar

We report on a remarkable class of exact solutions to force-free electrodynamics that has four-current along the light cones of an arbitrary timelike worldline in flat spacetime. No symmetry is assumed, and the solutions are given in terms of a free function of three variables. The field configuration should describe the outer magnetosphere of a pulsar moving on the worldline. The power radiated is the sum of an acceleration (Larmor-type) term and a pulsar-type term.

preprint2013arXiv

Exact Solutions to Force-Free Electrodynamics in Black Hole Backgrounds

A shared property of several of the known exact solutions to the equations of force-free electrodynamics is that their charge-current four-vector is \textit{null}. We examine the general properties of null-current solutions and then focus on the principal congruences of the Kerr black hole spacetime. We obtain a large class of exact solutions, which are in general time-dependent and non-axisymmetric. These solutions include waves that, surprisingly, propagate without scattering on the curvature of the black hole's background. They may be understood as generalizations to Robinson's solutions to vacuum electrodynamics associated with a shear-free congruence of null geodesics. When stationary and axisymmetric, our solutions reduce to those of Menon and Dermer, the only previously known solutions in Kerr. In Kerr, all of our solutions have null electromagnetic fields ($\vec{E} \cdot \vec{B} = 0$ and $E^2=B^2$). However, in Schwarzschild or flat spacetime there is freedom to add a magnetic monopole field, making the solutions magnetically dominated ($B^2>E^2$). This freedom may be used to reproduce the various flat-spacetime and Schwarzschild-spacetime (split) monopole solutions available in the literature (due to Michel and later authors), and to obtain a large class of time-dependent, non-axisymmetric generalizations. These generalizations may be used to model the magnetosphere of a conducting star that rotates with arbitrary prescribed time-dependent rotation axis and speed. We thus significantly enlarge the class of known exact solutions, while organizing and unifying previously discovered solutions in terms of their null structure.

preprint2013arXiv

Hidden Momentum and Black Hole Kicks

A stationary magnetic dipole immersed in an electric field carries "hidden" mechanical momentum. However, the fate of this momentum if the fields are turned off is unclear. We consider a charge-and-dipole hidden momentum configuration, and turn off the fields by collapsing a null shell onto the system, forming a black hole. In numerical calculations we find that the black hole receives a kick corresponding to 0.1% of the initial stored momentum. When extrapolated to apply to purely gravitational phenomena, this efficiency suggests a role for the hidden momentum kick mechanism in generating the binary black hole "superkicks" observed in numerical simulations of Einstein's equation.

preprint2013arXiv

Thermodynamics of a Black Hole with Moon

For a Kerr black hole perturbed by a particle on the "corotating" circular orbit (angular velocity equal to that of the unperturbed event horizon), the black hole remains in equilibrium in the sense that the perturbed event horizon is a Killing horizon of the helical Killing field. The associated surface gravity is constant over the horizon and should correspond to the physical Hawking temperature. We calculate the perturbation in surface gravity/temperature, finding it negative: the moon has a cooling effect on the black hole. We also compute the change in horizon angular frequency, which is positive, and the change in surface area/entropy, which vanishes.

preprint2012arXiv

Second Order Gravitational Self-Force

The second-order gravitational self-force on a small body is an important problem for gravitational-wave astronomy of extreme mass-ratio inspirals. We give a first-principles derivation of a prescription for computing the first and second perturbed metric and motion of a small body moving through a vacuum background spacetime. The procedure involves solving for a "regular field" with a specified (sufficiently smooth) "effective source", and may be applied in any gauge that produces a sufficiently smooth regular field.

preprint2011arXiv

A Note on the Coordinate Freedom in Describing the Motion of Particles in General Relativity

Our previous work developed a framework for treating the motion of a small body in general relativity, based on a one-parameter-family of solutions to Einstein's equation. Here we give an analysis of the coordinate freedom allowed within this framework, as is needed to determine the form of the equations of motion when they are expressed in general gauges.

preprint2011arXiv

A Rigorous Derivation of Gravitational Self-force

There is general agreement that the MiSaTaQuWa equations should describe the motion of a "small body" in general relativity, taking into account the leading order self-force effects. However, previous derivations of these equations have made a number of ad hoc assumptions and/or contain a number of unsatisfactory features. For example, all previous derivations have invoked, without proper justification, the step of "Lorenz gauge relaxation", wherein the linearized Einstein equation is written down in the form appropriate to the Lorenz gauge, but the Lorenz gauge condition is then not imposed--thereby making the resulting equations for the metric perturbation inequivalent to the linearized Einstein equations. In this paper, we analyze the issue of "particle motion" in general relativity in a systematic and rigorous way by considering a one-parameter family of metrics, $g_{ab} (λ)$, corresponding to having a body (or black hole) that is "scaled down" to zero size and mass in an appropriate manner. We prove that the limiting worldline of such a one-parameter family must be a geodesic of the background metric, $g_{ab} (λ=0)$. Gravitational self-force--as well as the force due to coupling of the spin of the body to curvature--then arises as a first-order perturbative correction in $λ$ to this worldline. No assumptions are made in our analysis apart from the smoothness and limit properties of the one-parameter family of metrics. Our approach should provide a framework for systematically calculating higher order corrections to gravitational self-force, including higher multipole effects, although we do not attempt to go beyond first order calculations here. The status of the MiSaTaQuWa equations is explained.

preprint2011arXiv

Dependence of Self-force on Central Object

For a particle in orbit about a static spherically symmetric body, we study the change in self-force that results when the central body type (i.e., the choice of interior metric for the Schwarzschild exterior) is changed. While a straight self-force is difficult to compute because of the need for regularization, such a "self-force difference" may be computed directly from the mode functions of the relevant wave equations. This technique gives a simple probe of the (non)locality of the force, as well as offers the practical benefit of an easy determination of the self-force on a body orbiting an arbitrary (static spherically symmetric) central body, once the corresponding result for a black hole (or some other reference interior) is known. We derive a general expression for the self-force difference at the level of a mode-sum in the case of a (possibly non-minimally coupled) scalar charge and indicate the generalization to the electromagnetic and gravitational cases. We then consider specific choices of orbit and/or central body. Our main findings are: (1) For charges held static at a large distance from the central body, the self-force is independent of the central body type in the minimally coupled scalar case and the electromagnetic case (but dependent in the nonminimally coupled scalar case); (2) For circular orbits about a thin-shell spacetime in the scalar case, the fractional change in self-force from a black hole spacetime is much larger for the radial (conservative) force than for the angular (dissipative) force; and (3) the radial self-force difference (between these spacetimes) agrees closely for a static charge and a circular orbit of the same radius.

preprint2011arXiv

Gauge and Averaging in Gravitational Self-force

A difficulty with previous treatments of the gravitational self-force is that an explicit formula for the force is available only in a particular gauge (Lorenz gauge), where the force in other gauges must be found through a transformation law once the Lorenz gauge force is known. For a class of gauges satisfying a ``parity condition'' ensuring that the Hamiltonian center of mass of the particle is well-defined, I show that the gravitational self-force is always given by the angle-average of the bare gravitational force. To derive this result I replace the computational strategy of previous work with a new approach, wherein the form of the force is first fixed up to a gauge-invariant piece by simple manipulations, and then that piece is determined by working in a gauge designed specifically to simplify the computation. This offers significant computational savings over the Lorenz gauge, since the Hadamard expansion is avoided entirely and the metric perturbation takes a very simple form. I also show that the rest mass of the particle does not evolve due to first-order self-force effects. Finally, I consider the ``mode sum regularization'' scheme for computing the self-force in black hole background spacetimes, and use the angle-average form of the force to show that the same mode-by-mode subtraction may be performed in all parity-regular gauges. It appears plausible that suitably modified versions of the Regge-Wheeler and radiation gauges (convenient to Schwarzschild and Kerr, respectively) are in this class.

preprint2010arXiv

Bobbing and Kicks in Electromagnetism and Gravity

We study systems analogous to binary black holes with spin in order to gain some insight into the origin and nature of "bobbing" motion and "kicks" that occur in this system. Our basic tool is a general formalism for describing the motion of extended test bodies in an external electromagnetic field in curved spacetime and possibly subject to other forces. We first show that bobbing of exactly the type as observed in numerical simulations of the binary black hole system occurs in a simple system consisting of two spinning balls connected by an elastic band in flat spacetime. This bobbing may be understood as arising from the difference between a spinning body's "lab frame centroid" and its true center of mass, and is purely "kinematical" in the sense that it will appear regardless of the forces holding two spinning bodies in orbit. Next, we develop precise rules for relating the motion of charged bodies in a stationary external electromagnetic field in flat spacetime with the motion of bodies in a weakly curved stationary spacetime. We then consider the system consisting of two orbiting charges with magnetic dipole moment and spin at a level of approximation corresponding to 1.5 post-Newtonian order. Here we find that considerable amounts of momentum are exchanged between the bodies and the electromagnetic field; however, the bodies store this momentum entirely as "hidden" mechanical momentum, so that the interchange does not give rise to any net bobbing. The net bobbing that does occur is due solely to the kinematical spin effect, and we therefore argue that the net bobbing of the electromagnetic binary is not associated with possible kicks. We believe that this conclusion holds in the gravitational case as well.

preprint2010arXiv

Motion of Small Bodies in Classical Field Theory

I show how prior work with R. Wald on geodesic motion in general relativity can be generalized to classical field theories of a metric and other tensor fields on four-dimensional spacetime that 1) are second-order and 2) follow from a diffeomorphism-covariant Lagrangian. The approach is to consider a one-parameter-family of solutions to the field equations satisfying certain assumptions designed to reflect the existence of a body whose size, mass, and various charges are simultaneously scaled to zero. (That such solutions exist places a further restriction on the class of theories to which our results apply.) Assumptions are made only on the spacetime region outside of the body, so that the results apply to ordinary bodies as well as black holes. The worldline "left behind" by the shrinking, disappearing body is interpreted as its lowest-order motion. An equation for this worldline follows from the "Bianchi identity" for the theory, without use of any properties of the field equations beyond their being second-order. The form of the force law for a theory therefore depends only on the ranks of its various tensor fields, with the field equations relevant only for determining the charges of a particular body from its exterior fields. I explicitly derive the force law (and mass-evolution law) in the case of scalar and vector fields, and give the recipe in the higher-rank case. Note that the vector force law is quite complicated, simplifying to the Lorentz force law only in the presence of the Maxwell gauge symmetry. Example applications of the results are the motion of "chameleon" bodies beyond the Newtonian limit, and the motion of bodies in (classical) non-Abelian gauge theory. I also make some comments on the role that scaling plays in the appearance of universality in the motion of bodies.