Search

We report on radio follow-up observations of the nearby Type II supernova, SN 2023ixf, spanning from 1.7 to 269.9 days after the explosion, conducted using three very long baseline interferometers (VLBIs), which are the Japanese VLBI Network (JVN), the VLBI Exploration of Radio Astrometry (VERA), and the Korean VLBI Network (KVN). In three observation epochs (152.3, 206.1, and 269.9 days), we detected emission at the 6.9 and 8.4 GHz bands, with a flux density of $\sim 5$ mJy. The flux density reached a peak at around 206.1 days, which is longer than the timescale to reach the peak observed in typical Type II supernovae. Based on the analytical model of radio emission, our late-time detections were inferred to be due to the decreasing optical depth. In this case, the mass-loss rate of the progenitor is estimated to have increased from $\sim 10^{-6} - 10^{-5}\, M_{\odot}\,{\rm yr^{-1}}$ to $\sim 10^{-4}\, M_{\odot}\,{\rm yr^{-1}}$ between 28 and 6 years before the explosion. Our radio constraints are also consistent with the mass-loss rate to produce a confined circumstellar medium proposed by previous studies, which suggest that the mass-loss rate increased from $\sim 10^{-4}\, M_{\odot}\,{\rm yr^{-1}}$ to $\gtrsim 10^{-2}\, M_{\odot}\,{\rm yr^{-1}}$ in the last few years before the explosion., Comment: 12 pages, 3 figures, 3 tables. Accepted for publication in ApJ

Black Hole Explorer (BHEX) is a space VLBI mission concept, which can probe the black hole spacetime and the plasma properties including the magnetic fields of the accretion flows and relativistic jets. We propose science cases anticipated to be addressed by BHEX mainly via the imaging of the target objects, whose observational features appear in several microarcsecond scale. An appearance of a crescent-shaped shadow in a bright state of the M87 will be able to constrain the magnitude of the black hole spin. A possible appearance of the plasma injection region in the vicinity of the black hole results in the formation of the multiple ring structure and may enable us to understand the jet formation processes. In addition, The separation of linear and circular polarization fluxes and reversal of circular polarization will constrain the magnetic field structure and the thermal properties of the electrons, respectively. Other topics including the test of the gravitational theory are also being discussed., Comment: 14 pages, 4figures, Proceedings of SPIE 2024

We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery will expose universal features of a black hole's spacetime that are distinct from the complex astrophysics of the emitting plasma, allowing the first direct measurements of a supermassive black hole's spin. In addition to studying the properties of the nearby supermassive black holes M87* and Sgr A*, BHEX will measure the properties of dozens of additional supermassive black holes, providing crucial insights into the processes that drive their creation and growth. BHEX will also connect these supermassive black holes to their relativistic jets, elucidating the power source for the brightest and most efficient engines in the universe. BHEX will address fundamental open questions in the physics and astrophysics of black holes that cannot be answered without submillimeter space VLBI. The mission is enabled by recent technological breakthroughs, including the development of ultra-high-speed downlink using laser communications, and it leverages billions of dollars of existing ground infrastructure. We present the motivation for BHEX, its science goals and associated requirements, and the pathway to launch within the next decade., Comment: Proceedings for SPIE Astronomical Telescopes and Instrumentation

The Black Hole Explorer (BHEX) is a next-generation space very long baseline interferometry (VLBI) mission concept that will extend the ground-based millimeter/submillimeter arrays into space. The mission, closely aligned with the science priorities of the Japanese VLBI community, involves an active engagement of this community in the development of the mission, resulting in the formation of the Black Hole Explorer Japan Consortium. Here we present the current Japanese vision for the mission, ranging from scientific objectives to instrumentation. The Consortium anticipates a wide range of scientific investigations, from diverse black hole physics and astrophysics studied through the primary VLBI mode, to the molecular universe explored via a potential single-dish observation mode in the previously unexplored 50-70\,GHz band that would make BHEX the highest-sensitivity explorer ever of molecular oxygen. A potential major contribution for the onboard instrument involves supplying essential elements for its high-sensitivity dual-band receiving system, which includes a broadband 300\,GHz SIS mixer and a space-certified multi-stage 4.5K cryocooler akin to those used in the Hitomi and XRISM satellites by the Japan Aerospace Exploration Agency. Additionally, the Consortium explores enhancing and supporting BHEX operations through the use of millimeter/submillimeter facilities developed by the National Astronomical Observatory of Japan, coupled with a network of laser communication stations operated by the National Institute of Information and Communication Technology., Comment: Submitted to SPIE Conference Series, 17 pages for the main text, 33 pages for the references, 4 figures, SPIE Astronomical Telescopes + Instrumentation 2024

