Showing 1–50 of 135 results for author: Matsushita, S

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

Metadata for the Flux Density Calibration of the April 2018 Event Horizon Telescope Data

Authors: J. Y. Koay, C. Romero-Cañizales, L. D. Matthews, M. Janssen, L. Blackburn, R. P. J. Tilanus, J. Park, K. Asada, S. Matsushita, A. -K. Baczko, N. La Bella, C. -K. Chan, G. B. Crew, V. Fish, N. Patel, V. Ramakrishnan, H. Rottmann, J. Wagner, K. Wiik, P. Friberg, C. Goddi, S. Issaoun, G. Keating, J. Kim, T. P. Krichbaum , et al. (7 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observations carried out in 2018 April at 1.3 mm wavelengths included 9 stations in the array, comprising 7 single-dish telescopes and 2 phased arrays. The metadata package for the 2018 EHT observing campaign contains calibration tables required for the a-priori amplitude calibration of the 2018 April visibility data. This memo is the official documentation accomp… ▽ More The Event Horizon Telescope (EHT) observations carried out in 2018 April at 1.3 mm wavelengths included 9 stations in the array, comprising 7 single-dish telescopes and 2 phased arrays. The metadata package for the 2018 EHT observing campaign contains calibration tables required for the a-priori amplitude calibration of the 2018 April visibility data. This memo is the official documentation accompanying the release of the 2018 EHT metadata package, providing an overview of the contents of the package. We describe how telescope sensitivities, gain curves and other relevant parameters for each station in the EHT array were collected, processed, and validated to produce the calibration tables. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 26 pages, 7 figures, EHT Memo Series 2023-L1-01

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

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

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

Submitted 5 December, 2023; originally announced December 2023.

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

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

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

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

Submitted 29 August, 2023; originally announced August 2023.

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

Milliarcsecond Core Size Dependence of the Radio Variability of Blazars

Authors: Po Chih Hsu, Jun Yi Koay, Satoki Matsushita, Chorng-Yuan Hwang, Talvikki Hovatta, Sebastian Kiehlmann, Anthony Readhead, Walter Max-Moerbeck, Rodrigo Reeves

Abstract: Studying the long-term radio variability (timescales of months to years) of blazars enables us to gain a better understanding of the physical structure of these objects on sub-parsec scales, and the physics of super massive black holes. In this study, we focus on the radio variability of 1157 blazars observed at 15 GHz through the Owens Valley Radio Observatory (OVRO) Blazar Monitoring Program. We… ▽ More Studying the long-term radio variability (timescales of months to years) of blazars enables us to gain a better understanding of the physical structure of these objects on sub-parsec scales, and the physics of super massive black holes. In this study, we focus on the radio variability of 1157 blazars observed at 15 GHz through the Owens Valley Radio Observatory (OVRO) Blazar Monitoring Program. We investigate the dependence of the variability amplitudes and timescales, characterized based on model fitting to the structure functions, on the milliarcsecond core sizes measured by Very Long Baseline Interferometry. We find that the most compact sources at milliarcsecond scales exhibit larger variability amplitudes and shorter variability timescales than more extended sources. Additionally, for sources with measured redshifts and Doppler boosting factors, the correlation between linear core sizes against variability amplitudes and intrinsic timescales are also significant. The observed relationship between variability timescales and core sizes is expected, based on light travel-time arguments. This variability vs core size relation extends beyond the core sizes measured at 15 GHz; we see significant correlation between the 15 GHz variability amplitudes (as well as timescales) and core sizes measured at other frequencies, which can be attributed to a frequency-source size relationship arising from the intrinsic jet structure. At low frequencies of 1 GHz where the core sizes are dominated by interstellar scattering, we find that the variability amplitudes have significant correlation with the 1 GHz intrinsic core angular sizes, once the scatter broadening effects are deconvoluted from the intrinsic core sizes. △ Less

Submitted 20 August, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 17 pages, 12 figures

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

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

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

ALMA Measurement of 10 kpc-scale Lensing Power Spectra towards the Lensed Quasar MG J0414+0534

Authors: Kaiki Taro Inoue, Takeo Minezaki, Satoki Matsushita, Koichiro Nakanishi

Abstract: The lensing power spectra for gravitational potential, astrometric shift, and convergence perturbations are powerful probes to investigate dark matter structures on small scales. We report the first lower and upper bounds of these lensing power spectra on angular scale $\sim 1$arcsec towards the anomalous quadruply lensed quasar MG$\,$J0414+0534 at a redshift $z=2.639$. To obtain the spectra, we c… ▽ More The lensing power spectra for gravitational potential, astrometric shift, and convergence perturbations are powerful probes to investigate dark matter structures on small scales. We report the first lower and upper bounds of these lensing power spectra on angular scale $\sim 1$arcsec towards the anomalous quadruply lensed quasar MG$\,$J0414+0534 at a redshift $z=2.639$. To obtain the spectra, we conducted observations of MG$\,$J0414+0534 using the Atacama Large Millimeter/submillimeter Array (ALMA) with high angular resolution (0.02 arcsec to 0.05 arcsec). We developed a new partially non-parametric method in which Fourier coefficients of potential perturbation are adjusted to minimize the difference between linear combinations of weighted mean de-lensed images. Using positions of radio jet components, extended dust emission on scales $>1\,$kpc, and mid-infrared flux ratios, the range of measured convergence, astrometric shift, and potential powers at an angular scale of $\sim 1.1$arcsec (corresponding to an angular wave number of $l=1.2\times 10^6$ or $\sim 9\,$kpc in the primary lens plane) within $1\,σ$ are $\varDelta_κ=0.021-0.028$, $\varDelta_α=7-9\,$mas, and $\varDelta_ψ=1.2-1.6\,$$\textrm{mas}^2$, respectively. Our result is consistent with the predicted abundance of halos in the line of sight and subhalos in cold dark matter models. Our partially non-parametric lens models suggest a presence of a clump in the vicinity of object Y, a possible dusty dwarf galaxy and some small clumps in the vicinity of other lensed quadruple images. Although much fainter than the previous report, we detected weak continuum emission possibly from object Y with a peak flux of $\sim 100\,μ\textrm{Jy}\, \textrm{beam} ^{-1}$ at the $\sim 4\,σ$ level. △ Less

Submitted 26 July, 2023; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 54 pages, 35 figures, version accepted in ApJ. Analyses of models with a possible dwarf galaxy, object Y and that of chi^2 fit in the visibility plane were added. A new section 7 about CDM consistency was added. The angular scale of the smallest angular bin in Figure 24 was corrected

Journal ref: The Astrophysical Journal (2023) 954, 197 (33pp) (Open Access)

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)

A Giant Molecular Cloud Catalog in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

Authors: E. R. Miura, D. Espada, A. Hirota, C. Henkel, S. Verley, M. Kobayashi, S. Matsushita, F. P. Israel, B. Vila-Vilaro, K. Morokuma-Matsui, J. Ott, C. Vlahakis, A. B. Peck, S. Aalto, M. Hogerheijde, N. Neumayer, D. Iono, K. Kohno, H. Takemura, S. Komugi

