Hung-Yi Pu

Energetic Gamma Radiation from Rapidly Rotating Black Holes

Supermassive black holes are believed to be the central power house of active galactic nuclei. Applying the pulsar outer-magnetospheric particle accelerator theory to black-hole magnetospheres, we demonstrate that an electric field is exerted along the magnetic field lines near the event horizon of a rotating black hole. In this particle accelerator (or a gap), electrons and positrons are created by photon-photon collisions and accelerated in the opposite directions by this electric field, efficiently emitting gamma-rays via curvature and inverse-Compton processes. It is shown that a gap arises around the null charge surface formed by the frame-dragging effect, provided that there is no current injection across the gap boundaries. The gap is dissipating a part of the holes rotational energy, and the resultant gamma-ray luminosity increases with decreasing plasma accretion from the surroundings. Considering an extremely rotating supermassive black hole, we show that such a gap reproduces the significant very-high-energy (VHE) gamma-ray flux observed from the radio galaxy IC 310, provided that the accretion rate becomes much less than the Eddington rate particularly during its flare phase. It is found that the curvature process dominates the inverse-Compton process in the magnetosphere of IC~310, and that the observed power-law-like spectrum in VHE gamma-rays can be explained to some extent by a superposition of the curvature emissions with varying curvature radius. It is predicted that the VHE spectrum extends into higher energies with increasing VHE photon flux.

Greenland telescope project: Direct confirmation of black hole with sub‐millimeter VLBI

A 12 m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline >9000 km to achieve an exceptional angular resolution of 20 µas at 350 GHz, which will enable us to resolve the shadow size of ~40 µas. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on 3 years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

STEADY GENERAL RELATIVISTIC MAGNETOHYDRODYNAMIC INFLOW/OUTFLOW SOLUTION ALONG LARGE-SCALE MAGNETIC FIELDS THAT THREAD A ROTATING BLACK HOLE

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.

LEPTON ACCELERATION IN THE VICINITY OF THE EVENT HORIZON: HIGH-ENERGY AND VERY-HIGH-ENERGY EMISSIONS FROM ROTATING BLACK HOLES WITH VARIOUS MASSES

We investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, applying the pulsar outer-gap model to black hole magnetospheres. During a low accretion phase, the radiatively inefficient accretion flow (RIAF) cannot emit enough MeV photons that are needed to sustain the force-free magnetosphere via two-photon collisions. In such a charge-starved region (or a gap), an electric field arises along the magnetic field lines to accelerate electrons and positrons into ultra-relativistic energies. These relativistic leptons emit copious gamma-rays via curvature and inverse-Compton (IC) processes. Some of such gamma-rays collide with the submillimeter-IR photons emitted from the RIAF to materialize as pairs, which polarize to partially screen the original acceleration electric field. It is found that the gap gamma-ray luminosity increases with decreasing accretion rate. However, if the accretion rate decreases too much, the diminished RIAF soft photon field can no longer sustain a stationary pair production within the gap. As long as a stationary gap is formed, the magnetosphere becomes force-free outside the gap by the cascaded pairs, irrespective of the BH mass. If a nearby stellar-mass black hole (BH) is in quiescence, or if a galactic intermediate-mass BH is in a very low accretion state, its curvature and IC emissions are found to be detectable with Fermi/LAT and imaging atmospheric Cherenkov telescopes (IACT). If a low-luminosity active galactic nucleus is located within a few tens of Mpc, the IC emission from its super-massive BH is marginally detectable with IACT.

INDICATION OF THE BLACK HOLE POWERED JET IN M87 BY VSOP OBSERVATIONS

