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.

Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, severely limiting our ability to infer the universal properties of molecular gas in galaxies. Hence, we conducted a legacy project with the 45 m telescope of the Nobeyama Radio Observatory, called the CO Multi-line Imaging of Nearby Galaxies (COMING), which simultaneously observed 147 galaxies with high far-infrared flux in $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ lines. The total molecular gas mass was derived using the standard CO-to-H$_2$ conversion factor and found to be positively correlated with the total stellar mass derived from the WISE $3.4 \mu$m band data. The fraction of the total molecular gas mass to the total stellar mass in galaxies does not depend on their Hubble types nor the existence of a galactic bar, although when galaxies in individual morphological types are investigated separately, the fraction seems to decrease with the total stellar mass in early-type galaxies and vice versa in late-type galaxies. No differences in the distribution of the total molecular gas mass, stellar mass, and the total molecular gas to stellar mass ratio was observed between barred and non-barred galaxies, which is likely the result of our sample selection criteria, in that we prioritized observing FIR bright (and thus molecular gas-rich) galaxies.

We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the volume densities of molecular gas $n(\rm{H}_2)$ in the eight different regions within the galactic disk with CO data combined with archival far-ultraviolet and 24 $\mu$m data. We confirmed that SFE in the bar region is lower by 39% than that in the spiral arms. Moreover, velocity-alignment stacking analysis was performed for the spectra in the individual regions. The integrated intensity ratios of $^{12}$CO to $^{13}$CO ($R_{12/13}$) range from 10 to 17 as the results of stacking. Fixing a kinetic temperature of molecular gas, $n(\rm{H}_2)$ was derived from $R_{12/13}$ via non-local thermodynamic equilibrium (non-LTE) analysis. The density $n(\rm{H}_2)$ in the bar is lower by 31-37% than that in the arms and there is a rather tight positive correlation between SFEs and $n(\rm{H}_2)$, with a correlation coefficient of $\sim 0.8$. Furthermore, we found a dependence of $n(\rm{H}_2)$ on the velocity dispersion of inter-molecular clouds ($\Delta V/ \sin i$). Specifically, $n(\rm{H}_2)$ increases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i < 100$ km s$^{-1}$. On the other hand, $n(\rm{H}_2)$ decreases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i > 100$ km s$^{-1}$. These relations indicate that the variations of SFE could be caused by the volume densities of molecular gas, and the volume densities could be governed by the dynamical influence such as cloud-cloud collisions, shear and enhanced inner-cloud turbulence.

We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galactocentric radius. We found the averaged SFE in 80 galaxies of $(1.69 \pm 1.1) \times 10^{-9}$ yr$^{-1}$, which is consistent with Leroy et al. 2008 if we consider the contribution of helium to the molecular gas mass evaluation and the difference in the assumed initial mass function between two studies. We compared SFE among different morphological (i.e., SA, SAB, and SB) types, and found that SFE within the inner radii ($r/r_{25} < 0.3$, where $r_{25}$ is $B$-band isophotal radius at 25 mag arcsec$^{-2}$) of SB galaxies is slightly higher than that of SA and SAB galaxies. This trend can be partly explained by the dependence of SFE on global stellar mass, which probably relates to the CO-to-H$_2$ conversion factor through the metallicity. For two representative SB galaxies in our sample, NGC 3367 and NGC 7479, the ellipse of $r/r_{25}$ = 0.3 seems to cover not only the central region but also the inner part of the disk, mainly the bar. These two galaxies show higher SFE in the bar than in spiral arms. However, we found an opposite trend in NGC 4303; SFE is lower in the bar than in spiral arms, which is consistent with earlier studies (e.g., Momose et al. 2010). These results suggest diversity of star formation activities in the bar.