Abstract: We present the first census of giant molecular clouds (GMCs) complete down to 10$^6 M_{\odot}$ and within the inner 4 kpc of the nearest giant elliptical and powerful radio galaxy, Centaurus A. We identified 689 GMCs using CO(1--0) data with 1" spatial resolution ($\sim 20$ pc) and 2 km/s velocity resolution obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). The $I$(CO)-… ▽ More We present the first census of giant molecular clouds (GMCs) complete down to 10$^6 M_{\odot}$ and within the inner 4 kpc of the nearest giant elliptical and powerful radio galaxy, Centaurus A. We identified 689 GMCs using CO(1--0) data with 1" spatial resolution ($\sim 20$ pc) and 2 km/s velocity resolution obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). The $I$(CO)-$N$(H$_2$) conversion factor based on the virial method is $X_{\rm CO}$ = $(2 \pm 1 )\times10^{20}$ cm$^{-2}$(K km/s)$^{-1}$ for the entire molecular disk, consistent with that of the disks of spiral galaxies including the Milky Way, and $X_{\rm CO}$ = $(5 \pm 2)\times10^{20}$ cm$^{-2}$(K km/s)$^{-1}$ for the circumnuclear disk (CND, within a galactocentric radius of 200 pc). We obtained the GMC mass spectrum distribution and find that the best-truncated power-law fit for the whole molecular disk, with index $γ\simeq -2.41 \pm 0.02$ and upper cutoff mass $\sim 1.3 \times 10^{7} M_{\odot}$, is also in agreement with that of nearby disk galaxies. A trend is found in the mass spectrum index from steep to shallow as we move to inner radii. Although the GMCs are in an elliptical galaxy, the general GMC properties in the molecular disk are as in spiral galaxies. However, in the CND, large offsets in the line-width-size scaling relations ($\sim$ 0.3 dex higher than those in the GMCs in the molecular disk), a different $X_{\rm CO}$ factor, and the shallowest GMC mass distribution shape ($γ= -1.1 \pm 0.2$) all suggest that there the GMCs are most strongly affected by the presence of the AGN and/or shear motions. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: 19 pages, 11 figures, 3 tables, accepted for publication in MNRAS

Black hole mass measurement using ALMA observations of [CI] and CO emissions in the Seyfert 1 galaxy NGC7469

Authors: Dieu D. Nguyen, Takuma Izumi, Sabine Thater, Masatoshi Imanishi, Taiki Kawamuro, Shunsuke Baba, Suzuka Nakano, Jean L. Turner, Kotaro Kohno, Satoki Matsushita, Sergio Martin, David S. Meier, Phuong M. Nguyen, Lam T. Nguyen

Abstract: We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-${\rm [CI]}$(1-0) and molecular-$^{12}$CO(1-0) emission lines at the spatial resolution of $\approx0.3$" (or $\approx$ 100 pc). These emissions reveal that NGC7469 hosts a circumnuclear gas disc (CND) with a ring-like stru… ▽ More We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-${\rm [CI]}$(1-0) and molecular-$^{12}$CO(1-0) emission lines at the spatial resolution of $\approx0.3$" (or $\approx$ 100 pc). These emissions reveal that NGC7469 hosts a circumnuclear gas disc (CND) with a ring-like structure and a two-arm/bi-symmetric spiral pattern within it, surrounded by a starbursting ring. The CND has a relatively low $σ/V\approx0.35$ ($r\sim0.5$") and $\sim0.19$ ($r>0.5"$), suggesting that the gas is dynamically settled and suitable for dynamically deriving the mass of its central source. As is expected from X-ray dominated region (XDR) effects that dramatically increase an atomic carbon abundance by dissociating CO molecules, we suggest that the atomic [CI](1-0) emission is a better probe of SMBH masses than CO emission in AGNs. Our dynamical model using the ${\rm [CI]}$(1-0) kinematics yields a $M_{\rm BH}=1.78^{+2.69}_{-1.10}\times10^7$M$_\odot$ and $M/L_{\rm F547M}=2.25^{+0.40}_{-0.43}$ (M$_\odot$/L$_\odot$). The model using the CO(1-0) kinematics also gives a consistent $M_{\rm BH}$ with a larger uncertainty, up to an order of magnitude, i.e.\ $M_{\rm BH}=1.60^{+11.52}_{-1.45}\times10^7$M$_\odot$. This newly dynamical $M_{\rm BH}$ is $\approx$ 2 times higher than the mass determined from the reverberation mapped (RM) method using emissions arising in the unresolved broad-line region (BLR). Given this new $M_{\rm BH}$, we are able to constrain the specific RM dimensionless scaling factor of $f=7.2^{+4.2}_{-3.4}$ for the AGN BLR in NGC7469. The gas within the unresolved BLR thus has a Keplerian virial velocity component and the inclination of $i\approx11.0^\circ$$_{-2.5}^{+2.2}$, confirming its face-on orientation in a Seyfert 1 AGN by assuming a geometrically thin BLR model. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 22 pages, 16 figures, 7 tables. Accepted for publication on MNRAS

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

ALMA High-frequency Long-baseline Campaign in 2017: A Comparison of the Band-to-band and In-band Phase Calibration Techniques and Phase-calibrator Separation Angles

Authors: Luke T. Maud, Yoshiharu Asaki, Edward B. Fomalont, William R. F. Dent, Akihiko Hirota, Satoki Matsushita, Neil M. Phillips, John M. Carpenter, Satoko Takahashi, Eric Villard, Tsuyoshi Sawada, Stuartt Corder

Abstract: The Atacama Large millimeter/submillimeter Array (ALMA) obtains spatial resolutions of 15 to 5 milli-arcsecond (mas) at 275-950GHz (0.87-0.32mm) with 16km baselines. Calibration at higher-frequencies is challenging as ALMA sensitivity and quasar density decrease. The Band-to-Band (B2B) technique observes a detectable quasar at lower frequency that is closer to the target, compared to one at the ta… ▽ More The Atacama Large millimeter/submillimeter Array (ALMA) obtains spatial resolutions of 15 to 5 milli-arcsecond (mas) at 275-950GHz (0.87-0.32mm) with 16km baselines. Calibration at higher-frequencies is challenging as ALMA sensitivity and quasar density decrease. The Band-to-Band (B2B) technique observes a detectable quasar at lower frequency that is closer to the target, compared to one at the target high-frequency. Calibration involves a nearly constant instrumental phase offset between the frequencies and the conversion of the temporal phases to the target frequency. The instrumental offsets are solved with a differential-gain-calibration (DGC) sequence, consisting of alternating low and high frequency scans of strong quasar. Here we compare B2B and in-band phase referencing for high-frequencies ($>$289GHz) using 2-15km baselines and calibrator separation angles between $\sim$0.68 and $\sim$11.65$^{\circ}$. The analysis shows that: (1) DGC for B2B produces a coherence loss $<$7% for DGC phase RMS residuals $<$30$^{\circ}$. (2) B2B images using close calibrators ( $<$1.67$^{\circ}$ ) are superior to in-band images using distant ones ( $>$2.42$^{\circ}$ ). (3) For more distant calibrators, B2B is preferred if it provides a calibrator $\sim$2$^{\circ}$ closer than the best in-band calibrator. (4) Decreasing image coherence and poorer image quality occur with increasing phase calibrator separation angle because of uncertainties in the antenna positions and sub-optimal phase referencing. (5) To achieve $>$70% coherence for long-baseline (16 km) band 7 (289GHz) observations, calibrators should be within $\sim$4$^{\circ}$ of the target. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 17 Figures, 10 Tables

