Showing 1–50 of 112 results for author: Nakamura, M

On the mechanism of black hole energy reduction in the Blandford-Znajek process

Authors: Kenji Toma, Fumio Takahara, Masanori Nakamura

Abstract: The Blandford-Znajek (BZ) process is electromagnetic energy release from rotating black holes (BHs) along magnetic field lines threading them and widely believed to drive relativistic jets. This process is successfully demonstrated in general relativistic magnetohydrodynamic (MHD) simulations with the coordinate system regular on the event horizon, by which one can estimate the outward Poynting fl… ▽ More The Blandford-Znajek (BZ) process is electromagnetic energy release from rotating black holes (BHs) along magnetic field lines threading them and widely believed to drive relativistic jets. This process is successfully demonstrated in general relativistic magnetohydrodynamic (MHD) simulations with the coordinate system regular on the event horizon, by which one can estimate the outward Poynting flux, although the direct energy release through the horizon shown in the simulations does not provide an intuitive picture. We revisit the mechanism of BH energy reduction by utilizing the coordinate system singular on the horizon, in which the falling membrane of past accreted matter should exist above the horizon. We find that the Poynting flux is produced at the boundary between the falling membrane and the magnetically-dominated inflow, and the front of the inflow creates the negative electromagnetic energy, which reduces the rotational energy of spacetime. We also clarify that the poloidal electric current does not form a closed circuit within the magnetically-dominated flow. Previous interpretations of the BZ process and possibilities of the ideal MHD violation and BH charging are also discussed. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: 14 pages, 2 figures; submitted to PTEP

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

The Interplay Between Collisionless Magnetic Reconnection and Turbulence

Authors: J. E. Stawarz, P. A. Muñoz, N. Bessho, R. Bandyopadhyay, T. K. M. Nakamura, S. Eriksson, D. Graham, J. Büchner, A. Chasapis, J. F. Drake, M. A. Shay, R. E. Ergun, H. Hasegawa, Yu. V. Khotyaintsev, M. Swisdak, F. Wilder

Abstract: Alongside magnetic reconnection, turbulence is another fundamental nonlinear plasma phenomenon that plays a key role in energy transport and conversion in space and astrophysical plasmas. From a numerical, theoretical, and observational point of view there is a long history of exploring the interplay between these two phenomena in space plasma environments; however, recent high-resolution, multi-s… ▽ More Alongside magnetic reconnection, turbulence is another fundamental nonlinear plasma phenomenon that plays a key role in energy transport and conversion in space and astrophysical plasmas. From a numerical, theoretical, and observational point of view there is a long history of exploring the interplay between these two phenomena in space plasma environments; however, recent high-resolution, multi-spacecraft observations have ushered in a new era of understanding this complex topic. The interplay between reconnection and turbulence is both complex and multifaceted, and can be viewed through a number of different interrelated lenses - including turbulence acting to generate current sheets that undergo magnetic reconnection (turbulence-driven reconnection), magnetic reconnection driving turbulent dynamics in an environment (reconnection-driven turbulence) or acting as an intermediate step in the excitation of turbulence, and the random diffusive/dispersive nature of magnetic field lines embedded in turbulent fluctuations enabling so-called stochastic reconnection. In this paper, we review the current state of knowledge on these different facets of the interplay between turbulence and reconnection in the context of collisionless plasmas, such as those found in many near-Earth astrophysical environments, from a theoretical, numerical, and observational perspective. Particular focus is given to several key regions in Earth's magnetosphere - Earth's magnetosheath, magnetotail, and Kelvin-Helmholtz vortices on the magnetopause flanks - where NASA's Magnetospheric Multiscale mission has been providing new insights on the topic. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: Chapter 2.2 of ISSI Book on Magnetic Reconnection, submitted to Space Science Reviews

New high-precision measurement system for electron-positron pairs from sub-GeV/GeV gamma-rays in the emulsion telescope

Authors: Yuya Nakamura, Shigeki Aoki, Tomohiro Hayakawa, Atsushi Iyono, Ayaka Karasuno, Kohichi Kodama, Ryosuke Komatani, Masahiro Komatsu, Masahiro Komiyama, Kenji Kuretsubo, Toshitsugu Marushima, Syota Matsuda, Kunihiro Morishima, Misaki Morishita, Naotaka Naganawa, Mitsuhiro Nakamura, Motoya Nakamura, Takafumi Nakamura, Noboru Nakano, Toshiyuki Nakano, Akira Nishio, Miyuki Oda, Hiroki Rokujo, Osamu Sato, Kou Sugimura , et al. (5 additional authors not shown)

Abstract: The GRAINE project observes cosmic gamma-rays, using a balloon-borne emulsion-film-based telescope in the sub-GeV/GeV energy band. We reported in our previous balloon experiment in 2018, GRAINE2018, the detection of the known brightest source, Vela pulsar, with the highest angular resolution ever reported in an energy range of $>$80 MeV. However, the emulsion scanning system used in the experiment… ▽ More The GRAINE project observes cosmic gamma-rays, using a balloon-borne emulsion-film-based telescope in the sub-GeV/GeV energy band. We reported in our previous balloon experiment in 2018, GRAINE2018, the detection of the known brightest source, Vela pulsar, with the highest angular resolution ever reported in an energy range of $>$80 MeV. However, the emulsion scanning system used in the experiment was designed to achieve a high-speed scanning, and it was not accurate enough to ensure the optimum spacial resolution of the emulsion film and limited the performance. Here, we report a new high-precision scanning system that can be used to greatly improve the observation result of GRAINE2018 and also be employed in future experiments. The system involves a new algorithm that recognizes each silver grain on an emulsion film and is capable of measuring tracks with a positional resolution for the passing points of tracks of almost the same as the intrinsic resolution of nuclear emulsion film ($\sim$70 nm). This resolution is approximately one order of magnitude smaller than that obtained with the high-speed scanning system. With this system, an angular resolution for gamma-rays of 0.1$^\circ$ at 1 GeV is expected to be achieved. Furthermore, we successfully combine the new high-precision system with the existing high-speed system, establishing the system to make a high-speed and high-precision measurement. Employing these systems, we reanalyze the gamma-ray events detected previously by only the high-speed system in GRAINE2018 and obtain an about three times higher angular resolution (0.22$^\circ$) in 500--700 MeV than the original value. The high-resolution observation may bring new insights into the gamma-ray emission from the Galactic center region and may realize polarization measurements of high-energy cosmic gamma-rays. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 20 pages, 11 figures

Jet Collimation and Acceleration in the Flat Spectrum Radio Quasar 1928+738

Authors: Kunwoo Yi, Jongho Park, Masanori Nakamura, Kazuhiro Hada, Sascha Trippe

Abstract: Using time-resolved multifrequency Very Long Baseline Array data and new KaVA (KVN and VERA Array) observations, we study the structure and kinematics of the jet of the flat spectrum radio quasar (FSRQ) 1928+738. We find two distinct jet geometries as function of distance from the central black hole, with the inner jet having a parabolic shape, indicating collimation, and the outer jet having a co… ▽ More Using time-resolved multifrequency Very Long Baseline Array data and new KaVA (KVN and VERA Array) observations, we study the structure and kinematics of the jet of the flat spectrum radio quasar (FSRQ) 1928+738. We find two distinct jet geometries as function of distance from the central black hole, with the inner jet having a parabolic shape, indicating collimation, and the outer jet having a conical shape, indicating free expansion of the jet plasma. Jet component speeds display a gradual outward acceleration up to a bulk Lorentz factor $Γ_{\rm max} \approx10$, followed by a deceleration further downstream. The location of the acceleration zone matches the region where the jet collimation occurs; this is the first direct observation of an acceleration and collimation zone (ACZ) in an FSRQ. The ACZ terminates approximately at a distance of 5.6$\times 10^6$ gravitational radii, which is in good agreement with the sphere of gravitational influence of the supermassive black hole, implying that the physical extent of the ACZ is controlled by the black hole gravity. Our results suggest that confinement by an external medium is responsible for the jet collimation and that the jet is accelerated by converting Poynting flux energy to kinetic energy. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 20 pages, 12 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 688, A94 (2024)

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

Observational Evidence to Support a Dense Ambient Medium Shaping the Jet in 3C 84

