Showing 1–47 of 47 results for author: Pu, H

Discovery of Limb-Brightening in the Parsec-Scale Jet of NGC 315 through Global VLBI Observations and Its Implications for Jet Models

Authors: Jongho Park, Guang-Yao Zhao, Masanori Nakamura, Yosuke Mizuno, Hung-Yi Pu, Keiichi Asada, Kazuya Takahashi, Kenji Toma, Motoki Kino, Ilje Cho, Kazuhiro Hada, Phil G. Edwards, Hyunwook Ro, Minchul Kam, Kunwoo Yi, Yunjeong Lee, Shoko Koyama, Do-Young Byun, Chris Phillips, Cormac Reynolds, Jeffrey A. Hodgson, Sang-Sung Lee

Abstract: We report the first observation of the nearby giant radio galaxy NGC 315 using a global VLBI array consisting of 22 radio antennas located across five continents, including high-sensitivity stations, at 22 GHz. Utilizing the extensive $(u,v)$-coverage provided by the array, coupled with the application of a recently developed super-resolution imaging technique based on the regularized maximum like… ▽ More We report the first observation of the nearby giant radio galaxy NGC 315 using a global VLBI array consisting of 22 radio antennas located across five continents, including high-sensitivity stations, at 22 GHz. Utilizing the extensive $(u,v)$-coverage provided by the array, coupled with the application of a recently developed super-resolution imaging technique based on the regularized maximum likelihood method, we were able to transversely resolve the NGC 315 jet at parsec scales for the first time. Previously known for its central ridge-brightened morphology at similar scales in former VLBI studies, the jet now clearly exhibits a limb-brightened structure. This finding suggests an inherent limb-brightening that was not observable before due to limited angular resolution. Considering that the jet is viewed at an angle of $\sim50^\circ$, the observed limb-brightening is challenging to reconcile with the magnetohydrodynamic models and simulations, which predict that the Doppler-boosted jet edges should dominate over the non-boosted central layer. The conventional jet model that proposes a fast spine and a slow sheath with uniform transverse emissivity may pertain to our observations. However, in this model, the relativistic spine would need to travel at speeds of $Γ\gtrsim6.0-12.9$ along the de-projected jet distance of (2.3-10.8) $\times 10^3$ gravitational radii from the black hole. We propose an alternative scenario that suggests higher emissivity at the jet boundary layer, resulting from more efficient particle acceleration or mass loading onto the jet edges, and consider prospects for future observations with even higher angular resolution. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: 25 pages, 12 figures, 1 table, accepted for publication in the Astrophysical Journal Letters

Observational features of reflection asymmetric black holes

Authors: Che-Yu Chen, Hung-Yi Pu

Abstract: The Kerr spacetime is symmetric with respect to a well-defined equatorial plane. When testing the equatorial reflection symmetry of an isolated black hole, one is at the same time testing the Kerr hypothesis in General Relativity. In this work, we investigate the possible observational features when a Keplerian disk is surrounding a rotating black hole without reflection symmetry. When such symmet… ▽ More The Kerr spacetime is symmetric with respect to a well-defined equatorial plane. When testing the equatorial reflection symmetry of an isolated black hole, one is at the same time testing the Kerr hypothesis in General Relativity. In this work, we investigate the possible observational features when a Keplerian disk is surrounding a rotating black hole without reflection symmetry. When such symmetry is broken, generically, the photon trajectories around the black hole and the Keplerian orbits on the accretion disk are distorted vertically away from the equatorial plane by an amount that depends on their distance to the black hole. In the reflection asymmetric spacetime we are considering, these two kinds of orbits are distorted in opposite directions. Interestingly, while the size and shape of black hole shadows closely resemble those of Kerr black holes, distinct observational characteristics can emerge in the disk image and emission line profiles. When observing the disk edge-on, a pronounced concave shape may appear along its innermost edge on the incoming side. Furthermore, distinctive horn-like features might be observed on the spectral line profile at the blue-shifted side. These special features can serve as compelling indicators of the reflection asymmetry present in rotating black holes. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 14 pages, 7 figures

Report number: RIKEN-iTHEMS-Report-24

Super-Eddington accretion in high-redshift black holes and the emergence of jetted AGN

Authors: Olmo Piana, Hung-Yi Pu, Kinwah Wu

Abstract: In this work we study the co-evolution of central black holes (BHs) and host galaxies by utilizing an advanced iteration of the DELPHI semi-analytical model of galaxy formation and evolution. Based on dark matter halo merger trees spanning the redshift range from $z=20$ to $z=4$, it now incorporates essential components such as gas heating and cooling, cold and hot BH accretion, jet and radiative… ▽ More In this work we study the co-evolution of central black holes (BHs) and host galaxies by utilizing an advanced iteration of the DELPHI semi-analytical model of galaxy formation and evolution. Based on dark matter halo merger trees spanning the redshift range from $z=20$ to $z=4$, it now incorporates essential components such as gas heating and cooling, cold and hot BH accretion, jet and radiative AGN feedback. We show how different BH growth models impact quasar and galaxy observables at $z \geq 5$, providing predictions that will help discriminate between super-Eddington and Eddington-limited accretion models: despite being both consistent with observed properties of SMBHs and their host galaxies at $z \sim 5-7$, they become very clearly distinguishable at higher redshift and in the intermediate mass regime. We find that the super-Eddington model, unlike the Eddington-limited scenario, predicts a gap in the BH mass function corresponding to the intermediate-mass range $10^4 \mathrm{M_\odot} < M_\mathrm{bh} < 10^6 \mathrm{M_\odot}$. Additionally, it predicts black holes up to two orders of magnitude more massive for the same stellar mass at $z=9$. The resulting velocity dispersion - BH mass relation at $z \geq 5$ is consistent with local measurements, suggesting that its slope and normalisation are independent of redshift. Depending on the Eddington ratio, we also model the emergence of AGN jets, predicting their duty cycle across as a function of BH mass and their potential impact on the observed number density distribution of high-redshift AGN in the hard X-ray band. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Absolute Flux Density Calibration of the Greenland Telescope Data for Event Horizon Telescope Observations

Authors: J. Y. Koay, K. Asada, S. Matsushita, C. -Y. Kuo, C. -W. L. Huang, C. Romero-Cañizales, S. Koyama, J. Park, W. -P. Lo, G. Bower, M. -T. Chen, S. -H. Chang, C. -C. Chen, R. Chilson, C. C. Han, P. T. P. Ho, Y. -D. Huang, M. Inoue, B. Jeter, H. Jiang, P. M. Koch, D. Kubo, C. -T. Li, C. -T. Liu, K. -Y. Liu , et al. (13 additional authors not shown)