Journal ref: 2020ApJS..250...18M

ALMA Observations of Multiple-CO and C Lines Toward the Active Galactic Nucleus of NGC 7469: X-Ray-dominated Region Caught in the Act

Authors: Takuma Izumi, Dieu D. Nguyen, Masatoshi Imanishi, Taiki Kawamuro, Shunsuke Baba, Suzuka Nakano, Kotaro Kohno, Satoki Matsushita, David S. Meier, Jean L. Turner, Tomonari Michiyama, Nanase Harada, Sergio Martín, Kouichiro Nakanishi, Shuro Takano, Tommy Wiklind, Naomasa Nakai, Pei-Ying Hsieh

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surround… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $α_{\rm CO} = 4.1$ and $α_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments. △ Less

Submitted 3 July, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 25 pages, 12 figures, 8 tables. Accepted for publication in ApJ. Minor updates in a replacement

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

The MALATANG Survey: Dense Gas and Star Formation from High Transition HCN and HCO+ maps of NGC253

Authors: Xue-Jian Jiang, Thomas R. Greve, Yu Gao, Zhi-Yu Zhang, Qinghua Tan, Richard de Grijs, Luis C. Ho, Michal J. Michalowski, Malcolm J. Currie, Christine D. Wilson, Elias Brinks, Yiping Ao, Yinghe Zhao, Jinhua He, Nanase Harada, Chentao Yang, Qian Jiao, Aeree Chung, Bumhyun Lee, Matthew W. L. Smith, Daizhong Liu, Satoki Matsushita, Yong Shi, Masatoshi Imanishi, Mark G. Rawlings , et al. (4 additional authors not shown)

Abstract: To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission feat… ▽ More To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $Σ_\text{stellar}$, and SFE$_\text{dense}$ versus $Σ_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $Σ_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars. △ Less

Submitted 14 March, 2020; originally announced March 2020.

Comments: accepted to MNRAS

ALMA 50-parsec resolution imaging of jet-ISM interaction in the lensed quasar MGJ0414+0534

Authors: Kaiki Taro Inoue, Satoki Matsushita, Kouichiro Nakanishi, Takeo Minezaki

Abstract: We report our high-resolution (0.03 arcsec to 0.07 arcsec) ALMA imaging of the quadruply lensed radio-loud quasar MG\,J0414+0534 at redshift $z=2.639$ in the continuum and the broad CO(11-10) line at $\sim 340\,$GHz. With the help of strong lensing magnification and ALMA's high-resolution, we succeeded in resolving the jet/dust and CO gas in the quasar host galaxy both extending up to $\sim 1\,$kp… ▽ More We report our high-resolution (0.03 arcsec to 0.07 arcsec) ALMA imaging of the quadruply lensed radio-loud quasar MG\,J0414+0534 at redshift $z=2.639$ in the continuum and the broad CO(11-10) line at $\sim 340\,$GHz. With the help of strong lensing magnification and ALMA's high-resolution, we succeeded in resolving the jet/dust and CO gas in the quasar host galaxy both extending up to $\sim 1\,$kpc, with a resolution of $\sim 50\,$pc for the first time. Both the continuum emission and the CO(11-10) line have a similar bimodal structure aligned with the quasar jets ($\sim 200\,$pc) observed by Very Long Baseline Interferometry (VLBI) at $5\,$GHz and $8.4\,$GHz. The CO gas in the vicinity of both the eastern and western jet components at the location of $\sim 80\,$pc from the quasar core are moving at high velocities up to $\pm 600\,$$\textrm{km}\,\textrm{s}^{-1}$ relative to the core. The observed features show clear evidence of strong interaction between the jets and interstellar medium (ISM). High temperature and high density environments in the ISM of the quasar host galaxy, as suggested from CO line ratios, also support this result. The small scale of the jets, the jet-ISM interaction, and the continuum spectral energy distribution of this source indicate that we are watching the infancy stage of quasar radio activity. △ Less

Submitted 28 March, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 7 pages, 4 figures, typos corrected

Journal ref: The Astrophysical Journal Letters, 892, Issue 2, id.L18 (online publication on March 27, 2020)

A More Efficient Search for H2O Megamaser Galaxies : The Power of the X-ray and Mid-infrared Photometry

Authors: C. Y. Kuo, J. Y. Hsiang, H. H. Chung, A. Constantin, Y. -Y. Chang, E. da Cunha, D. Pesce, W. T. Chien, B. Y. Chen, J. A. Braatz, Ingyin Zaw, S. Matsushita, J. C. Lin

Abstract: We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk… ▽ More We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk maser) detection rate increases sharply for galaxies with 12 micron AGN luminosity L^{AGN}_{12 micron} greater than 10^42 erg/s, from ~<3%(~<2%) to ~12%(~5%). By using the ratio of the observed X-ray to mid-IR AGN luminosity as an indicator of N_H, we also find that the megamaser (disk maser) detection rates are boosted to 15%(7%) and 20%(9%) for galaxies with N_H >= 10^23 cm^{-2} and N_H >= 10^{24} cm^{-2}, respectively. Combining these column density cuts with a constraint for high L^{AGN}_{12 micron} (>=10^42 erg/s) predicts further increases in the megamaser (disk maser) detection rates to 19%(8%) and 27%(14%), revealing unprecedented potential boosts of the megamaser and disk maser detection rates by a factor of 7-15 relative to the current rates, depending on the chosen sample selection criteria. A noteworthy aspect of these new predictions is that the completeness rates are only compromised mildly, with the rates remaining at the level of ~95%(~50%) for sources with N_H >= 10^{23} cm^{-2} (N_H >= 10^24 cm^-2). Applying these selection methods to the current X-ray AGN surveys predicts the detection of >~15 new megamaser disks. △ Less

Submitted 15 February, 2020; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 22 pages, 11 figures, 4 tables, Accepted by ApJ

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

Star Formation Efficiencies at Giant Molecular Cloud Scales in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

Authors: D. Espada, S. Verley, R. E. Miura, F. P. Israel, C. Henkel, S. Matsushita, B. Vila-Vilaro, J. Ott, K. Morokuma-Matsui, A. B. Peck, A. Hirota, S. Aalto, A. C. Quillen, M. R. Hogerheijde, N. Neumayer, C. Vlahakis, D. Iono, K. Kohno

Abstract: We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previ… ▽ More We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $τ=$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas. △ Less

Submitted 16 December, 2019; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: Minor update. 20 pages, 11 figures, accepted for publication in ApJ

