Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Boosting jet power in black hole spacetimes

David Neilsen, Luis Lehner, Carlos Palenzuela, Eric W. Hirschmann, Steven L. Liebling, Patrick M. Motl, and Travis Garrett
PNAS August 2, 2011 108 (31) 12641-12646; https://doi.org/10.1073/pnas.1019618108
David Neilsen
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Luis Lehner
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: llehner@perimeterinstitute.ca
Carlos Palenzuela
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric W. Hirschmann
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Steven L. Liebling
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick M. Motl
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Travis Garrett
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Roger D. Blandford, Stanford University, Menlo Park, CA, and approved June 21, 2011 (received for review January 7, 2011)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

The extraction of rotational energy from a spinning black hole via the Blandford–Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the interaction of the black hole with its surrounding plasma. We study the resulting Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter case, we find that increasing the orbital angular momenta of the system and/or the spins of the individual black holes results in an enhanced Poynting flux.

  • electromagnetic radiation
  • gravitational waves
  • multimessenger astronomy

Enormously powerful events illuminate the universe and challenge our understanding of the cosmos. Both gamma ray bursts and active galactic nuclei are among the most distant objects seen—a testament to the tremendous amounts of energy these sources can tap and radiate. Intense observational and theoretical efforts are ongoing in order to unravel these fascinating phenomena. Although the full details remain elusive, one of the natural ingredients in theoretical models is the inclusion of a rotating black hole which serves to convert binding and rotational energy of the system to electromagnetic radiation in a highly efficient manner. The starting point for these theoretical models can be traced back to ideas laid out by Penrose (1) and Blandford and Znajek (BZ) (2), which explain the extraction of energy from a rotating black hole. These seminal studies, along with subsequent work (see, e.g., refs. 3–5), have provided a basic understanding of highly energetic emissions from single black hole systems interacting with their surroundings. However, despite important theoretical and observational advances, we still lack a thorough understanding of these systems (e.g., refs. 6 and 7).

Recent work has indicated that related systems can also tap kinetic energy and lead to powerful jets (8). This work concentrated on nonspinning black holes moving through a plasma and highlighted that relative black hole motion alone (with respect to a stationary electromagnetic field topology at far distances from it) can induce the production of jets. Furthermore, subsequent work (9) demonstrated that even black holes with spins misaligned with respect to the asymptotic magnetic field direction induce strong emissions with power comparable to the aligned case. These studies suggested that, independent of their inclination, astrophysical jets might be powered by the efficient extraction of both rotational and translational kinetic energy of black holes, becoming even more powerful than the standard BZ mechanism would suggest.

Galactic mergers provide a likely scenario for yielding both the binary black hole and single black hole systems considered here (10, 11). In such a merger, the supermassive black holes associated with each galaxy will ultimately form a binary surrounded by a circumbinary disk in the resulting merged galaxy. A variety of interactions will tighten the black hole binary and eventually its dynamics will be governed by gravitational radiation reaction which drives the binary to merge. The circumbinary disk will likely be magnetized and thereby anchor magnetic field lines, some of which will traverse the central region containing the binary and eventual final black hole (see Fig. 1). Furthermore, as a result of the merger process, this final black hole will likely acquire a recoil velocity and will move relative to the circumbinary disk until the disk relaxes to the black hole’s motion. Thus, the system will have, at different stages, both binary and single black holes moving through an ambient magnetic field sourced by the surrounding disk. Preliminary observational evidence for supermassive black hole binaries resulting from galactic mergers has already been presented (12–14). Another example of a system where plasma could tap kinetic energy from a moving black hole corresponds to black-hole–neutron star binary (15). In such a system, the black hole would move through the star’s magnetic field. As discussed by BZ (2), an ambient magnetic field threading the black hole region populates a low density plasma surrounding the black hole. Even for black holes with no spin, it was recently shown that the orbiting binary interacting with the surrounding plasma can lead to a collimated Poynting flux (8). In this work, we consider this basic paradigm of energy extraction from black holes with the additional complexity of intrinsic black hole spin. Binaries consisting of spinning black holes demonstrate similar, albeit energetically enhanced, phenomenology. Furthermore, we investigate the dependence of the energy flux on the black hole velocity and highlight a resulting strong boost in the emitted power.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Diagram of the initial physical conditions. Two black holes (black circles) orbit inside a central region of a circumbinary disk (red ellipse) which sources the ambient magnetic field (blue dashed lines). The computational domain is restricted to a small section enclosing the black holes.