We present high-resolution images of the submillimeter SiO line emissions of a massive young stellar object Orion Source I using the Atacama Large Millimeter/ Submillimeter Array (ALMA) at band 8. We detected the 464 GHz SiO v=4 J=11-10 line in Source I, which is the first detection of the SiO v=4 line in star-forming regions, together with the 465 GHz 29SiO v=2 J=11-10 and the 428 GHz SiO v=2 J=10-9 lines with a resolution of 50 AU. The 29SiO v=2 J=11-10 and SiO v=4 J=11-10 lines have compact structures with the diameter of <80 AU. The spatial and velocity distribution suggest that the line emissions are associated with the base of the outflow and the surface of the edge-on disk. In contrast, SiO v=2 J=10-9 emission shows a bipolar structure in the direction of northeast-southwest low-velocity outflow with ~200 AU scale. The emission line exhibits a velocity gradient along the direction of the disk elongation. With the assumption of the ring structure with Keplerian rotation, we estimated the lower limit of the central mass to be 7 solar mass and the radius of 12 AU< r <26 AU.

The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies exhibit a ratio of noncircular to circular velocities of molecular gas larger by a factor of 1.5-2 than non-barred (SA) spiral galaxies at radii within the bar semimajor axis $a_\mathrm{b}$ at 1 kpc resolution, with a maximum at a radius of $R/a_\mathrm{b}\sim0.3$. Residual velocity field images, created by subtracting model velocity fields from the data, reveal that this trend is caused by kpc-scale streaming motions of molecular gas in the bar region. Applying a new method based on radial velocity reversal, we estimated the corotation radius $R_\mathrm{CR}$ and bar pattern speed $\Omega_\mathrm{b}$ in seven SAB and SB systems. The ratio of the corotation to bar radius is found to be in a range of $\mathcal{R}\equiv R_\mathrm{CR}/a_\mathrm{b}\sim0.8\mathrm{-}1.6$, suggesting that intermediate (SBb-SBc), luminous barred spiral galaxies host fast and slow rotator bars. Tentative negative correlations are found for $\Omega_\mathrm{b}$ vs. $a_\mathrm{b}$ and $\Omega_\mathrm{b}$ vs. total stellar mass $M_\ast$, indicating that bars in massive disks are larger and rotate slower, possibly a consequence of angular momentum transfer. The kinematic properties of SAB and SB galaxies, derived from Fourier decomposition, are compared with recent numerical simulations that incorporate various rotation curve models and galaxy interactions.