The first bird's-eye view of a gravitationally unstable accretion disk in high-mass star formation

Authors: Kazuhito Motogi, Tomoya Hirota, Masahiro N. Machida, Yoshinori Yonekura, Mareki Honma, Shigehisa Takakuwa, Satoki Matsushita

Abstract: We report on the first bird's-eye view of the innermost accretion disk around the high-mass protostellar object G353.273+0.641, taken by Atacama Large Millimter/submillimeter Array long-baselines. The disk traced by dust continuum emission has a radius of 250 au, surrounded by the infalling rotating envelope traced by thermal CH$_3$OH lines. This disk radius is consistent with the centrifugal radi… ▽ More We report on the first bird's-eye view of the innermost accretion disk around the high-mass protostellar object G353.273+0.641, taken by Atacama Large Millimter/submillimeter Array long-baselines. The disk traced by dust continuum emission has a radius of 250 au, surrounded by the infalling rotating envelope traced by thermal CH$_3$OH lines. This disk radius is consistent with the centrifugal radius estimated from the specific angular momentum in the envelope. The lower-limit envelope mass is $\sim$5-7 M$_{\odot}$ and accretion rate onto the stellar surface is 3 $\times$ 10$^{-3}$ M$_{\odot}$ yr$^{-1}$ or higher. The expected stellar age is well younger than 10$^{4}$ yr, indicating that the host object is one of the youngest high-mass objects at present. The disk mass is 2-7 M$_{\odot}$, depending on the dust opacity index. The estimated Toomre's $Q$ parameter is typically 1-2 and can reach 0.4 at the minimum. These $Q$ values clearly satisfy the classical criteria for the gravitational instability, and are consistent with the recent numerical studies. Observed asymmetric and clumpy structures could trace a spiral arm and/or disk fragmentation. We found that 70$\%$ of the angular momentum in the accretion flow could be removed via the gravitational torque in the disk. Our study has indicated that the dynamical nature of a self-gravitating disk could dominate the early phase of high-mass star formation. This is remarkably consistent with the early evolutionary scenario of a low-mass protostar. △ Less

Submitted 30 May, 2019; originally announced May 2019.

Comments: Accepted to ApJL, 11 pages, 3 figures, 2 tables

Precipitable Water Vapor, Temperature, and Wind Statistics At Sites Suitable for mm and Submm Wavelength Astronomy in Northern Chile

Authors: Angel Otarola, Carlos De Breuck, Tony Travouillon, Satoki Matsushita, Lars-Åke Nyman, Al Wootten, Simon J. E. Radford, Marc Sarazin, Florian Kerber, Juan P. Pérez-Beaupuits

Abstract: Atmospheric water vapor is the main limiting factor of atmospheric transparency in the mm and submm wavelength spectral windows. Thus, dry sites are needed for the installation and successful operation of radio astronomy observatories exploiting those spectral windows. Temperature and wind are variables of special consideration when planning the installation and operation of large-aperture radio t… ▽ More Atmospheric water vapor is the main limiting factor of atmospheric transparency in the mm and submm wavelength spectral windows. Thus, dry sites are needed for the installation and successful operation of radio astronomy observatories exploiting those spectral windows. Temperature and wind are variables of special consideration when planning the installation and operation of large-aperture radio telescopes, as these parameters affect the mechanical response of radio telescopes exposed to the environmental conditions. Temperature, and in particular temperature gradients, induce thermal deformation of mechanical structures, while high wind speeds and gusts induce pointing jitter affecting the tracking accuracy of astronomical sources during the observations. This work summarizes the statistics of precipitable water vapor (PWV), temperature, and wind monitored at sites by the coastal mountain range, as well as on the west slope of the Andes mountain range in the region of Antofagasta, Chile. This information could prove useful for the planning of extended baselines for the Atacama Large Millimeter and Submillimeter Array (ALMA), and/or new radio telescope projects, such as the Atacama Large Aperture Submm/ mm Telescope (AtLAST) initiative. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 28 pages, 13 Figures. Manuscript accepted under Open Access terms to the Publications of the Astronomical Society of the Pacific (PASP) Journal

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

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

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

Submitted 25 September, 2018; originally announced September 2018.

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

Molecular gas and Star Formation Properties in Early Stage Mergers: SMA CO(2-1) Observations of the LIRGs NGC 3110 and NGC 232

Authors: D. Espada, S. Martin, S. Verley, A. R. Pettitt, S. Matsushita, M. Argudo-Fernandez, Z. Randriamanakoto, P. -Y. Hsieh, T. Saito, R. E. Miura, Y. Kawana, J. Sabater, L. Verdes-Montenegro, P. T. P. Ho, R. Kawabe, D. Iono

Abstract: Mergers of galaxies are an important mode for galaxy evolution because they serve as an efficient trigger of powerful starbursts. However, observational studies of the molecular gas properties during their early stages are scarce. We present interferometric CO(2-1) maps of two luminous infrared galaxies (LIRGs), NGC 3110 and NGC 232, obtained with the Submillimeter Array (SMA) with ~ 1 kpc resolut… ▽ More Mergers of galaxies are an important mode for galaxy evolution because they serve as an efficient trigger of powerful starbursts. However, observational studies of the molecular gas properties during their early stages are scarce. We present interferometric CO(2-1) maps of two luminous infrared galaxies (LIRGs), NGC 3110 and NGC 232, obtained with the Submillimeter Array (SMA) with ~ 1 kpc resolution. While NGC 3110 is a spiral galaxy interacting with a minor (14:1 stellar mass) companion, NGC 232 is interacting with a similarly sized object. We find that such interactions have likely induced in these galaxies enhancements in the molecular gas content and central concentrations, partly at the expense of atomic gas. The obtained molecular gas surface densities in their circumnuclear regions are $Σ_{\rm mol}~\gtrsim10^{2.5}$ M$_\odot$ pc$^{-2}$, higher than in non-interacting objects by an order of magnitude. Gas depletion times of ~ 0.5 - 1 Gyr are found for the different regions, lying in between non-interacting disk galaxies and the starburst sequence. In the case of NGC 3110, the spiral arms show on average 0.5 dex shorter depletion times than in the circumnuclear regions if we assume a similar H$_2$-CO conversion factor. We show that even in the early stages of the interaction with a minor companion, a starburst is formed along the circumnuclear region and spiral arms, where a large population of SSCs is found (~350), and at the same time a large central gas concentration is building up which might be the fuel for an active galactic nucleus. The main morphological properties of the NGC 3110 system are reproduced by our numerical simulations and allow us to estimate that the current epoch of the interaction is at ~ 150 Myrs after closest approach. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 30 pages. Accepted for publication in ApJ

Lightning black holes as unidentified TeV sources

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

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

Submitted 9 August, 2018; originally announced August 2018.

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

Electronics Instrumentation for the Greenland Telescope