In addition to exquisite and powerful electromagnetic detectors, soon gravitational waves will be added to the arsenal of phenomena employed to understand our cosmos. These studies suggest excellent prospects for the coincident detection of both electromagnetic and gravitational signals from binary black hole systems. Certainly dual detection of electromagnetic and gravitational wave signals would transform our understanding of these systems and lead to the refinement of theoretical models (e.g., refs. 16–19).

Results

In this work, we discuss the basic phenomenology arising from the interaction of black hole systems immersed within a plasma environment. We study both boosted and orbiting black holes moving through an initially constant magnetic field and find that the interaction with the magnetic field induces jets of strong electromagnetic radiation. We find that, remarkably, both the spin and the translational motion of the black hole contribute to this electromagnetic luminosity. Indeed, the luminosity L of the jets is naturally decomposed in terms of each contribution as Embedded Image[1]Here the contribution from the spin, Lspin, is consistent with the expected BZ luminosity, and Lspeed represents the energy extracted from the black hole’s translational kinetic energy. In what follows, we conduct numerical experiments to detail how each component scales.

To fix ideas, we envision a magnetic field anchored in the disk with its associated magnetic dipole aligned, in the binary case, with the orbital angular momentum (chosen along the Embedded Image direction). In the single black hole case, we consider propagation velocities orthogonal to Embedded Image. Our resulting initial configuration thus has the magnetic field perpendicular to the velocity of the black holes and the electric field is set to zero. Because the electromagnetic field is affected by the spacetime curvature, it is dynamically distorted from its initial configuration and generates a transient burst of electromagnetic waves as it settles into a physically relevant and dynamical configuration. The initial magnitude of the magnetic field is B0 = 104 G. This value is chosen to be consistent with possible astrophysical magnetic field strengths (20, 21). We present our results for a field strength which is bounded by the Eddington magnetic field strength B ≃ 6 × 104(M/108 M⊙)-1/2 G (22). For these values, the plasma’s energy is several orders of magnitude smaller than that of the gravitational field. Thus, although the plasma is profoundly affected by the black holes, it has a negligible influence on the black holes.

For simplicity, we introduce the notation L43 = 1043 erg/s and compute fully dimensional quantities with respect to a representative system with total mass 108 M⊙ immersed in a magnetic field with strength B0 = 104 G. The luminosities that we calculate scale as L ∝ M2B2, and so luminosities for other masses and magnetic field strengths can be easily obtained from our results. We provide proportionality factors for other quantities of interest. Quantities in geometric units are calculated by setting Newton’s gravitation constant and the speed of light to unity, G = c = 1. For example, the Schwarzschild radius in geometric units, RS = 2M, is converted to other unit systems using the constants G and c as RS = 2GM/c2. Thus, RS ≃ 2.95 km for a solar mass black hole, and 2.95 × 108 km for the 108 M⊙ black hole considered here. Time in geometric units is measured in terms of a mass, and the conversion factor is G/c3. A time interval of 1 M⊙ corresponds to about 4.9 × 10-6 s, or 4.9 × 102 s for an interval of 108 M⊙.

The result captured in Eq. 1 has important consequences for black hole binaries because the luminosity that results from the motion, Lspeed, can be significantly higher than that resulting from the spin, Lspin. In addition, eccentric orbits and spin–orbit interactions driving orbital plane precession can induce important variabilities in the luminosity, which can also aid in the detection of these systems.

Single Black Holes.

We first consider the radiation from a single black hole immersed in a magnetic field anchored by currents in a distant circumbinary disk. Such a black hole can be spinning and can also have some constant velocity with respect to the initially uniform magnetic field. We calculate the resulting electromagnetic collimated flux energy and its dependence on both velocity and black hole spin. For these runs, we adopt a computational domain defined by [-160RS,160RS]3 (RS is the Schwarzschild radius of a black hole). We set the initial linear velocity of the black hole to Embedded Image (with n an integer that we vary between 0 and 4), and we set the intrinsic angular momentum parameter Embedded Image to either 0 or Embedded Image (recall that black hole spins must satisfy |a| ≤ 1 according to the Kerr bound).

Fig. 2 illustrates the qualitative features of the electromagnetic emission for a boosted black hole with and without spin. A collimated emission is clearly induced along the asymptotic magnetic field direction. An even stronger emission is obtained for the spinning black hole, as expected.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Zoom-in of the electromagnetic emission |rΦ2|2 from a boosted black hole (with Embedded Image). The left frame shows the nonspinning, single black hole case (a = 0), and the right frame shows the spinning case (a = 0.6). The color scale is the same for both frames and the black hole is boosted to the right. The stronger emission in the vertically collimated region of the right frame is apparent.