The nearby radio galaxy M87 offers a unique opportunity to explore the connections between the central supermassive black hole and relativistic jets. Previous studies of the inner region of M87 revealed a wide opening angle for the jet originating near the black hole. The Event Horizon Telescope resolved the central radio source and found an asymmetric ring structure consistent with expectations from General Relativity. With a baseline of 17 years of observations, there was a shift in the jet's transverse position, possibly arising from an eight to ten-year quasi-periodicity. However, the origin of this sideways shift remains unclear. Here we report an analysis of radio observations over 22 years that suggests a period of about 11 years in the position angle variation of the jet. We infer that we are seeing a spinning black hole that induces the Lense-Thirring precession of a misaligned accretion disk. Similar jet precession may commonly occur in other active galactic nuclei but has been challenging to detect owing to the small magnitude and long period of the variation., Comment: 41 pages, 7 figures, 7 tables

We report dual-frequency very long baseline interferometry (VLBI) observations at 22 and 43 GHz toward the nucleus of a nearby radio galaxy NGC 4261. In particular, we present a VLBI image of the 22 GHz H2O maser line and its location in the circumnuclear region of NGC 4261. H2O maser emission is marginally detected above the three times the rms level within a velocity range of approximately 2250-2450 km/s, slightly red-shifted with respect to the systemic velocity. H2O maser emission is located approximately 1 milliarcsecond (mas) east of the brightest continuum component at 22 GHz, where the continuum spectrum is optically thick, that is at the free-free absorbed receding jet by ionized gas. A positional coincidence between H2O maser emission and an ionized gas disk implies that the H2O maser emission arises from the near side of the disk, amplifying continuum emission from the background receding jet. If the disk axis is oriented 64 degree relative to the line of sight, the H2O maser emission is expected to be at a mean radius of 0.3 pc in the disk. The broad line width of the H2O maser emission can be attributed to complex kinematics in the immediate vicinity of the supermassive black hole (SMBH), including ongoing gas infall onto the SMBH, turbulence, and outflow. This is analogous to the multi-phase circumnuclear torus model in the nearest radio-loud H2O megamaser source NGC 1052. An alternative explanation for H2O maser association is the shock region between the jet and the ambient molecular clouds. However, this explanation fails to describe the explicit association of H2O maser emission only with the free-free absorbed receding jet., Comment: 10 pages, 5 figures, submitted to PASJ

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

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

We present a detection of a bright burst from FRB 20201124A, which is one of the most active repeating FRBs, based on S-band observations with the 64-m radio telescope at the Usuda Deep Space Center/JAXA. This is the first FRB observed by using a Japanese facility. Our detection at 2 GHz in February 2022 is the highest frequency for this FRB and the fluence of $>$ 189 Jy ms is one of the brightest bursts from this FRB source. We place an upper limit on the spectral index $\alpha$ = -2.14 from the detection of the S band and non-detection of the X band at the same time. We compare an event rate of the detected burst with ones of the previous research, and suggest that the power-law of the luminosity function might be broken at lower fluence, and the fluences of bright FRBs distribute up to over 2 GHz with the power-law against frequency. In addition, we show the energy density of the burst detected in this work was comparable to the bright population of one-off FRBs. We propose that repeating FRBs can be as bright as one-off FRBs, and only their brightest bursts could be detected so some of repeating FRBs intrinsically might have been classified as one-off FRBs., Comment: 8 pages, 5 figures, accepted for publication in Publications of the Astronomical Society of Japan (PASJ)

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

We report multi-epoch Very Long Baseline Interferometric (VLBI) observations of the compact radio lobe in the radio galaxy 3C84 (NGC1275) during 2016 - 2020. The image sequence of 3C84 reveals that the hotspot in the radio lobe showed the one-year long frustration in 2017 within a compact region of about 0.07 parsec, suggesting a strong collision between the jet and a compact dense cloud with the estimated average density about $(4\-- 6)\times 10^{5}~{\rm cm^{-3}}$. Although the hotspot and the radio lobe began to move south again after its breakout, the radio lobe showed a morphological transition from FR II- to FR I-class radio lobe and its radio flux became fainter. This is the first detection of the dynamical feedback from the cloud to the jet where the cloud located on the jet axis significantly interferes with the jet propagation and evolution at the central one-parsec region in active galactic nucleus., Comment: The Astrophysical Journal Letters, Volume 920, Issue 1, id.L24, 7 pp