Abstract: Starting from the observing campaign in April 2018, the Greenland Telescope (GLT) has been added as a new station of the Event Horizon Telescope (EHT) array. Visibilities on baselines to the GLT, particularly in the North-South direction, potentially provide valuable new constraints for the modeling and imaging of sources such as M87*. The GLT's location at high Northern latitudes adds unique chal… ▽ More Starting from the observing campaign in April 2018, the Greenland Telescope (GLT) has been added as a new station of the Event Horizon Telescope (EHT) array. Visibilities on baselines to the GLT, particularly in the North-South direction, potentially provide valuable new constraints for the modeling and imaging of sources such as M87*. The GLT's location at high Northern latitudes adds unique challenges to its calibration strategies. Additionally, the performance of the GLT was not optimal during the 2018 observations due to it being only partially commissioned at the time. This document describes the steps taken to estimate the various parameters (and their uncertainties) required for the absolute flux calibration of the GLT data as part of the EHT. In particular, we consider the non-optimized status of the GLT in 2018, as well as its improved performance during the 2021 EHT campaign. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 17 pages, 4 figures, EHT Memo Series 2023-L1-02

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

The Greenland Telescope: Construction, Commissioning, and Operations in Pituffik

Authors: Ming-Tang Chen, Keiichi Asada, Satoki Matsushita, Philippe Raffin, Makoto Inoue, Paul T. P. Ho, Chih-Chiang Han, Derek Kubo, Timothy Norton, Nimesh A. Patel, George Nystrom, Chih-Wei L. Huang, Pierre Martin-Cocher, Jun Yi Koay, Cristina Romero-Cañizales, Ching-Tang Liu, Teddy Huang, Kuan-Yu Liu, Tashun Wei, Shu-Hao Chang, Ryan Chilson, Peter Oshiro, Homin Jiang, Chao-Te Li, Geoffrey Bower , et al. (29 additional authors not shown)

Abstract: In 2018, the Greenland Telescope (GLT) started scientific observation in Greenland. Since then, we have completed several significant improvements and added new capabilities to the telescope system. This paper presents a full review of the GLT system, a summary of our observation activities since 2018, the lessons learned from the operations in the Arctic regions, and the prospect of the telescope… ▽ More In 2018, the Greenland Telescope (GLT) started scientific observation in Greenland. Since then, we have completed several significant improvements and added new capabilities to the telescope system. This paper presents a full review of the GLT system, a summary of our observation activities since 2018, the lessons learned from the operations in the Arctic regions, and the prospect of the telescope. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 26 pages, 11 figures, and 8 tables. This is the version of the article before publication editing, as submitted by an author to Publications of the Astronomical Society of the Pacific. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record will be added when it becomes available

i-SPin 2: An integrator for general spin-s Gross-Pitaevskii systems

Authors: Mudit Jain, Mustafa A. Amin, Han Pu

Abstract: We provide an algorithm for evolving general spin-$s$ Gross-Pitaevskii / non-linear Schrödinger systems carrying a variety of interactions, where the $2s+1$ components of the `spinor' field represent the different spin-multiplicity states. We consider many nonrelativistic interactions up to quartic order in the Schrödinger field (both short and long-range, and spin-dependent and spin-independent i… ▽ More We provide an algorithm for evolving general spin-$s$ Gross-Pitaevskii / non-linear Schrödinger systems carrying a variety of interactions, where the $2s+1$ components of the `spinor' field represent the different spin-multiplicity states. We consider many nonrelativistic interactions up to quartic order in the Schrödinger field (both short and long-range, and spin-dependent and spin-independent interactions), including explicit spin-orbit couplings. The algorithm allows for spatially varying external and/or self-generated vector potentials that couple to the spin density of the field. Our work can be used for scenarios ranging from laboratory systems such as spinor Bose-Einstein condensates (BECs), to cosmological/astrophysical systems such as self-interacting bosonic dark matter. As examples, we provide results for two different setups of spin-$1$ BECs that employ a varying magnetic field and spin-orbit coupling, respectively, and also collisions of spin-$1$ solitons in dark matter. Our symplectic algorithm is second-order accurate in time, and is extensible to the known higher-order accurate methods. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: 13 pages, 3 figures, 2 appendices

Journal ref: Phys. Rev. E 108, 055305, 15 November 2023

A ring-like accretion structure in M87 connecting its black hole and jet

Authors: Ru-Sen Lu, Keiichi Asada, Thomas P. Krichbaum, Jongho Park, Fumie Tazaki, Hung-Yi Pu, Masanori Nakamura, Andrei Lobanov, Kazuhiro Hada, Kazunori Akiyama, Jae-Young Kim, Ivan Marti-Vidal, José L. Gómez, Tomohisa Kawashima, Feng Yuan, Eduardo Ros, Walter Alef, Silke Britzen, Michael Bremer, Avery E. Broderick, Akihiro Doi, Gabriele Giovannini, Marcello Giroletti, Paul T. P. Ho, Mareki Honma , et al. (96 additional authors not shown)

Abstract: The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the comp… ▽ More The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of 8.4_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 50 pages, 18 figures, 3 tables, author's version of the paper published in Nature

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

The Photon Ring in M87*

Authors: Avery E. Broderick, Dominic W. Pesce, Paul Tiede, Hung-Yi Pu, Roman Gold, Richard Anantua, Silke Britzen, Chiara Ceccobello, Koushik Chatterjee, Yongjun Chen, Nicholas S. Conroy, Geoffrey B. Crew, Alejandro Cruz-Osorio, Yuzhu Cui, Sheperd S. Doeleman, Razieh Emami, Joseph Farah, Christian M. Fromm, Peter Galison, Boris Georgiev, Luis C. Ho, David J. James, Britton Jeter, Alejandra Jimenez-Rosales, Jun Yi Koay , et al. (26 additional authors not shown)

Abstract: We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical… ▽ More We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Journal ref: ApJ, 935, 61 (2022)

Hybrid Very Long Baseline Interferometry Imaging and Modeling with Themis

Authors: Avery E. Broderick, Dominic W. Pesce, Paul Tiede, Hung-Yi Pu, Roman Gold

Abstract: Generating images from very long baseline interferometric observations poses a difficult, and generally not unique, inversion problem. This problem is simplified by the introduction of constraints, some generic (e.g., positivity of the intensity) and others motivated by physical considerations (e.g., smoothness, instrument resolution). It is further complicated by the need to simultaneously addres… ▽ More Generating images from very long baseline interferometric observations poses a difficult, and generally not unique, inversion problem. This problem is simplified by the introduction of constraints, some generic (e.g., positivity of the intensity) and others motivated by physical considerations (e.g., smoothness, instrument resolution). It is further complicated by the need to simultaneously address instrumental systematic uncertainties and sparse coverage in the u-v plane. We report a new Bayesian image reconstruction technique in the parameter estimation framework Themis that has been developed for the Event Horizon Telescope. This has two key features: first, the full Bayesian treatment of the image reconstruction makes it possible to generate a full posterior for the images, permitting a rigorous and quantitative investigation into the statistical significance of image features. Second, it is possible to seamlessly incorporate directly modeled features simultaneously with image reconstruction. We demonstrate this second capability by incorporating a narrow, slashed ring in reconstructions of simulated M87 data in an attempt to detect and characterize the photon ring. We show that it is possible to obtain high-fidelity photon ring sizes, enabling mass measurements with accuracies of 2%-5% that are essentially insensitive to astrophysical uncertainties, and creating opportunities for precision tests of general relativity. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Journal ref: ApJ, 898, 9 (2020)