To achieve a more quantitative understanding, we compute a measure of the collimated electromagnetic luminosity, Lcollimated, by integrating the flux over a 15° cone of a spherical surface centered along the moving black hole with a radius r = 40RS. We present our results with respect to the velocity, v (the measured coordinate velocity of the black hole), and plot in Fig. 3 the collimated luminosity achieved once the system reaches a quasi-stationary state (in other words, after the initial transient stage) for the spinning and nonspinning cases.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Collimated luminosity as a function of the boost velocity of the black hole for two different spin magnitudes. The dashed lines are obtained from a quadratic fit of the form Lcollimated = Lspin + L2v2 where Lspin and L2 are two constants. The constant L2 is the same for both fits, supporting the argument that the luminosity does contain two different components, one for spin and one for boost velocity.

Several key observations are evident from the figure. For v = 0, the electromagnetic energy luminosity does not vanish for the spinning black hole, although it does vanish for the nonspinning black hole. That there is radiation in the spinning case is expected, as the spinning black hole interacts with a surrounding plasma and radiates by the BZ mechanism (2). This luminosity results from the plasma’s ability to extract rotational energy from the black hole and power a jet with an energy luminosity scaling as Embedded Image (9, 23) [with ΩH ≡ a/(2RH) the rotation frequency associated with the black hole].

For v ≠ 0, both the nonspinning and spinning black holes have a nonzero energy flux. In the former case, this flux arises solely from the ability of the system to tap translational kinetic energy from the black hole, whereas the latter results from the extraction of both translational and rotational kinetic energies.

The nonspinning black hole demonstrates an electromagnetic luminosity that increases with v2 (see Fig. 3). This behavior is consistent with treating the black hole within the membrane paradigm of ref. 24. A spinless black hole moving with speed v through an ambient magnetic field behaves as a conductor and acquires an induced charge proportional to the speed (8, 9, 25, 26). In this regard, the charge separation on the surface of the black hole is completely analogous to the classical Hall effect. With this observation and the induction equation, it is straightforward to show that the electromagnetic energy flux will increase as ∝ v2B2. This quadratic dependence on speed is apparent in the figure. These results are also consistent with the work of Drell et al. (27), who studied the Poynting flux associated with a moving conductor in a magnetized plasma as applied to the motion of artificial satellites in orbit. They found that the flux obeys Lv ≈ (v/valf)2B2 (27), where valf, the propagation speed of the Alfvén modes, is the speed of light in a force-free environment. Furthermore, a misalignment is expected between the collimated energy flux and the original magnetic field orientation such that tan(α) = v/valf. With the cautionary note that measuring this angle is ambiguous in the curved spacetime around the black hole, such a relation is indeed manifested in our results. For instance, for a = 0, vx = 0.10, the measured angle is tan(α) = 0.07, in good agreement with the predicted value of 0.08 (27).

When we allow the black hole to spin, we once again find the same quadratic dependence on velocity, with an additional contribution due to the hole’s spin. Notice that the difference between the obtained luminosities in the spinning and spinless cases remains fairly constant for the different values of v. The roughly constant offset between spinning and nonspinning cases supports the decomposition given in Eq. 1 in which the luminosity receives two contributions, one from the black hole spin and another from its velocity.

We have seen that Embedded Image, and that Lspeed ∝ v2. In order to be more definite, we introduce two fitting constants, L1 and L2, such that the general luminosity dependence of Eq. 1 can be expressed as Embedded Image[2]A fit of the luminosity for the nonspinning black hole (a = 0) gives L2 = 127. Similarly, a fit for a spinning black hole with a = 0.6 provides the constant L1 = 0.87. In the above equation, we express the spin contribution in a general way so that, to obtain the expected luminosity for an arbitrary black hole, one need only provide the ratio of the rotational frequency to that of our fiducial a = 0.6 black hole.

Notice that, for a = 0.6, the nonrotational contribution Lspeed to the emitted power becomes larger than the rotational component for speeds approximately v ≳ 0.08c. This relationship, for example, would predict a luminosity for an a = 0.95, v = 0.5c black hole to be ≃36L43(M8B4)2. This phenomenology strongly suggests that the Poynting flux can tap both rotational and translational kinetic energies from the black hole and that faster and more rapidly spinning black holes have a stronger associated power output.

Binary Black Holes.