We present the results of 3 GHz radio continuum observations of 23 superluminous supernovae (SLSNe) and their host galaxies by using the Karl G. Jansky Very Large Array conducted 5-21 years after the explosions. The sample consists of 15 Type I and 8 Type II SLSNe at z < 0.3, providing one of the largest sample of SLSNe with late-time radio data. We detected radio emission from one SLSN (PTF10hgi) and 5 hosts with a significance of >5$\sigma$. No time variability is found in late-time radio light curves of the radio-detected sources in a timescale of years except for PTF10hgi, whose variability is reported in a separate study. Comparison of star-formation rates (SFRs) derived from the 3 GHz flux densities with those derived from SED modeling based on UV-NIR data shows that four hosts have an excess of radio SFRs, suggesting obscured star formation. Upper limits for undetected hosts and stacked results show that the majority of the SLSN hosts do not have a significant obscured star formation. By using the 3 GHz upper limits, we constrain the parameters for afterglows arising from interaction between initially off-axis jets and circumstellar medium (CSM). We found that the models with higher energies ($E_{\rm iso} \gtrsim$ several $\times 10^{53}$ erg) and CSM densities ($n \gtrsim 0.01$ cm$^{-3}$) are excluded, but lower energies or CSM densities are not excluded with the current data. We also constrained the models of pulsar wind nebulae powered by a newly born magnetar for a subsample of SLSNe with model predictions in the literature., Comment: Accepted for publication in ApJ, 16 pages, 7 figures, 6 tables

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

We report the time variability of the late-time radio emission in a Type-I superluminous supernova (SLSN), PTF10hgi, at z = 0.0987. The Karl G. Jansky Very Large Array 3 GHz observations at 8.6 and 10 years after the explosion both detected radio emission with a ~40% decrease in flux density in the second epoch. This is the first report of a significant variability of the late-time radio light curve in a SLSN. Through combination with previous measurements in two other epochs, we constrained both the rise and decay phases of the radio light curve over three years, peaking at approximately 8-9 years after the explosion with a peak luminosity of L(3GHz) = 2 x 10^21 W/Hz. Possible scenarios for the origin of the variability are an active galactic nucleus (AGN) in the host galaxy, an afterglow caused by the interaction between an off-axis jet and circumstellar medium, and a wind nebula powered by a newly-born magnetar. Comparisons with models show that the radio light curve can be reproduced by both the afterglow model and magnetar wind nebula model. Considering the flat radio spectrum at 1-15 GHz and an upper limit at 0.6 GHz obtained in previous studies, plausible scenarios are a low-luminosity flat-spectrum AGN or a magnetar wind nebula with a shallow injection spectral index., Comment: 9 pages, 4 figures, accepted for publication in ApJ Letters

We observed the gravitationally lensed blazar JVAS B0218+357 with the KVN and VERA Array (KaVA) at 22, 43, and 86 GHz. The source has recently been identified as an active gamma-ray source up to GeV/TeV energy bands, rendering a unique target for studying relativistic jets through gravitational lensing. Here we report the first robust VLBI detection and imaging of the lensed images up to 86 GHz. The detected mas-scale/parsec-scale morphology of the individual lensed images (A and B) is consistent with that previously seen at 22 and 15 GHz, showing the core-jet morphology with the jet direction being the same as at the low frequencies. The radio spectral energy distributions of the lensed images become steeper at higher frequencies, indicating that the innermost jet regions become optically thin to synchrotron emission. Our findings confirm that the absorption effects due to the intervening lensing galaxy become negligible at millimeter wavelengths. These results indicate that high-frequency VLBI observations are a powerful tool to better recover the intrinsic properties of lensed active galactic nucleus jets, which therefore allow us to study the interplay between the low- and high-energy emission., Comment: 11pages, 7 pages, 4 tables, accepted for publication in ApJ