Black Hole Shadow with Soft Hair

Authors: Feng-Li Lin, Avani Patel, Hung-Yi Pu

Abstract: Light bending by the strong gravity around the black hole will form the so-called black hole shadow, the shape of which can shed light on the structure of the near-horizon geometry to possibly reveal novel physics of strong gravity and black hole. In this work, we adopt both analytical and ray-tracing methods to study the black hole shadow in the presence of the infrared structure of gravity theor… ▽ More Light bending by the strong gravity around the black hole will form the so-called black hole shadow, the shape of which can shed light on the structure of the near-horizon geometry to possibly reveal novel physics of strong gravity and black hole. In this work, we adopt both analytical and ray-tracing methods to study the black hole shadow in the presence of the infrared structure of gravity theory, which manifests the asymptotic symmetries of spacetime as the supertranslation soft hairs of the black hole. Though the black hole metrics with and without the soft hair are related by large gauge transformations, the near horizon geometries relevant for the shape of the shadow are quite different. Moreover, the Hamiltonian for the geodesic seems intrinsically different, i.e., the loss of separability due to the breaking of spherical symmetry by soft hair. By applying ray-tracing computations, we find that the soft hair, although not affecting the shape of the shadow, may change the average size and position of the shadow. Images resulting from soft hair black holes with surrounding accretion flows are also discussed. △ Less

Submitted 16 August, 2022; v1 submitted 28 February, 2022; originally announced February 2022.

Comments: 28 pages, 11 figures, including analytical calculation of shadow without any change in final conclusion, matches with the version accepted by JHEP

Collimation of the relativistic jet in the quasar 3C 273

Authors: Hiroki Okino, Kazunori Akiyama, Keiichi Asada, José L. Gómez, Kazuhiro Hada, Mareki Honma, Thomas P. Krichbaum, Motoki Kino, Hiroshi Nagai, Uwe Bach, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Geoffrey B. Crew, Sheperd S. Doeleman, Vincent L. Fish, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Svetlana Jorstad, Shoko Koyama, Colin J. Lonsdale, Ru-sen Lu, Ivan Martí-Vidal, Lynn D. Matthews , et al. (10 additional authors not shown)

Abstract: The collimation of relativistic jets launched from the vicinity of supermassive black holes (SMBHs) at the centers of active galactic nuclei (AGN) is one of the key questions to understand the nature of AGN jets. However, little is known about the detailed jet structure for AGN like quasars since very high angular resolutions are required to resolve these objects. We present very long baseline int… ▽ More The collimation of relativistic jets launched from the vicinity of supermassive black holes (SMBHs) at the centers of active galactic nuclei (AGN) is one of the key questions to understand the nature of AGN jets. However, little is known about the detailed jet structure for AGN like quasars since very high angular resolutions are required to resolve these objects. We present very long baseline interferometry (VLBI) observations of the archetypical quasar 3C 273 at 86 GHz, performed with the Global Millimeter VLBI Array, for the first time including the Atacama Large Millimeter/submillimeter Array. Our observations achieve a high angular resolution down to $\sim$60 ${\rm μ}$as, resolving the innermost part of the jet ever on scales of $\sim 10^5$ Schwarzschild radii. Our observations, including close-in-time High Sensitivity Array observations of 3C 273 at 15, 22, and 43 GHz, suggest that the inner jet collimates parabolically, while the outer jet expands conically, similar to jets from other nearby low luminosity AGN. We discovered the jet collimation break around $10^{7}$ Schwarzschild radii, providing the first compelling evidence for structural transition in a quasar jet. The location of the collimation break for 3C 273 is farther downstream the sphere of gravitational influence (SGI) from the central SMBH. With the results for other AGN jets, our results show that the end of the collimation zone in AGN jets is governed not only by the SGI of the SMBH but also by the more diverse properties of the central nuclei. △ Less

Submitted 7 October, 2022; v1 submitted 22 December, 2021; originally announced December 2021.

Comments: 22 pages, 13 figures, accepted for publication in ApJ

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

Relativistic jet acceleration region in a black hole magnetosphere

Authors: Masaaki Takahashi, Motoki Kino, Hung-Yi Pu

Abstract: We discuss stationary and axisymmetric trans-magnetosonic outflows in the magnetosphere of a rotating black hole (BH). Ejected plasma from the plasma source located near the BH is accelerated far away to form a relativistic jet. In this study, the plasma acceleration efficiency and conversion of fluid energy from electromagnetic energy are considered by employing the trans-fast magnetosonic flow s… ▽ More We discuss stationary and axisymmetric trans-magnetosonic outflows in the magnetosphere of a rotating black hole (BH). Ejected plasma from the plasma source located near the BH is accelerated far away to form a relativistic jet. In this study, the plasma acceleration efficiency and conversion of fluid energy from electromagnetic energy are considered by employing the trans-fast magnetosonic flow solution derived by Takahashi & Tomimatsu (2008). Considering the parameter dependence of magnetohydrodynamical flows, we search for the parameters of the trans-magnetosonic outflow solution to the recent M87 jet observations and obtain the angular velocity values of the magnetic field line and angular momentum of the outflow in the magnetized jet flow. Therefore, we estimate the locations of the outer light surface, Alfvén surface, and separation surface of the flow. We also discuss the electromagnetic energy flux from the rotating BH (i.e., the Blandford-Znajek process), which suggests that the energy extraction mechanism is effective for the M87 relativistic jet. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 16 pages, 7 figures

A revised view of the linear polarization in the subparsec core of M87 at 7 mm

Authors: Jongho Park, Keiichi Asada, Masanori Nakamura, Motoki Kino, Hung-Yi Pu, Kazuhiro Hada, Evgeniya V. Kravchenko, Marcello Giroletti

