Search

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.06 + 0.12 μ as and a corresponding black hole mass of (7.13 ± 0.39) × 109 M ⊙, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

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.

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., 16 pages, 7 figures

Full list of authors: Okino, Hiroki; Akiyama, Kazunori; Asada, Keiichi; Gomez, Jose L.; Hada, Kazuhiro; Honma, Mareki; Krichbaum, Thomas P.; Kino, Motoki; Nagai, Hiroshi; Bach, Uwe; Blackburn, Lindy; Bouman, Katherine L.; Chael, Andrew; Crew, Geoffrey B.; Doeleman, Sheperd S.; Fish, Vincent L.; Goddi, Ciriaco; Issaoun, Sara; Johnson, Michael D.; Jorstad, Svetlana; Koyama, Shoko; Lonsdale, Colin J.; Lu, Ru-Sen; Marti-Vidal, Ivan; Matthews, Lynn D.; Mizuno, Yosuke; Moriyama, Kotaro; Nakamura, Masanori; Pu, Hung-Yi; Ros, Eduardo; Savolainen, Tuomas; Tazaki, Fumie; Wagner, Jan; Wielgus, Maciek; Zensus, Anton.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited., The collimation of relativistic jets launched from the vicinity of supermassive black holes (SMBHs) at the centers of active galactic nuclei (AGNs) 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 ∼60 μas, resolving the innermost part of the jet ever on scales of ∼105 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 AGNs. We discovered the jet collimation break around 107 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 from 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. © 2022. The Author(s). Published by the American Astronomical Society., This work was financially supported by grants from the National Science Foundation (NSF; AST-1440254; AST-1614868; AST-2034306) and the Japan Society for Promotion of Science (JSPS; JP21H01137; JP20H01951; JP19KK0081; JP18H03721.). This work has been partially supported by the Generalitat Valenciana GenT Project CIDEGENT/2018/021 and by the MICINN Research Project PID2019-108995GB-C22. This work was partially supported by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) under grant 2021/01183-8. K. Akiyama is also financially supported by the following NSF grants (AST-1935980, AST-2107681, AST-2132700, OMA-2029670). A.C. was supported by Hubble Fellowship grant HST-HF2-51431.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. A.C. also gratefully acknowledges support from the Princeton Gravity Initiative. S.I. is supported by NASA Hubble Fellowship grant HST-HF2-51482.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. R.-S. L. is supported by the Max Planck Partner Group of the MPG and the CAS and acknowledges support by the Key Program of the National Natural Science Foundation of China (grant No. 11933007), the Key Research Program of Frontier Sciences, CAS (grant No. ZDBS-LY-SLH011), and the Shanghai Pilot Program for Basic Research—Chinese Academy of Science, Shanghai Branch (JCYJ-SHFY-2022-013). Y.M. acknowledges the support of the National Natural Science Foundation of China (grant No. 12273022). The black hole Initiative at Harvard University is financially supported by a grant from the John Templeton Foundation., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

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., Comment: 22 pages, 13 figures, 4 tables, accepted for publication in ApJ

