Stability of relativistic jets from rotating, accreting black holes via fully three-dimensional magnetohydrodynamic simulations

GRMHD Simulations and Modeling for Jet Formation and Acceleration Region in AGNs

Relativistic jets are collimated plasma outflows with relativistic speeds. Astrophysical objects involving relativistic jets are a system comprising a compact object such as a black hole, surrounded by rotating accretion flows, with the relativistic jets produced near the central compact object. The most accepted models explaining the origin of relativistic jets involve magnetohydrodynamic (MHD) processes. Over the past few decades, many general relativistic MHD (GRMHD) codes have been developed and applied to model relativistic jet formation in various conditions. This short review provides an overview of the recent progress of GRMHD simulations in generating relativistic jets and their modeling for observations.

The Origin of Matter at the Base of Relativistic Jets in Active Galactic Nuclei

The generation of relativistic jets in active sources such as blazars is a complex problem with many aspects, most of them still not fully understood. Relativistic jets are likely produced by the accretion of matter and magnetic fields onto spinning black holes. Ergospheric dragging effects launch a Poynting-dominated outflow in the polar directions of these systems. Observations with very high resolution of the jet in the nearby radio galaxy M87 and evidence of extremely fast variability in the non-thermal radiation of several other objects indicate that charged particles produce synchrotron emission and gamma rays very close to the base of the jet. How these particles are injected into the magnetically shielded outflow is a mystery. Here we explore the effects of various processes in the hot accretion inflow close to the black hole that might result in the copious production of neutral particles which, through annihilation and decay in the jet’s funnel, might load the outflow with mass and charged particles on scales of a few Schwarzschild radii.

PKS 2155-304: A Case Study of Blazar Variability Power Spectrum at the Highest Energies and on the Longest Timescales

We present the results of our Power Spectral Density (PSD) analysis for the BL Lac object PKS 2155−304, utilizing the nightly-binned long-term light curve from the decade-long monitoring, as well as the minute-binned intra-night light curve from the High Energy Stereoscopic Survey (H.E.S.S.; >200 GeV). The source is unique for exhibiting the shortest flux-doubling timescale at Very High Energy (VHE) among its class and thus provides a rare opportunity to study the particle acceleration on the smallest spatial scales in blazar jets. The light curves are modeled in terms of the Continuous-Time Auto-Regressive Moving Average (CARMA) process. The combined long-term and intra-night PSD extends up to ∼6 decades in the temporal frequency range; unprecedented at the TeV energies for a blazar source. Our systematic approach reveals that PKS 2155−304 shows, on average, a complex shape of variability power spectrum, with more variability power on longer timescales. The long-term variability is best modeled by the CARMA(2,1) process, while the intra-night variability is modeled by a CARMA(1,0) process. We note that the CARMA(1,0) process refers to an Ornstein–Uhlenbeck process where the power-law PSD slope (PSD varies as a function of variability frequency to the power of the negative slope) changes from two to zero, above a certain “characteristic/relaxation” timescale. Even though the derived power spectrum of the intra-night light curve did not reveal a flattening, we speculate such relaxation must occur on timescales longer than a few hours for the source.

Approaching the Black Hole by Numerical Simulations

Black holes represent extreme conditions of physical laws. Predicted about a century ago, they are now accepted as astrophysical reality by most of the scientific community. Only recently has more direct evidence of their existence been found—the detection of gravitational waves from black hole mergers and of the shadow of a supermassive black hole in the center of a galaxy. Astrophysical black holes are typically embedded in an active environment which is affected by the strong gravity. When the environmental material emits radiation, this radiation may carry imprints of the black hole that is hosting the radiation source. In order to understand the physical processes that take place in the close neighborhood of astrophysical black holes, numerical methods and simulations play an essential role. This is simply because the dynamical evolution and the radiative interaction are far too complex in order to allow for an analytic solution of the physical equations. A huge progress has been made over the last decade(s) in the numerical code development, as well as in the computer power that is needed to run these codes. This review tries to summarize the basic questions and methods that are involved in the undertaking of investigating the astrophysics of black holes by numerical means. It is intended for a non-expert audience interested in an overview over this broad field. The review comes along without equations and thus without a detailed expert discussion of the underlying physical processes or numerical specifics. Instead, it intends to illustrate the richness of the field and to motivate further reading. The review puts some emphasis on magneto-hydrodynamic simulations but also touches radiation transfer and merger simulations, in particular pointing out differences in these approaches.

Numerical Simulations of Jets from Active Galactic Nuclei

Numerical simulations have been playing a crucial role in the understanding of jets from active galactic nuclei (AGN) since the advent of the first theoretical models for the inflation of giant double radio galaxies by continuous injection in the late 1970s. In the almost four decades of numerical jet research, the complexity and physical detail of simulations, based mainly on a hydrodynamical/magneto-hydrodynamical description of the jet plasma, have been increasing with the pace of the advance in theoretical models, computational tools and numerical methods. The present review summarizes the status of the numerical simulations of jets from AGNs, from the formation region in the neighborhood of the supermassive central black hole up to the impact point well beyond the galactic scales. Special attention is paid to discuss the achievements of present simulations in interpreting the phenomenology of jets as well as their current limitations and challenges.

Magnetorotational collapse of supermassive stars: Black hole formation, gravitational waves, and jets