Abstract: The linear polarization images of the jet in the giant elliptical galaxy M87 have previously been observed with Very Long Baseline Array at 7 mm. They exhibit a complex polarization structure surrounding the optically thick and compact subparsec-scale core. However, given the low level of linear polarization in the core, it is required to verify that this complex structure does not originate from… ▽ More The linear polarization images of the jet in the giant elliptical galaxy M87 have previously been observed with Very Long Baseline Array at 7 mm. They exhibit a complex polarization structure surrounding the optically thick and compact subparsec-scale core. However, given the low level of linear polarization in the core, it is required to verify that this complex structure does not originate from residual instrumental polarization signals in the data. We have performed a new analysis of the same data sets observed in four epochs by using the Generalized Polarization CALibration pipeline (GPCAL). This novel instrumental polarization calibration pipeline overcomes the limitations of LPCAL, a conventional calibration tool used in the previous M87 studies. The resulting images show a compact linear polarization structure with its peak nearly coincident with the total intensity peak, which is significantly different from the results of previous studies. The core linear polarization is characterized as fractional polarization of $\sim0.2-0.6$% and polarization angles of $\sim66-92^\circ$, showing moderate variability. We demonstrate that, based on tests with synthetic data sets, LPCAL using calibrators having complex polarization structures cannot achieve sufficient calibration accuracy to obtain the true polarization image of M87 due to a breakdown of the "similarity approximation". We find that GPCAL obtains more accurate D-terms than LPCAL by using observed closure traces of calibrators that are insensitive to both antenna gain and polarization leakage corruptions. This study suggests that polarization imaging of very weakly polarized sources has become possible with the advanced instrumental polarization calibration techniques. △ Less

Submitted 14 September, 2021; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: 22 pages, 13 figures, 4 tables, accepted for publication in ApJ

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

Properties of Trans-fast Magnetosonic Jets in Black Hole Magnetospheres

Authors: Hung-Yi Pu, Masaaki Takahashi

Abstract: Traveling across several order of magnitude in distance, relativistic jets from strong gravity region to asymptotic flat spacetime region are believed to consist of several general relativistic magnetohydrodynamic (GRMHD) processes. We present a semi-analytical approach for modeling the global structures of a trans-fast magnetosonic relativistic jet, which should be ejected from a plasma source ne… ▽ More Traveling across several order of magnitude in distance, relativistic jets from strong gravity region to asymptotic flat spacetime region are believed to consist of several general relativistic magnetohydrodynamic (GRMHD) processes. We present a semi-analytical approach for modeling the global structures of a trans-fast magnetosonic relativistic jet, which should be ejected from a plasma source nearby a black hole in a funnel region enclosed by dense accreting flow and also disk corona around the black hole. Our model consistently includes the inflow and outflow part of the GRMHD solution along the magnetic field lines penetrating the black hole horizon. After the rotational energy of the black hole is extracted electromagnetically by the negative energy GRMHD inflow, the huge electromagnetic energy flux then propagates from the inflow to the outflow region across the plasma source, and in the outflow region the electromagnetic energy converts to the fluid kinetic energy. Eventually, the accelerated outflow must exceed the fast-magnetosonic wave speed. We apply the semi-analytical trans-fast magnetosonic flow model to the black hole magnetosphere for both parabolic and split-monopole magnetic field configurations, and discuss the general flow properties; that is, jet acceleration, jet magnetization, and the locations of some characteristic surfaces of the black hole magnetosphere. We have confirmed that, at large distance, the GRMHD jet solutions are in good agreement with the previously known trans-fast special relativistic magnetohydrodynamic (SRMHD) jet properties, as expected. The flexibility of the model provides a prompt and heuristic way to approximate the global GRMHD trans-fast magnetosonic jet properties. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: accpeted for publication in the Astrophysical Journal

Spacetime Tomography Using The Event Horizon Telescope

Authors: Paul Tiede, Hung-Yi Pu, Avery E. Broderick, Roman Gold, Mansour Karami, Jorge A. Preciado-López

Abstract: We have now entered the new era of high-resolution imaging astronomy with the beginning of the Event Horizon Telescope (EHT). The EHT can resolve the dynamics of matter in the immediate vicinity around black holes at and below the horizon scale. One of the candidate black holes, Sagittarius A* flares 1\-4 times a day depending on the wavelength. A possible interpretation of these flares could be h… ▽ More We have now entered the new era of high-resolution imaging astronomy with the beginning of the Event Horizon Telescope (EHT). The EHT can resolve the dynamics of matter in the immediate vicinity around black holes at and below the horizon scale. One of the candidate black holes, Sagittarius A* flares 1\-4 times a day depending on the wavelength. A possible interpretation of these flares could be hotspots generated through magnetic reconnection events in the accretion flow. In this paper, we construct a semi-analytical model for hotspots that include the effects of shearing as a spot moves along the accretion flow. We then explore the ability of the EHT to recover these hotspots. Even including significant systematic uncertainties, such as thermal noise, diffractive scattering, and background emission due to an accretion disk, we were able to recover the hotspots and spacetime structure to sub-percent precision. Moreover, by observing multiple flaring events we show how the EHT could be used to tomographically map spacetime. This provides new avenues for testing relativistic fluid dynamics and general relativity near the event horizon of supermassive black holes. △ Less

Submitted 17 February, 2020; v1 submitted 13 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ (17 pages, 15 figures) (revised edition corrected some affiliations)

Parabolic Jets from the Spinning Black Hole in M87

Authors: Masanori Nakamura, Keiichi Asada, Kazuhiro Hada, Hung-Yi Pu, Scott Noble, Chihyin Tseng, Kenji Toma, Motoki Kino, Hiroshi Nagai, Kazuya Takahashi, Juan-Carlos Algaba, Monica Orienti, Kazunori Akiyama, Akihiro Doi, Gabriele Giovannini, Marcello Giroletti, Mareki Honma, Shoko Koyama, Rocco Lico, Kotaro Niinuma, Fumie Tazaki

Abstract: The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determin… ▽ More The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 29 pages, 23 figures, accepted for publication in ApJ

High-energy and very-high-energy emission from stellar-mass black holes moving in gaseous clouds

Authors: Kouichi Hirotani, Hung-Yi Pu, Sabrina Outmani, Hsinhao Huang, Dawoon Kim, Yoogeun Song, Satoki Matsushita, Albert K. H Kong

Abstract: We investigate the electron-positron pair cascade taking place in the magnetosphere of a rapidly rotating black hole. Because of the spacetime frame dragging, the Goldreich-Julian charge density changes sign in the vicinity of the event horizon, which leads to an occurrence of a magnetic-field aligned electric field, in the same way as the pulsar outer-magnetospheric accelerator. In this lepton ac… ▽ More We investigate the electron-positron pair cascade taking place in the magnetosphere of a rapidly rotating black hole. Because of the spacetime frame dragging, the Goldreich-Julian charge density changes sign in the vicinity of the event horizon, which leads to an occurrence of a magnetic-field aligned electric field, in the same way as the pulsar outer-magnetospheric accelerator. In this lepton accelerator, electrons and positrons are accelerated in the opposite directions, to emit copious gamma-rays via the curvature and inverse-Compton processes. We examine a stationary pair cascade, and show that a stellar-mass black hole moving in a gaseous cloud can emit a detectable very-high-energy flux, provided that the black hole is extremely rotating and that the distance is less than about 1 kpc. We argue that the gamma-ray image will have a point-like morphology, and demonstrate that their gamma-ray spectra have a broad peak around 0.01-1 GeV and a sharp peak around 0.1 TeV, that the accelerators become most luminous when the mass accretion rate becomes about 0.01% of the Eddington rate, and that the predicted gamma-ray flux little changes in a wide range of magnetospheric currents. An implication of the stability of such a stationary gap is discussed. △ Less