Waseda University Nasu telescope array is a spatial fast Fourier transform (FFT) interferometer consisting of eight linearly aligned antennas with 20 m spherical dishes. This type of interferometer was developed to survey transient radio sources with an angular resolution as high as that of a 160 m dish with a field of view as wide as that of a 20 m dish. We have been performing drift-scan-mode observations, in which the telescope scans the sky around a selected declination as the earth rotates. The black hole X-ray binary V404 Cygni underwent a new outburst in 2015 June after a quiescent period of 26 years. Because of the interest in black hole binaries, a considerable amount of data on this outburst at all wavelengths was accumulated. Using the above telescope, we had been monitoring V404 Cygni daily from one month before the X-ray outburst, and two radio flares at 1.4 GHz were detected on June 21.73 and June 26.71. The flux density and time-scale of each flare were 313+/-30 mJy and 1.50+/-0.49 days, 364+/-30 mJy and 1.70+/-0.16 days, respectively. We have also confirmed the extreme variation of radio spectra within a short period by collecting other radio data observed with several radio telescopes. Such spectral behaviors are considered to reflect the change in the opacity of the ejected blobs associated with these extreme activities in radio and X-ray. Our 1.4 GHz radio data are expected to be helpful for studying the physics of the accretion and ejection phenomena around black holes., Comment: 9 pages, 10 figures, accepted for publication in the Publications of the Astronomical Society of Japan

We present the results of CO(1-0) observations of the host galaxy of a Type I superluminous supernova (SLSN-I), SN2017egm, one of the closest SLSNe-I at z = 0.03063, by using the Atacama Large Millimeter/submillimeter Array. The molecular gas mass of the host galaxy is $M_{\rm gas} = (4.8 \pm 0.3) \times 10^9$ $M_{\odot}$, placing it on the sequence of normal star-forming galaxies in an $M_{\rm gas}$-star-formation rate (SFR) plane. The molecular hydrogen column density at the location of SN2017egm is higher than that of the Type II SN PTF10bgl, which is also located in the same host galaxy, and those of other Type II and Ia SNe located in different galaxies, suggesting that SLSNe-I have a preference for a dense molecular gas environment. On the other hand, the column density at the location of SN2017egm is comparable to those of Type Ibc SNe. The surface densities of molecular gas and the SFR at the location of SN2017egm are consistent with those of spatially resolved local star-forming galaxies and follow the Schmidt-Kennicutt relation. These facts suggest that SLSNe-I can occur in environments with the same star-formation mechanism as in normal star-forming galaxies., Comment: 6 pages, 5 figures, accepted for publication in PASJ

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

We investigate the molecular gas in, and star-formation properties of, the host galaxy (CGCG 137-068) of a mysterious transient, AT2018cow, at kpc and larger scales, using archival band-3 data from the Atacama Large Millimeter/submillimeter Array (ALMA). AT2018cow is the nearest Fast-Evolving Luminous Transient (FELT), and this is the very first study unveiling molecular-gas properties of FELTs. The achieved rms and beam size are 0.21 mJy beam$^{-1}$ at a velocity resolution of $40$ km s$^{-1}$ and $3".66\times2".71$ ($1.1~{\rm kpc} \times 0.8~{\rm kpc}$), respectively. CO($J$=1-0) emission is successfully detected. The total molecular gas mass inferred from the CO data is $(1.85\pm0.04)\times10^8$ M$_\odot$ with the Milky Way CO-to-H$_2$ conversion factor. The H$_2$ column density at the AT2018cow site is estimated to be $8.6\times10^{20}$ cm$^{-2}$. The ALMA data reveal that (1) CGCG 137-068 is a normal star-forming (SF) dwarf galaxy in terms of its molecular gas and star-formation properties and (2) AT2018cow is located between a CO peak and a blue star cluster. These properties suggest on-going star formation and favor the explosion of a massive star as the progenitor of AT2018cow. We also find that CGCG 137-068 has a solar or super-solar metallicity. If the metallicity of the other FELT hosts is not higher than average, then some property of SF dwarf galaxies other than metallicity may be related to FELTs., Comment: 8 pages, 3 figures, 2 tables, accepted for publication in ApJL

An observational result of a radio continuum and H92$\alpha$ 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.2^{\circ}$ 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 as a part of the GCL, but coincides with a ridge of the radio continuum at a galactic longitude of $0^{\circ}$. 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 the one 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 an explosive activity in the galactic center is unlikely, but favor that the GCL is a giant HII region.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c., Comment: 11 pages, 9 figures, 3 tables. To appear in MNRAS (submitted 2018 June 7; accepted 2019 March 28)

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets., Comment: 8 pages, 6 figures, Published in JKAS

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