The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20 GHz as the first step. This paper describes the system design and development of a series of observing instruments and bus systems. In addition to the advantages of avoiding the atmospheric effects of absorption and fluctuation in high frequency radio observation, the mobility of a station can improve the sampling coverage (`uv-coverage') by increasing the number of baselines by the number of ground-based counterparts for each observation day. This benefit cannot be obtained with conventional arrays that solely comprise ground-based stations. The balloon-borne VLBI can contribute to a future progress of research fields such as black holes by direct imaging.

One of the outstanding problems in star-formation theory concerns the transfer of angular momentum such that mass can accrete onto a newly born young stellar object (YSO). From a theoretical standpoint, outflows and jets are predicted to play an essential role in angular momentum transfer and their rotation motions have been reported for both low- and high-mass YSOs. However, little quantitative discussion on outflow launching mechanisms have been presented for high-mass YSOs due to a lack of observational data. Here we present a clear signature of rotation in the bipolar outflow driven by Orion Source I, a high-mass YSO candidate, using the Atacama Large Millimeter/Submillimeter Array (ALMA). A rotational transition of silicon monoxide (Si18O) reveals a velocity gradient perpendicular to the outflow axis which is consistent with that of the circumstellar disk traced by a high-excitation water (H2O) line. The launching radii and outward velocity of the outflow are estimated to be >10 au and 10 km s-1, respectively. These parameters rule out a possibility that the observed outflow is produced by entrainment of a high-velocity jet, and that contribution from stellar-wind or X-wind which have smaller launching radii are significant in the case of Source I. Thus, present results provide a convincing evidence of a rotating outflow directly driven by the magneto-centrifugal disk wind launched by a high-mass YSO candidate.

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

We present observational results of the submillimeter H2O and SiO lines toward a candidate high-mass young stellar object Orion Source I using ALMA. The spatial structures of the high excitation lines at lower-state energies of >2500 K show compact structures consistent with the circumstellar disk and/or base of the northeast-southwest bipolar outflow with a 100 au scale. The highest excitation transition, the SiO (v=4) line at band 8, has the most compact structure. In contrast, lower-excitation transitions are more extended than 200 au tracing the outflow. Almost all the line show velocity gradients perpendicular to the outflow axis suggesting rotation motions of the circumstellar disk and outflow. While some of the detected lines show broad line profiles and spatially extended emission components indicative of thermal excitation, the strong H2O lines at 321 GHz, 474 GHz, and 658 GHz with brightness temperatures of >1000 K show clear signatures of maser action.

The FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope (FUGIN) project is one of the legacy projects using the new multi-beam FOREST receiver installed on the Nobeyama 45-m telescope. This project aims to investigate the distribution, kinematics, and physical properties of both diffuse and dense molecular gas in the Galaxy at once by observing 12CO, 13CO, and C18O J=1-0 lines simultaneously. The mapping regions are a part of the 1st quadrant (10d < l < 50d, |b| < 1d) and the 3rd quadrant (198d < l <236d, |b| < 1d) of the Galaxy, where spiral arms, bar structure, and the molecular gas ring are included. This survey achieves the highest angular resolution to date (~20") for the Galactic plane survey in the CO J=1-0 lines, which makes it possible to find dense clumps located farther away than the previous surveys. FUGIN will provide us with an invaluable dataset for investigating the physics of the galactic interstellar medium (ISM), particularly the evolution of interstellar gas covering galactic scale structures to the internal structures of giant molecular clouds, such as small filament/clump/core. We present an overview of the FUGIN project, observation plan, and initial results, which reveal wide-field and detailed structures of molecular clouds, such as entangled filaments that have not been obvious in previous surveys, and large-scale kinematics of molecular gas such as spiral arms.

Recently a new method to discover obscured active galactic nuclei (AGNs) by utilizing X-ray and Infrared data has been developed. We carried out a survey of H2O maser emission toward ten obscured AGNs with the Nobeyama 45-m telescope. We newly detected the maser emission with the signal-noise-ratio (SNR) of above 4 from two AGNs, NGC 1402 and NGC 7738. We also found a tentative detection with its SNR > 3 in NGC 5037. The detection rate of 20% is higher than those of previous surveys (usually several percents).

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.

We have measured continuum flux densities of a high-mass protostar candidate, a radio source I in the Orion KL region (Orion Source I) using the Atacama Large Millimeter/Submillimeter Array (ALMA) at band 8 with an angular resolution of 0.1". The continuum emission at 430, 460, and 490 GHz associated with Source I shows an elongated structure along the northwest-southeast direction perpendicular to the so-called low-velocity bipolar outflow. The deconvolved size of the continuum source, 90 au times 20 au, is consistent with those reported previously at other millimeter/submillimeter wavelength. The flux density can be well fitted to the optically thick black-body spectral energy distribution (SED), and the brightness temperature is evaluated to be 700-800 K. It is much lower than that in the case of proton-electron or H- free-free radiations. Our data are consistent with the latest ALMA results by Plambeck & Wright (2016), in which the continuum emission have been proposed to arise from the edge-on circumstellar disk via thermal dust emission, unless the continuum source consists of an unresolved structure with the smaller beam filling factor.

We present simultaneous mappings of J=1-0 emission of 12CO, 13CO, and C18O molecules toward the whole disk (8' x 5' or 20.8 kpc x 13.0 kpc) of the nearby barred spiral galaxy NGC 2903 with the Nobeyama Radio Observatory 45-m telescope at an effective angular resolution of 20" (or 870 pc). We detected 12CO(J=1-0) emission over the disk of NGC 2903. In addition, significant 13CO(J=1-0) emission was found at the center and bar-ends, whereas we could not detect any significant C18O(J=1-0) emission. In order to improve the signal-to-noise ratio of CO emission and to obtain accurate line ratios of 12CO(J=2-1)/12CO(J=1-0) ($R_{2-1/1-0}$) and 13CO(J=1-0)/12CO(J=1-0) ($R_{13/12}$), we performed the stacking analysis for our 12CO(J=1-0), 13CO(J=1-0), and archival 12CO(J=2-1) spectra with velocity-axis alignment in nine representative regions of NGC 2903. We successfully obtained the stacked spectra of the three CO lines, and could measure averaged $R_{2-1/1-0}$ and $R_{13/12}$ with high significance for all the regions. We found that both $R_{2-1/1-0}$ and $R_{13/12}$ differ according to the regions, which reflects the difference in the physical properties of molecular gas; i.e., density ($n_{\rm H_2}$) and kinetic temperature ($T_K$). We determined $n_{\rm H_2}$ and $T_K$ using $R_{2-1/1-0}$ and $R_{13/12}$ based on the large velocity gradient approximation. The derived $n_{\rm H_2}$ ranges from ~ 1000 cm$^{-3}$ (in the bar, bar-ends, and spiral arms) to 3700 cm$^{-3}$ (at the center) and the derived $T_K$ ranges from 10 K (in the bar and spiral arms) to 30 K (at the center). We examined the dependence of star formation efficiencies (SFEs) on $n_{\rm H_2}$ and $T_K$, and found the positive correlation between SFE and $n_{\rm H_2}$ with the correlation coefficient for the least-square power-law fit $R^2$ of 0.50. This suggests that molecular gas density governs the spatial variations in SFEs.

We conducted a high-sensitivity radio detection survey for forty narrow-line Seyfert 1 (NLS1) galaxies using very-long-baseline interferometry (VLBI) at 22 GHz through phase-referencing long-time integration and using a newly developing recorder with a data rate of 8 Gbps, which is a candidate of the next generation VLBI data recording systems for the Japanese VLBI Network. The baseline sensitivity was typically a few mJy. The observations resulted in a detection rate of 12/40 for our radio-selected NLS1 sample; 11 out of the detected 12 NLS1s showed inverted radio spectra between 1.4 and 22 GHz on the basis of the Very Large Array flux densities and the VLBI detections. These high fractions suggest that a compact radio core with a high brightness temperature is frequently associated with NLS1 nuclei. On the other hand, at least half of the sample indicated apparently steep spectra even with the limited VLBI sensitivity. Both the inverted and steep spectrum radio sources are included in the NLS1 population.

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.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

Multi-epoch phase-referencing VLBI (Very Long Baseline Interferometry) observations with VERA (VLBI Exploration of Radio Astrometry) were performed for the symbiotic star R Aquarii (R Aqr) from September 2005 to Oct 2006. Tracing one of the $v=2$, $J=1-0$ SiO maser spots, we measured an annual parallax of $\pi = 4.59\pm0.24$ mas, corresponding to a distance of $218_{-11}^{+12}$ pc. Our result is consistent with earlier distance measurements, but yields the highest accuracy of about $5\%$ level. Applying our distance, we derived an absolute K-band magnitude of $M_{\mathrm{K}} = -7.71 \pm 0.11$, which is consistent with the recent Period-Luminosity relation by VLBI parallax measurements for 5 OH-Mira variables. In addition, the expansion age of an inner nebulae around R Aqr is found to be about 240 years, corresponds to about the year 1773.

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.

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

In this Letter, we report detections of SiO v=3 J=1--0 maser emission in very long baseline interferometric (VLBI) observations towards 4 out of 12 long-period variable stars: WX Psc, R Leo, W Hya, and T Cep. The detections towards WX Psc and T Cep are new ones. We also present successful astrometric observations of SiO v=2 and v=3 J=1--0 maser emissions associated with two stars: WX Psc and W Hya and their position-reference continuum sources: J010746.0+131205 and J135146.8-291218 with the VLBI Exploration of Radio Astrometry (VERA). The relative coordinates of the position-reference continuum source and SiO v=3 maser spots were measured with respect to those of an SiO v=2 maser spot adopted as fringe-phase reference. Thus the faint continuum sources were inversely phase-referenced to the bright maser sources. It implies possible registration of multiple SiO maser line maps onto a common coordinate system with 10 microarcsecond-level accuracy.

We present the first upper limit on gravitational wave (GW) backgrounds at an unexplored frequency of 0.2 Hz using a torsion-bar antenna (TOBA). A TOBA was proposed to search for low-frequency GWs. We have developed a small-scaled TOBA and successfully found {\Omega}gw(f) < 4.3 \times 1017 at 0.2 Hz as demonstration of the TOBA's capabilities, where {\Omega}gw (f) is the GW energy density per logarithmic frequency interval in units of the closure density. Our result is the first nonintegrated limit to bridge the gap between the LIGO band (around 100 Hz) and the Cassini band (10-6 - 10-4 Hz).

Multi-epoch Very Long Baseline Interferometry (VLBI) study of the sub-pc scale jet of 3C 84 is presented. We carried out 14-epoch VLBI observations during 2006-2009 with the Japanese VLBI Network (JVN) and the VLBI Exploration of Radio Astrometry (VERA), immediately following the radio outburst that began in 2005. We confirmed that the outburst was associated with the central ~1 pc core, accompanying the emergence of a new component. This is striking evidence of the recurrence of jet activity. The new component became brighter during 2008, in contrast to the constant gamma-ray emission that was observed with the Fermi Gamma-ray Space Telescope during the same time. We found that the projected speed of the new component is 0.23c from 2007/297 (2007 October 24) to 2009/114 (2009 April 24). The direction of movement of this component differs from that of the pre-existing component by ~40 degree. This is the first measurement of kinematics of a sub-pc jet in a gamma-ray active phase. Possible detection of jet deceleration and the jet kinematics in connection with the gamma-ray emission is discussed.

A non-LTE analysis of K I resonance lines at 7664.91 and 7698.97 A was carried out for 15 red giants belonging to three globular clusters of different metallicity (M 4, M 13, and M 15) along with two reference early-K giants (rho Boo and alpha Boo), in order to check whether the K abundances are uniform within a cluster and to investigate the behavior of [K/Fe] ratio at the relevant metallicity range of -2.5 <[Fe/H] < -1. We confirmed that [K/H] (as well as [Fe/H]) is almost homogeneous within each cluster to a precision of < ~0.1 dex, though dubiously large deviations are exceptionally seen for two peculiar stars showing signs of considerably increased turbulence in the upper atmosphere. The resulting [K/Fe] ratios are mildly supersolar by a few tenths of dex for three clusters, tending to gradually increase from ~+0.1-0.2 at [Fe/H] ~-1 to ~+0.3 at [Fe/H] ~ -2.5. This result connects reasonably well with the [K/Fe] trend of disk stars (-1 < [Fe/H]) and that of extremely metal-poor stars (-4 <[Fe/H] < -2.5). That is, [K/Fe] appears to continue a gradual increase from [Fe/H]~0 toward a lower metallicity regime down to [Fe/H]~-3, where a broad maximum of [K/Fe]~+0.3-0.4 is attained, possibly followed by a slight downturn at [Fe/H]<~-3.

We analyze the sidereal daily variations observed between 1985 and 2006 at Matsushiro, Japan (MAT) and between 1993 and 2005 at Liapootah, Tasmania (LPT). These stations comprise the two hemisphere network (THN) of underground muon detectors in Japan and Australia. Yearly mean harmonic vectors at MAT and LPT are more or less stable without any significant change in phase and amplitude in correlation with either the solar activity- or magnetic-cycles. In this paper, therefore, we analyze the average anisotropy over the entire observation periods, i.e. 1985-2006 for MAT and 1993-2005 for LPT. We apply to the THN data a best-fitting analysis based on a model anisotropy in space identical to that adopted by Amenomori et al. (2007) for Tibet III data. The median energies of primary cosmic rays recorded are ~0.5 TeV for THN and ~5 TeV for the Tibet III experiment. It is shown that the intensity distribution of the best-fit anisotropy is quite similar to that derived from Tibet III data, regardless of the order of magnitude difference in energies of primary particles. This, together with the THN observations, confirms that the analysis by Amenomori et al. (2007) based on the Tibet III experiment in the northern hemisphere is not seriously biased. The best-fit amplitudes of the anisotropy, on the other hand, are only one third or less of those reported by the Tibet III experiment, indicating attenuation due to solar modulation. The rigidity dependence of the anisotropy amplitude in the sub-TeV region is consistent with the spectrum reported by Hall et al. (1999), smoothly extending to the Tibet III result in the multi-TeV region. The amplitude at higher energies appears almost constant or gradually decreasing with increasing rigidity.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.