Submitted 25 September, 2018; originally announced September 2018.

Comments: 26 pages, 9 figures; accepted for publication in the Astrophysical Journal

Lightning black holes as unidentified TeV sources

Authors: Kouichi Hirotani, Hung-Yi Pu, Satoki Matsushita

Abstract: Imaging Atmospheric Cherenkov Telescopes have revealed more than 100 TeV sources along the Galactic Plane, around 45% of them remain unidentified. However, radio observations revealed that dense molecular clumps are associated with 67% of 18 unidentified TeV sources. In this paper, we propose that an electron-positron magnetospheric accelerator emits detectable TeV gamma-rays when a rapidly rotati… ▽ More Imaging Atmospheric Cherenkov Telescopes have revealed more than 100 TeV sources along the Galactic Plane, around 45% of them remain unidentified. However, radio observations revealed that dense molecular clumps are associated with 67% of 18 unidentified TeV sources. In this paper, we propose that an electron-positron magnetospheric accelerator emits detectable TeV gamma-rays when a rapidly rotating black hole enters a gaseous cloud. Since the general-relativistic effect plays an essential role in this magnetospheric lepton accelerator scenario, the emissions take place in the direct vicinity of the event horizon, resulting in a point-like gamma-ray image. We demonstrate that their gamma-ray spectra have two peaks around 0.1 GeV and 0.1 TeV and that the accelerators become most luminous when the mass accretion rate becomes about 0.01% of the Eddington accretion rate. We compare the results with alternative scenarios such as the cosmic-ray hadron scenario, which predicts an extended morphology of the gamma-ray image with a single power-law photon spectrum from GeV to 100 TeV. △ Less

Submitted 9 August, 2018; originally announced August 2018.

Comments: Accepted for publication in Journal of Astrophysics and Astronomy, 12 pages, 4 figures

Probing the innermost accretion flow geometry of Sgr A* with Event Horizon Telescope

Authors: Hung-Yi Pu, Avery E. Broderick

Abstract: Upcoming Event Horizon Telescope (EHT) observations will provide a unique opportunity to reveal the innermost region of the radiative inefficient accretion flow (RIAF) around the Galactic black hole, Sgr A*. Depending on the flow dynamics and accumulated magnetic flux, the innermost region of an RIAF could have a quasi-spherical or disk-like geometry. Here we present a phenomenological model to in… ▽ More Upcoming Event Horizon Telescope (EHT) observations will provide a unique opportunity to reveal the innermost region of the radiative inefficient accretion flow (RIAF) around the Galactic black hole, Sgr A*. Depending on the flow dynamics and accumulated magnetic flux, the innermost region of an RIAF could have a quasi-spherical or disk-like geometry. Here we present a phenomenological model to investigate the characteristics of the black hole shadow images with different flow geometries, together with the effect of black hole spin and flow dynamics. The resulting image consists in general of two major components: a crescent, which may surround the funnel region of the black hole or the black hole itself, and a photon ring, which may be partially luminous and overlapped with the crescent component. Compared to a quasi-spherical flow case, a disk-like flow in the vicinity of a black hole exhibits the following image features: (i) due to less material near the funnel region, the crescent structure has a smaller size, and (ii) due to the combination of emission from the flow beside and behind the black hole, the crescent structure has a more irregular shape, and a less smooth brightness distribution. How these features can result in different observables for EHT observations is discussed. △ Less

Submitted 4 September, 2018; v1 submitted 4 July, 2018; originally announced July 2018.

Comments: 11 pages, 8 figures; published in the ApJ

Superresolution Interferometric Imaging with Sparse Modeling Using Total Squared Variation --- Application to Imaging the Black Hole Shadow

Authors: Kazuki Kuramochi, Kazunori Akiyama, Shiro Ikeda, Fumie Tazaki, Vincent L. Fish, Hung-Yi Pu, Keiichi Asada, Mareki Honma

Abstract: We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super… ▽ More We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super-resolution of $\sim 30$% compared to the traditional CLEAN beam size using synthetic observations of two point sources. Second, we present simulated observations of three physically motivated static models of Sgr A* with the Event Horizon Telescope (EHT) to show the performance of proposed techniques in greater detail. We find that $\ell_1$+TSV regularization outperforms $\ell_1$+TV regularization with the popular isotropic TV term and the Cotton-Schwab CLEAN algorithm, demonstrating that TSV is well-matched to the expected physical properties of the astronomical images, which are often nebulous. Remarkably, in both the image and gradient domains, the optimal beam size minimizing root-mean-squared errors is $\lesssim 10$% of the traditional CLEAN beam size for $\ell_1$+TSV regularization, and non-convolved reconstructed images have smaller errors than beam-convolved reconstructed images. This indicates that the traditional post-processing technique of Gaussian convolution in interferometric imaging may not be required for the $\ell_1$+TSV regularization. We also propose a feature extraction method to detect circular features from the image of a black hole shadow with the circle Hough transform (CHT) and use it to evaluate the performance of the image reconstruction. With our imaging technique and the CHT, the EHT can constrain the radius of the black hole shadow with an accuracy of $\sim 10-20$% in present simulations for Sgr A*. △ Less

Submitted 13 March, 2018; v1 submitted 15 February, 2018; originally announced February 2018.

Comments: 18 pages, 7 figures, ApJ in Press

Searching for High Energy, Horizon-scale Emissions from Galactic Black Hole Transients during Quiescence

Authors: Lupin Chun-Che Lin, Hung-Yi Pu, Kouichi Hirotani, Albert K. H Kong, Satoki Matsushita, Hsiang-Kuang Chang, Makoto Inoue, Pak-Hin T. Tam

Abstract: We search for the gamma-ray counterparts of stellar-mass black holes using long-term Fermi archive to investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, by applying the pulsar outer-gap model to their magnetosphere. When a black hole transient (BHT) is in a low-hard or quiescent state, the radiatively inefficient accretion flow cannot emit e… ▽ More We search for the gamma-ray counterparts of stellar-mass black holes using long-term Fermi archive to investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, by applying the pulsar outer-gap model to their magnetosphere. When a black hole transient (BHT) is in a low-hard or quiescent state, the radiatively inefficient accretion flow cannot emit enough MeV photons that are required to sustain the force-free magnetosphere in the polar funnel via two-photon collisions. In this charge-starved gap region, 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 the curvature and inverse-Compton (IC) processes. It is found that these gamma-ray emissions exhibit a flaring activity when the plasma accretion rate stays typically between 0.01 and 0.005 percent of the Eddington value for rapidly rotating, stellar-mass black holes. By analyzing the detection limit determined from archival Fermi/LAT data, we find that the 7-year averaged duty cycle of such flaring activities should be less than 5% and 10% for XTE J1118+480 and 1A 0620-00, respectively, and that the detection limit is comparable to the theoretical prediction for V404 Cyg. It is predicted that the gap emission can be discriminated from the jet emission, if we investigate the high-energy spectral behaviour or observe nearby BHTs during deep quiescence simultaneously in infrared wavelength and very-high energies. △ Less