We present measurements of the properties of the central radio source in M87 using Event Horizon Telescope data obtained during the 2017 campaign. We develop and fit geometric crescent models (asymmetric rings with interior brightness depressions) using two independent sampling algorithms that consider distinct representations of the visibility data. We show that the crescent family of models is statistically preferred over other comparably complex geometric models that we explore. We calibrate the geometric model parameters using general relativistic magnetohydrodynamic (GRMHD) models of the emission region and estimate physical properties of the source. We further fit images generated from GRMHD models directly to the data. We compare the derived emission region and black hole parameters from these analyses with those recovered from reconstructed images. There is a remarkable consistency among all methods and data sets. We find that >50% of the total flux at arcsecond scales comes from near the horizon, and that the emission is dramatically suppressed interior to this region by a factor >10, providing direct evidence of the predicted shadow of a black hole. Across all methods, we measure a crescent diameter of 42 ± 3 μas and constrain its fractional width to be, The authors of this Letter thank the following organizations and programs: the Academy of Finland (projects 274477, 284495, 312496); the Advanced European Network of E-infrastructures for Astronomy with the SKA (AENEAS) project, supported by the European Commission Framework Programme Horizon 2020 Research and Innovation action under grant agreement 731016; the Alexander von Humboldt Stiftung; the Black Hole Initiative at Harvard University, through a grant (60477) from the John Templeton Foundation; the China Scholarship Council; Comisión Nacional de Investigación Científica y Tecnológica (CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB170002, ALMA-conicyt 31140007); Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico, projects 104497, 275201, 279006, 281692); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Dirección General de Asuntos del Personal Académico-Universidad Nacional 9 The Astrophysical Journal Letters, 875:L1 (17pp), 2019 April 10 The EHT Collaboration et al. Autónoma de México (DGAPA-UNAM, project IN112417); the European Research Council (ERC) Synergy Grant “BlackHoleCam: Imaging the Event Horizon of Black Holes” (grant 610058); the Generalitat Valenciana postdoctoral grant APOSTD/2018/177; the Gordon and Betty Moore Foundation (grants GBMF-3561, GBMF-5278); the Istituto Nazionale di Fisica Nucleare (INFN) sezione di Napoli, iniziative specifiche TEONGRAV; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; the Jansky Fellowship program of the National Radio Astronomy Observatory (NRAO); the Japanese Government (Monbukagakusho: MEXT) Scholarship; the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for JSPS Research Fellowship (JP17J08829); JSPS Overseas Research Fellowships; the Key Research Program of Frontier Sciences, Chinese, Academy of Sciences (CAS, grants QYZDJ-SSW-SLH057, QYZDJ-SSW-SYS008); the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft (MPG); the Max Planck Partner Group of the MPG and the CAS; the MEXT/JSPS KAKENHI (grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, 25120007); the MIT International Science and Technology Initiatives (MISTI) Funds; the Ministry of Science and Technology (MOST) of Taiwan (105-2112-M-001-025-MY3, 106-2112-M001-011, 106-2119-M-001-027, 107-2119-M-001-017, 107- 2119-M-001-020, and 107-2119-M-110-005); the National Aeronautics and Space Administration (NASA, Fermi Guest Investigator grant 80NSSC17K0649); the National Institute of Natural Sciences (NINS) of Japan; the National Key Research and Development Program of China (grant 2016YFA0400704, 2016YFA0400702); the National Science Foundation (NSF, grants AST-0096454, AST-0352953, AST-0521233, AST0705062, AST-0905844, AST-0922984, AST-1126433, AST1140030, DGE-1144085, AST-1207704, AST-1207730, AST1207752, MRI-1228509, OPP-1248097, AST-1310896, AST1312651, AST-1337663, AST-1440254, AST-1555365, AST1715061, AST-1615796, AST-1614868, AST-1716327, OISE1743747, AST-1816420); the Natural Science Foundation of China (grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312, 11873028, 11873073, U1531245, 11473010); the Natural Sciences and Engineering Research Council of Canada (NSERC, including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships-Doctoral Program); the National Youth Thousand Talents Program of China; the National Research Foundation of Korea (grant 2015-R1D1A1A01056807, the Global PhD Fellowship Grant: NRF-2015H1A2A1033752, and the Korea Research Fellowship Program: NRF-2015H1D3A1066561); the Netherlands Organization for Scientific Research (NWO) VICI award (grant 639.043.513) and Spinoza Prize (SPI 78-409); the New Scientific Frontiers with Precision, Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology (DST) of South Africa; the Onsala Space Observatory (OSO) national infrastructure, for the provisioning of its facilities/observational support (OSO receives funding through the Swedish Research Council under grant 2017-00648); the Perimeter Institute for Theoretical Physics (research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Economic Development, Job Creation and Trade); the Russian Science Foundation (grant 17-12-01029); the Spanish Ministerio de Economía y Competitividad (grants AYA2015-63939-C2-1-P, AYA2016-80889-P); the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709); the Toray Science Foundation; the US Department of Energy (USDOE) through the Los Alamos National Laboratory (operated by Triad National Security, LLC, for the National Nuclear Security Administration of the USDOE (Contract 89233218CNA000001)); the Italian Ministero dell’Istruzione Università e Ricerca through the grant Progetti Premiali 2012-iALMA (CUP C52I13000140001); the European Unionʼs Horizon 2020 research and innovation programme under grant agreement No 730562 RadioNet; ALMA North America Development Fund; Chandra TM6-17006X

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., Comment: 15 pages, 16 figures; Accepted for publication in the Astrophysical Journal. arXiv admin note: text overlap with arXiv:1610.07819

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.

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., Comment: 12 pages, 4 figures, 3 tables

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., Comment: 12 pages, 12 figures. Websites are ready for downloading the source code (Odyssey) and the educational software (Odyssey_Edu). Accepted by ApJ

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., Comment: 13 Pages, incliding 2 tables and 5 Figures; accepted by ApJ

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 ($\chi$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $\chi$ 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., Comment: 11 pages, 4 figures, accepted by ApJL