Search

In 2019 September, a sudden flare of the 6.7 GHz methanol maser was observed toward the high-mass young stellar object (HMYSO) G24.33+0.14. This may represent the fourth detection of a transient mass accretion event in an HMYSO after S255IR NIRS3, NGC 6334I-MM1, and G358.93−0.03-MM1. G24.33+0.14 is unique among these sources as it clearly shows a repeating flare with an 8 yr interval. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the millimeter continuum and molecular lines toward G24.33+0.14 in the pre-flare phase in 2016 August (ALMA Cycle 3) and the mid-flare phase in 2019 September (ALMA Cycle 6). We identified three continuum sources in G24.33+0.14, and the brightest source, C1, which is closely associated with the 6.7 GHz maser emission, shows only a marginal increase in flux density with a flux ratio (Cycle 6$/$Cycle 3) of 1.16 ± 0.01, considering an additional absolute flux calibration uncertainty of $10\%$. We identified 26 transitions from 13 molecular species other than methanol, and they exhibit similar levels of flux differences with an average flux ratio of 1.12 ± 0.15. In contrast, eight methanol lines observed in Cycle 6 are brighter than those in Cycle 3 with an average flux ratio of 1.23 ± 0.13, and the higher excitation lines tend to show a larger flux increase. If this systematic increasing trend is real, it would suggest radiative heating close to the central HMYSO due to an accretion event which could expand the size of the emission region and/or change the excitation conditions. Given the low brightness temperatures and small flux changes, most of the methanol emission is likely to be predominantly thermal, except for the 229.759 GHz (8−1–70 E) line known as a class I methanol maser. The flux change in the millimeter continuum of G24.33+0.14 is smaller than in S255IR NIRS3 and NGC 6334I-MM1 but is comparable with that in G358.93−0.03-MM1, suggesting different amounts of accreted mass in these events.

Powerful relativistic jets in active galactic nuclei (AGN) are believed to be originated from the accretion of material onto the supermassive black hole (SMBH)1. The nearby radio galaxy M87 is one of the best examples of this phenomenon, and the recent detection of a black-hole shadow with the Event Horizon Telescope provided compelling evidence for this paradigm2. However, whether the central SMBH in M87 has a spin or not, a fundamental parameter of a BH along with the mass, remains unconstrained by observations since the appearance of the photon ring is rather insensitive to the spin2-4. A clue to this challenge is to trace the long-term evolution of the innermost jet base where the central BH and accretion disk tightly regulate its flow dynamics. Here we report an extensive analysis of the morphological evolution of the M87 jet nozzle based on 170 high-resolution radio images spanning more than two decades. The ensemble of the data reveals a periodic (a period of 11 years) variation of the jet direction with a peak-to-peak amplitude of ~10 deg on the sky. A successful fit of a precession model to the data as well as general-relativistic magnetohydrodynamic simulations strongly supports the presence of a spinning BH in M87, inducing a Lense-Thirring precession of a misaligned accretion disk by the frame-dragging effect5. Our results suggest that jet precession with a spinning BH can commonly exist in other more distant AGN but had evaded detection due to the small magnitude and long period of position angle variation.