Submitted 27 July, 2017; originally announced July 2017.

Comments: 7 pages, 3 figures; accepted for publication in ApJ

Enhanced gamma radiation toward the rotation axis from the immediate vicinity of extremely rotating black holes

Authors: Yoogeun Song, Hung-Yi Pu, Kouichi Hirotani, Satoki Matsushita, Albert K. H. Kong, Hsiang-Kuang Chang

Abstract: We investigate the acceleration of electrons and positrons by magnetic-field-aligned electric fields in the polar funnel of an accreting black hole (BH). Applying the pulsar outer-gap theory to BH magnetospheres, we find that such a lepton accelerator arises in the immediate vicinity of the event horizon due to frame-dragging, and that their gamma-ray luminosity increases with decreasing accretion… ▽ More We investigate the acceleration of electrons and positrons by magnetic-field-aligned electric fields in the polar funnel of an accreting black hole (BH). Applying the pulsar outer-gap theory to BH magnetospheres, we find that such a lepton accelerator arises in the immediate vicinity of the event horizon due to frame-dragging, and that their gamma-ray luminosity increases with decreasing accretion rate. Furthermore, we demonstrate that the gamma-ray flux is enhanced along the rotation axis by more than an order of magnitude if the BH spin increases from $a=0.90M$ to $a=0.9999M$. As a result, if a ten-solar-mass, almost-maximally rotating BH is located within 3 kpc, when its accretion rate is between 0.005% and 0.01% of the Eddington rate, its high-energy flare becomes detectable with the Fermi/Large Area Telescope, provided that the flare lasts longer than 1.2 months and that we view the source nearly along the rotation axis. In addition, its very-high-energy flux is marginally detectable with the Cherenkov Telescope Array, provided that the flare lasts longer than a night and that our viewing angle is about 45 degrees with respect to the rotation axis. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: 5 figures

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

Authors: Hung-Yi Pu, Kinwah Wu, Ziri Younsi, Keiichi Asada, Yosuke Mizuno, Masanori Nakamura

Abstract: 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… ▽ More 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. △ Less

Submitted 21 July, 2017; originally announced July 2017.

Comments: 19 pages, 12 figure; accepted for publication in ApJ

Lepton acceleration in the vicinity of the event horizon: Very-high-energy emissions from super-massive black holes

Authors: Kouichi Hirotani, Hung-Yi Pu, Lupin Chun-Che Lin, Albert K. H Kong, Satoki Matsushita, Keiichi Asada, Hsiang-Kuang Chang, Pak-Hin T. Tam

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: 15 pages, 16 figures; Accepted for publication in the Astrophysical Journal. arXiv admin note: text overlap with arXiv:1610.07819

White Paper on East Asian Vision for mm/submm VLBI: Toward Black Hole Astrophysics down to Angular Resolution of 1~R$_{S}$

Authors: K. Asada, M. Kino, M. Honma, T. Hirota, R. -S. Lu, M. Inoue, B. -W. Sohn, Z. -Q. Shen, P. T. P. Ho, K. Akiyama, J-C. Algaba, T. An, G. Bower, D-Y. Byun, R. Dodson, A. Doi, P. G. Edwards, K. Fujisawa, M-F. Gu, K. Hada, Y. Hagiwara, P. Jaroenjittichai, T. Jung, T. Kawashima, S. Koyama , et al. (13 additional authors not shown)

Abstract: This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in Nor… ▽ More This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 93 pages, 54 figures. White Paper from East Asia VLBI community

The Greenland Telescope: Antenna Retrofit Status and Future Plans

Authors: Philippe Raffin, Paul T. P. Ho, Keiichi Asada, Raymond Blundell, Geoffrey C. Bower, Roberto Burgos, Chih-Cheng Chang, Ming-Tang Chen, You-Hua Chu, Paul K. Grimes, C. C. Han, Chih-Wei L. Huang, Yau-De Huang, Fang-Chia Hsieh, Makoto Inoue, Patrick M. Koch, Derek Kubo, Steve Leiker, Lupin Lin, Ching-Tang Liu, Shih-Hsiang Lo, Pierre Martin-Cocher, Satoki Matsushita, Masanori Nakamura, Zheng Meyer-Zhao , et al. (10 additional authors not shown)

Abstract: Since the ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO), SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA) are working jointly to relocate the antenna to Greenland. This paper shows the status of the antenna retrofit and the work carried out after the recommissioning and subsequent disassembly of the antenna at the VLA h… ▽ More Since the ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO), SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA) are working jointly to relocate the antenna to Greenland. This paper shows the status of the antenna retrofit and the work carried out after the recommissioning and subsequent disassembly of the antenna at the VLA has taken place. The next coming months will see the start of the antenna reassembly at Thule Air Base. These activities are expected to last until the fall of 2017 when commissioning should take place. In parallel, design, fabrication and testing of the last components are taking place in Taiwan. △ Less

Submitted 9 December, 2016; originally announced December 2016.

Lepton acceleration in the vicinity of the event horizon: High-energy and Very-high-energy emissions from rotating black holes with various masses

Authors: Kouichi Hirotani, Hung-Yi Pu, Lupin Chun-Che Lin, Hsiang-Kuang Chang, Makoto Inoue, Albert K. H Kong, Satoki Matsushita, Pak-Hin T. Tam

Abstract: 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 regi… ▽ More 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. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 17 pages, 28 figures; accepted for publication in ApJ

Structural Transition in the NGC 6251 Jet: An Interplay with the Supermassive Black Hole and Its Host Galaxy

Authors: Chih-Yin Tseng, Keiichi Asada, Masanori Nakamura, Hung-Yi Pu, Juan-Carlos Algaba, Wen-Ping Lo

Abstract: 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 supermas… ▽ More 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. △ Less

Submitted 24 January, 2017; v1 submitted 20 October, 2016; originally announced October 2016.

Comments: 12 pages, 4 figures, 3 tables

Journal ref: ApJ 833 (2016) 1-8

The Effects of Accretion Flow Dynamics on the Black Hole Shadow of Sagittarius A$^{*}$

Authors: Hung-Yi Pu, Kazunori Akiyama, Keiichi Asada