Authors: Derek Kubo, Chih-Chiang Han, Hiroaki Nishioka, Ryan Chilson, Ranjani Srinivasan, Sheng- Feng Yen, Kuo-Chieh Fu, Homin Jiang, Kuan-Yu Liu, Ta-Shun Wei, Chih-Wei Huang, Chen- Yu Yu, Peter Oshiro, Shu-Hao Chang, Chung-Cheng Chen, Philippe Raffin, Yau-De Huang, Pierre Martin-Cocher, Ming-Tang Chen, Makoto Inoue, Satoki Matsushita, Keiichi Asada, Shoko Koyama, Patrick Koch, Paul T. P. Ho , et al. (15 additional authors not shown)

Abstract: The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope… ▽ More The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope is currently at Thule Air Base in Greenland with the long-term goal to move the telescope over the Greenland ice sheet to Summit Station. The GLT currently has a single cryostat which houses three dual polarization receivers that cover 84-96 GHz, 213-243 GHz and 271-377 GHz bands. A hydrogen maser frequency source in conjunction with high frequency synthesizers are used to generate the local oscillator references for the receivers. The intermediate frequency outputs of each receiver cover 4-8 GHz and are heterodyned to baseband for digitization within a set of ROACH-2 units then formatted for recording onto Mark-6 data recorders. A separate set of ROACH-2 units operating in parallel provides the function of auto-correlation for real-time spectral analysis. Due to the stringent instrumental stability requirements for interferometry a diagnostic test system was incorporated into the design. Tying all of the above equipment together is the fiber optic system designed to operate in a low temperature environment and scalable to accommodate a larger distance between the control module and telescope for Summit Station. A report on the progress of the above electronics instrumentation system will be provided. △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: 19 pages, 37 figures, oral presentation at the June 2018 SPIE conference: Astronomical Telescopes + Instrumentation, in the section on Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy IX

The MALATANG Survey: the L_gas-L_IR correlation on sub-kiloparsec scale in six nearby star-forming galaxies as traced by HCN J=4-3 and HCO^+ J=4-3

Authors: Qing-Hua Tan, Yu Gao, Zhi-Yu Zhang, Thomas R. Greve, Xue-Jian Jiang, Christine D. Wilson, Chen-Tao Yang, Ashley Bemis, Aeree Chung, Satoki Matsushita, Yong Shi, Yi-Ping Ao, Elias Brinks, Malcolm J. Currie, Timothy A. Davis, Richard de Grijs, Luis C. Ho, Masatoshi Imanishi, Kotaro Kohno, Bumhyun Lee, Harriet Parsons, Mark G. Rawlings, Dimitra Rigopoulou, Erik Rosolowsky, Joanna Bulger , et al. (21 additional authors not shown)

Abstract: We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where… ▽ More We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 μm/100 μm flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies. △ Less

Submitted 15 May, 2018; originally announced May 2018.

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

An ALMA view of star formation efficiency suppression in early-type galaxies after gas-rich minor mergers

Authors: Freeke van de Voort, Timothy A. Davis, Satoki Matsushita, Kate Rowlands, Stanislav S. Shabala, James R. Allison, Yuan-Sen Ting, Anne E. Sansom, Paul P. van der Werf

Abstract: Gas-rich minor mergers contribute significantly to the gas reservoir of early-type galaxies (ETGs) at low redshift, yet the star formation efficiency (SFE; the star formation rate divided by the molecular gas mass) appears to be strongly suppressed following some of these events, in contrast to the more well-known merger-driven starbursts. We present observations with the Atacama Large Millimeter/… ▽ More Gas-rich minor mergers contribute significantly to the gas reservoir of early-type galaxies (ETGs) at low redshift, yet the star formation efficiency (SFE; the star formation rate divided by the molecular gas mass) appears to be strongly suppressed following some of these events, in contrast to the more well-known merger-driven starbursts. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of six ETGs, which have each recently undergone a gas-rich minor merger, as evidenced by their disturbed stellar morphologies. These galaxies were selected because they exhibit extremely low SFEs. We use the resolving power of ALMA to study the morphology and kinematics of the molecular gas. The majority of our galaxies exhibit spatial and kinematical irregularities, such as detached gas clouds, warps, and other asymmetries. These asymmetries support the interpretation that the suppression of the SFE is caused by dynamical effects stabilizing the gas against gravitational collapse. Through kinematic modelling we derive high velocity dispersions and Toomre Q stability parameters for the gas, but caution that such measurements in edge-on galaxies suffer from degeneracies. We estimate merger ages to be about 100~Myr based on the observed disturbances in the gas distribution. Furthermore, we determine that these galaxies lie, on average, two orders of magnitude below the Kennicutt-Schmidt relation for star-forming galaxies as well as below the relation for relaxed ETGs. We discuss potential dynamical processes responsible for this strong suppression of star formation surface density at fixed molecular gas surface density. △ Less

Submitted 24 January, 2018; originally announced January 2018.

Comments: Accepted for publication in MNRAS

Journal ref: MNRAS, 2018, Volume 476, Issue 1, p.122-132

1000 AU Exterior Arcs Connected to the Protoplanetary Disk around HL Tau

Authors: Hsi-Wei Yen, Shigehisa Takakuwa, You-Hua Chu, Naomi Hirano, Paul T. P. Ho, Kazuhiro D. Kanagawa, Chin-Fei Lee, Hauyu Baobab Liu, Sheng-Yuan Liu, Tomoaki Matsumoto, Satoki Matsushita, Takayuki Muto, Kazuya Saigo, Ya-Wen Tang, Alfonso Trejo, Chun-Ju Wu

Abstract: The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resoluti… ▽ More The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: Accepted by A&A

Journal ref: A&A 608, A134 (2017)

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

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

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

Submitted 27 July, 2017; originally announced July 2017.

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

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

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

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

Submitted 25 July, 2017; originally announced July 2017.

Comments: 5 figures

On the disappearance of a cold molecular torus around the low-luminosity active galactic nucleus of NGC 1097

Authors: T. Izumi, K. Kohno, K. Fathi, E. Hatziminaoglou, R. I. Davies, S. Martín, S. Matsushita, E. Schinnerer, D. Espada, S. Aalto, K. Onishi, J. L. Turner, M. Imanishi, K. Nakanishi, D. S. Meier, K. Wada, N. Kawakatu, T. Nakajima

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $μ$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: 9 pages, 5 figures. Accepted for publication in ApJL

Phase correction for ALMA - Investigating water vapour radiometer scaling:The long-baseline science verification data case study

Authors: L. T. Maud, R. P. J. Tilanus, T. A. van Kempen, M. R. Hogerheijde, M. Schmalzl, I. Yoon, Y. Contreras, M. C. Toribio, Y. Asaki, W. R. F. Dent, E. Fomalont, S. Matsushita

Abstract: The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in… ▽ More The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data. △ Less

Submitted 11 July, 2017; originally announced July 2017.

Comments: 27 pages, 11 figures, 5 tables - Accepted to A&A

Journal ref: A&A 605, A121 (2017)

ALMA Observations of the Gravitational Lens SDP.9

Authors: Kenneth C. Wong, Tsuyoshi Ishida, Yoichi Tamura, Sherry H. Suyu, Masamune Oguri, Satoki Matsushita