We turn our attention now to orbiting binary black holes. We consider equal-mass black holes with orbital diameter D ≈ 8RS. We adopt a computational domain of [-320RS,320RS]3 and consider black holes with either no spin (abbreviated “0” below), or with a spin aligned (spin up or “u”) or antialigned (spin down or “d”) with the z direction. In particular, we concentrate on three cases: a binary with no individual spin; a binary with zero net spin (by considering equal, but antialigned, individual spins); and a binary with equal, aligned spins. For the spinning black hole cases, all spinning holes have spin parameter |a| = 0.515. Where required, we denote these cases by 0/0, u/d, and u/u, respectively, and associated measured quantities have this notation as subindices. These configurations are summarized in Table 1.

View this table:
  • View inline
  • View popup
Table 1.

Binary black hole configurations considered here

Before discussing the results from these binaries, we first describe our analysis that enables a quantitative comparison between the evolutions. First, we evaluate the luminosities as functions of gravitational wave frequency, as this is an observable and allows for a direct comparison of the different cases. We obtain frequencies from the dominant l = 2, m = 2 gravitational mode. Second, we compute three different luminosities for each case: (i) the collimated luminosity Lcollimated obtained by integrating the electromagnetic flux over a cone of points within 15° of vertical from the center of mass of the system; (ii) the noncollimated, or “isotropic,” luminosity Lisotropic obtained from the integral over an encompassing sphere minus the collimated luminosity of (i); and (iii) the gravitational wave luminosity, LGW. These different luminosities are displayed for the three binary configurations in Fig. 4.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Luminosity as a function of the orbital frequency for the binary black hole configurations described in Table 1, assuming M = 108 M⊙ and B = 104 G. (Upper Left) The collimated luminosity associated with the jets. (Upper Right) The isotropic luminosity representing electromagnetic flux not associated with the jets. (Lower) The gravitational wave output.

We first consider the 0/0 and u/d binaries, which have essentially the same total angular momentum. Their qualitative behavior is illustrated in Fig. 5. Extensive numerical simulations (see, for instance, ref. 28) and simple estimates (29) indicate that both these binaries will merge into a final black hole with essentially the same spin (a ≃ 0.67), and so for late times after the merger, the expected jet structure should be quite similar: namely, that described by the standard BZ mechanism.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Electromagnetic flux corresponding to the 0/0 (top row) and u/d (bottom row) cases at times (-16.7,-2.5) h (10-8 M/M⊙) and (-10.5,3.6) h (10-8 M/M⊙) (with respect to the merger time as marked by the peak in strength of gravitational wave emission), respectively. The orbiting stage leaves its imprint as twisted tubes.

The binary with individual spins aligned reaches a higher orbital velocity before merger than the previous two cases; thus, the expected maximum power should be higher. Moreover, the resulting final black hole spins faster (a ≃ 0.8) and thus its BZ associated power will be higher than that of the previous cases. Fig. 4 illustrates this expected behavior by presenting the Poynting flux energy vs. gravitational wave frequency for the three cases.

As evident from the figure, at low frequencies where the orbital dynamics are the same in all cases, both spinning cases have a higher output than the nonspinning one and the difference is provided by the spin contribution to the jet emission. Furthermore, both spinning cases have equal collimated power output because the spin contribution to the luminosity depends only on the spin magnitude.

We can further examine the basic relation given in Eq. 1 by comparing the fluxes corresponding to the three cases studied. We estimate the black hole coordinate velocities v and (collimated) luminosities at three representative frequencies Ωi = {1; 1.5; 2}10-5 [Hz/M8] (i = 1,…,3) before the strong nonlinear interaction starts. Because for these frequencies the measured speeds for the 0/0 and u/d cases are essentially the same, differences in their corresponding fluxes should be due to the contribution from the spin. This contribution can be estimated to be Embedded Image[3]where Lu/d and L0/0 are the measured luminosities for both cases. The obtained values for the frequencies considered are as follows: Notice that the Embedded Image remains fairly constant through these frequencies, which is evident also in Fig. 4. That this value is essentially the same for all cases further supports the basic relation in Eq. 1. Furthermore, we can employ this value to estimate the luminosity for the u/u case (illustrated in Fig. 6) exploiting the one measured for the 0/0 case. To do so, we observe that for a given frequency Ω the difference between the two cases is the contribution due to the spin (absent in the 0/0 case) together with a suitable scaling of Lspeed as Embedded Image[4]The obtained values can be compared with those computed directly from the numerical evolution, Lu/u. Using our measured orbital velocities, Eq. 4 provides the following estimates which we show next to our computed values: The estimates are within ≃10% of the measured values, when the black holes are sufficiently separated that their jets do not strongly interact. Finally, notice that, at even higher frequencies, Fig. 4 illustrates that the aligned (u/u) case indeed has a higher associated power.

View this table:
  • View inline
  • View popup