Abstract: A radiatively inefficient accretion flow (RIAF), which is commonly characterized by its sub-Keplerian nature, is a favored accretion model for the supermassive black hole at Galactic center, Sagittarius A$^{*}$. To investigate the observable features of a RIAF, we compare the modeled shadow images, visibilities, and spectra of three flow models with dynamics characterized by (i) a Keplerian shell… ▽ More A radiatively inefficient accretion flow (RIAF), which is commonly characterized by its sub-Keplerian nature, is a favored accretion model for the supermassive black hole at Galactic center, Sagittarius A$^{*}$. To investigate the observable features of a RIAF, we compare the modeled shadow images, visibilities, and spectra of three flow models with dynamics characterized by (i) a Keplerian shell which is rigidly-rotating outside the innermost stable circular orbit (ISCO) and infalling with a constant angular momentum inside ISCO, (ii) a sub-Keplerian motion, and (iii) a free-falling motion with zero angular momentum at infinity. At near-mm wavelengths the emission is dominated by the flow within several Schwarzschild radii. The energy shift due to the flow dynamics becomes important and distinguishable, suggesting that the flow dynamics are an important model parameter for interpreting the mm/submillimeter very long baseline interferometric observations with the forthcoming, fully assembled Event Horizon Telescope (EHT). As an example, we demonstrate that structural variations of Sagittarius A$^{*}$ on event horizon-scales detected in previous EHT observations can be explained by the non-stationary dynamics of a RIAF. △ Less

Submitted 9 August, 2016; originally announced August 2016.

Comments: 8 pages, 4 figures. Accepted by ApJ

Odyssey: A Public GPU-Based Code for General-Relativistic Radiative Transfer in Kerr Spacetime

Authors: Hung-Yi Pu, Kiyun Yun, Ziri Younsi, Suk-Jin Yoon

Abstract: General-relativistic radiative transfer (GRRT) calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra and light curves from matter in the vicinity of black holes. Such studies are especially important for ongoing and upcoming millimeter/submillimeter (mm/sub-mm) Very Long Baseline Interferometry (VLBI) observations of… ▽ More General-relativistic radiative transfer (GRRT) calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra and light curves from matter in the vicinity of black holes. Such studies are especially important for ongoing and upcoming millimeter/submillimeter (mm/sub-mm) Very Long Baseline Interferometry (VLBI) observations of the supermassive black holes at the centres of Sgr A^{*} and M87. To this end we introduce Odyssey, a Graphics Processing Unit(GPU)-based code for ray tracing and radiative transfer in the Kerr spacetime. On a single GPU, the performance of Odyssey can exceed 1 nanosecond per photon, per Runge-Kutta integration step. Odyssey is publicly available, fast, accurate, and flexible enough to be modified to suit the specific needs of new users. Along with a Graphical User Interface (GUI) powered by a video-accelerated display architecture, we also present an educational software tool, Odyssey_Edu, for showing in real time how null geodesics around a Kerr black hole vary as a function of black hole spin and angle of incidence onto the black hole. △ Less

Submitted 19 January, 2016; v1 submitted 8 January, 2016; originally announced January 2016.

Comments: 12 pages, 12 figures. Websites are ready for downloading the source code (Odyssey) and the educational software (Odyssey_Edu). Accepted by ApJ

Energetic Gamma Radiation from Rapidly Rotating Black Holes

Authors: Kouichi Hirotani, Hung-Yi Pu

Abstract: 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… ▽ More 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 hole's 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. △ Less

Submitted 15 December, 2015; originally announced December 2015.

Comments: 21 pages, 20 figures. Accepted for publication in ApJ, December 14, 2015

First-generation Science Cases for Ground-based Terahertz Telescopes

Authors: Hiroyuki Hirashita, Patrick M. Koch, Satoki Matsushita, Shigehisa Takakuwa, Masanori Nakamura, Keiichi Asada, Hauyu Baobab Liu, Yuji Urata, Ming-Jye Wang, Wei-Hao Wang, Satoko Takahashi, Ya-Wen Tang, Hsian-Hong Chang, Kuiyun Huang, Oscar Morata, Masaaki Otsuka, Kai-Yang Lin, An-Li Tsai, Yen-Ting Lin, Sundar Srinivasan, Pierre Martin-Cocher, Hung-Yi Pu, Francisca Kemper, Nimesh Patel, Paul Grimes , et al. (11 additional authors not shown)

Abstract: Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is the… ▽ More Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: 39 pages, 17 figures, accepted for publication in PASJ (as a review paper)

Steady General Relativistic Magnetohydrodynamic Inflow/Outflow Solution along Large-Scale Magnetic Fields that Thread a Rotating Black Hole

Authors: Hung-Yi Pu, Masanori Nakamura, Kouichi Hirotani, Yosuke Mizuno, Kinwah Wu, Keiichi Asada

Abstract: 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… ▽ More 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. △ Less

Submitted 18 April, 2015; v1 submitted 9 January, 2015; originally announced January 2015.

Comments: 13 Pages, incliding 2 tables and 5 Figures; accepted by ApJ

The Greenland Telescope (GLT): Antenna status and future plans

Authors: Philippe Raffin, Juan Carlos Algaba-Marcos, Keichi Asada, Raymond Blundell, Roberto Burgos, Chih-Cheng Chang, Ming-Tang Chen, Robert Christensen, Paul K. Grimes, C. C. Han, Paul T. P. Ho, Yau-De Huang, Makoto Inoue, Patrick M. Koch, Derek Kubo, Steve Leiker, Ching-Tang Liu, Pierre Martin-Cocher, Satoki Matsushita, Masanori Nakamura, Hiroaki Nishioka, George Nystrom, Scott N. Paine, Nimesh A. Patel, Nicolas Pradel , et al. (7 additional authors not shown)

Abstract: The ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO) in 2011. SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA), SAO's main partner for this project, are working jointly to relocate the antenna to Greenland to carry out millimeter and submillimeter VLBI observations. This paper presents the work carried out on upgrading the… ▽ More The ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO) in 2011. SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA), SAO's main partner for this project, are working jointly to relocate the antenna to Greenland to carry out millimeter and submillimeter VLBI observations. This paper presents the work carried out on upgrading the antenna to enable operation in the Arctic climate by the GLT Team to make this challenging project possible, with an emphasis on the unexpected telescope components that had to be either redesigned or changed. Five-years of inactivity, with the antenna laying idle in the desert of New Mexico, coupled with the extreme weather conditions of the selected site in Greenland have it necessary to significantly refurbish the antenna. We found that many components did need to be replaced, such as the antenna support cone, the azimuth bearing, the carbon fiber quadrupod, the hexapod, the HVAC, the tiltmeters, the antenna electronic enclosures housing servo and other drive components, and the cables. We selected Vertex, the original antenna manufacturer, for the main design work, which is in progress. The next coming months will see the major antenna components and subsystems shipped to a site of the US East Coast for test-fitting the major antenna components, which have been retrofitted. The following step will be to ship the components to Greenland to carry out VLBI and single dish observations. Antenna reassembly at Summit Station should take place during the summer of 2018. △ Less