Authors: Jongho Park, Motoki Kino, Hiroshi Nagai, Masanori Nakamura, Keiichi Asada, Minchul Kam, Jeffrey A. Hodgson

Abstract: Highly collimated relativistic jets are a defining feature of certain active galactic nuclei (AGN), yet their formation mechanism remains elusive. Previous observations and theoretical models have proposed that the ambient medium surrounding the jets could exert pressure, playing a crucial role in shaping the jets. However, direct observational confirmation of such a medium has been lacking. In th… ▽ More Highly collimated relativistic jets are a defining feature of certain active galactic nuclei (AGN), yet their formation mechanism remains elusive. Previous observations and theoretical models have proposed that the ambient medium surrounding the jets could exert pressure, playing a crucial role in shaping the jets. However, direct observational confirmation of such a medium has been lacking. In this study, we present very long baseline interferometric (VLBI) observations of 3C 84 (NGC 1275), located at the center of the Perseus Cluster. Through monitoring observations with the Very Long Baseline Array (VLBA) at 43 GHz, a jet knot was detected to have been ejected from the sub-parsec scale core in the late 2010s. Intriguingly, this knot propagated in a direction significantly offset from the parsec-scale jet direction. To delve deeper into the matter, we employ follow-up VLBA 43 GHz observations, tracing the knot's trajectory until the end of 2022. We discovered that the knot abruptly changed its trajectory in the early 2020s, realigning itself with the parsec-scale jet direction. Additionally, we present results from an observation of 3C 84 with the Global VLBI Alliance (GVA) at 22 GHz, conducted near the monitoring period. By jointly analyzing the GVA 22 GHz image with a VLBA 43 GHz image observed about one week apart, we generated a spectral index map, revealing an inverted spectrum region near the edge of the jet where the knot experienced deflection. These findings suggest the presence of a dense, cold ambient medium characterized by an electron density exceeding $\sim10^5\ {\rm cm^{-3}}$, which guides the jet's propagation on parsec scales and significantly contributes to the overall shaping of the jet. △ Less

Submitted 28 December, 2023; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 18 pages, 18 figures, accepted for publication in Astronomy and Astrophysics

A kinematic study of the disc-outflow system around a high-mass protostar G59.783+0.065 probed by methanol and water masers

Authors: M. Nakamura, K. Motogi, H. Nakamura, Y. Yonekura, K. Fujisawa

Abstract: Class II CH3OH masers are used as a convenient tracer of disc-like structures in high-mass star formation. However, more than half of them show a complex distribution in Very Long Baseline Interferometry (VLBI) maps. The origin of such a complex distribution is still unknown. We conducted VLBI monitoring observations to unveil the origin of a complex class II CH3OH maser in the high-mass star-form… ▽ More Class II CH3OH masers are used as a convenient tracer of disc-like structures in high-mass star formation. However, more than half of them show a complex distribution in Very Long Baseline Interferometry (VLBI) maps. The origin of such a complex distribution is still unknown. We conducted VLBI monitoring observations to unveil the origin of a complex class II CH3OH maser in the high-mass star-forming region G59.783+0.065. We observed the CH3OH maser at 6.7 GHz and the H2O maser at 22 GHz to probe detailed circumstellar kinematics and structures by the Japanese VLBI network and the VLBI Exploration of Radio Astrometry. We found similar bipolar distributions in both masers, specifically two clusters located 2000 au apart along the East-West direction. We detected a linear distribution of CH3OH masers in the Western cluster. A position-velocity diagram shows that the Western CH3OH masers trace a rotating disc-wind or infalling component inside an edge-on disc-like structure. In contrast to the simple bipolar expanding motions of the H2O masers, the CH3OH masers exhibited complex motions despite their spatial coincidence. Some of the Eastern CH3OH masers showed bipolar expansions similar to the H2O masers, while others displayed random or even inward motions. Such complex kinematics and their close association with the H2O maser could occur at the boundary between outflow and inflow. We suggest that the complex distribution of class II CH3OH masers, like G59.783+0.065 arises from several distinct circumstellar structures that simultaneously achieve maser excitation. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: Accepted to MNRAS, 22 pages, 17 figures, 13 tables

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

Advanced methods for analyzing in-situ observations of magnetic reconnection

Authors: H. Hasegawa, M. R. Argall, N. Aunai, R. Bandyopadhyay, N. Bessho, I. J. Cohen, R. E. Denton, J. C. Dorelli, J. Egedal, S. A. Fuselier, P. Garnier, V. Genot, D. B. Graham, K. J. Hwang, Y. V. Khotyaintsev, D. B. Korovinskiy, B. Lavraud, Q. Lenouvel, T. C. Li, Y. -H. Liu, B. Michotte de Welle, T. K. M. Nakamura, D. S. Payne, S. M. Petrinec, Y. Qi , et al. (11 additional authors not shown)

Abstract: There is ample evidence for magnetic reconnection in the solar system, but it is a nontrivial task to visualize, to determine the proper approaches and frames to study, and in turn to elucidate the physical processes at work in reconnection regions from in-situ measurements of plasma particles and electromagnetic fields. Here an overview is given of a variety of single- and multi-spacecraft data a… ▽ More There is ample evidence for magnetic reconnection in the solar system, but it is a nontrivial task to visualize, to determine the proper approaches and frames to study, and in turn to elucidate the physical processes at work in reconnection regions from in-situ measurements of plasma particles and electromagnetic fields. Here an overview is given of a variety of single- and multi-spacecraft data analysis techniques that are key to revealing the context of in-situ observations of magnetic reconnection in space and for detecting and analyzing the diffusion regions where ions and/or electrons are demagnetized. We focus on recent advances in the era of the Magnetospheric Multiscale mission, which has made electron-scale, multi-point measurements of magnetic reconnection in and around Earth's magnetosphere. △ Less

Submitted 24 June, 2024; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: submitted to Space Science Reviews (116 pages, incl. 31 figures, 7 tables)

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

Spectral analysis of a parsec-scale jet in M87: Observational constraint on the magnetic field strengths in the jet

Authors: Hyunwook Ro, Motoki Kino, Bong Won Sohn, Kazuhiro Hada, Jongho Park, Masanori Nakamura, Yuzhu Cui, Kunwoo Yi, Aeree Chung, Jeffrey Hodgson, Tomohisa Kawashima, Tao An, Sascha Trippe, Juan-Carlos Algaba, Jae-Young Kim, Satoko Sawada-Satoh, Kiyoaki Wajima, Zhiqiang Shen, Xiaopeng Cheng, Ilje Cho, Wu Jiang, Taehyun Jung, Jee-Won Lee, Kotaro Niinuma, Junghwan Oh , et al. (27 additional authors not shown)

Abstract: Because of its proximity and the large size of its black hole, M87 is one of the best targets for studying the launching mechanism of active galactic nucleus jets. Currently, magnetic fields are considered to be an essential factor in the launching and accelerating of the jet. However, current observational estimates of the magnetic field strength of the M87 jet are limited to the innermost part o… ▽ More Because of its proximity and the large size of its black hole, M87 is one of the best targets for studying the launching mechanism of active galactic nucleus jets. Currently, magnetic fields are considered to be an essential factor in the launching and accelerating of the jet. However, current observational estimates of the magnetic field strength of the M87 jet are limited to the innermost part of the jet or to HST-1. No attempt has yet been made to measure the magnetic field strength in between. We aim to infer the magnetic field strength of the M87 jet out to a distance of several thousand $r_s$ by tracking the distance-dependent changes in the synchrotron spectrum of the jet from high-resolution very long baseline interferometry observations. In order to obtain high-quality spectral index maps, quasi-simultaneous observations at 22 and 43 GHz were conducted using the KVN and VERA Array (KaVA) and the VLBA. We compared the spectral index distributions obtained from the observations with a model and placed limits on the magnetic field strengths as a function of distance. The overall spectral morphology is broadly consistent over the course of these observations. The observed synchrotron spectrum rapidly steepens from $α_{22-43 GHz}$ ~ -0.7 at ~ 2 mas to $α_{22-43 GHz}$ ~ -2.5 at ~ 6 mas. A spectral index model in which nonthermal electron injections inside the jet decrease with distance can adequately reproduce the observed trend. This suggests the magnetic field strength of the jet at a distance of 2 - 10 mas (~ 900 $r_s$ - ~ 4500 $r_s$ in the deprojected distance) has a range of $B=(0.3 - 1.0 G)(z/2 mas)^{-0.73}$. Extrapolating to the EHT scale yields consistent results, suggesting that the majority of the magnetic flux of the jet near the black hole is preserved out to ~ 4500 $r_s$ without significant dissipation. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 16 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 673, A159 (2023)