View this table:
  • View inline
  • View popup
Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Electromagnetic flux corresponding to the u/u case at times (-2.5,11.6,40,68) h (10-8M/M⊙) (with respect to the merger time as marked by the peak in strength of gravitational wave emission). The orbiting behavior leaves a clear imprint in the resulting flux of energy and, as the merger proceeds, emission along all directions is evident. At late times, the system settles to a single jet as dictated by the BZ mechanism.

In all cases, a significant noncollimated emission is induced during the merger phase (evident in the upper-right plot of Fig. 4 and illustrated in the second frames of Figs. 5 and 6). Clearly, the simple-minded picture of a jet produced by the superposition of the orbital and spinning effect cannot fully capture the complete behavior at the merger epoch, although it serves to understand the main qualitative features and provides a means to estimate the power of the electromagnetic emission.

Discussion

We have studied the impact of black hole motion through a plasma and indicated how the interaction can induce powerful electromagnetic emissions even for nonspinning black holes. Despite having examined a very small subset of the binary black hole parameter space, the results presented both here and in refs. 8 and 9 suggest broad applicability to general black hole binaries. Moreover, as the plasma generally has a negligible effect on the dynamics of the black holes, one needs only to know the dynamics of the black holes, say by numerical solution or other approximate methods, in order to estimate the expected electromagnetic luminosity. A recent example is the work of McWilliams (30), which uses the expected ≃B2v2 scaling for nonspinning black hole binaries coupled with the known distance dependence in time to obtain excellent agreement with the luminosity from fully nonlinear numerical evolutions.

Naturally, spinning black holes produce stronger jets, and these jets, as shown earlier (9), will be aligned with the asymptotic magnetic field direction.* The resulting electromagnetic luminosity can be estimated to be L ≃ Lspin + Lspeed, which can be significant and have associated time variabilities tied to the dynamical behavior of the system. In particular, the luminosity tied to the motion can become significantly higher than that tied to the spin. In addition, eccentric orbits and spin–orbit interactions driving orbital plane precessions can induce important variabilities that can aid in the detection of these systems. Furthermore, a significant pulse of nearly isotropic radiation is emitted during merger, thereby allowing observations of the system along directions not aligned with the jet.

Consequently, binary black hole interactions with surrounding plasmas can yield powerful electromagnetic outputs and allow for observing these systems through both gravitational and electromagnetic radiation. Gravitational waves from these systems corresponding to the last year before the merger could be observed to large distances with the Laser Interferometer Space Antenna (up to redshifts of 5–10, ref. 31, for masses ≃104–7 M⊙) or earlier in the orbiting phase (and possibly through merger) via Pulsar Timing Array observations (32), targeting binaries with masses of M ≃ 107–10 M⊙. As we have indicated here, both scenarios can have strong associated electromagnetic emissions. Our inferred luminosities of several 1043(M8B4)2 erg/s correspond to an isotropic bolometric flux of Fx ≃ 10-15 erg (M8B4)2/(cm2 s) that could be detected to redshifts of z ≈ 1 and beyond depending on anisotropies and on the efficiency of processes which tap this available energy and produce observable signals.

Materials and Methods

The combined gravitational and electromagnetic systems that we consider consist of black hole spacetimes in which the black holes can be regarded as immersed in an external magnetic field. Such fields, as mentioned earlier, will be anchored to a disk. We consider this disk to be outside our computational domain but its influence is realized through the imposition of suitable boundary conditions on the incoming electric and magnetic field modes on the boundaries of our domain (8, 9). These conditions essentially correspond to setting to zero the electric field while constant values for the magnetic field (given by Embedded Image.) With regard to the magnetosphere around the black holes, we assume that the energy density of the magnetic field dominates over its tenuous density such that the inertia of the plasma can be neglected. The magnetosphere is therefore treated within the force-free approximation (2, 33). We note that the contribution of the energy density of the plasma to the dynamics of the spacetime is negligible and we can ignore its back reaction on the spacetime.

We use the Baumgarte–Shapiro–Shibata–Nakamura formulation (34, 35) of the Einstein equations and the force-free equations as described in refs. 8, 9, and 36. We discretize the equations using finite difference techniques on a regular Cartesian grid and use adaptive mesh refinement (AMR) to ensure that sufficient resolution is available where required in an efficient manner. We use the Hybrid Adaptive computational infrastructure, which provides distributed, Berger–Oliger style AMR (37, 38) with full subcycling in time, together with an improved treatment of artificial boundaries (39). Refinement regions are determined using truncation error estimation provided by a shadow hierarchy (40), which adapts dynamically to ensure the estimated error is bounded by a prespecified tolerance. Typically our adopted values result in a grid hierarchy yielding a resolution such that 40 grid points in each direction cover each black hole. We use a fourth-order accurate spatial discretization and a third-order accurate in time Runge–Kutta integration scheme (41). We adopt a Courant parameter of λ = 0.4, so that Δtℓ = 0.4Δxℓ on each refinement level ℓ. In tests performed here for the coupled system (and in our previous works for the force-free Maxwell equations), the code demonstrates convergence while maintaining small constraint residuals for orbiting black holes. Furthermore, we obtain for orbiting black hole evolutions agreement with runs from other codes for the same initial data.