We present a kinematic study of the subparsec-scale radio jet of the radio galaxy 3C 84/NGC 1275 with the VLBI Exploration of Radio Astrometry (VERA) array at 22 GHz for 80 epochs from 2007 October to 2013 December. The averaged radial velocity of the bright component "C3" with reference to the radio core is found to be $0.27 pm 0.02c$ between 2007 October and 2013 December. This constant velocity of C3 is naturally explained by the advancing motion of the head of the mini-radio lobe. We also find a non-linear component in the motion of C3 with respect to the radio core. We briefly discuss possible origins of this non-linear motion., Comment: 11 pages, 7 figures, 8 tables (table 1 - 5 are supplementaries), accepted for publication on PASJ

We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Array) and three additional telescopes in Japan (Takahagi 32m, Yamaguchi 32m and Nobeyama 45m). By comparing the independently obtained gain values based on these two methods, we found that the gain values from each method were consistent within 10% at elevations higher than 10 degree. We also found that the total flux densities of two images produced from the different amplitude calibrations were in agreement within 10% at both 22 and 43 GHz. By using the template spectrum method, furthermore, the additional radio telescopes can participate in the KaVA (i.e. EAVN) so that it can give a notable sensitivity increase. Therefore, our results will constrain the detailed conditions to reliably measure the VLBI amplitude using EAVN and give a potential to extend possible telescopes comprising EAVN., Comment: in press, 26 pages, 8 figures, 4 tables

Mid-infrared properties are reported of the west hot spot of the radio galaxy Pictor A with the Wide-field Infrared Survey Explorer (WISE). The mid-infrared counterpart to the hot spot, WISE J051926.26-454554.1, is listed in the AllWISE source catalog. The source was detected in all the four WISE photometric bands. A comparison between the WISE and radio images reinforces the physical association of the \wise\ source to the hot spot. The WISE flux density of the source was carefully evaluated. A close investigation of the multi-wavelength synchrotron spectral energy distribution from the object reveals a mid-infrared excess at the wavelength of $\lambda=22$ $\mu$m with a statistical significance of $4.8 \sigma$ over the simple power-law extrapolation from the synchrotron radio spectrum. The excess is reinforced by single and double cutoff power-law modeling of the radio-to-optical spectral energy distribution. The synchrotron cutoff frequency of the main and excess components was evaluated as $7.1 \times 10^{14}$ Hz and $5.5 \times 10^{13}$ Hz, respectively. From the cutoff frequency, the magnetic field of the emission region was constrained as a function of the region size. In order to interpret the excess component, an electron population different from the main one dominating the observed radio spectrum is necessary. The excess emission is proposed to originate in a sub structure within the hot spot, in which the magnetic field is by a factor of a few stronger than that in the minimum energy condition. The relation of the mid-infrared excess to the X-ray emission is briefly discussed., Comment: 11 pages, 4 figures, ApJ accepted

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., Comment: 17 pages, 11 figures

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations., Comment: 22 pages, 7 figures, 2 tables. Accepted for publication in PASJ

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets., Comment: 10 pages, 9 figures, accepted for publication in PASJ

This chapter summarizes radio astrometry in relation to current very long baseline interferometry (VLBI) projects and describes its perspectives with the SKA. The scientific goals of the astrometry with the SKA have been discussed in the international and Japanese communities of researchers, whose major issues are shown here. We have demonstrated some of the issues, such as censuses of possible targets and the technical feasibility of astrometry in the SKA frequency bands. The preliminary results of our case studies on SKA astrometry are also presented. In addition, possible synergy and commensality of the SKA astrometric projects with those in the optical and infrared astrometric missions, especially JASMINE (Japan Astrometry Satellite Mission for INfrared Exploration) are discussed., Comment: 16 pages, 13 figures, Chapter in the SKA Japan Science Book, provided by SKA Japan Astrometry Science Working Group