Submitted 18 July, 2014; originally announced July 2014.

Comments: 10 pages, 8 figures

Journal ref: Proc. SPIE 9145-15, 2014

Greenland Telescope Project --- Direct Confirmation of Black Hole with Sub-millimeter VLBI

Authors: M. Inoue, J. C. Algaba-Marcos, K. Asada, C. -C. Chang, M. -T. Chen, J. Han, H. Hirashita, P. T. P. Ho, S. -N. Hsieh, T. Huang, H. Jiang, P. M. Koch, D. Y. Kubo, C. -Y. Kuo, B. Liu, P. Martin-Cocher, S. Matsushita, Z. Meyer-Zhao, M. Nakamura, H. Nishioka, G. Nystrom, N. Pradel, H. -Y. Pu, P. A. Raffin, H. -Y. Shen , et al. (14 additional authors not shown)

Abstract: 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 > 9,000 km t… ▽ More 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 > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. 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 three 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. △ Less

Submitted 9 July, 2014; originally announced July 2014.

Comments: 15 pages, 6 figures, accepted by Radio Science journal

Measuring Mass Accretion Rate onto the Supermassive Black Hole in M 87 Using Faraday Rotation Measure with the Submillimeter Array

Authors: C. Y. Kuo, K. Asada, R. Rao, M. Nakamura, J. C. Algaba, H. B. Liu, M. Inoue, P. M. Koch, P. T. P. Ho, S. Matsushita, H. -Y. Pu, K. Akiyama, H. Nishioka, N. Pradel

Abstract: We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($χ$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $χ$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the… ▽ More We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($χ$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $χ$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the M 87 core to be between $-$7.5$\times$10$^{5}$ and 3.4$\times$10$^{5}$ rad/m$^{2}$. Assuming a density profile of the accretion flow that follows a power-law distribution and a magnetic field that is ordered, radial, and has equipartition strength, the limit on the rotation measure constrains the mass accretion rate $\dot{M}$ to be below 9.2$\times$10$^{-4}$ M$_{\odot}$~yr$^{-1}$ at a distance of 21 Schwarzchild radii from the central black hole. This value is at least two orders of magnitude smaller than the Bondi accretion rate, suggesting significant suppression of the accretion rate in the inner region of the accretion flow. Consequently, our result disfavors the classical \emph{advection dominated accretion flow} (ADAF) and prefers the \emph{adiabatic inflow-outflow solution} (ADIOS) or \emph{convection-dominated accretion flow} (CDAF) for the hot accretion flow in M 87. △ Less

Submitted 21 February, 2014; originally announced February 2014.

Comments: 11 pages, 4 figures, accepted by ApJL

Greenland Telescope (GLT) Project: "A Direct Confirmation of Black Hole with Submillimeter VLBI"

Authors: M. Nakamura, J. -C. Algaba, K. Asada, B. Chen, M. -T. Chen, J. Han, P. H. P. Ho, S. -N. Hsieh, T. Huang, M. Inoue, P. Koch, C. -Y. Kuo, P. Martin-Cocher, S. Matsushita, Z. Meyer-Zhao, H. Nishioka, G. Nystom, N. Pradel, H. -Y. Pu, P. Raffin, H. -Y. Shen, C. -Y. Tseng, the Greenland Telescope Project Team

Abstract: The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $μ$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to ob… ▽ More The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $μ$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $μ$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016. △ Less

Submitted 6 October, 2013; originally announced October 2013.

Comments: 4 pages, 5 figures; to appear in the conference proceedings of "Innermost Regions of Relativistic Jets and Their Magnetic Fields", June 10-14, 2013, Granada, Spain

Two Types of Ergospheric Jets from Accreting Black Holes: The Dichotomy of Fanaroff-Riley Galaxies

Authors: Hung-Yi Pu, Kouichi Hirotani, Yosuke Mizuno, Hsiang-Kuang Chang

Abstract: We investigate the extraction of the rotational energy of a black hole under different accreting environment. When the accretion rate is moderate, the accretion disk consists of an outer thin disk and an inner advection-dominated accretion flow. In such a combined disk, the outer thin disk can sustain a magnetic field with moderate strength at the event horizon, leading to the formation of relativ… ▽ More We investigate the extraction of the rotational energy of a black hole under different accreting environment. When the accretion rate is moderate, the accretion disk consists of an outer thin disk and an inner advection-dominated accretion flow. In such a combined disk, the outer thin disk can sustain a magnetic field with moderate strength at the event horizon, leading to the formation of relativistic jets with moderate luminosity and speed via the magnetohrodynamic Penrose process. When the accretion rate increases enough, on the other hand, the disk becomes geometrically thin near the horizon. In this slim disk, the denser plasmas can sustain a stronger magnetic field than that in a combined disk, leading to the formation of jets with greater luminosity and speed via the Blandford-Znajek processs. It is discussed that the former jets are associated with the Fanaroff-Riley (FR) I galaxies and the latter with FR II galaxies. △ Less

Submitted 7 November, 2012; originally announced November 2012.

Comments: 6 pages, including 1 table and 4 figures. Submitted

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

Authors: Hung-Yi Pu, Kouichi Hirotani, Hsiang-Kuang Chang

Abstract: 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, $\dot{m}$, 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 t… ▽ More 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, $\dot{m}$, 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 $\dot{m}\lesssim0.1$), 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 brings more rest-mass energy than what is extracted from the hole electromagnetically to stop the extraction. Several other 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. △ Less

Submitted 21 September, 2012; originally announced September 2012.

Comments: Accepted for publication in the Astrophysical Journal. 17pages, 9 figures

Journal ref: The Astrophysical Journal, 758(2012), 113

On Spin Dependence of Relativistic Acoustic Geometry

Authors: Hung-Yi Pu, Ishita Maiti, Tapas Kumar Das, Hsiang-Kuang Chang

Abstract: This work makes the first ever attempt to understand the influence of the black hole background space-time in determining the fundamental properties of the embedded relativistic acoustic geometry. To accomplish such task, 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 v… ▽ More This work makes the first ever attempt to understand the influence of the black hole background space-time in determining the fundamental properties of the embedded relativistic acoustic geometry. To accomplish such task, 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 space time) has been investigated for axially symmetric inflow of hydrodynamic fluid onto a rotating black hole. 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. For certain values of the initial boundary conditions describing the aforementioned flow, more than one acoustic horizons, namely two black hole type and one white hole type, may form, where the surface gravity may become formally infinite at the acoustic white hole. The connection between the corresponding analogue Hawking temperature with astrophysically relevant observables associated with the spectral signature has been discussed. △ Less

Submitted 5 April, 2012; originally announced April 2012.

Comments: 22 pages, 11 figures, Comments welcome