High-mass stars are thought to accumulate much of their mass via short, infrequent bursts of disk-aided accretion. Such accretion events are rare and difficult to observe directly but are known to drive enhanced maser emission. In this Letter we report high-resolution, multi-epoch methanol maser observations toward G358.93-0.03 which reveal an interesting phenomenon; the sub-luminal propagation of a thermal radiation "heat-wave" emanating from an accreting high-mass proto-star. The extreme transformation of the maser emission implies a sudden intensification of thermal infrared radiation from within the inner (40 mas, 270 au) region. Subsequently, methanol masers trace the radial passage of thermal radiation through the environment at $\geq$ 4-8\% the speed of light. Such a high translocation rate contrasts with the $\leq$ 10 km s$^{-1}$ physical gas motions of methanol masers typically observed using very long baseline interferometry (VLBI). The observed scenario can readily be attributed to an accretion event in the high-mass proto-star G358.93-0.03-MM1. While being the third case in its class, G358.93-0.03-MM1 exhibits unique attributes hinting at a possible `zoo' of accretion burst types. These results promote the advantages of maser observations in understanding high-mass star formation, both through single-dish maser monitoring campaigns and via their international cooperation as VLBI arrays., Comment: Published in Nature Astronomy in 2020

In the near future, two radio telescopes will be installed in the Indonesian region. These telescopes are proposed to be integrated into the existing Very Long Baseline Interferometry (VLBI) network both for astronomical and geodetical purposes. Here we simulate the impact of the inclusion of the future Indonesian antennas to the estimates of Earth Orientation Parameters and the station position. The simulation was performed in two types of VLBI sessions. First, we analyse the contribution of Indonesian antennas to the existing intensive session INT3 (IN320-314), which focuses on the estimation of dUT1. We found that the addition of Indonesian antennas has reduced the estimated dUT1 repeatability value by about 25%. Next, we simulate the 24-hour session by considering two existing network configurations, which are R4 (R4934) and AOV (AOV049). Overall we found that the addition of Indonesian antennas to each network configuration decreases the repeatability value of the Earth Orientation Parameter by about 20%. Meanwhile, the repeatability value of station position is reduced up to 12%. This reduction is already achieved even when we include only one Indonesian antenna.

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

We present the discovery of the simultaneous flux variation of a massive young stellar object (MYSO) G036.70+00.09 (G036.70) both in the maser emission and mid-infrared (MIR;$\lambda=3$--$5$~$\mu$m) bands. Utilizing the ALLWISE and NEOWISE archival databases covering a long time span of approximately 10 years with a cadence of 6 months, we confirmed that G036.70 indicates a stochastic year-long MIR variability with no signs of the WISE band color change of W1 (3.4~$\mu$m) $-$W2 (4.6~$\mu$m). Cross-matching the MIR data set with the high-cadence 6.7~GHz class II methanol maser flux using a Hitachi 32-m radio telescope that discovered its periodicity in the methanol maser of 53.0--53.2 days, we also determine the flux correlations between the two bands at two different timescales, year-long and day-long, both of which have never been reported in MYSOs except when they are in a state of the accretion burst phase. The results of our study support the scenario that a class II methanol maser is pumped up by infrared emission from accreting disks of MYSOs. We also discuss the possible origins of MIR and maser variability. To explain the two observed phenomena, a stochastic year-long MIR variability with no signs of significant color change and maser-MIR variability correlation, change in mass accretion rate and line-of-sight extinction because of nonaxisymmetric dust density distribution in a rotating accretion disk are possible origins. Observations through spectroscopic monitoring of accretion-related emission lines are essential for determining the origin of the observed variability in G036.70., Comment: 10 pages, 6 figures, accepted for publication in ApJ

Numerical simulations have predicted that substructures such as spiral arms can be produced through a gravitationally unstable disk around high-mass young stellar objects (HMYSOs). Recent high-resolution observations from the Atacama Large Millimeter/submillimeter Array have investigated these substructures at a spatial resolution of ∼100 au. An accretion burst, which is a manifestation of an increase in the accretion rate caused by a gravitational instability in the disk, can result in luminosity outbursting phenomena - as has been seen in several HMYSOs. However, no clear relationship between the accretion bursts and disk substructures has been established. Here we report the detections of three new molecular maser species, HDO, HNCO and 13CH3OH, from the direction of the HMYSO G358.93-0.03 during a 6.7 GHz methanol maser flaring event. High-quality imaging of the three new maser species exhibits consistent observational evidence that these masers closely trace the spiral-arm substructures around this HMYSO. The rapid decay of the spectral lines emitted from these molecules suggests that these are transient phenomena (for only ∼1 month), probably associated with rapid changes in radiation field due to an accretion burst. Therefore, these new maser species provide evidence linking the spiral-arm substructure with an accretion burst, both expected from massive disk instabilities.

We report the earliest detection of the 19.967-GHz [transition 21–30E (t = 0)] methanol maser associated with the massive star-forming region G358.93–0.03. The flare was detectable from 2019 January 23 to March 5, for only 44 d. It turned out to be the most powerful 19.967-GHz maser in the Galaxy in the entire history of observations, taking into account the 104-Jy flux from it on the Earth’s surface and the distance to the source, 6.75 kpc. The 19.967-GHz maser flared contemporaneously with the first of two flares detected in associated 20.971-GHz methanol masers. We estimated that the ratio of flux densities between these two transitions is F20.971/F19.967 = 14 ± 4, increasing to >520 in the second flare. We discuss the differences between the two flares in the 20.971-GHz methanol masers and the consequence thereof.

Fermi J1544-0649 is a transient GeV source first detected during its GeV flares in 2017. Multiwavelength observations during the flaring time demonstrate variability and spectral energy distributions that are typical of a blazar. Other than the flare time, Fermi J1544-0649 is quiet in the GeV band and has looked rather like a quiet galaxy (2MASX J15441967-0649156) for a decade. Together with the broad absorption-lines-like feature we further explore the “misaligned blazar scenario.” We analyzed the Very Long Baseline Array (VLBA) and East Asian VLBI Network (EAVN) data from 2018 to 2020 and discovered the four jet components from Fermi J1544-0649. We found a viewing angle around 3.°7 to 7.°4. The lower limit of the viewing angle indicates a blazar with an extreme low duty cycle of gamma-ray emission; the upper limit of it supports the “misaligned blazar scenario.” Follow-up multiwavelength observations after 2018 show Fermi J1544-0649 remains quiet in GeV, X-ray, and optical bands. A multimessenger search of neutrinos is also performed, and an excess of 3.1σ significance is found for this source.

We report the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93-0.03, which are flaring on similarly short timescales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded., 9 pages, 3 figures, 2 tables, Published in MNRAS

Class II methanol masers are signs of massive young stellar objects (MYSOs). Recent findings show that MYSO accretion bursts cause flares of these masers. Thus, maser monitoring can be used to identify such bursts. Burst-induced SED changes provide valuable information on a very intense phase of high-mass star formation. In mid-January 2019, a maser flare of the MYSO G358.93-0.03 was reported. ALMA and SMA imaging resolved the core of the star forming region and proved the association of the masers with the brightest continuum source MM1. However, no significant flux rise of the (sub)mm dust continuum was found. Thus, we performed NIR imaging with GROND and IFU spectroscopy with FIFI-LS aboard SOFIA to detect possible counterparts to the (sub)mm sources, and compare their photometry to archival measurements. The comparison of pre-burst and burst SEDs is of crucial importance to judge whether a luminosity increase due to the burst is present and if it triggered the maser flare. The FIR fluxes of MM1 measured with FIFI-LS exceed those from Herschel significantly, which clearly confirms the presence of an accretion burst. The second epoch data, taken about 16 months later, still show increased fluxes. Our RT modeling yielded major burst parameters and suggests that the MYSO features a circumstellar disk which might be transient. From the multi-epoch SEDs, conclusions on heating and cooling time-scales could be drawn. Circumstances of the burst-induced maser relocation have been explored. The verification of the accretion burst from G358 is another confirmation that Class II methanol maser flares represent an alert for such events. The few events known to date already indicate that there is a broad range in burst strength and duration as well as environmental characteristics. The G358 event is the shortest and least luminous MYSO accretion burst so far., Comment: 22 pages, 11 figures, accepted by A&A. Abstract abridged for arxiv submission

We carried out observations of the 22 GHz H2O masers in a high-mass protostar S255 NIRS 3 by using VERA. We measured the proper motions of the 22 GHz H2O masers associated with a bipolar outflow. The expansion velocity of the blueshifted bow shock traced by the 22 GHz H2O masers was 28 km s-1 corresponding to a dynamical timescale of 60 years. The direction of the maser outflow is slightly tilted compared with the radio jet, which could suggest a more recent ejection episode during the accretion burst event. The total flux density of the 22 GHz H2O masers has gradually increased from early 2017 and has become almost constant in 2018. For the first time, we reveal extended H2O maser emission at 22 GHz in a star forming region, which is partly resolved out by VERA and even by the most extended VLA configurations. We find that the flux variation of such an extended component is similar to that of the unresolved maser emission. We also conducted observations of the submillimeter continuum and the 321 GHz H2O masers with ALMA at Band 7. The continuum emission does not show significant variations compared with the previous observations performed 5 months before. We mapped the 321 GHz H2O masers in S255 NIRS 3 providing the fourth example, for this maser, of the spatial distribution in a high-mass star-forming region. The lower ratio of the 22 GHz/321 GHz maser luminosity in the blueshifted bow shock suggests a temperature (>1000 K), higher than for the other maser features in this region. We conclude that the bow shock structure traced by the 22 GHz H2O maser features is unlikely to be originating at the interface between the radio jet powered by the recent accretion outburst and the surrounding medium. The brightening of the 22 GHz H2O masers could be due to radiative excitation by photons form the infrared outburst escaping along the cavity created by the newly ejected material., 19 pages, 28 figures, accepted for publication in Astronomy & Astrophysics

Full list of authors: Cho, Ilje; Zhao, Guang-Yao; Kawashima, Tomohisa; Kino, Motoki; Akiyama, Kazunori; Johnson, Michael D.; Issaoun, Sara; Moriyama, Kotaro; Cheng, Xiaopeng; Algaba, Juan-Carlos; Jung, Taehyun; Sohn, Bong Won; Krichbaum, Thomas P.; Wielgus, Maciek; Hada, Kazuhiro; Lu, Ru-Sen; Cui, Yuzhu; Sawada-Satoh, Satoko; Shen, Zhiqiang; Park, Jongho; Jiang, Wu; Ro, Hyunwook; Yi, Kunwoo; Wajima, Kiyoaki; Lee, Jee Won; Hodgson, Jeffrey; Tazaki, Fumie; Honma, Mareki; Niinuma, Kotaro; Trippe, Sascha; An, Tao; Zhang, Yingkang; Lee, Jeong Ae; Oh, Se-Jin; Byun, Do-Young; Lee, Sang-Sung; Kim, Jae-Young; Oh, Junghwan; Koyama, Shoko; Asada, Keiichi; Wang, Xuezheng; Cui, Lang; Hagiwara, Yoshiaki; Nakamura, Masanori; Takamura, Mieko; Hirota, Tomoya; Sugiyama, Koichiro; Kawaguchi, Noriyuki; Kobayashi, Hideyuki; Oyama, Tomoaki; Yonekura, Yoshinori; Kim, Jongsoo; Hwang, Ju-Yeon; Jung, Dong-Kyu; Kim, Hyo-Ryoung; Kim, Jeong-Sook; Oh, Chung-Sik; Roh, Duk-Gyoo; Yeom, Jae-Hwan; Xia, Bo; Zhong, Weiye; Li, Bin; Zhao, Rongbing; Wang, Jinqing; Liu, Qinghui; Chen, Zhong.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited., Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets in which to resolve the innermost region of an SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multiwavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ∼704 ± 102 μas (axial ratio ∼${1.19}_{-0.19}^{+0.24}$) and ∼300 ± 25 μas (axial ratio ∼1.28 ± 0.2) at 1.349 cm and 6.95 mm, respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power law, with an index ∼1.2 ± 0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes. © 2022. The Author(s). Published by the American Astronomical Society., I.C. and J.P. are supported by the National Research Foundation of Korea (NRF) via a Global PhD Fellowship grant (NRF-2015H1A2A1033752 and NRF-2014H1A2A1018695, respectively). I.C. acknowledges financial support by the Spanish Ministerio de Economía y Competitividad (grants PID2019-108995GB-C21) and the Junta de Andalucía (grant P18-FR-1769). G.-Y.Z. and T.J. are supported by the Korea Research Fellowship Program through the NRF (NRF-2015H1D3A1066561). I.C. and G.-Y.Z. acknowledge financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). K.A. is supported by the National Science Foundation (NSF) through grants AST-1440254, AST-1614868, and AST-2034306. This work is partially supported by JSPS KAKENHI grant Nos. JP18K03656 (M.K.), JP18H03721 (K.N., M.K., K.H.), JP18K13594 (T.K.), JP18KK0090 (K.H.), JP19H01943 (K.H.), JP21H01137 (K.H., M.K.), and the Mitsubishi Foundation grant No. 201911019 (K.H.). J.-C.A. acknowledges support from the Malaysian Fundamental Research grant Scheme (FRGS) FRGS/1/2019/STG02 /UM/02/6. M.D.J. acknowledges support from the National Science Foundation (AST-1716536, AST-1935980) and the Gordon and Betty Moore Foundation (GBMF-5278). X.-P.C. and B.-W.S. are supported by Brain Pool Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (2019H1D3A1A01102564). R.-S.L. is supported by the Max Planck Partner Group of the MPG and the CAS and acknowledges the support by the Key Program of the National Natural Science Foundation of China (grant No. 11933007) and Research Program of Fundamental and Frontier Sciences, CAS (grant No. ZDBS-LY-SLH011). J.P. acknowledges financial support through the EACOA Fellowship awarded by the East Asia Core Observatories Association, which consists of the Academia Sinica Institute of Astronomy and Astrophysics, the National Astronomical Observatory of Japan, Center for Astronomical Mega-Science, Chinese Academy of Sciences, and the Korea Astronomy and Space Science Institute. S.T. acknowledges financial support from the NRF via Basic Research grant 2019R1F1A1059721. This work was also supported in part by MEXT SPIRE, MEXT as "Priority Issue on post-K computer" (Elucidation of the Fundamental Laws and Evolution of the Universe) and as "Program for Promoting Researches on the Supercomputer Fugaku" (High-energy astrophysical phenomena in black holes and supernovae). J.S.K. (grant Nos. 2016R1A5A1013277 and 2020R1A2C1007219) and J.O. (grant No. 2021R1A6A3A01086420) has been supported by the Basic Science Research Program through the NRF funded by the Ministry of Education.

We report results of 0.05"-resolution observations toward the O-type proto-binary system IRAS 16547-4247 with the Atacama Large Millimeter/submillimeter Array (ALMA). We present dynamical and chemical structures of the circumbinary disk, circumstellar disks, outflows and jets, illustrated by multi-wavelength continuum and various molecular lines. In particular, we detect sodium chloride, silicon compounds, and vibrationally-excited water lines as probes of the individual protostellar disks at a scale of 100 au. These are complementary to typical hot-core molecules tracing the circumbinary structures on a 1000-au scale. The H2O line tracing inner-disks has an upper-state energy of Eu/k>3000K, indicating a high temperature of the disks. On the other hand, despite the detected transitions of NaCl, SiO, and SiS not necessarily having high upper-state energies, they are enhanced only in the vicinity of the protostars. We interpret that these molecules are the products of dust destruction, which only happens in the inner disks. This is the second detection of alkali metal halide in protostellar systems after the case of the disk of Orion Source I, and also one of few massive protostellar disks associated with high-energy transition water and silicon compounds. These new results suggest these "hot-disk" lines may be common in innermost disks around massive protostars, and have great potential for future research of massive star formation. We also tentatively find that the twin disks are counter-rotating, which might give a hint of the origin of the massive proto-binary system IRAS 16547-4247., 15 pages, 5 figures, 2 appendix figures. Published in ApJ Letters

We present results of continuum and spectral line observations with ALMA and 22 GHz water (H$_2$O) maser observations using KaVA and VERA toward a high-mass star-forming region, G25.82-0.17. Multiple 1.3 mm continuum sources are revealed, indicating the presence of young stellar objects (YSOs) at different evolutionary stages, namely an ultra-compact HII region, G25.82-E, a high-mass young stellar object (HM-YSO), G25.82-W1, and starless cores, G25.82-W2 and G25.82-W3. Two SiO outflows, at N-S and SE-NW orientations, are identified. The CH$_3$OH 8$_{-1}$-7$_{0}$ E line, known to be a class I CH$_3$OH maser at 229 GHz is also detected showing a mixture of thermal and maser emission. Moreover, the H$_2$O masers are distributed in a region ~0.25" shifted from G25.82-W1. The CH$_3$OH 22$_{4}$-21$_{5}$ E line shows a compact ring-like structure at the position of G25.82-W1 with a velocity gradient, indicating a rotating disk or envelope. Assuming Keplerian rotation, the dynamical mass of G25.82-W1 is estimated to be $>$25 M$_{\odot}$ and the total mass of 20 M$_\odot$-84 M$_\odot$ is derived from the 1.3 mm continuum emission. The driving source of the N-S SiO outflow is G25.82-W1 while that of the SE-NW SiO outflow is uncertain. Detection of multiple high-mass starless$/$protostellar cores and candidates without low-mass cores implies that HM-YSOs could form in individual high-mass cores as predicted by the turbulent core accretion model. If this is the case, the high-mass star formation process in G25.82 would be consistent with a scaled-up version of low-mass star formation., 28 pages, 9 figures, accepted to ApJ

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively., Comment: 35 pages, 7 figures, accepted for publication in the Publications of the Astronomical Society of Japan

We followed up the massive young stellar object S255-NIRS3 (= S255-IRS1b) during its recent accretion outburst event in the $K_{\rm s}$ band with Kanata/HONIR for four years after its burst and obtained a long-term light curve. This is the most complete near-infrared light curve of the S255-NIRS3 burst event that has ever been presented. The light curve showed a steep increase reaching a peak flux that was 3.4 mag brighter than the quiescent phase and then a relatively moderate year-scale fading until the last observation, similar to that of the accretion burst events such as EXors found in lower-mass young stellar objects. The behavior of the $K_{\rm s}$-band light curve is similar to that observed in 6.7 GHz class II methanol maser emission, with a sudden increase followed by moderate year-scale fading. However, the maser emission peaks appear 30–50 d earlier than that of the $K_{\rm s}$ band emission. The similarities confirmed that the origins of the maser emission and the $K_{\rm s}$-band continuum emission are common, as previously shown from other infrared and radio observations by Stecklum et al. (2016, Astronomer’s Telegram, 8732), Caratti o Garatti et al. (2017b, Nature Phys., 13, 276), and Moscadelli et al. (2017, A&A, 600, L8). However, the differences in energy transfer paths, such as the exciting/emitting/scattering structures, may cause the delay in the flux-peak dates.

In this proceeding paper, we introduce the recent results of Galactic maser astrometry by mainly focusing on those obtained with Japanese VLBI array VERA. So far we have obtained parallaxes for 86 sources including preliminary results, and combination with the data obtained with VLBA/BeSSeL provides astrometric results for 159 sources. With these most updated results we conduct preliminary determinations of Galactic fundamental parameters, obtaining R0 = 8.16 ± 0.26 kpc and Θ0 = 237 ± 8 km/s. We also derive the rotation curve of the Milky Way Galaxy and confirm the previous results that the rotation curve is fairly flat between 5 kpc and 16 kpc, while a remarkable deviation is seen toward the Galactic center region. In addition to the results on the Galactic structure, we also present brief overviews on other science topics related to masers conducted with VERA, and also discuss the future prospect of the project.

We report on interferometric observations of a face-on accretion system around the high mass young stellar object, G353.273+0.641. The innermost accretion system of 100-au radius was resolved in a 45-GHz continuum image taken with the Jansky Very Large Array. Our SED analysis indicated that the continuum could be explained by optically-thick dust emission. 6.7 GHz CH3OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution thermal continuum and line observations with ALMA long baselines.

Despite their importance in the formation and evolution of stellar clusters and galaxies, the formation of high-mass stars remains poorly understood. We recently started a systematic observational study of the 22 GHz water and 44 GHz class I methanol masers in high-mass star-forming regions as a four-year KaVA large program. Our sample consists of 87 high-mass young stellar objects (HM-YSOs) in various evolutionary phases, many of which are associated with two or more different maser species. The primary scientific goals are to measure the spatial distributions and 3-dimensional velocity fields of multiple maser species, and understand the dynamical evolution of HM-YSOs and their circumstellar structures, in conjunction with follow-up observations with JVN/EAVN (6.7 GHz class II methanol masers), VERA, and ALMA. In this paper we present details of our KaVA large program, including the first-year results and observing/data analysis plans for the second year and beyond.

An observational result of a radio continuum and H92α radio recombination line of the Galactic center lobe (GCL), using the Yamaguchi 32 m radio telescope, is reported. The obtained spatial intensity distribution of the radio recombination line shows two distinctive ridge-like structures extending from the Galactic plane vertically to the north at the eastern and western sides of the Galactic center, which are connected to each other at a latitude of ${1{^{\circ}_{.}}2}$ to form a loop-like structure as a whole. This suggests that most of the radio continuum emission of the GCL is free–free emission, and that the GCL is filled with thermal plasma. The east ridge of the GCL observed with the radio recombination line separates 30 pc from the radio arc, which has been considered a part of the GCL, but coincides with a ridge of the radio continuum at a Galactic longitude of 0°. The radial velocity of the radio recombination line is found to be between −4 and +10 km s−1 across the GCL. This velocity is much smaller than expected from the Galactic rotation, and hence indicates that the GCL exists apart from the Galactic center. These characteristics of the GCL suggest that the long-standing hypothesis that the GCL was created by explosive activity in the Galactic center is unlikely, but favor that the GCL is a giant H ii region.

We report the first detection of isotopic methanol (13CH3OH) maser emission in interstellar space. The emission was detected toward the high-mass young stellar object G358.93-0.03 during monitoring of a flare in the 6.7 GHz methanol (CH3OH) maser emission in this source. We find that the spectral and spatial distribution of the 13CH3OH masers differs from the CH3OH masers imaged at the same epoch, contrary to expectations from similarity of their pumping. This conclusively demonstrates that isotopic methanol masers are bright under different physical conditions and suggests that they can provide additional, complementary information to the CH3OH masers from the same source. We detect a rapid decay of the 13CH3OH maser lines suggesting that they are transient phenomena (masing for only a few months), likely associated with rapid changes in radiation field due to an accretion burst induced by massive disk fragmentation. Changes in the line flux density are faster than required to achieve equilibrium in the energy level populations, indicating that the pumping of these masers is likely variable. © 2020. The American Astronomical Society. All rights reserved..

We initiated a long-term and highly frequent monitoring project toward 442 methanol masers at 6.7 GHz (Dec >−30 deg) using the Hitachi 32-m radio telescope in December 2012. The observations have been carried out daily, monitoring a spectrum of each source with intervals of 9–10 days. In September 2015, the number of the target sources and intervals were redesigned into 143 and 4–5 days, respectively. This monitoring provides us complete information on how many sources show periodic flux variations in high-mass star-forming regions, which have been detected in 20 sources with periods of 29.5–668 days so far (e.g., Goedhart et al. 2004). We have already obtained new detections of periodic flux variations in 31 methanol sources with periods of 22–409 days. These periodic flux variations must be a unique tool to investigate high-mass protostars themselves and their circumstellar structure on a very tiny spatial scale of 0.1–1 au.

We report on interferometric observations of a face-on accretion system around the High-Mass young stellar object, G353.273+0.641. The innermost accretion system of 100 au radius was resolved in a 45 GHz continuum image taken with the Jansky-Very Large Array. Our spectral energy distribution analysis indicated that the continuum could be explained by optically thick dust emission. The total mass of the dusty system is $\sim$ 0.2 $M_{\sun}$ at minimum and up to a few $M_{\sun}$ depending on the dust parameters. 6.7 GHz CH$_{3}$OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. This is clearly smaller than that of typical accretion disks in High-Mass star formation, indicating that the initial angular momentum was very small, or the CH$_{3}$OH masers selectively trace accreting material that has small angular momentum. In the former case, the initial specific angular momentum is estimated to be 8 $\times$ 10$^{20}$ ($M_{*}$$/$10 $M_{\sun}$)$^{0.5}$ cm$^{2}$ s$^{-1}$, or a significant fraction of the initial angular momentum was removed outside of 100 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution ($\sim$ 10 mas) thermal continuum and line observations with ALMA long baselines., 17 pages, 11 figures

We detected flaring flux variability that regularly occurred with the period of 23.9 days on a 6.7 GHz methanol maser emission at Vlsr = 25.30 km/s in G 014.23-00.50 through highly frequent monitoring using the Hitachi 32-m radio telescope. By analyzing data from 05 January 2013 to 21 January 2016, the periodic variability has persisted in at least 47 cycles, corresponding to approximately 1,100 days. The period of 23.9 days is the shortest one observed in masers at around high-mass young stellar objects so far. The flaring component normally falls below the detection limit (3 sigma) of 0.9 Jy. In the flaring periods, the component rises above the detection limit with the ratio of the peak flux density more than 180 in comparison with a quiescent phase, showing intermittent periodic variability. The time-scale of the flux rise was typically two days or shorter, and both symmetric and asymmetric profiles of flux variability were observed through intraday monitoring. These characteristics might be explained by a change in the flux of seed photons by a colliding-wind binary (CWB) system or a variation of the dust temperature by an extra heating source of a shock formed by the CWB system within a gap region in a circumbinary disk, in which the orbital semi-major axes of the binary are 0.26-0.34 au., 12 pages, 9 figures, 2 tables, accepted for publication in PASJ

For the sake of high-sensitivity 6.7 GHz methanol maser observations, we developed a new technology for coherently combining the two signals from the Hitachi 32 m radio telescope and the Takahagi 32 m radio telescope of the Japanese Very long baseline interferometer Network (JVN), where the two telescopes were separated by about 260 m. After the two telescopes were phased as a twofold larger single telescope, the mean signal-to-noise ratio (SNR) of the 6.7 GHz methanol masers observed by the phased telescopes was improved to 1.254-fold higher than that of the single dish, through a Very Long Baseline Interferometry (VLBI) experiment on the 50 km baseline of the Kashima 34 m telescope and the 1000 km baseline of the Yamaguchi 32 m telescope. Furthermore, we compared the SNRs of the 6.7 GHz maser spectra for two methods. One is a VLBI method and the other is the newly developed digital position switching, which is a similar technology to that used in noise-cancelling headphones. Finally, we confirmed that the mean SNR of method of the digital position switching (ON-OFF) was 1.597-fold higher than that of the VLBI method., Comment: 17 pages, 10 figures, accepted for publication in Publications of the Astronomical Society of the Pacific (PASP)

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements., 10 pages, 2 figures, 3 tables, accepted for publication in The Astrophysical Journal Letters. All maser transitions detected in this object can be displayed via an option on the advanced interface of http://splatalogue.net

We report the first detection of a radio-continuum and molecular jet associated with a dominant blue-shifted maser source, G353.273+0.641. A radio jet is extended 3000 au along NW-SE direction. H$_{2}$O masers are found to be clustered in the root of a bipolar radio jet. A molecular jet is detected by thermal SiO ($\upsilon$ = 0, $J$ = 2-1) emission. The SiO spectrum is extremely wide (-120 -- +87 km s$^{-1}$) and significantly blue-shift dominated, similar to the maser emission. The observed geometry and remarkable spectral similarity between H$_{2}$O maser and SiO strongly suggests the existence of a maser-scale ($\sim$ 340 au) molecular jet that is enclosed by the extended radio jet. We propose a "disc-masking" scenario as the origin of the strong blue-shift dominance, where an optically thick disc obscures a red-shifted lobe of a compact jet., 5 pages, 3 figures, 2 tables, accepted for publication in MNRAS

The water maser site associated with G353.273+0.641 is classified as a dominant blueshifted H2O maser, which shows an extremely wide velocity range (± 100 km s−1) with almost all flux concentrated in the highly blueshifted emission. The previous study has proposed that this peculiar H2O maser site is excited by a pole-on jet from high mass protostellar object. We report on the monitoring of 22-GHz H2O maser emission from G353.273+0.641 with the VLBI Exploration of Radio Astrometry (VERA) and the Tomakamai 11-m radio telescope. Our VLBI imaging has shown that all maser features are distributed within a very small area of 200 × 200 au2, in spite of the wide velocity range (> 100 km s−1). The light curve obtained by weekly single-dish monitoring shows notably intermittent variation. We have detected three maser flares during three years. Frequent VLBI monitoring has revealed that these flare activities have been accompanied by a significant change of the maser alignments. We have also detected synchronized linear acceleration (−5 km s−1yr−1) of two isolated velocity components, suggesting a lower-limit momentum rate of 10−3 M⊙ km s−1yr−1 for the maser acceleration. All our results support the previously proposed pole-on jet scenario, and finally, a radio jet itself has been detected in our follow-up ATCA observation. If highly intermittent maser flares directly reflect episodic jet-launchings, G353.273+0.641 and similar dominant blueshifted water maser sources can be suitable targets for a time-resolved study of high mass protostellar jet.

We have started survey observations of the 22 GHz water maser sources associated with high-mass young stellar objects (HM-YSOs) as a part of the KaVA (KVN and VERA Array) large program (LP). The aim of our LP is to understand dynamical evolution of jets/outflows from HM-YSOs by analyzing 3D velocity structures of water maser features. In the first year (2016-2017), an imaging survey toward 25 HM-YSOs has been conducted and the 22 GHz water masers are detected toward 21 sources. Spatial distributions of maser features for individual sources are mapped. To complement physical properties in the vicinity of HM-YSOs, we have carried out ALMA cycle 3 observations of thermal molecular lines and continuum emissions toward 11 selected samples. Summary of the KaVA first year observations and the initial results from the ALMA toward one of our targets, G25.82-0.17, are reported.

Accepted: 2010-12-16, 資料番号: SA1002907000

We have observed 13 methanol maser sources associated with massive star-forming regions; W3(OH), Mon R2, S 255, W 33A, IRAS 18151-1208, G 24.78+0.08, G 29.95-0.02, IRAS 18556+0136, W 48, OH 43.8-0.1, ON 1, Cep A and NGC 7538 at 6.7 GHz using the Japanese VLBI Network (JVN). Twelve of the thirteen sources were detected at our longest baseline of $\sim$50 M$\lambda$, and their images are presented. Seven of them are the first VLBI images at 6.7 GHz. This high detection rate and the small fringe spacing of $\sim$4 milli-arcsecond suggest that most of the methanol maser sources have compact structure. Given this compactness as well as the known properties of long-life and small internal-motion, this methanol maser line is suitable for astrometry with VLBI., Comment: 15 pages, 13 figures, accepted for publication in PASJ

We have initiated single-dish monitoring observations of ∼400 methanol maser sources at 6.7 GHz using the Hitachi 32-m radio telescope from December 2012 to systematically research periodic flux variations, which are observed in some methanol maser sources associated with high-mass (proto-)stars. In our monitoring, we have made daily monitoring, so that each source has been observed every nine days with an integration time of 5 min (typical 3 σ detection sensitivities of 0.9 Jy). The monitoring observations help us statistically understand periodic flux variations with a period longer than 50 days. As an initial result, we present a new detection of periodic flux variations in the 6.7 GHz methanol maser source G 036.70+00.09. The period of the flux variations is ∼53 days (∼0.019 cycles day−1), and seems to be stable over 9 cycles, at least until the middle of August 2014.

We report on a long-term single-dish and VLBI monitoring for intermittent flare activities of a Dominant Blue-Shifted H$_{2}$O Maser (DBSM) associated with a southern high mass young stellar object, G353.273+0.641. Bi-weekly single-dish monitoring using Hokkaido University Tomakomai 11-m radio telescope has shown that a systematic acceleration continues over four years beyond a lifetime of individual maser features. This fact suggests that the H$_{2}$O maser traces a region where molecular gas is steadily accelerated. There were five maser flares during five-years monitoring, and maser distributions in four of them were densely monitored by the VLBI Exploration of Radio Astrometry (VERA). The overall distribution of the maser features suggests the presence of a bipolar jet, with the 3D kinematics indicating that it is almost face-on (inclination angle of $\sim$ 8$^{\fdg}$--17$^{\fdg}$ from the line-of-sight). Most of maser features were recurrently excited within a region of 100$\times$100 AU$^{2}$ around the radio continuum peak, while their spatial distributions significantly varied between each flare. This confirms that episodic propagations of outflow shocks recurrently invoke intermittent flare activities. We also measured annual parallax, deriving the source distance of 1.70 $^{+0.19}_{-0.16}$ kpc that is consistent with the commonly-used photometric distance., Accepted to PASJ, 17 pages, 12 figure

We have observed bursting variability of the 6.7 GHz methanol maser of G33.641-0.228. Five bursts were detected in the observation period of 294 days from 2009 to 2012. The typical burst is a large flux density rise in about one day followed by a slow fall. A non-typical burst observed in 2010 showed a large and rapid flux density enhancement from the stable state, but the rise and fall of the flux density were temporally symmetric and a fast fluctuation continued 12 days. On average, the bursts occurred once every 59 days, although bursting was not periodic. Since the average power required for causing the burst of order of 10^21 Js^-1 is far smaller than the luminosity of G33.641-0.228, a very small fraction of the source's power would be sufficient to cause the burst occasionally. The burst can be explained as a solar-flare like event in which the energy is accumulated in the magnetic field of the circumstellar disk, and is released for a short time. However, the mechanism of the energy release and the dust heating process are still unknown., Comment: 11 pages, 4 figures, accepted to PASJ

We present the first internal motion measurement of the 6.7-GHz methanol maser within S269, a small HII region in the outer Galaxy, which was carried out in 2006 and 2011 using the Japanese VLBI Network (JVN). Several maser groups and weak isolated spots were detected in an area spanning by ~200 mas (1000 AU). Three remarkable maser groups are aligned at a position angle of 80 degree. Two of three maser groups were also detected by a previous observation in 1998, which allowed us to study a long-term position variation of maser spots from 1998 to 2011. The angular separation between the two groups increased ~10 mas, which corresponds to an expansion velocity of ~10 km s^{-1}. Some velocity gradient (~10^{-2} km s^{-1} mas^{-1}) in the overall distribution was found. The internal motion between the maser groups support the hypothesis that the methanol masers in S269 could trace a bipolar outflow., 8 pages, 6 figures, accepted for publication in PASJ

The 6.7 GHz methanol maser emission, a tracer of forming massive stars, sometimes shows enigmatic periodic flux variations over several 10-100 days. In this Letter, we propose that this periodic variations could be explained by the pulsation of massive protostars growing under rapid mass accretion with rates of Mdot > 10^-3 Msun/yr. Our stellar evolution calculations predict that the massive protostars have very large radius exceeding 100 Rsun at maximum, and we here study the pulsational stability of such the bloated protostars by way of the linear stability analysis. We show that the protostar becomes pulsationally unstable with various periods of several 10-100 days, depending on different accretion rates. With the fact that the stellar luminosity when the star is pulsationally unstable also depends on the accretion rate, we derive the period-luminosity relation log (L/Lsun) = 4.62 + 0.98log(P/100 day), which is testable with future observations. Our models further show that the radius and mass of the pulsating massive protostar should also depend on the period. It would be possible to infer such protostellar properties and the accretion rate with the observed period. Measuring the maser periods enables a direct diagnosis of the structure of accreting massive protostars, which are deeply embedded in dense gas and inaccessible with other observations., 5 pages, 3 figures, 1 table, accepted for publication in ApJL

We present the results of daily monitoring of 6.7 GHz methanol maser in Cepheus A (Cep A) using Yamaguchi 32-m radio telescope as well as the results of imaging observations conducted with the JVN (Japanese VLBI Network). We indentified five spectral features, which are grouped into red-shifted (-1.9 and -2.7 km/s) and blue-shifted (-3.8, -4.2, and -4.9 km/s), and we detected rapid variabilities of these maser features within a monitoring period of 81 days. The red-shifted features decreased in flux density to 50% of its initial value, while the flux density of the blue-shifted features rapidly increased within a 30 days. The time variation of these maser features showed two remarkable properties; synchronization and anti-correlation between the red-shifted and the blue-shifted. The spatial distribution of the maser spots obtained by the JVN observation showed an arclike structure with a scale of $\sim$1400 AU, and separations of the five maser features were found to be larger than 100 AU. The absolute position of the methanol maser was also obtained based on the phase-referencing observations, and the arclike structure were found to be associated with the Cep A-HW2 object, with the elongation of the arclike structure nearly perpendicularly to the radio continuum jet from the Cep A-HW2 object. These properties of the masers, namely, the synchronization of flux variation, and the spectral and spatial isolation of features, suggest that the collisional excitation by shock wave from a common exciting source is unlikely. Instead, the synchronized time variation of the masers can be explained if all the maser features are excited by infrared radiation from dust which is heated by a common exciting source with a rapid variability., 7 pages, 3 figures, 2 tables, accepted for publication in PASJ

We present the basic features and the activities of Japanese VLBI network (JVN), a newly-established VLBI network with baselines ranging from 50 to 2560 km spreading across the Japanese islands, and capable of observing at 6.7, 8.4, and 22 GHz. We show a number of results of JVN observations: 8.4-GHz continuum images of a Giga-hertz Peaked Spectrum (GPS) source and radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s), the spatial and velocity structures of water masers in NML Cygni as well as methanol masers in Cep A, and demonstrative observations with the bigradient phase referencing., Comment: 3 pages, 2 figures, to appear in the proceedings of the 8th EVN symposium held in Torun, Poland, September 2006