We perform magnetohydrodynamic simulations in full general relativity of uniformly rotating stars that are marginally unstable to collapse. These simulations model the direct collapse of supermassive stars (SMSs) to seed black holes that can grow to become the supermassive black holes at the centers of quasars and active galactic nuclei. They also crudely model the collapse of massive Population III stars to black holes, which could power a fraction of distant, long gamma-ray bursts. The initial stellar models we adopt are Γ = 4/3 polytropes initially with a dynamically unimportant dipole magnetic field. We treat initial magnetic-field configurations either confined to the stellar interior or extending out from the stellar interior into the exterior. We find that the black hole formed following collapse has mass MBH ≃ 0.9M (where M is the mass of the initial star) and dimensionless spin parameter aBH/MBH ≃ 0.7. A massive, hot, magnetized torus surrounds the remnant black hole. At Δt ~ 400-550M ≈ 2000 - 2700(M/106M⊙)s following the gravitational wave peak amplitude, an incipient jet is launched. The disk lifetime is Δt ~ 105(M/106M⊙)s, and the outgoing Poynting luminosity is LEM ~ 1051-52 ergs/s. If > ˜ 1 % - 10 % of this power is converted into gamma rays, Swift and Fermi could potentially detect these events out to large redshifts z ~ 20. Thus, SMSs could be sources of ultra-long gamma-ray bursts (ULGRBs), and massive Population III stars could be the progenitors that power a fraction of the long GRBs observed at redshift z ~ 5-8. Gravitational waves are copiously emitted during the collapse and peak at ~15(106M⊙/M) mHz [~0.15(104 M⊙/M) Hz], i.e., in the LISA (DECIGO/BBO) band; optimally oriented SMSs could be detectable by LISA (DECIGO/BBO) at z < ˜ 3 ( z < ˜ 11 ) .Hence, 104M⊙ SMSs collapsing at z ~ 10 are promising multimessenger sources of coincident gravitational and electromagnetic waves.

Riemann solvers and Alfven waves in black hole magnetospheres

In the magnetosphere of a rotating black hole, an inner Alfven critical surface (IACS) must be crossed by inflowing plasma. Inside the IACS, Alfven waves are inward directed toward the black hole. The majority of the proper volume of the active region of spacetime (the ergosphere) is inside of the IACS. The charge and the totally transverse momentum flux (the momentum flux transverse to both the wave normal and the unperturbed magnetic field) are both determined exclusively by the Alfven polarization. Thus, it is important for numerical simulations of black hole magnetospheres to minimize the dissipation of Alfven waves. Elements of the dissipated wave emerge in adjacent cells regardless of the IACS, there is no mechanism to prevent Alfvenic information from crossing outward. Thus, numerical dissipation can affect how simulated magnetospheres attain the substantial Goldreich-Julian charge density associated with the rotating magnetic field. In order to help minimize dissipation of Alfven waves in relativistic numerical simulations we have formulated a one-dimensional Riemann solver, called HLLI, which incorporates the Alfven discontinuity and the contact discontinuity. We have also formulated a multidimensional Riemann solver, called MuSIC, that enables low dissipation propagation of Alfven waves in multiple dimensions. The importance of higher order schemes in lowering the numerical dissipation of Alfven waves is also catalogued.

Grid-based Methods in Relativistic Hydrodynamics and Magnetohydrodynamics

An overview of grid-based numerical methods used in relativistic hydrodynamics (RHD) and magnetohydrodynamics (RMHD) is presented. Special emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods. Results of a set of demanding test bench simulations obtained with different numerical methods are compared in an attempt to assess the present capabilities and limits of the various numerical strategies. Applications to three astrophysical phenomena are briefly discussed to motivate the need for and to demonstrate the success of RHD and RMHD simulations in their understanding. The review further provides FORTRAN programs to compute the exact solution of the Riemann problem in RMHD, and to simulate 1D RMHD flows in Cartesian coordinates.

High-energy transients

We present an overview of high-energy transients in astrophysics, highlighting important advances over the past 50 years. We begin with early discoveries of γ-ray transients, and then delve into physical details associated with a variety of phenomena. We discuss some of the unexpected transients found by Fermi and Swift, many of which are not easily classifiable or in some way challenge conventional wisdom. These objects are important insofar as they underscore the necessity of future, more detailed studies.

Foundations of Black Hole Accretion Disk Theory

This review covers the main aspects of black hole accretion disk theory. We begin with the view that one of the main goals of the theory is to better understand the nature of black holes themselves. In this light we discuss how accretion disks might reveal some of the unique signatures of strong gravity: the event horizon, the innermost stable circular orbit, and the ergosphere. We then review, from a first-principles perspective, the physical processes at play in accretion disks. This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks, and advection-dominated accretion flows (ADAFs). After presenting the models we discuss issues of stability, oscillations, and jets. Following our review of the analytic work, we take a parallel approach in reviewing numerical studies of black hole accretion disks. We finish with a few select applications that highlight particular astrophysical applications: measurements of black hole mass and spin, black hole vs. neutron star accretion disks, black hole accretion disk spectral states, and quasi-periodic oscillations (QPOs).

Alignment of magnetized accretion disks and relativistic jets with spinning black holes

Accreting black holes (BHs) produce intense radiation and powerful relativistic jets, which are affected by the BH's spin magnitude and direction. Although thin disks might align with the BH spin axis via the Bardeen-Petterson effect, this does not apply to jet systems with thick disks. We used fully three-dimensional general relativistic magnetohydrodynamical simulations to study accreting BHs with various spin vectors and disk thicknesses and with magnetic flux reaching saturation. Our simulations reveal a "magneto-spin alignment" mechanism that causes magnetized disks and jets to align with the BH spin near BHs and to reorient with the outer disk farther away. This mechanism has implications for the evolution of BH mass and spin, BH feedback on host galaxies, and resolved BH images for the accreting BHs in SgrA* and M87.