In order to study the collimation and acceleration mechanism of relativistic jets, the jet streamline of M87 at milliarcsecond scale is extensively investigated with images from VSOP observations at 1.6 and 5 GHz. Thanks to the higher angular resolution of VSOP, especially in the direction transverse to the jet, we resolved the jet streamline into three ridgelines at the scale of milli arcseconds. While the properties of the outer two ridgelines are in good agreement with those measured in previous observations and can be expressed by one power-law line with a power law index of 1.7, an inner ridgeline is clearly observed for the first time. We compared the measured size with the outermost streamline expected by Blandford & Znajeks parabolic solutions, which are anchored at the event horizon, with different black hole spin parameters. We revealed that the observed inner ridgeline is narrower than the prediction, suggesting the origin of the inner ridgeline to be part of a spine originating from the spinning black hole. The inner ridgeline becomes very dim at large distances from the central engine at 5 GHz. We considered two possible cases for this; Doppler beaming and/or radiative cooling. Either case seems to be reasonable for its explanation, and future multi-frequency observations will discriminate those two possibilities.

On spin dependence of relativistic acoustic geometry

This work makes the first ever attempt to understand the influence of the black hole background spacetime in determining the fundamental properties of the embedded relativistic acoustic geometry. To accomplish such task, we investigate the role of the spin angular momentum of the astrophysical black hole (the Kerr parameter a—a representative feature of the background black hole metric) in estimating the value of the acoustic surface gravity (the representative feature of the corresponding analogue spacetime). Since almost all astrophysical black holes are supposed to posses some degree of intrinsic rotation, the influence of the Kerr parameter on classical analogue models is very important to understand. We study the general relativistic, axially symmetric, non-self-gravitating inflow of the hydrodynamic fluid onto a rotating astrophysical black hole from the dynamical systems point of view. In this work the location of the acoustic horizon inside such fluid flow is identified and the associated acoustic surface gravity is estimated. We study the dependence of such surface gravity as a function of the Kerr parameter as well as with other dynamical and thermodynamic variables governing the fluid flow under strong gravity, and demonstrate that for retrograde flow, the surface gravity (and hence the associated analogue Hawking temperature) correlates with the black hole spin in general, whereas for the prograde flow, the surface gravity as well as the analogue temperature correlates with the black hole spin for slow to moderately rotating holes, but anti-correlates with the spin for fast to extremely rotating holes. We found that for certain values of the initial boundary conditions, more than one acoustic horizons, namely two black hole types and one white hole type, may form, and the surface gravity may become formally infinite at the acoustic white hole. We discuss the possible connection between the corresponding analogue Hawking temperature and astrophysically relevant observables associated with the spectral signature of the black hole candidates. Our result indicates that the modified dispersion relation evaluated at the close proximity of the acoustic horizon (and hence the nonuniversal feature of Hawking-like effects) is a sensitive function of the spin angular momentum of the astrophysical black hole. We propose that the black hole spin dependence of such dispersion relation may be used to distinguish a corotating flow from a counter rotating flow for axisymmetric accretion onto a Kerr black hole.

Observable Emission Features of Black Hole GRMHD Jets on Event Horizon Scales

The general-relativistic magnetohydrodynamical (GRMHD) formulation for black hole-powered jets naturally gives rise to a stagnation surface, wherefrom inflows and outflows along magnetic field lines that thread the black hole event horizon originate. We derive a conservative formulation for the transport of energetic electrons which are initially injected at the stagnation surface and subsequently transported along flow streamlines. With this formulation the energy spectra evolution of the electrons along the flow in the presence of radiative and adiabatic cooling is determined. For flows regulated by synchrotron radiative losses and adiabatic cooling, the effective radio emission region is found to be finite, and geometrically it is more extended along the jet central axis. Moreover, the emission from regions adjacent to the stagnation surface is expected to be the most luminous as this is where the freshly injected energetic electrons concentrate. An observable stagnation surface is thus a strong prediction of the GRMHD jet model with the prescribed non-thermal electron injection. Future millimeter/sub-millimeter (mm/sub-mm) very-long-baseline interferometric (VLBI) observations of supermassive black hole candidates, such as the one at the center of M87, can verify this GRMHD jet model and its associated non-thermal electron injection mechanism.