We present results from 43 GHz (VLBA, six epochs from 2012.2 to 2013.2) and 86 GHz (GMVA, one epoch in 2012.4) observations toward the basis of the jet in the TeV Blazar Mrk 501. The 43-GHz data analysis reveals a new feature located northeast of the radio core, with a flux density of several tens of mJy, perpendicularly to the jet axis. The 86-GHz image shows the jet feature located 0.75 mas southeast of the radio core, which is consistent with the previous result. The location of Gaussian model for 0.75 mas feature does not coincide with those for the jets in the 43-GHz image, however, a distribution of emission is found. We also discuss the spectral indices of the core, the northeast feature, and the jet feature between 43 GHz and 86 GHz, which show flat-to-steep, steep, and flat-to-invert, respectively., Comment: 4 pages, 2 figures; proceedings of the 12th European VLBI Network Symposium and Users Meeting, 7-10 October 2014, Cagliari, Italy

KaVA (KVN and VERA Array) is a new combined VLBI array with KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). First, we briefly review the imaging capabilities of KaVA array which actually achieves more than three times better dynamic range than that achieved by VERA alone. The KaVA images clearly show detailed structures of extended radio jets in AGNs. Next, we represent the key science program to be led by KaVA AGN sub working group. We will conduct the monitoring observations of Sgr A* and M87 because of the largeness of their central super-massive black hole angular sizes. The main science goals of the program are (i) testing magnetically-driven-jet paradigm by mapping velocity fields of the M87 jet, and (ii) obtaining tight constraints on physical properties of radio emitting region in Sgr A*., Comment: Proceedings of the 12th European VLBI Network Symposium and Users Meeting, 7-10 October 2014, Cagliari, Italy

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., Comment: 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

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays., Comment: 15 pages, 10 figures, 6 tables, accepted for publication in PASJ

We have detected periodic flares of the 6.7 GHz methanol maser from an intermediate-mass star-forming region IRAS22198+6336. The maser was monitored daily in 2011, 2012, and 2013. Six flares were observed with a period of 34.6 days. The variation pattern is intermittent, and the flux ratio of the flaring and the quiescent states exceeds 30. Such intermittent variation with the short period uniquely characterizes the variation of the IRAS22198+6336 maser. At least five spectral components were identified. The spectral components varied almost synchronously, but their peak times differed by 1.8 days. These characteristics can be explained by the colliding-wind binary model., Comment: 17 pages, 6 figures, accepted to PASJ

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion., Comment: 77 pages, 77 figures

We introduce the GENJI program (Gamma-ray Emitting Notable AGN Monitoring by Japanese VLBI), which is a monitoring program of gamma-ray bright AGNs with the VERA array (VLBI Exploration of Radio Astrometry). The GENJI programme aims a dense monitoring at 22 GHz towards the $\gamma$-ray emitting active galactic nuclei (AGNs) to investigate the radio time variation of the core and possible ejection of new radio component, motion of jets, and their relation with the emission at other wavelengths especially in $\gamma$-rays. Currently we are monitoring 8 $\gamma$-ray-emitting notable AGNs (DA 55, 3C 84, M 87, PKS 1510-089, DA 406, NRAO 530, BL Lac, 3C 454.3) about once every two weeks. This programme is promising to trace the trend of radio time variation on shorter timescale than conventional VLBI monitoring programme and to provide complimentary data with them (e.g., MOJAVE, Boston University Blazar Project). In particular, we successfully coordinated quick follow-up observations after the GeV $\gamma$-ray flare in NRAO 530 and 3C 454.3 reported by the Fermi Gamma-ray Space Telescope. Here we present the initial results of morphology and light curves for the first 7-month operation., Comment: 15 pages, 13 figures, accepted for publication in PASJ

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., Comment: 5 pages, 3 figures, 2 tables, accepted for publication in MNRAS

We measured the trigonometric annual parallax of H$_2$O maser source associated with the massive star-forming regions IRAS 06061+2151 with VERA. The annual parallax of $0.496\pm0.031$ mas corresponding to a distance of $2.02^{+0.13}_{-0.12}$ kpc was obtained by 10 epochs of observations from 2007 October to 2009 September. This distance was obtained with a higher accuracy than the photometric distance previously measured, and places IRAS 06061+2151 in the Perseus spiral arm. We found that IRAS 06061+2151 also has a peculiar motion of larger than 15 km s$^{-1}$ counter to the Galactic rotation. That is similar to five sources in the Perseus spiral arm, whose parallaxes and proper motions have already been measured with higher accuracy. Moreover, these sources move at on average 27 km s$^{-1}$ toward the Galactic center and counter to the Galactic rotation., Comment: 8 pages, 4 figures, accepted to PASJ