Overview of the Observing System and Initial Scientific Accomplishments of the East Asian VLBI Network (EAVN)

Authors: Kazunori Akiyama, Juan-Carlos Algaba, Tao An, Keiichi Asada, Kitiyanee Asanok, Do-Young Byun, Thanapol Chanapote, Wen Chen, Zhong Chen, Xiaopeng Cheng, James O. Chibueze, Ilje Cho, Se-Hyung Cho, Hyun-Soo Chung, Lang Cui, Yuzhu Cui, Akihiro Doi, Jian Dong, Kenta Fujisawa, Wei Gou, Wen Guo, Kazuhiro Hada, Yoshiaki Hagiwara, Tomoya Hirota, Jeffrey A. Hodgson , et al. (79 additional authors not shown)

Abstract: The East Asian VLBI Network (EAVN) is an international VLBI facility in East Asia and is operated under mutual collaboration between East Asian countries, as well as part of Southeast Asian and European countries. EAVN currently consists of 16 radio telescopes and three correlators located in China, Japan, and Korea, and is operated mainly at three frequency bands, 6.7, 22, and 43 GHz with the lon… ▽ More The East Asian VLBI Network (EAVN) is an international VLBI facility in East Asia and is operated under mutual collaboration between East Asian countries, as well as part of Southeast Asian and European countries. EAVN currently consists of 16 radio telescopes and three correlators located in China, Japan, and Korea, and is operated mainly at three frequency bands, 6.7, 22, and 43 GHz with the longest baseline length of 5078 km, resulting in the highest angular resolution of 0.28 milliarcseconds at 43 GHz. One of distinct capabilities of EAVN is multi-frequency simultaneous data reception at nine telescopes, which enable us to employ the frequency phase transfer technique to obtain better sensitivity at higher observing frequencies. EAVN started its open-use program in the second half of 2018, providing a total observing time of more than 1100 hours in a year. EAVN fills geographical gap in global VLBI array, resulting in enabling us to conduct contiguous high-resolution VLBI observations. EAVN has produced various scientific accomplishments especially in observations toward active galactic nuclei, evolved stars, and star-forming regions. These activities motivate us to initiate launch of the 'Global VLBI Alliance' to provide an opportunity of VLBI observation with the longest baselines on the earth. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 27 pages, appeared in Galaxies special issue 'Challenges in Understanding Black Hole Powered Jets with VLBI' as an invited review

Journal ref: Galaxies 2022, 10(6), 113

Spectral evolution of dark asteroid surfaces induced by space weathering over a decade

Authors: Sunao Hasegawa, Francesca E. DeMeo, Michael Marsset, Josef Hanus, Chrysa Avdellidou, Marco Delbo, Schelte J. Bus, Hidekazu Hanayama, Takashi Horiuchi, Driss Takir, Emmanuel Jehin, Marin Ferrais, Jooyeon Geem, Myungshin Im, Jinguk Seo, Yoonsoo P. Bach, Sunho Jin, Masateru Ishiguro, Daisuke Kuroda, Richard P. Binzel, Akiko M. Nakamura, Bin Yang, Pierre Vernazza

Abstract: The surface of airless bodies like asteroids in the Solar System are known to be affected by space weathering. Experiments simulating space weathering are essential for studying the effects of this process on meteorite samples, but the problem is that the time spent to reproduce space weathering in these experiments is billions of times shorter than the actual phenomenon. In December 2010, the T-t… ▽ More The surface of airless bodies like asteroids in the Solar System are known to be affected by space weathering. Experiments simulating space weathering are essential for studying the effects of this process on meteorite samples, but the problem is that the time spent to reproduce space weathering in these experiments is billions of times shorter than the actual phenomenon. In December 2010, the T-type asteroid 596 Scheila underwent a collision with a few-tens-of-meters impactor. A decade later, there is an opportunity to study how the surface layer of this asteroid is being altered by space weathering after the impact. To do so, we performed visible spectrophotometric and near-infrared spectroscopic observations of 596 Scheila. The acquired spectrum is consistent with those observed shortly after the 2010 impact event within the observational uncertainty range. This indicates that the surface color of dark asteroids is not noticeably changed by space weathering over a 10-year period. This study is the first to investigate color changes due to space weathering on an actual asteroid surface in the Solar System. Considering that fresh layers are regularly created on asteroid surfaces by collisions, we suggest a genetic link between D/T-type and dark (low albedo) X-complex asteroids and very red objects such as 269 Justitia, 732 Tjilaki (and 203 Pompeja). New observations show that 203 Pompeja has a X-type-like surface, with some local surface areas exhibiting a very red spectrum. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 16 pages, 9 figures, 2 tables, Accepted for publication in ApJ Letters

Reflectivity of Venus' dayside disk during the 2020 observation campaign: outcomes and future perspectives

Authors: Yeon Joo Lee, Antonio García Muñoz, Atsushi Yamazaki, Eric Quémerais, Stefano Mottola, Stephan Hellmich, Thomas Granzer, Gilles Bergond, Martin Roth, Eulalia Gallego-Cano, Jean-Yves Chaufray, Rozenn Robidel, Go Murakami, Kei Masunaga, Murat Kaplan, Orhan Erece, Ricardo Hueso, Petr Kabáth, Magdaléna Špoková, Agustín Sánchez-Lavega, Myung-Jin Kim, Valeria Mangano, Kandis-Lea Jessup, Thomas Widemann, Ko-ichiro Sugiyama , et al. (6 additional authors not shown)

Abstract: We performed a unique Venus observation campaign to measure the disk brightness of Venus over a broad range of wavelengths in August and September 2020. The primary goal of the campaign is to investigate the absorption properties of the unknown absorber in the clouds. The secondary goal is to extract a disk mean SO$_2$ gas abundance, whose absorption spectral feature is entangled with that of the… ▽ More We performed a unique Venus observation campaign to measure the disk brightness of Venus over a broad range of wavelengths in August and September 2020. The primary goal of the campaign is to investigate the absorption properties of the unknown absorber in the clouds. The secondary goal is to extract a disk mean SO$_2$ gas abundance, whose absorption spectral feature is entangled with that of the unknown absorber at the ultraviolet (UV) wavelengths. A total of 3 spacecraft and 6 ground-based telescopes participated in this campaign, covering the 52 to 1700~nm wavelength range. After careful evaluation of the observational data, we focused on the data sets acquired by 4 facilities. We accomplished our primary goal by analyzing the reflectivity spectrum of the Venus disk over the 283-800 nm wavelengths. Considerable absorption is present in the 350-450 nm range, for which we retrieved the corresponding optical depth by the unknown absorber. The result shows a consistent wavelength dependence of the relative optical depth with that at low latitudes during the Venus flyby by MESSENGER in 2007 (Pérez-Hoyos et al. 2018), which was expected because the overall disk reflectivity is dominated by low latitudes. Last, we summarize the experience obtained during this first campaign that should enable us to accomplish our second goal in future campaigns. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: Accepted for publication in PSJ, 45 pages, 13 figures

Disruption of Saturn's ring particles by thermal stress

Authors: Naoyuki Hirata, Ryuji Morishima, Keiji Ohtsuki, Akiko M. Nakamura

Abstract: Spacecraft and ground-based observations show that the main rings of Saturn lack particles larger than 10 m. Tidal or collisional destruction of satellites/comets have been proposed as the origin of the main rings; however, Saturn's tide alone cannot grind km-sized fragments into submeter-sized particles because of the high mechanical strength of water ice and rock. The question arises as to why s… ▽ More Spacecraft and ground-based observations show that the main rings of Saturn lack particles larger than 10 m. Tidal or collisional destruction of satellites/comets have been proposed as the origin of the main rings; however, Saturn's tide alone cannot grind km-sized fragments into submeter-sized particles because of the high mechanical strength of water ice and rock. The question arises as to why such large particles are not left in the current ring. It is known that thermal stress induced by diurnal and seasonal temperature variations can cause weathering and fragmentation of boulders and contribute to dust and regolith production on the Moon and terrestrial planets, and then such thermal stress can break particles larger than a critical radius while cannot smaller than the critical radius. In this study, we examined the role of thermal stress acting on Saturn's ring particles. We found that thermal stress can grind porous ring particles larger than 10-20 m, which explains the lack of particles larger than 10 m in Saturn's ring. Also, fragmentation by thermal stress can be adoptable for the Epsilon rings of Uranus. Furthermore, thermal stress caused by diurnal or seasonal temperature variation acting on boulders on surfaces of icy satellites and asteroids may play an important role in the evolution of their sizes. Our calculations explain the lack of boulders on icy satellites, except in the geologically active provinces such as the tiger stripes of Enceladus, where boulders are supplied by recent geological activity. We predict that future observations can find numerous boulders around Europa's geologically active cracks. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: 35 pages, 8 figures