Launching and Quenching of Black Hole Relativistic Jets at Low Accretion Rate

Relativistic jets are launched from black hole (BH) X-ray binaries and active galactic nuclei when the disk accretion rate is below a certain limit (i.e., when the ratio of the accretion rate to the Eddingtion accretion rate, , is below about 0.01) but quenched when above. We propose a new paradigm to explain this observed coupling between the jet and the accretion disk by investigating the extraction of the rotational energy of a BH when it is surrounded by different types of accretion disk. At low accretion rates (e.g., when ), the accretion near the event horizon is quasi-spherical. The accreting plasmas fall onto the event horizon in a wide range of latitudes, breaking down the force-free approximation near the horizon. To incorporate the plasma inertia effect, we consider the magnetohydrodynamical (MHD) extraction of the rotational energy from BHs by the accreting MHD fluid, as described by the MHD Penrose process. It is found that the energy extraction operates, and hence a relativistic jet is launched, preferentially when the accretion disk consists of an outer Shakura-Sunyaev disk (SSD) and an inner advection-dominated accretion flow. When the entire accretion disk type changes into an SSD, the jet is quenched because the plasmas bring more rest-mass energy than what is extracted from the hole electromagnetically to stop the extraction. Several other issues related to observed BH disk-jet couplings, such as why the radio luminosity increases with increasing X-ray luminosity until the radio emission drops, are also explained.

Lepton Acceleration in the Vicinity of the Event Horizon: Very High Energy Emissions from Supermassive Black Holes

Around a rapidly rotating black hole (BH), when the plasma accretion rate is much less than the Eddington rate, the radiatively inefficient accretion flow (RIAF) cannot supply enough MeV photons that are capable of materializing as pairs. In such a charge-starved BH magnetosphere, the force-free condition breaks down in the polar funnels. Applying the pulsar outer-magnetospheric lepton accelerator theory to super-massive BHs, we demonstrate that a strong electric field arises along the magnetic field lines in the direct vicinity of the event horizon in the funnels, that the electrons and positrons are accelerated up to 100~TeV in this vacuum gap, and that these leptons emit copious photons via inverse-Compton (IC) process between 0.1~TeV and 30~TeV for a distant observer. It is found that these IC fluxes will be detectable with Imaging Atmospheric Cherenkov Telescopes, provided that a low-luminosity active galactic nucleus is located within 1~Mpc for a million-solar-mass central BH or within 30~Mpc for a billion-solar-mass central BH. These very-high-energy fluxes are beamed in a relatively small solid angle around the rotation axis because of the inhomogeneous and anisotropic distribution of the RIAF photon field, and show an anti-correlation with the RIAF submillimeter fluxes. The gap luminosity little depends on the three-dimensional magnetic-field configuration, because the Goldreich-Julian charge density, and hence the exerted electric field is essentially governed by the frame-dragging effect, not by the magnetic field configuration.

STRUCTURAL TRANSITION IN THE NGC 6251 JET: AN INTERPLAY WITH THE SUPERMASSIVE BLACK HOLE AND ITS HOST GALAXY

The structure of the NGC 6251 jet on the milliarcsecond scale is investigated using images taken with the European VLBI Network and the Very Long Baseline Array. We detect a structural transition of the jet from a parabolic to a conical shape at a distance of (1-2) x 10^5 times the Schwarzschild radius from the central engine, which is close to the sphere of gravitational influence of the supermassive black hole (SMBH). We also examine the jet pressure profiles with the synchrotron minimum energy assumption to discuss the physical origin of the structural transition. The NGC 6251 jet, together with the M87 jet, suggests a fundamental process of structural transition in the jets of active galactic nuclei (AGNs). Collimated AGN jets are characterized by their external galactic medium, showing that AGN jets interplay with the SMBH and its host galaxy.