Abstract: We present long-baseline ALMA observations of the strong gravitational lens H-ATLAS J090740.0-004200 (SDP.9), which consists of an elliptical galaxy at $z_{\mathrm{L}}=0.6129$ lensing a background submillimeter galaxy into two extended arcs. The data include Band 6 continuum observations, as well as CO $J$=6$-$5 molecular line observations, from which we measure an updated source redshift of… ▽ More We present long-baseline ALMA observations of the strong gravitational lens H-ATLAS J090740.0-004200 (SDP.9), which consists of an elliptical galaxy at $z_{\mathrm{L}}=0.6129$ lensing a background submillimeter galaxy into two extended arcs. The data include Band 6 continuum observations, as well as CO $J$=6$-$5 molecular line observations, from which we measure an updated source redshift of $z_{\mathrm{S}}=1.5747$. The image morphology in the ALMA data is different from that of the HST data, indicating a spatial offset between the stellar, gas, and dust component of the source galaxy. We model the lens as an elliptical power law density profile with external shear using a combination of archival HST data and conjugate points identified in the ALMA data. Our best model has an Einstein radius of $θ_{\mathrm{E}}=0.66\pm0.01$ and a slightly steeper than isothermal mass profile slope. We search for the central image of the lens, which can be used constrain the inner mass distribution of the lens galaxy including the central supermassive black hole, but do not detect it in the integrated CO image at a 3$σ$ rms level of 0.0471 Jy km s$^{-1}$. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: Accepted for publication in ApJL; 6 pages, 2 figures, 3 tables

Disentangling the Circumnuclear Environs of Centaurus A: III. An Inner Molecular Ring, Nuclear Shocks and the CO to warm H2 interface

Authors: D. Espada, S. Matsushita, R. E. Miura, F. P. Israel, N. Neumayer, S. Martin, C. Henkel, T. Izumi, D. Iono, S. Aalto, J. Ott, A. B. Peck, A. C. Quillen, K. Kohno

Abstract: We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist with… ▽ More We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist within the CND, which form a ring-like structure with an unprojected diameter of 9 x 6 arcsec (162pc x 108pc) and a position angle PA = 155deg. Inside the nuclear ring, there are two leading and straight filamentary structures with lengths of about 30-60pc at PA = 120deg on opposite sides of the AGN, with a rotational symmetry of 180deg and steeper position-velocity diagrams, which are interpreted as nuclear shocks due to non-circular motions. Along the filaments, and unlike other nearby AGNs, several dense molecular clumps present low HCN/HCO+(4-3) ratios (~0.5). The filaments abruptly end in the probed transitions at r = 20pc from the AGN, but previous near-IR H2 (J=1-0) S(1) maps show that they continue in an even ~1000 K), winding up in the form of nuclear spirals, and forming an inner ring structure with another set of symmetric filaments along the N-S direction and within r = 10pc. The molecular gas is governed primarily by non-circular motions, being the successive shock fronts at different scales where loss of angular momentum occurs, a mechanism which may feed efficiently powerful radio galaxies down to parsec scales. △ Less

Submitted 18 June, 2017; originally announced June 2017.

Comments: 46 pages. Accepted for publication in ApJ

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

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

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

Submitted 9 June, 2017; originally announced June 2017.

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

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

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

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

Submitted 12 May, 2017; originally announced May 2017.

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

Luminous Infrared Galaxies with the Submillimeter Array V: Molecular Gas in Intermediate to Late-stage Mergers

Authors: Kazimierz Sliwa, Christine D. Wilson, Satoki Matsushita, Alison B. Peck, Glen R. Petitpas, Toshiki Saito, Min Yun

Abstract: We present new high-resolution ALMA (13CO J=1-0 and J= 2-1) and CARMA (12CO and 13CO J=1-0) observations of two Luminous Infrared Galaxies (LIRGs): Arp 55 and NGC 2623. The new data are complementary to published and archival Submillimeter Array observations of 12CO J=2-1 and J=3-2. We perform a Bayesian likelihood non-local thermodynamic equilibrium analysis to constrain the molecular gas physica… ▽ More We present new high-resolution ALMA (13CO J=1-0 and J= 2-1) and CARMA (12CO and 13CO J=1-0) observations of two Luminous Infrared Galaxies (LIRGs): Arp 55 and NGC 2623. The new data are complementary to published and archival Submillimeter Array observations of 12CO J=2-1 and J=3-2. We perform a Bayesian likelihood non-local thermodynamic equilibrium analysis to constrain the molecular gas physical conditions such as temperature, column and volume densities and the [12CO]/[13CO] abundance ratio. For Arp 55, an early/intermediate staged merger, the line measurements are consistent with cold (~10-20 K), dense (>10$^{3.5}$ cm$^{-3}$) molecular gas. For NGC 2623, the molecular gas is warmer (~110 K) and less dense (~10$^{2.7}$ cm$^{-3}$). Since Arp 55 is an early/intermediate stage merger while NGC 2623 is a merger remnant, the difference in physical conditions may be an indicator of merger stage. Comparing the temperature and volume density of several LIRGs shows that the molecular gas, averaged over ~kpc scale, of advanced mergers is in general warmer and less dense than early/intermediate stage mergers. We also find that the [12CO]/[13CO] abundance ratio of NGC 2623 is unusually high (>250) when compared to the Milky Way; however, it follows a trend seen with other LIRGs in literature. This high [12CO]/[13CO] value is very likely due to stellar nucleosynthesis enrichment of the interstellar medium. On the other hand, Arp 55 has a more Galactic [12CO]/[13CO] value with the most probable [12CO]/[13CO] value being 20-30. We measure the CO-to-H2 conversion factor, $α_{CO}$, to be ~0.1 and ~0.7 (3x10$^{-4}$/x$_{CO}$) M$_{\odot}$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$ for Arp 55 and NGC 2623, respectively. Since Arp 55 is an early/intermediate , this suggests that the transition from a Galactic conversion factor to a LIRG value happens at an even earlier merger stage. △ Less

Submitted 21 March, 2017; originally announced March 2017.

Comments: 19 pages, 9 figures, Accepted for publication in ApJ

Evidence for a Dusty Dark Dwarf Galaxy in the Quadruple Lens MG0414+0534

Authors: Kaiki Taro Inoue, Satoki Matsushita, Takeo Minezaki, Masashi Chiba