Journal ref: Icarus Volume 378, 15 May 2022, 114919

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

The intrinsic structure of Sagittarius A* at 1.3 cm and 7 mm

Authors: Ilje Cho, Guang-Yao Zhao, Tomohisa Kawashima, Motoki Kino, Kazunori Akiyama, Michael D. Johnson, Sara Issaoun, Kotaro Moriyama, Xiaopeng Cheng, Juan-Carlos Algaba, Taehyun Jung, Bong Won Sohn, Thomas P. Krichbaum, Maciek Wielgus, Kazuhiro Hada, Ru-Sen Lu, Yuzhu Cui, Satoko Sawada-Satoh, Zhiqiang Shen, Jongho Park, Wu Jiang, Hyunwook Ro, Kunwoo Yi, Kiyoaki Wajima, Jee Won Lee , et al. (41 additional authors not shown)

Abstract: Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets to resolve the innermost region of SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multi-wavelength campaign of the Event Horizo… ▽ More Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets to resolve the innermost region of SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multi-wavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ~704$\pm$102 $μ$as (axial ratio $\sim$1.19$^{+0.24}_{-0.19}$) and $\sim$300$\pm$25 $μ$as (axial ratio $\sim$1.28$\pm$0.2) at 1.349 cm and 6.95 mm respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power-law, with an index $\sim$1.2$\pm$0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: accepted to ApJ

The appearance of a 'fresh' surface on 596 Scheila as a consequence of the 2010 impact event

Authors: Sunao Hasegawa, Michael Marsset, Francesca E. DeMeo, Schelte J. Bus, Masateru Ishiguro, Daisuke Kuroda, Richard P. Binzel, Josef Hanus, Akiko M. Nakamura, Bin Yang, Pierre Vernazza

Abstract: Dust emission was detected on main-belt asteroid 596 Scheila in December 2010, and attributed to the collision of a few-tens-of-meters projectile on the surface of the asteroid. In such impact, the ejected material from the collided body is expected to mainly comes from its fresh, unweathered subsurface. Therefore, it is expected that the surface of 596 was partially or entirely refreshed during t… ▽ More Dust emission was detected on main-belt asteroid 596 Scheila in December 2010, and attributed to the collision of a few-tens-of-meters projectile on the surface of the asteroid. In such impact, the ejected material from the collided body is expected to mainly comes from its fresh, unweathered subsurface. Therefore, it is expected that the surface of 596 was partially or entirely refreshed during the 2010 impact. By combining spectra of 596 from the literature and our own observations, we show that the 2010 impact event resulted in a significant slope change in the near-infrared (0.8 to 2.5 μm) spectrum of the asteroid, from moderately red (T-type) before the impact to red (D-type) after the impact. This provides evidence that red carbonaceous asteroids become less red with time due to space weathering, in agreement with predictions derived from laboratory experiments on the primitive Tagish Lake meteorite, which is spectrally similar to 596. This discovery provides the very first telescopic confirmation of the expected weathering trend of asteroids spectrally analog to Tagish Lake and/or anhydrous chondritic porous interplanetary dust particles. Our results also suggest that the population of implanted objects from the outer solar system is much larger than previously estimated in the main-belt, but many of these objects are hidden below their space-weathered surface. △ Less

Submitted 8 December, 2021; originally announced December 2021.

Comments: 11 pages, 5 figures, 1 tables, Accepted for publication in ApJ Letters

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

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)

East Asian VLBI Network Observations of Active Galactic Nuclei Jets: Imaging with KaVA+Tianma+Nanshan

Authors: Yuzhu Cui, Kazuhiro Hada, Motoki Kino, Bong Won Sohn, Jongho Park, Hyun Wook Ro, Satoko Sawada-Satoh, Wu Jiang, Lang Cui, Mareki Honma, Zhi Qiang Shen, Fumie Tazaki, Tao An, Ilje Cho, Guang Yao Zhao, Xiao Peng Cheng, Kotaro Niinuma, Kiyoaki Wajima, Ying Kang Zhang, Noriyuki Kawaguchi, Juan Carlos Algaba, Shoko Koyama, Tomoya Hirota, Yoshinori Yonekura, Nobuyuki Sakai , et al. (52 additional authors not shown)

Abstract: The East Asian very-long-baseline interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan, and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Explorati… ▽ More The East Asian very-long-baseline interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan, and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Exploration of Radio Astrometry (VERA) into KaVA, further expansion with the joint array in East Asia is actively promoted. Here we report the first imaging results (at 22 and 43 GHz) of bright radio sources obtained with KaVA connected to Tianma 65-m and Nanshan 26-m Radio Telescopes in China. To test the EAVN imaging performance for different sources, we observed four active galactic nuclei (AGN) having different brightness and morphology. As a result, we confirmed that Tianma 65-m Radio Telescope (TMRT) significantly enhances the overall array sensitivity, a factor of 4 improvement in baseline sensitivity and 2 in image dynamic range compared to the case of KaVA only. The addition of Nanshan 26-m Radio Telescope (NSRT) further doubled the east-west angular resolution. With the resulting high-dynamic-range, high-resolution images with EAVN (KaVA+TMRT+NSRT), various fine-scale structures in our targets, such as the counter-jet in M87, a kink-like morphology of the 3C273 jet and the weak emission in other sources, are successfully detected. This demonstrates the powerful capability of EAVN to study AGN jets and to achieve other science goals in general. Ongoing expansion of EAVN will further enhance the angular resolution, detection sensitivity and frequency coverage of the network. △ Less

Submitted 14 April, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: 19 pages, 9 figures, accepted by Research in Astronomy and Astrophysics (RAA)

Packing fraction of clusters formed in free-falling granular streams based on flash X-ray radiography

Authors: Yuuya Nagaashi, Akiko M. Nakamura, Sunao Hasegawa, Koji Wada

Abstract: We study the packing fraction of clusters in free-falling streams of spherical and irregularly shaped particles using flash X-ray radiography. The estimated packing fraction of clusters is low enough to correspond to coordination numbers less than 6. Such coordination numbers in numerical simulations correspond to aggregates that collide and grow without bouncing. Moreover, the streams of irregula… ▽ More We study the packing fraction of clusters in free-falling streams of spherical and irregularly shaped particles using flash X-ray radiography. The estimated packing fraction of clusters is low enough to correspond to coordination numbers less than 6. Such coordination numbers in numerical simulations correspond to aggregates that collide and grow without bouncing. Moreover, the streams of irregular particles evolved faster and formed clusters of larger sizes with lower packing fraction. This result on granular streams suggests that particle shape has a significant effect on the agglomeration process of granular materials. △ Less

Submitted 8 March, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: Accepted for publication in Physical Review E

Crater shape as a possible record of the impact environment of metallic bodies: Effects of temperature, impact velocity and impactor density

Authors: Ryo Ogawa, Akiko M. Nakamura, Ayako I. Suzuki, Sunao Hasegawa