To extract physical information, we monitor the Newman–Penrose electromagnetic (Φ2) and gravitational (Ψ4) radiative scalars (42). These scalars are computed by contracting the Maxwell and the Weyl tensors, respectively, with a suitably defined null tetrad (as discussed in ref. 43), Embedded Image[5]and they allow us to account for the energy carried off by outgoing waves at infinity. The scalar Φ2 (whose modulus squared is essentially the radial component of the Poynting vector) provides a measure of the electromagnetic radiation at large distances from an isolated system. However, as the system studied here has a ubiquitous magnetic field, special care must be taken to compute the energy flux. We account for this difficulty by subtracting the scalar Φ0 = -Fablamb from Φ2. Hence, by Φ2, we mean the difference Φ2 → Φ2 - Φ0. The luminosities in electromagnetic and gravitational waves are given by the integrals of the fluxes Embedded Image[6]Embedded Image[7]

Acknowledgments

The authors thank J. Arons, P. Chang, B. MacNamara, K. Menou, E. Quataert, and C. Thompson, as well as our long time collaborators Matthew Anderson, Miguel Megevand, and Oscar Reula for useful discussions and comments. We acknowledge support from National Science Foundation Grants PHY-0803629 (to Louisiana State University), PHY-0969811 (to Brigham Young University), PHY-0969827 (to Long Island University), as well as the Natural Sciences and Engineering Research Council through a Discovery Grant. Research at Perimeter Institute is supported through Industry Canada and by the Province of Ontario through the Ministry of Research and Innovation. Computations were performed at Louisiana Optical Network Initiative, Teragrid, and Scinet.

Footnotes

  • ↵1To whom correspondence should be addressed. E-mail: llehner{at}perimeterinstitute.ca.
  • Author contributions: L.L. designed research; D.N., L.L., C.P., and T.G. performed research; L.L., C.P., E.W.H., S.L.L., P.M.M., and T.G. contributed new analytic tools; D.N., C.P., and E.W.H. code development; D.N., L.L., C.P., S.L.L., and P.M.M. analyzed data; and D.N., L.L., C.P., E.W.H., and S.L.L. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • ↵*The magnitude of the BZ-associated emission diminishes for nonaligned cases but it is nevertheless significant: The orthogonal case is only half as powerful as the aligned case.