Abstract: We report the $4 \, σ$ detection of a faint object with a flux of ~ 0.3 mJy, in the vicinity of the quadruply lensed QSO MG0414+0534 using the Atacama Large Millimeter/submillimeter array (ALMA) Band 7. The object is most probably a dusty dark dwarf galaxy, which has not been detected in either the optical, near-infrared (NIR) or radio (cm) bands. An anomaly in the flux ratio of the lensed images… ▽ More We report the $4 \, σ$ detection of a faint object with a flux of ~ 0.3 mJy, in the vicinity of the quadruply lensed QSO MG0414+0534 using the Atacama Large Millimeter/submillimeter array (ALMA) Band 7. The object is most probably a dusty dark dwarf galaxy, which has not been detected in either the optical, near-infrared (NIR) or radio (cm) bands. An anomaly in the flux ratio of the lensed images observed in Band 7 and the mid-infrared (MIR) band and the reddening of the QSO light color can be simultaneously explained if we consider the object as a lensing substructure with an ellipticity ~ 0.7 at a redshift of $0.5 \lesssim z \lesssim 1$. Using the best-fit lens models with three lenses, we find that the dark matter plus baryon mass associated with the object is $\sim 10^9\, M_{\odot}$, the dust mass is $\sim 10^7\,M_{\odot}$ and the linear size is $\gtrsim 5\,$kpc. Thus our findings suggest that the object is a dusty dark dwarf galaxy. A substantial portion of faint submillimeter galaxies (SMGs) in the universe may be attributed to such dark objects. △ Less

Submitted 18 January, 2017; originally announced January 2017.

Comments: 8 pages, 4 figures, accepted in ApJL

Journal ref: The Astrophysical Journal Letters, 835, Issue 2, article id. L23, 6 pp. (published on Jan. 30, 2017)

ALMA Long Baseline Campaigns: Phase Characteristics of Atmosphere at Long Baselines in the Millimeter and Submillimeter Wavelengths

Authors: Satoki Matsushita, Yoshiharu Asaki, Edward B. Fomalont, Koh-Ichiro Morita, Denis Barkats, Richard E. Hills, Ryohei Kawabe, Luke T. Maud, Bojan Nikolic, Remo P. J. Tilanus, Catherine Vlahakis, Nicholas D. Whyborn

Abstract: This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallowe… ▽ More This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallower (average of ~0.2-0.3) at baseline lengths longer than ~1 km, namely showing a turn-over in SSF. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV >1 mm, which reduces the degree of phase fluctuations by a factor of two in many cases. However, phase fluctuations still remain after the WVR phase correction, suggesting the existence of other turbulent constituent that cause the phase fluctuation. This is supported by occasional SSFs that do not exhibit any turn-over; these are only seen when the PWV is low or after WVR phase correction. This means that the phase fluctuation caused by this turbulent constituent is inherently smaller than that caused by water vapor. Since there is no turn-over in the SSF up to the maximum baseline length of ~15 km, this turbulent constituent must have scale height of 10 km or more, and thus cannot be water vapor, whose scale height is around 1 km. This large scale height turbulent constituent is likely to be water ice or a dry component. Excess path length fluctuation after the WVR phase correction at a baseline length of 10 km is large (>200 micron), which is significant for high frequency (>450 GHz or <700 micron) observations. These results suggest the need for an additional phase correction method, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching, and the result suggests that it works well, with shorter cycle times linearly improving the coherence. △ Less

Submitted 11 January, 2017; originally announced January 2017.

Comments: 12 pages, 7 figures. PASP, in press

The Greenland Telescope: Antenna Retrofit Status and Future Plans

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

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

Submitted 9 December, 2016; originally announced December 2016.

Ground-Based Mid-Infrared Study of the Compton-Thick AGN in M51 at 10-100 pc Scale

Authors: Youichi Ohyama, Satoki Matsushita, Nagisa Oi, Ai-Lei Sun

Abstract: We performed near-diffraction-limited (~0.4 FWHM) N-band imaging of one of the nearest Active Galactic Nucleus (AGN) in M51 with 8.2m Subaru telescope to study the nuclear structure and spectral energy distribution (SED) at 8-13 um. We found that the nucleus is composed of an unresolved core (at ~13 pc resolution, or intrinsic size corrected for the instrumental effect of <6 pc) and an extended ha… ▽ More We performed near-diffraction-limited (~0.4 FWHM) N-band imaging of one of the nearest Active Galactic Nucleus (AGN) in M51 with 8.2m Subaru telescope to study the nuclear structure and spectral energy distribution (SED) at 8-13 um. We found that the nucleus is composed of an unresolved core (at ~13 pc resolution, or intrinsic size corrected for the instrumental effect of <6 pc) and an extended halo (at a few tens pc scale), and each of their SEDs is almost flat. We examined the SED by comparing with the archival Spitzer IRS spectrum processed to mimic our chopping observation of the nucleus, and the published radiative-transfer model SEDs of the AGN clumpy dusty torus. The halo SED is likely due to circumnuclear star formation showing little Polycyclic Aromatic Hydrocarbon (PAH) emission due to the AGN. The core SED is likely dominated by the AGN because of the following two reasons. Firstly, the clumpy torus model SEDs can reproduce the red mid-infrared continuum with apparently moderate silicate 9.7 um absorption. Secondly, the core 12 um luminosity and the absorption-corrected X-ray luminosity at 2-10 keV in the literature follow the mid-infrared-X-ray luminosity correlation known for the nearby AGNs including the Compton-thick ones. △ Less

Submitted 9 December, 2016; originally announced December 2016.

Comments: 16 pages, 4 figures, accepted for publication in ApJ

3.5-Year Monitoring of 225 GHz Opacity at the Summit of Greenland

Authors: Satoki Matsushita, Keiichi Asada, Pierre L. Martin-Cocher, Ming-Tang Chen, Paul T. P. Ho, Makoto Inoue, Patrick M. Koch, Scott N. Paine, David D. Turner

Abstract: We present the 3.5-yr monitoring results of 225 GHz opacity at the summit of the Greenland ice sheet (Greenland Summit Camp) at an altitude of 3200 m using a tipping radiometer. We chose this site as our submillimeter telescope (Greenland Telescope; GLT) site, because its location offers favorable baselines to existing submillimeter telescopes for global-scale VLBI. The site shows a clear seasonal… ▽ More We present the 3.5-yr monitoring results of 225 GHz opacity at the summit of the Greenland ice sheet (Greenland Summit Camp) at an altitude of 3200 m using a tipping radiometer. We chose this site as our submillimeter telescope (Greenland Telescope; GLT) site, because its location offers favorable baselines to existing submillimeter telescopes for global-scale VLBI. The site shows a clear seasonal variation with the average opacity lower by a factor of two during winter. For the winter quartiles of 25% and 50%, the Greenland site is about 10%-30% worse than the ALMA or the South Pole sites. Estimated atmospheric transmission spectra in winter season are similar to the ALMA site at lower frequencies (<450 GHz), which are transparent enough to perform astronomical observations almost all of the winter time with opacities <0.5, but 10%-25% higher opacities at higher frequencies (>450 GHz) than those at the ALMA site. This is due to the lower altitude of the Greenland site. Nevertheless, half of the winter time at the Greenland site can be used for astronomical observations at frequencies between 450 GHz and 1000 GHz with opacities <1.2, and 10% of the time show >10% transmittance in the THz (1035 GHz, 1350 GHz, and 1500 GHz) windows. One major advantage of the Greenland site in winter is that there is no diurnal variation due to the polar night condition, and therefore the durations of low-opacity conditions are significantly longer than at the ALMA site. Opacities lower than 0.05 or 0.04 can continue for more than 100 hours. Such long stable opacity conditions do not occur as often even at the South Pole; it happens only for the opacity lower than 0.05. Since the opacity variation is directly related to the sky temperature (background) variation, the Greenland site is suitable for astronomical observations that need unusually stable sky background. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Comments: 9 pages, 8 figures. Accepted in PASP