Abstract: Metallic bodies that were the cores of differentiated bodies are sources of iron meteorites and are considered to have formed early in the terrestrial planet region before migrating to the main asteroid belt. Surface temperatures and mutual collision velocities differ between the terrestrial planet region and the main asteroid belt. To investigate the dependence of crater shape on temperature, vel… ▽ More Metallic bodies that were the cores of differentiated bodies are sources of iron meteorites and are considered to have formed early in the terrestrial planet region before migrating to the main asteroid belt. Surface temperatures and mutual collision velocities differ between the terrestrial planet region and the main asteroid belt. To investigate the dependence of crater shape on temperature, velocity and impactor density, we conducted impact experiments on room- and low-temperature iron meteorite and iron alloy targets (carbon steel SS400 and iron-nickel alloy) with velocities of 0.8-7 km/s. The projectiles were rock cylinders and metal spheres and cylinders. Oblique impact experiments were also conducted using stainless steel projectiles and SS400 steel targets which produced more prominent radial patterns downrange at room temperature than at low temperature. Crater diameters and depths were measured and compiled using non-dimensional parameter sets based on the $π$-group crater scaling relations. Two-dimensional numerical simulations were conducted using iSALE-2D code with the Johnson-Cook strength model. Both experimental and numerical results showed that the crater depth and diameter decreased with decreasing temperature, which strengthened the target, and with decreasing impact velocity. The decreasing tendency was more prominent for depth than for diameter, i.e., the depth/diameter ratio was smaller for the low temperature and low velocity conditions. The depth/diameter ratios of craters formed by rock projectiles were shallower than those of craters formed by metallic projectiles. Our results imply that the frequency distribution of the depth/diameter ratio for craters on the surface of metallic bodies may be used as a probe of the past impact environment of metallic bodies. △ Less

Submitted 15 March, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: Icarus, accepted (corrected version)

Cohesion of regolith: Measurements of meteorite powders

Authors: Yuuya Nagaashi, Takanobu Aoki, Akiko M. Nakamura

Abstract: The cohesion of particles has a significant effect on the properties of small bodies. In this study, we measured in open air, the cohesive forces of tens of micron-sized irregularly shaped meteorite, silica sand, glass powder, and spherical glass particles, using a centrifugal method. In addition, we estimated the amount of water vapor adsorbed on the particles under the measurement conditions. Th… ▽ More The cohesion of particles has a significant effect on the properties of small bodies. In this study, we measured in open air, the cohesive forces of tens of micron-sized irregularly shaped meteorite, silica sand, glass powder, and spherical glass particles, using a centrifugal method. In addition, we estimated the amount of water vapor adsorbed on the particles under the measurement conditions. The measured cohesive forces of the meteorite particles are tens of times smaller than those of an ideally spherical silica particle and correspond to the submicron-scale effective (or equivalent) curvature radius of the particle surface. Moreover, based on the estimated amount of water vapor adsorbed on the particles, we expect the cohesive forces of the particles in airless bodies to be approximately 10 times larger than those measured in open air. Based on the measurement results, we estimate that the cohesive forces of the particles on asteroids are typically in the sub-micro-Newton range, and that the particles on fast-rotating asteroids are tens of microns in size. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 32 pages, 9 figures, 7 tables. Accepted for publication in Icarus

Primordial porous structure of chondrite parent bodies due to self-gravity

Authors: Tomomi Omura, Akiko M. Nakamura

Abstract: The porosity of an asteroid is important when studying the evolution of our solar system through small bodies and for planning mitigation strategies to avoid disasters due to asteroid impacts. Our knowledge of asteroid porosity largely relies on meteorites sampled on Earth. However, chondrites sampled on Earth are suggested to be sorted by strength. In this study, we obtained an estimate of the mo… ▽ More The porosity of an asteroid is important when studying the evolution of our solar system through small bodies and for planning mitigation strategies to avoid disasters due to asteroid impacts. Our knowledge of asteroid porosity largely relies on meteorites sampled on Earth. However, chondrites sampled on Earth are suggested to be sorted by strength. In this study, we obtained an estimate of the most porous structure of primordial "granular" chondrite parent bodies based on measurements of the compaction behavior of chondrite component analogs. We measured compaction curves of dust and dust-bead mixture samples. The dust sample consisted of various spherical and irregular particles with diameters on the order of 10^0-10^1 $μ$m. The mixture sample consisted of dust and beads with different dust volume fractions (~0.2-1). We used 1.5 and 4.8 $μ$m particles as dust as a first step, although the typical size of materials in matrix may be much smaller. We approximated the compaction curve of each sample with a power-law form and calculated the porosity structure of the primordial chondrite parent bodies using the experimental results. Our results show that the primordial parent bodies are likely to have higher porosity than the chondrites. Moreover, the relatively higher volume fraction of the matrix may be one of the reasons why most meteorites with high porosity are carbonaceous chondrites. △ Less

Submitted 1 March, 2021; v1 submitted 23 January, 2021; originally announced January 2021.

Comments: 20 pages, 4 tables, 5 figures; typos corrected

Journal ref: Planet. Sci. J. 2 (2021) 41

Jet collimation and acceleration in the giant radio galaxy NGC 315

Authors: Jongho Park, Kazuhiro Hada, Masanori Nakamura, Keiichi Asada, Guang-Yao Zhao, Motoki Kino

Abstract: We study the collimation and acceleration of the jets in the nearby giant radio galaxy NGC 315, using multifrequency Very Long Baseline Array observations and archival High Sensitivity Array and Very Large Array data. We find that the jet geometry transitions from a semi-parabolic shape into a conical/hyperbolic shape at a distance of $\approx10^5$ gravitational radii. We constrain the frequency-d… ▽ More We study the collimation and acceleration of the jets in the nearby giant radio galaxy NGC 315, using multifrequency Very Long Baseline Array observations and archival High Sensitivity Array and Very Large Array data. We find that the jet geometry transitions from a semi-parabolic shape into a conical/hyperbolic shape at a distance of $\approx10^5$ gravitational radii. We constrain the frequency-dependent position of the core, from which we locate the jet base. The jet collimation profile in the parabolic region is in good agreement with the steady axisymmetric force-free electrodynamic solution for the outermost poloidal magnetic field line anchored to the black hole event horizon on the equatorial plane, similar to the nearby radio galaxies M87 and NGC 6251. The velocity field derived from the asymmetry in brightness between the jet and counterjet shows gradual acceleration up to the bulk Lorentz factor of $Γ\sim 3$ in the region where the jet collimation occurs, confirming the presence of the jet acceleration and collimation zone. These results suggest that the jets are collimated by the pressure of the surrounding medium and accelerated by converting Poynting flux to kinetic energy flux. We discover limb-brightening of the jet in a limited distance range where the angular resolution of our data is sufficient to resolve the jet transverse structure. This indicates that either the jet has a stratified velocity field of fast-inner and slow-outer layers or the particle acceleration process is more efficient in the outer layer due to the interaction with the surroundings on pc-scales. △ Less

Submitted 28 December, 2020; originally announced December 2020.

Comments: 32 pages, 18 figures, 5 tables, accepted for publication in ApJ

The dual nature of blazar fast variability. Space and ground observations of S5 0716+714

Authors: C. M. Raiteri, M. Villata, D. Carosati, E. Benítez, S. O. Kurtanidze, A. C. Gupta, D. O. Mirzaqulov, F. D'Ammando, V. M. Larionov, T. Pursimo, J. A. Acosta-Pulido, G. V. Baida, B. Balmaverde, G. Bonnoli, G. A. Borman, M. I. Carnerero, W. -P. Chen, V. Dhiman, A. Di Maggio, S. A. Ehgamberdiev, D. Hiriart, G. N. Kimeridze, O. M. Kurtanidze, C. S. Lin, J. M. Lopez , et al. (17 additional authors not shown)

Abstract: Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to i… ▽ More Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to investigate the spectral variability during the TESS pointing. The spectral analysis is further extended in frequency to the UV and X-ray bands with data from the Neil Gehrels Swift Observatory. We develop a new method to unveil the shortest optical variability time-scales. This is based on progressive de-trending of the TESS light curve by means of cubic spline interpolations through the binned fluxes, with decreasing time bins. The de-trended light curves are then analysed with classical tools for time-series analysis (periodogram, auto-correlation and structure functions). The results show that below 3 d there are significant characteristic variability time-scales of about 1.7, 0.5, and 0.2 d. Variability on time-scales < 0.2 d is strongly chromatic and must be ascribed to intrinsic energetic processes involving emitting regions, likely jet sub-structures, with dimension less than about milliparsec. In contrast, flux changes on time-scales > 0.5 d are quasi-achromatic and are probably due to Doppler factor changes of geometric origin. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 17 pages, 18 figures, accepted by MNRAS

A Jet-Bases Emission Model of the EHT 2017 Image of M87*