References

  1. ↵
    1. Penrose R
    (1969) Gravitational collapse: The role of general relativity. Riv Nuovo Cimento Soc Ital Fis 1:252–276.
    OpenUrlCrossRef
  2. ↵
    1. Blandford RD,
    2. Znajek RL
    (1977) Electromagnetic extractions of energy from Kerr black holes. Mon Not R Astron Soc, 179–433, 456.
  3. ↵
    1. Lee HK,
    2. Wijers RAMJ,
    3. Brown GE
    (2000) The Blandford-Znajek process as a central engine for a gamma ray burst. Phys Rep 325:83–114.
    OpenUrlCrossRef
  4. ↵
    1. Gilfanov M,
    2. Sunyeav R,
    3. Churazov E
    1. Blandford RD
    (2002) in To the Lighthouse, Lighthouses of the Universe: The Most Luminous Celestial Objects and Their Use for Cosmology, eds Gilfanov M, Sunyeav R, Churazov E (Springer, Berlin), pp 381–404.
  5. ↵
    1. Punsly B
    , ed (2009) Black Hole Gravitohydromagnetics, Astrophysics and Space Science Library (Springer, Berline), 355.
  6. ↵
    1. Lijun G,
    2. et al.
    (2010) The spin of the black hole in the soft X-ray transient A0620-00. Astrophys J Lett 718:L122–L126.
    OpenUrlCrossRef
  7. ↵
    1. Fender R,
    2. Gallo E,
    3. Russell D
    (2010) No evidence for black hole spin powering of jets in X-ray binaries. Mon Not R Astron Soc 406:1425–1434.
    OpenUrl
  8. ↵
    1. Palenzuela C,
    2. Lehner L,
    3. Liebling SL
    (2010) Dual jets from binary black holes. Science 329:927–930.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Palenzuela C,
    2. Garrett T,
    3. Lehner L,
    4. Liebling SL
    (2010) Magnetospheres of black hole systems in force-free plasma. Phys Rev D 82:044045.
    OpenUrlCrossRef
  10. ↵
    1. Begelman MC,
    2. Blandford RD,
    3. Rees MJ
    (1980) Massive black hole binaries in active galactic nuclei. Nature 287:307–309.
    OpenUrlCrossRef
  11. ↵
    1. Milosavljevic M,
    2. Phinney ES
    (2005) The afterglow of massive black hole coalescence. Astrophys J 622:L93–L96.
    OpenUrlCrossRef
  12. ↵
    1. Komossa S,
    2. et al.
    (2003) Discovery of a binary AGN in the ultraluminous infrared galaxy NGC 6240 using Chandra. Astrophys J 582:L15–L20.
    OpenUrlCrossRef
  13. ↵
    1. Rodriguez C,
    2. et al.
    (2006) A compact supermassive binary black hole system. Astrophys J 646:49–60.
    OpenUrlCrossRef
  14. ↵
    1. Comerford JM,
    2. et al.
    (2009) Inspiralling supermassive black holes: A new signpost for galaxy mergers. Astrophys J 698:956–965.
    OpenUrlCrossRef
  15. ↵
    1. McWilliams ST,
    2. Levin J
    (2011) Electromagnetic extraction of energy from black hole-neutron star binaries. arXiv:1101.1969.
  16. ↵
    1. Gilfanov M,
    2. Sunyeav R,
    3. Churazov E
    1. Schutz BF
    (2002) in Lighthouses of Gravitational Wave Astronomy Lighthouses of the Universe: The Most Luminous Celestial Objects and Their Use for Cosmology, eds Gilfanov M, Sunyeav R, Churazov E (Springer, Berlin), pp 207–224.
  17. ↵
    1. Sylvestre J
    (2003) Prospects for the detection of electromagnetic counterparts to gravitational wave events. Astrophys J 591:1152–1156.
    OpenUrlCrossRef
  18. ↵
    1. Bloom JS,
    2. et al.
    (2009) Coordinated science in the gravitational and electromagnetic skies. Astro2010 Decadal Survey Whitepaper, arXiv:0902.1527.
  19. ↵
    1. Phinney ES
    (2009) Finding and using electromagnetic counterparts of gravitational wave sources. Astro2010 Decadal Survey Whitepaper, arXiv:0903.0098.
  20. ↵
    1. Massi M,
    2. Kaufman M
    (2008) Magnetic field upper limits for jet formation. Astron Astrophys 477:1–7.
    OpenUrlCrossRef
  21. ↵
    1. Field GB,
    2. Rogers RD
    (1993) Radiation from magnetized accretion disks in active galactic nuclei. Astrophys J 403:94.
    OpenUrlCrossRef
  22. ↵
    1. Dermer CD,
    2. Finke JD,
    3. Menon G
    (2008) Black-hole engine kinematics, flares from PKS 2155-304, and multiwavelength blazar analysis. arXiv:0810.1055.
  23. ↵
    1. Tchekhovskoy A,
    2. Narayan R,
    3. McKinney JC
    (2010) Black hole spin and the radio loud/quiet dichotomy of active galactic nuclei. Astrophys J 711:50–63.
    OpenUrlCrossRef
  24. ↵
    1. Thorne KS,
    2. Price RH,
    3. MacDonald DA
    (1986) Black Holes: The Membrane Paradigm (Yale Univ Press, New Haven, CT).
  25. ↵
    1. Palenzuela C,
    2. Anderson M,
    3. Lehner L,
    4. Liebling SL,
    5. Neilsen D
    (2009) Binary black holes’ effects on electromagnetic fields. Phys Rev Lett 103:081101.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Palenzuela C,
    2. Lehner L,
    3. Yoshida S
    (2010) Understanding possible electromagnetic counterparts to loud gravitational wave events. Phys Rev D 81:084007.
    OpenUrlCrossRef
  27. ↵
    1. Drell SD,
    2. Foley HM,
    3. Ruderman MA
    (1965) Drag and propulsion of large satellites in the ionosphere: An Alfvén propulsion engine in space. Phys Rev Lett 14:171–175.
    OpenUrlCrossRef
  28. ↵
    1. Hannam M,
    2. et al.
    (2009) The Samurai Project: Verifying the consistency of black- hole-binary waveforms for gravitational-wave detection. Phys Rev D 79:084025.
    OpenUrlCrossRef
  29. ↵
    1. Buonanno A,
    2. Kidder LE,
    3. Lehner L
    (2008) Estimating the final spin of a binary black hole coalescence. Phys Rev D 77:026004.
    OpenUrlCrossRef
  30. ↵
    1. McWilliams ST
    (2011) The status of black-hole binary merger simulations with numerical relativity. Classical Quantum Gravity 28:134001.
    OpenUrlCrossRef
  31. ↵
    1. Vecchio A
    (2004) LISA observations of rapidly spinning massive black hole binary systems. Phys Rev D 70:042001.
    OpenUrlCrossRef
  32. ↵
    1. Demorest P,
    2. Lazio J,
    3. Lommen A,
    4. NANOGRAV collaboration
    (2009) Gravitational Wave Astronomy Using Pulsars: Massive Black Hole Mergers & The Early Universe. Astro2010 Decadal Survey Whitepaper, arXiv:0902.2968.
  33. ↵
    1. Goldreich P,
    2. Julian WH
    (1969) Pulsar electrodynamics. Astrophys J 157:869–880.
    OpenUrlCrossRef
  34. ↵
    1. Baumgarte TW,
    2. Shapiro SL
    (1999) On the numerical integration of Einstein’s field equations. Phys Rev D 59:024007.
    OpenUrl
  35. ↵
    1. Shibata M,
    2. Nakamura T
    (1995) Evolution of three-dimensional gravitational waves: Harmonic slicing case. Phys Rev D 52:5428–5444.
    OpenUrlCrossRef
  36. ↵
    1. Palenzuela C,
    2. Bona C,
    3. Lehner L,
    4. Reula O
    (2011) Robustness of the Blandford-Znajek mechanism. Classical Quantum Gravity 28:134007.
    OpenUrlCrossRef
  37. ↵
    1. Liebling SL,
    2. et al.
    (2008) HAD. http://www.had.liu.edu/.
  38. ↵
    1. Liebling SL
    (2002) Singularity threshold of the nonlinear sigma model using 3D adaptive mesh refinement. Phys Rev D 66:041703.
    OpenUrlCrossRef
  39. ↵
    1. Lehner L,
    2. Liebling SL,
    3. Reula O
    (2006) AMR, stability and higher accuracy. Classical Quantum Gravity 23:S421–S445.
    OpenUrlCrossRef
  40. ↵
    1. Pretorius F
    (2002) Numerical simulations of gravitational collapse. PhD thesis (University of British Columbia, Vancouver).
  41. ↵
    1. Anderson M,
    2. et al.
    (2008) Simulating binary neutron stars: Dynamics and gravitational waves. Phys Rev D 77:024006.
    OpenUrlCrossRef
  42. ↵
    1. Newman E,
    2. Penrose R
    (1962) An approach to gravitational radiation by a method of spin coefficients. J Math Phys 3:566–578.
    OpenUrlCrossRef
  43. ↵
    1. Lehner L,
    2. Moreschi OM
    (2007) Dealing with delicate issues in waveforms calculations. Phys Rev D 76:124040.
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Boosting jet power in black hole spacetimes
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Boosting jet power in black hole spacetimes
David Neilsen, Luis Lehner, Carlos Palenzuela, Eric W. Hirschmann, Steven L. Liebling, Patrick M. Motl, Travis Garrett
Proceedings of the National Academy of Sciences Aug 2011, 108 (31) 12641-12646; DOI: 10.1073/pnas.1019618108

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Boosting jet power in black hole spacetimes
David Neilsen, Luis Lehner, Carlos Palenzuela, Eric W. Hirschmann, Steven L. Liebling, Patrick M. Motl, Travis Garrett
Proceedings of the National Academy of Sciences Aug 2011, 108 (31) 12641-12646; DOI: 10.1073/pnas.1019618108
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 108 (31)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Physics

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Surgeons hands during surgery
Inner Workings: Advances in infectious disease treatment promise to expand the pool of donor organs
Despite myriad challenges, clinicians see room for progress.
Image credit: Shutterstock/David Tadevosian.
Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Double helix
Journal Club: Noncoding DNA shown to underlie function, cause limb malformations
Using CRISPR, researchers showed that a region some used to label “junk DNA” has a major role in a rare genetic disorder.
Image credit: Nathan Devery.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Multi-color molecular model
Enzymatic breakdown of PET plastic
A study demonstrates how two enzymes—MHETase and PETase—work synergistically to depolymerize the plastic pollutant PET.
Image credit: Aaron McGeehan (artist).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490