Journal ref: PASP, 129, 025001 (2017)

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

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

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

Submitted 25 October, 2016; originally announced October 2016.

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

The Molecular Baryon Cycle of M82

Authors: John Chisholm, Satoki Matsushita

Abstract: Baryons cycle into galaxies from the inter-galactic medium, are converted into stars, and a fraction of the baryons are ejected out of galaxies by stellar feedback. Here we present new high resolution (3.9"; 68 pc) CO(2-1) and CO(3-2) images that probe these three stages of the baryon cycle in the nearby starburst M 82. We combine these new observations with previous CO(1-0) and [Fe II] images to… ▽ More Baryons cycle into galaxies from the inter-galactic medium, are converted into stars, and a fraction of the baryons are ejected out of galaxies by stellar feedback. Here we present new high resolution (3.9"; 68 pc) CO(2-1) and CO(3-2) images that probe these three stages of the baryon cycle in the nearby starburst M 82. We combine these new observations with previous CO(1-0) and [Fe II] images to study the physical conditions within the molecular gas. Using a Bayesian analysis and the radiative transfer code RADEX, we model molecular Hydrogen temperatures and densities, as well as CO column densities. Besides the disc, we concentrate on two regions within the galaxy: an expanding super-bubble and the base of a molecular streamer. Shock diagnostics, kinematics, and optical extinction suggest that the streamer is an inflowing filament, with a molecular gas mass inflow rate of 3.5 M$_\odot$ yr$^{-1}$. We measure the molecular gas mass outflow rate of the expanding super-bubble to be 17 M$_\odot$ yr$^{-1}$, 5 times higher than the inferred inflow rate, and 1.3 times the star formation rate of the galaxy. The high mass outflow rate and large star formation rate will deplete the galaxy of molecular gas within eight million years, unless there are additional sources of molecular gas. △ Less

Submitted 2 August, 2016; originally announced August 2016.

Comments: 14 pages, 15 figures, accepted for publication in ApJ

The Fossil Nuclear Outflow in the Central 30 pc of the Galactic Center

Authors: Pei-Ying Hsieh, Paul T. P. Ho, Chorng-Yuan Hwang, Yoshito Shimajiri, Satoki Matsushita, Patrick M. Koch, Daisuke Iono

Abstract: We report a new 1-pc (30") resolution CS($J=2-1$) line map of the central 30 pc of the Galactic Center (GC), made with the Nobeyama 45m telescope. We revisit our previous study of the extraplanar feature called polar arc (PA), which is a molecular cloud located above SgrA* with a velocity gradient perpendicular to the Galactic plane. We find that the PA can be traced back to the Galactic disk. Thi… ▽ More We report a new 1-pc (30") resolution CS($J=2-1$) line map of the central 30 pc of the Galactic Center (GC), made with the Nobeyama 45m telescope. We revisit our previous study of the extraplanar feature called polar arc (PA), which is a molecular cloud located above SgrA* with a velocity gradient perpendicular to the Galactic plane. We find that the PA can be traced back to the Galactic disk. This provides clues of the launching point of the PA , roughly $6\times10^{6}$ years ago. Implications of the dynamical time scale of the PA might be related to the Galactic Center Lobe (GCL) at parsec scale. Our results suggest that in the central 30 pc of the GC, the feedback from past explosions could alter the orbital path of the molecular gas down to the central tenth of parsec. In the follow-up work of our new CS($J=2-1$) map, we also find that near the systemic velocity, the molecular gas shows an extraplanar hourglass-shaped feature (HG-feature) with a size of $\sim$13 pc. The latitude-velocity diagrams show that the eastern edge of the HG-feature is associated with an expanding bubble B1, $\sim$7 pc away from SgrA*. The dynamical time scale of this bubble is $\sim3\times10^{5}$ years. This bubble is interacting with the 50 km s$^{-1}$ cloud. Part of the molecular gas from the 50 km s$^{-1}$ cloud was swept away by the bubble to $b=-0.2deg$. The western edge of the HG-feature seems to be the molecular gas entrained from the 20 km s$^{-1}$ cloud towards the north of the Galactic disk. Our results suggest a fossil explosion in the central 30 pc of the GC a few 10$^{5}$ years ago. △ Less

Submitted 13 July, 2016; originally announced July 2016.

Comments: 38 pages,12 figures, accepted for ApJ

Linearly Polarized Millimeter and Submillimeter Continuum Emission of Sgr A* Constrained by ALMA

Authors: Hauyu Baobab Liu, Melvyn C. H. Wright, Jun-Hui Zhao, Christiaan D. Brinkerink, Paul T. P. Ho, Elisabeth A. C. Mills, Sergio Martín, Heino Falcke, Satoki Matsushita, Ivan Martí-Vidal

Abstract: Our aim is to characterize the polarized continuum emission properties including intensity, polarization position angle, and polarization percentage of Sgr A* at $\sim$100 (3.0 mm), $\sim$230 (1.3 mm), $\sim$345 (0.87 mm), $\sim$500 (0.6 mm), and $\sim$700 GHz (0.43 mm). We report continuum emission properties of Sgr A* at the above frequency bands, based on the Atacama Large Millimeter Array (ALM… ▽ More Our aim is to characterize the polarized continuum emission properties including intensity, polarization position angle, and polarization percentage of Sgr A* at $\sim$100 (3.0 mm), $\sim$230 (1.3 mm), $\sim$345 (0.87 mm), $\sim$500 (0.6 mm), and $\sim$700 GHz (0.43 mm). We report continuum emission properties of Sgr A* at the above frequency bands, based on the Atacama Large Millimeter Array (ALMA) observations. We measured flux densities of Sgr A* from ALMA single pointing and mosaic observations. We performed sinusoidal fittings to the observed (XX-YY)/I intensity ratios, to derive the polarization position angles and polarization percentages. We successfully detect polarized continuum emission from all observed frequency bands. We observed lower Stokes I intensity at $\sim$700 GHz than that at $\sim$500 GHz, which suggests that emission at $\gtrsim$500 GHz is from optically thin part of a synchrotron emission spectrum. Both the Stokes I intensity and the polarization position angle at our highest observing frequency of $\sim$700 GHz, may be varying with time. However, we do not yet detect variation in the polarization percentage at $>$500 GHz. The polarization percentage at $\sim$700 GHz is likely lower than that at $\sim$500 GHz. By comparing the $\sim$500 GHz and $\sim$700 GHz observations with the observations at lower frequency bands, we suggest that the intrinsic polarization position angle of Sgr A* is varying with time. This paper also reports the measurable polarization properties from the observed calibration quasars. The future simultaneous multi-frequency polarization observations are required for clarifying the time and frequency variation of polarization position angle and polarization percentage. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: 11 pages, 10 figures, submitted to A&A; first referee report received. arXiv admin note: text overlap with arXiv:1604.00599