Authors: Tomohisa Kawashima, Kenji Toma, Motoki Kino, Kazunori Akiyama, Masanori Nakamura, Kotaro Moriyama

Abstract: We carry out general relativistic ray-tracing radiative-transfer calculations to study whether a localized emission from $e^{\pm}$ plasma rings created at the stagnation surface in the jet funnel, to which we refer as stagnation rings, can explain the ring image of M87 observed by Event Horizon Telescope (EHT) 2017. The resultant images consist of the direct image of the stagnation rings and the r… ▽ More We carry out general relativistic ray-tracing radiative-transfer calculations to study whether a localized emission from $e^{\pm}$ plasma rings created at the stagnation surface in the jet funnel, to which we refer as stagnation rings, can explain the ring image of M87 observed by Event Horizon Telescope (EHT) 2017. The resultant images consist of the direct image of the stagnation rings and the ring images formed via the strong deflection by the black-hole (BH) gravity, to which we refer as "quasi photon-ring". For the model with the BH spin $a_* = 0.99$, the direct image of the counter-jet and quasi photon-ring are almost coincident to the photon ring with diameter $\sim 40 μ{\rm as}$, while the approaching jet shows the small ring-image inside them. The synthetic observation image assuming the EHT 2017 array is consistent with that observed in M87, because the array is a bit sparse to detect the inner ring image. This indicates that the ring image in M87 might contain the important feature of the jet bases in addition to the photon ring. We find that forthcoming EHT observations can resolve the stagnation-ring image and may enable us to explore the plasma-injection mechanism into the jet funnel. △ Less

Submitted 21 December, 2020; v1 submitted 18 September, 2020; originally announced September 2020.

Comments: 8 pages, 3 figures, accepted for publication in ApJ

Linear polarization in the nucleus of M87 at 7 mm and 1.3 cm

Authors: E. V. Kravchenko, M. Giroletti, K. Hada, D. L. Meier, M. Nakamura, J. Park, R. C. Walker

Abstract: We report on high angular resolution polarimetric observations of the nearby radio galaxy M87 using the Very Long Baseline Array at 24 GHz ($λ=$1.3 cm) and 43 GHz ($λ=$7 mm) in 2017-2018. New images of the linear polarization substructure in the nuclear region are presented, characterized by a two-component pattern of polarized intensity and smooth rotation of the polarization plane around the 43… ▽ More We report on high angular resolution polarimetric observations of the nearby radio galaxy M87 using the Very Long Baseline Array at 24 GHz ($λ=$1.3 cm) and 43 GHz ($λ=$7 mm) in 2017-2018. New images of the linear polarization substructure in the nuclear region are presented, characterized by a two-component pattern of polarized intensity and smooth rotation of the polarization plane around the 43 GHz core. From a comparison with an analogous dataset from 2007, we find that this global polarization pattern remains stable on a time interval of 11 yr, while showing smaller month-scale variability. We discuss the possible Faraday rotation toward the M87 nucleus at centimeter to millimeter wavelengths. These results can be interpreted in a scenario where the observed polarimetric pattern is associated with the magnetic structure in the confining magnetohydrodynamic wind, which also serves as the source of the observed Faraday rotation. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: 9 pages, including 4 pages of Appendix; published in A&A Letters

A Long-lived Sharp Disruption on the Lower Clouds of Venus

Authors: J. Peralta, T. Navarro, C. W. Vun, A. Sánchez-Lavega, K. McGouldrick, T. Horinouchi, T. Imamura, R. Hueso, J. P. Boyd, G. Schubert, T. Kouyama, T. Satoh, N. Iwagami, E. F. Young, M. A. Bullock, P. Machado, Y. J. Lee, S. S. Limaye, M. Nakamura, S. Tellmann, A. Wesley, P. Miles

Abstract: Planetary-scale waves are thought to play a role in powering the yet-unexplained atmospheric superrotation of Venus. Puzzlingly, while Kelvin, Rossby and stationary waves manifest at the upper clouds (65--70 km), no planetary-scale waves or stationary patterns have been reported in the intervening level of the lower clouds (48--55 km), although the latter are probably Lee waves. Using observations… ▽ More Planetary-scale waves are thought to play a role in powering the yet-unexplained atmospheric superrotation of Venus. Puzzlingly, while Kelvin, Rossby and stationary waves manifest at the upper clouds (65--70 km), no planetary-scale waves or stationary patterns have been reported in the intervening level of the lower clouds (48--55 km), although the latter are probably Lee waves. Using observations by the Akatsuki orbiter and ground-based telescopes, we show that the lower clouds follow a regular cycle punctuated between 30$^{\circ}$N--40$^{\circ}$S by a sharp discontinuity or disruption with potential implications to Venus's general circulation and thermal structure. This disruption exhibits a westward rotation period of $\sim$4.9 days faster than winds at this level ($\sim$6-day period), alters clouds' properties and aerosols, and remains coherent during weeks. Past observations reveal its recurrent nature since at least 1983, and numerical simulations show that a nonlinear Kelvin wave reproduces many of its properties. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Comments: 21 pages, 10 figures, 2 animated figures and 2 tables

Journal ref: Geophysical Research Letters, Volume 47, Issue 11, 16 June 2020, e2020GL087221

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

Collisional disruption of highly porous targets in the strength regime: Effects of mixture

Authors: Yuichi Murakami, Akiko M. Nakamura, Koki Yokoyama, Yusuke Seto, Sunao Hasegawa

Abstract: Highly porous small bodies are thought to have been ubiquitous in the early solar system. Therefore, it is essential to understand the collision process of highly porous objects when considering the collisional evolution of primitive small bodies in the solar system. To date, impact disruption experiments have been conducted using high-porosity targets made of ice, pumice, and glass, and numerical… ▽ More Highly porous small bodies are thought to have been ubiquitous in the early solar system. Therefore, it is essential to understand the collision process of highly porous objects when considering the collisional evolution of primitive small bodies in the solar system. To date, impact disruption experiments have been conducted using high-porosity targets made of ice, pumice, and glass, and numerical simulations of impact fracture of porous bodies have also been conducted. However, a variety of internal structures of high-porosity bodies are possible. Therefore,laboratory experiments and numerical simulations in the wide parameter space are necessary. In this study, high-porosity targets of sintered hollow glass beads and targets made by mixing perlite with hollow beads were used in a collision disruption experiment to investigate the effects of the mixture on collisional destruction of high-porosity bodies. Among the targets prepared under the same sintering conditions, it was found that the targets with more impurities tend to have lower compressive strength and lower resistance against impact disruption. Further, destruction of the mixture targets required more impact energy density than would have been expected from compressive strength. It is likely that the perlite grains in the target matrix inhibit crack growth through the glass framework. The mass fraction of the largest fragment collapsed to a single function of a scaling parameter of energy density in the strength regime (Π_s) when assuming ratios of tensile strength to compressive strength based on a relationship obtained for ice-silicate mixtures. △ Less

Submitted 8 December, 2019; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: Planetary and Space Science,in press

Kinematics of the M87 jet in the collimation zone: gradual acceleration and velocity stratification

Authors: Jongho Park, Kazuhiro Hada, Motoki Kino, Masanori Nakamura, Jeffrey Hodgson, Hyunwook Ro, Yuzhu Cui, Keiichi Asada, Juan-Carlos Algaba, Satoko Sawada-Satoh, Sang-Sung Lee, Ilje Cho, Zhiqiang Shen, Wu Jiang, Sascha Trippe, Kotaro Niinuma, Bong Won Sohn, Taehyun Jung, Guang-Yao Zhao, Kiyoaki Wajima, Fumie Tazaki, Mareki Honma, Tao An, Kazunori Akiyama, Do-Young Byun , et al. (14 additional authors not shown)

Abstract: We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal s… ▽ More We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($Γ\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification. △ Less

Submitted 7 November, 2019; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: Accepted for publication in ApJ; This is a new version after updating Figure 11 in the manuscript

The Two-sided Jet Structures of NGC 1052 at Scales from 300 to $4 \times 10^7$ Schwarzschild Radii

Authors: Satomi Nakahara, Akihiro Doi, Yasuhiro Murata, Masanori Nakamura, Kazuhiro Hada, Keiichi Asada, Satoko Sawada-Satoh, Seiji Kameno

Abstract: We investigated the jet width profile with distance along the jet in the nearby radio galaxy NGC 1052 at radial distances between $\sim300$ to $4 \times 10^7$ Schwarzschild Radii($R_{\rm S}$) from the central engine on both their approaching and receding jet sides. The width of jets was measured in images obtained with the VLBI Space Observatory Programme (VSOP), the Very Long Baseline Array (VLBA… ▽ More We investigated the jet width profile with distance along the jet in the nearby radio galaxy NGC 1052 at radial distances between $\sim300$ to $4 \times 10^7$ Schwarzschild Radii($R_{\rm S}$) from the central engine on both their approaching and receding jet sides. The width of jets was measured in images obtained with the VLBI Space Observatory Programme (VSOP), the Very Long Baseline Array (VLBA), and the Very Large Array (VLA). The jet-width profile of receding jets are apparently consistent with that of approaching jets throughout the measuring distance ranges, indicating symmetry at least up to the sphere of gravitational influence of the central black hole. The power-law index $a$ of the jet-width profile ($w_{\rm{jet}} \propto r^{a}$, where $w_{\rm jet}$ is the jet width, $r$ is the distance from the central engine in the unit of $R_{\rm S}$) apparently shows a transition from $a \sim 0$ to $a \sim 1$, i.e., the cylindrical-to-conical jet structures, at a distance of $\sim1\times10^{4} \ R_{\mathrm{S}}$. The cylindrical jet shape at the small distances is reminiscent of the innermost jets in 3C 84. Both the central engines of NGC 1052 and 3C 84 are surrounded by dense material, part of which is ionized and causes heavy free-free absorption. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 8 figures, 2 tables, accepted for publication in AJ

Angle of repose of Martian wet sand using discrete element method: Implication for the seasonal cycle of recurring slope lineae(RSL) by relative humidity

Authors: Daigo Shoji, Shoko Imamura, Maya Nakamura, Rina Noguchi

Abstract: Seasonal cycle of Martian relative humidity shows that, in the southern and the northern RSL sites, humidity decreases when RSL increment while humidity is high on the fading season. We suggest that the recurrence of RSL is caused by the different angle of response between wet and dry grains. 2D DEM simulation indicates that, when the ratio of a liquid bridge volume to a particle volume is the ord… ▽ More Seasonal cycle of Martian relative humidity shows that, in the southern and the northern RSL sites, humidity decreases when RSL increment while humidity is high on the fading season. We suggest that the recurrence of RSL is caused by the different angle of response between wet and dry grains. 2D DEM simulation indicates that, when the ratio of a liquid bridge volume to a particle volume is the order of 10^{-7}, the repose angle of Martian sand can increase to ~40 degrees due to cohesion, which is near the slope around RSL starting area. Thus, if this amount of liquid bridge is generated by deliquescence or capillary condensation of vapor and if the water evaporates following the humidity and temperature cycles, Martian sand can induce avalanche on RSL slope. When grains are wet again, sand and dust deposit on RSL, which may be the mechanism of RSL disappearance and recharge of grains for the next season. In our model, because liquid induces RSL indirectly (RSL itself is the granular flow), the amount of liquid can be drastically reduced compared with the wet mechanism. In this case, as well as temperature, the seasonal cycle of relative humidity plays an important role in the recurrence of RSL. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: 22 pages, 5 figures, 1 table. Comments are welcome

Planetary-scale variations in winds and UV brightness at the Venusian cloud top: Periodicity and temporal evolution

Authors: Masataka Imai, Toru Kouyama, Yukihiro Takahashi, Atsushi Yamazaki, Shigeto Watanabe, Manabu Yamada, Takeshi Imamura, Takehiko Satoh, Masato Nakamura, Shin-ya Murakami, Kazunori Ogohara, Takeshi Horinouchi

Abstract: Planetary-scale waves at the Venusian cloud-top cause periodic variations in both winds and ultraviolet (UV) brightness. While the wave candidates are the 4-day Kelvin wave and 5-day Rossby wave with zonal wavenumber 1, their temporal evolutions are poorly understood. Here we conducted a time series analysis of the 365-nm brightness and cloud-tracking wind variations, obtained by the UV Imager onb… ▽ More Planetary-scale waves at the Venusian cloud-top cause periodic variations in both winds and ultraviolet (UV) brightness. While the wave candidates are the 4-day Kelvin wave and 5-day Rossby wave with zonal wavenumber 1, their temporal evolutions are poorly understood. Here we conducted a time series analysis of the 365-nm brightness and cloud-tracking wind variations, obtained by the UV Imager onboard the Japanese Venus Climate Orbiter Akatsuki from June to October 2017, revealing a dramatic evolution of planetary-scale waves and corresponding changes in planetary-scale UV features. We identified a prominent 5-day periodicity in both the winds and brightness variations, whose phase velocities were slower than the dayside mean zonal winds (or the super-rotation) by >35 m s$^{-1}$. The reconstructed planetary-scale vortices were nearly equatorially symmetric and centered at ~35° latitude in both hemispheres, which indicated that they were part of a Rossby wave. The amplitude of winds variation associated with the observed Rossby wave packet were amplified gradually over ~20 days and attenuated over ~50 days. Following the formation of the Rossby wave vortices, brightness variation emerges to form rippling white cloud belts in the 45°-60° latitudes of both hemispheres. ~3.8-day periodic signals were observed in the zonal wind and brightness variations in the equatorial region before the Rossby wave amplification. Although the amplitude and significance of the 3.8-day mode were relatively low in the observation season, this feature is consistent with a Kelvin wave, which may be the cause of the dark clusters in the equatorial region. △ Less

Submitted 18 September, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 44 pages, 19 figures

ALMA Observations of the Terahertz Spectrum of Sagittarius A*

Authors: Geoffrey C. Bower, Jason Dexter, Keiichi Asada, Christiaan D. Brinkerink, Heino Falcke, Paul Ho, Makoto Inoue, Sera Markoff, Daniel P. Marrone, Satoki Matsushita, Monika Moscibrodzka, Masanori Nakamura, Alison Peck, Ramprasad Rao

Abstract: We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $α\approx -0.3$, where the flux density $S \propto ν^α$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities… ▽ More We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $α\approx -0.3$, where the flux density $S \propto ν^α$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities $n = (2-5) \times 10^6$ cm$^{-3}$, electron temperatures $T_e = (1-3) \times 10^{11}$ K, and magnetic field strength $B = 10-50$ G. The parameter range can be further constrained using the observed quiescent X-ray luminosity. The flat submm spectrum results in a high electron temperature and implies that the emitting electrons are efficiently heated. We also find that the emission is most likely optically thin at 233 GHz. These results indicate that millimeter and submillimeter wavelength very long baseline interferometry of Sgr A* including those of the Event Horizon Telescope should see a transparent emission region down to event horizon scales. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Comments: Accepted for publication in ApJL

Global structure of thermal tides in the upper cloud layer of Venus revealed by LIR onboard Akatsuki

Authors: T. Kouyama, M. Taguchi, T. Fukuhara, T. Imamura, T. Horinouchi, T. M. Sato, S. Murakami, G. L. Hashimoto, Y. J. Lee, M. Futaguchi, T. Yamada, M. Akiba, T. Satoh, M. Nakamura

Abstract: Longwave Infrared Camera (LIR) onboard Akatsuki first revealed the global structure of the thermal tides in the upper cloud layer of Venus. The data were acquired over three Venusian years, and the analysis was done over the areas from the equator to the mid-latitudes in both hemispheres and over the whole local time. Thermal tides at two vertical levels were analyzed by comparing data at two diff… ▽ More Longwave Infrared Camera (LIR) onboard Akatsuki first revealed the global structure of the thermal tides in the upper cloud layer of Venus. The data were acquired over three Venusian years, and the analysis was done over the areas from the equator to the mid-latitudes in both hemispheres and over the whole local time. Thermal tides at two vertical levels were analyzed by comparing data at two different emission angles. Dynamical wave modes consisting of tides were identified; the diurnal tide consisted mainly of Rossby-wave and gravity-wave modes, while the semidiurnal tide predominantly consisted of a gravity-wave mode. The revealed vertical structures were roughly consistent with the above wave modes, but some discrepancy remained if the waves were supposed to be monochromatic. In turn, the heating profile that excites the tidal waves can be constrained to match this discrepancy, which would greatly advance the understanding of the Venusian atmosphere. △ Less

Submitted 18 August, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: 27 pages, 8 figures

Resolving the Innermost Geometry of Relativistic Jets in Active Galactic Nuclei

Authors: Juan Carlos Algaba, Masanori Nakamura, Keiichi Asada, Sang Sung Lee

Abstract: In the current paradigm, it is believed that the compact VLBI radio core of radio-loud active galactic nuclei (AGNi) represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with a roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to… ▽ More In the current paradigm, it is believed that the compact VLBI radio core of radio-loud active galactic nuclei (AGNi) represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with a roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral distribution properties and it has been recently found that, at least in some nearby radio galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNi. Data combined from the sources considered here are not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNi needs to be re-examined by millimeter/sub-millimeter VLBI observations. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Proceedings for the 14th European VLBI Network Symposium & Users Meeting (EVN 2018) 8-11 October 2018, Granada, Spain

Radio jet structures at ~100 parsec and larger scales of the $γ$-ray-emitting narrow-line Seyfert 1 galaxy PMN J0948+0022

Authors: Akihiro Doi, Satomi Nakahara, Masanori Nakamura, Motoki Kino, Nozomu Kawakatu, Hiroshi Nagai

Abstract: The narrow-line Seyfert 1 (NLS1) galaxy PMN J0948+0022 is an archetype of gamma-ray-emitting NLS1s in active galactic nuclei (AGNs). In this study, we investigate its radio structures using archival data obtained using the Karl G.~Jansky very large array (VLA) and the very long baseline array (VLBA). The new VLA images reveal an outermost radio emission separated by 9.1 arcsec. Its resolved struct… ▽ More The narrow-line Seyfert 1 (NLS1) galaxy PMN J0948+0022 is an archetype of gamma-ray-emitting NLS1s in active galactic nuclei (AGNs). In this study, we investigate its radio structures using archival data obtained using the Karl G.~Jansky very large array (VLA) and the very long baseline array (VLBA). The new VLA images reveal an outermost radio emission separated by 9.1 arcsec. Its resolved structure and steep spectrum suggest a terminal shock in a radio lobe energized by the jet from the PMN J0948+0022 nucleus. This large-scale radio component is found at almost the same position angle as that of the pc-scale jet, indicating a stable jet direction up to ~1 Mpc. Its apparent one-sidedness implies a moderate advancing speed ($β>0.27$). The kinematic age of 1 x 10^7 year of the jet activity is consistent with the expected NLS1 phase of ~10^7 year in the AGN lifetime. The VLBA image reveals the jet structure at distances ranging from r=0.82 milliarcsec to 3.5 milliarcsec, corresponding to approximately 100 pc, where superluminal motions were found. The jet width profile ($\propto r^{1.12}$) and flux-density profile ($\propto r^{-1.44}$) depending on the distance $r$ along the jet suggest that the jet kinetic energy is converted into internal energy in this region. The jet is causally connected in a nearly conical streamline, which is supported by ambient pressure at 100 pc scales in the host galaxy of PMN J0948+0022. △ Less

Submitted 14 May, 2019; originally announced May 2019.

Comments: 10 pages, 5 figures, accepted, accepted for publication in MNRAS

M87 black hole mass and spin estimate through the position of the jet boundary shape break

Authors: E. E. Nokhrina, L. I. Gurvits, V. S. Beskin, M. Nakamura, K. Asada, K. Hada

Abstract: We propose a new method of estimating a mass of a super massive black hole residing in the center of an active galaxy. The active galaxy M87 offers a convenient test case for the method due to the existence of a large amount of observational data on the jet and ambient environment properties in the central area of the object. We suggest that the observed transition of a jet boundary shape from a p… ▽ More We propose a new method of estimating a mass of a super massive black hole residing in the center of an active galaxy. The active galaxy M87 offers a convenient test case for the method due to the existence of a large amount of observational data on the jet and ambient environment properties in the central area of the object. We suggest that the observed transition of a jet boundary shape from a parabolic to a conical form is associated with the flow transiting from the magnetically dominated regime to the energy equipartition between plasma bulk motion and magnetic field. By coupling the unique set of observations available for the jet kinematics, environment and boundary profile with our MHD modelling under assumption on the presence of a dynamically important magnetic field in the M87 jet, we estimate the central black hole mass and spin. The method leads us to believe that the M87 super massive black hole has a mass somewhat larger than typically accepted so far. △ Less

Submitted 30 July, 2019; v1 submitted 11 April, 2019; originally announced April 2019.

Comments: 10 pages, 1 figure, 3 tables, accepted for publication by MNRAS

Relativistic Jets in the Accretion & Collimation Zone: New Challenges Enabled by New Instruments

Authors: Eric S. Perlman, Mark Birkinshaw, Matthias Kadler, Serguei Komissarov, Matthew Lister, David Meier, Eileen Meyer, Masanori Nakamura, Kristina Nyland, Christopher ODea, Diana Worrall, Andrzej Zdziarski

Abstract: Jets are a ubiquitous part of the accretion process, seen in a wide variety of objects ranging from active galaxies (AGN) to X-ray binary stars and even newly formed stars. AGN jets are accelerated by the supermassive black hole of their host galaxy by a coupling between the magnetic field and inflowing material. They are the source for many exciting phenomena and can profoundly influence the larg… ▽ More Jets are a ubiquitous part of the accretion process, seen in a wide variety of objects ranging from active galaxies (AGN) to X-ray binary stars and even newly formed stars. AGN jets are accelerated by the supermassive black hole of their host galaxy by a coupling between the magnetic field and inflowing material. They are the source for many exciting phenomena and can profoundly influence the larger galaxy and surrounding cluster. This White Paper points out what advances can be achieved in the field by new technologies, concentrating on the zone where jets are accelerated to relativistic speeds and collimated. The ngVLA and new space VLBI missions will give higher angular resolution, sensitivity and fidelity in the radio, penetrating this zone for additional objects and allowing us to resolve fundamental questions over the physics of jet acceleration and collimation. Interferometry in other bands would allow us to probe directly flaring components. We also emphasize the need for polarimetry, which is essential to revealing the role and configuration of magnetic fields. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 8 pages, 5 figures, submitted to Astro2020, the Decadal Survey of Astronomy & Astrophysics

Faraday rotation in the jet of M87 inside the Bondi radius: indication of winds from hot accretion flows confining the relativistic jet

Authors: Jongho Park, Kazuhiro Hada, Motoki Kino, Masanori Nakamura, Hyunwook Ro, Sascha Trippe

Abstract: We study Faraday rotation in the jet of M87 inside the Bondi radius using eight Very Long Baseline Array data sets, one at 8 GHz, four at 5 GHz, and three at 2 GHz. We obtain Faraday rotation measures (RMs) measured across the bandwidth of each data set. We find that the magnitude of RM systematically decreases with increasing distance from the black hole from 5,000 to 200,000 Schwarzschild radii.… ▽ More We study Faraday rotation in the jet of M87 inside the Bondi radius using eight Very Long Baseline Array data sets, one at 8 GHz, four at 5 GHz, and three at 2 GHz. We obtain Faraday rotation measures (RMs) measured across the bandwidth of each data set. We find that the magnitude of RM systematically decreases with increasing distance from the black hole from 5,000 to 200,000 Schwarzschild radii. The data, showing predominantly negative RM sign without significant difference of the RMs on the northern and southern jet edges, suggest that the spatial extent of the Faraday screen is much larger than the jet. We apply models of hot accretion flows, thought to be prevalent in active galactic nuclei having relatively low luminosity such as M87, and find that the decrease of RM is described well by a gas density profile $ρ\propto r^{-1}$. This behavior matches the theoretically expected signature of substantial winds, nonrelativistic un-collimated gas outflows from hot accretion flows, which is consistent with the results of various numerical simulations. The pressure profile inferred from the density profile is flat enough to collimate the jet, which can result in gradual acceleration of the jet in a magneto-hydrodynamical process. This picture is in good agreement with the observed gradual collimation and acceleration of the M87 jet inside the Bondi radius. The dominance of negative RMs suggests that jet and wind axis are misaligned such that the jet emission exposes only one side of the toroidal magnetic fields permeating the winds. △ Less

Submitted 20 December, 2018; originally announced December 2018.

Comments: 24 pages, 14 figures, 5 tables. accepted for publication in ApJ

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