We report ALMA observations of NGC 1052 to quest mass accretion in a gas-poor active galactic nucleus (AGN). We detected CO emission representing a rotating ring-like circumnuclear disk (CND) seen edge-on with the gas mass of $5.3 \times 10^{5}$ M$_{\odot}$. The CND has smaller gas mass than that in typical Seyfert galaxies with circumnuclear star formation and is too gas-poor to drive mass accretion onto the central engine. The continuum emission casts molecular absorption features of CO, HCN, HCO$^+$, SO, SO$_2$, CS, CN, and H$_2$O, with H$^{13}$CN and HC$^{15}$N and vibrationally-excited (v$_2 = 1$) HCN and HCO$^+$. Broader absorption line widths than CND emission line widths imply presence of a geometrically thick molecular torus with a radius of $2.4 \pm 1.3$ pc and a thickness ratio of $0.7 \pm 0.3$. We estimate the H$_2$ column density of $(3.3 \pm 0.7) \times 10^{25}$ cm$^{-2}$ using H$^{12}$CN, H$^{13}$CN, and HCO$^{+}$ absorption features and adopting abundance ratio of $^{12}$C-to-$^{13}$C and a HCO$^{+}$-to-H$_2$, and derived the torus gas mass of $(1.3 \pm 0.3) \times 10^7$ M$_{\odot}$, which is $\sim 9$\% of the central black-hole mass. The molecular gas in the torus is clumpy with the estimated covering factor of $0.17^{+0.06}_{-0.03}$. The gas density of clumps inside the torus is inferred to be $(6.4 \pm 1.3) \times 10^7$ cm$^{-3}$, which meets the excitation conditions of H$_2$O maser. Specific angular momentum in the torus exceeds a flat-rotation curve extrapolated from that of the CND, indicating a Keplerian rotation inside a 14.4-pc sphere of influence.

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.

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.

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

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.

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.

In the circumnuclear region of the nearby radio galaxy NGC~1052. Our $\lambda$3mm VLBI observations with the Korean VLBI Network have spatially resolved the broad HCO$^{+}$ absorption at a velocity range of 1350--1850 km~s$^{-1}$ against a double-sided nuclear jet, and have revealed that the HCO$^{+}$ absorption is concentrated on the receding jet and the nuclear components. The distribution of the HCO$^{+}$ absorbing gas strongly supports the circumnuclear torus surrounding the supermassive black hole. From the estimations of the column density and the volume density of molecular hydrogen, the size of the molecular gas region in the torus is at least 1~pc. The broad spectral profile of HCO$^{+}$ is likely to be a blend of multiple gas clumps with various velocities. The HCO$^{+}$ absorption of NGC~1052 could trace complex kinematics in the vicinity of the active galactic nucleus, such as inflow, outflow, turbulence, and so on.

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.

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.

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.

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.

We present the first VLBI detection of HCN molecular absorption in the nearby active galactic nucleus NGC 1052. Utilizing the 1 milliarcsecond resolution achieved by the Korean VLBI Network, we have spatially resolved the HCN absorption against a double-sided nuclear jet structure. Two velocity features of HCN absorption are detected significantly at the radial velocity of 1656 and 1719 km/s, redshifted by 149 and 212 km/s with respect to the systemic velocity of the galaxy. The column density of the HCN molecule is estimated to be 10^{15}-10^{16} cm^{-2}, assuming the excitation temperature of 100-230 K. The absorption features show high optical depth localized on the receding jet side, where the free-free absorption occurred due to the circumnuclear torus. The size of the foreground absorbing molecular gas is estimated to be on approximately one-parsec scales, which agrees well with the approximate size of the circumnuclear torus. HCN absorbing gas is likely to be several clumps smaller than 0.1 parsec inside the circumnuclear torus. The redshifted velocities of the HCN absorption features imply that HCN absorbing gas traces ongoing infall motion inside the circumnuclear torus onto the central engine.

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at 120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi Large Area Telescope. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure in the radio band. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central black hole.

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at >120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure at radio. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central supermassive black hole.

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.

We report our intensive radio monitoring observations of the jet in M87 with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi-LAT. During this period, an elevated level of the M87 flux is reported at VHE gamma rays. We detected a remarkable increase of the radio flux density from the unresolved jet base (radio core) with VERA at 22 and 43GHz coincident with the VHE activity. Meanwhile, we confirmed with EVN at 5GHz that HST-1 (an alternative gamma-ray production candidate site) remained quiescent in terms of its flux density and structure. These results in the radio bands strongly suggest that the VHE gamma-ray activity in 2012 originates in the jet base within 0.03pc or 56 Schwarzschild radii from the central supermassive black hole. We further conducted VERA astrometry for the M87 core during the flaring period, and detected core shifts between 22 and 43GHz. We also discovered a clear frequency-dependent evolution of the radio core flare at 43, 22 and 5GHz; the radio flux density increased more rapidly at higher frequencies with a larger amplitude, and the light curves clearly showed a time-lag between the peaks at 22 and 43GHz. This indicates that a new radio-emitting component was created near the black hole in the period of the VHE event, and then propagated outward with progressively decreasing synchrotron opacity. By combining these results, we estimated an apparent speed of the newborn component, and derived a sub-luminal speed of less than ~0.2c. This value is significantly slower than the super-luminal (~1.1c) features that appeared from the core during the prominent VHE flaring event in 2008, suggesting that the stronger VHE activity can be associated with the production of the higher Lorentz factor jet.

We have measured the internal proper motions of the 6.7 GHz methanol masers associated with Cepheus A (Cep A) HW2 using Very Long Baseline Interferometery (VLBI) observations. We conducted three epochs of VLBI monitoring observations of the 6.7 GHz methanol masers in Cep A-HW2 with the Japanese VLBI Network (JVN) over the period between 2006-2008. In 2006, we were able to use phase-referencing to measure the absolute coordinates of the maser emission with an accuracy of a few milliarcseconds. We compared the maser distribution with other molecular line observations that trace the rotating disk. We measured the internal proper motions for 29 methanol maser spots, of which 19 were identified at all three epochs and the remaining ten at only two epochs. The magnitude of proper motions ranged from 0.2 to 7.4 km/s, with an average of 3.1 km/s. Although there are large uncertainties in the observed internal proper motions of the methanol maser spots in Cep A, they are well fitted by a disk that includes both rotation and infall velocity components. The derived rotation and infall velocities at the disk radius of 680 au are 0.5 +- 0.7 and 1.8 +- 0.7 km/s, respectively. Assuming that the modeled disk motion accurately represents the accretion disk around the Cep A-HW2 high-mass YSO, we estimated the mass infall rate to be 3 x 10^{-4} n_8 Msun/yr (n_8 is the gas volume density in units of 10^{8} cm^{-3}). The combination of the estimated mass infall rate and the magnitude of the fitted infall velocity suggests that Cep A-HW2 is at an evolutionary phase of active gas accretion from the disk onto the central high-mass YSO. The infall momentum rate is estimated to be 5 x 10^{-4} n_8 Msun/yr km/s, which is larger than the estimated stellar radiation pressure of the HW2 object, supporting the hypothesis that this object is in an active gas accretion phase.

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.

The KVN+VERA array is a joint VLBI project of seven VLBI stations spread throughout Korea and Japan. Since the first fringe detection in 2008, the early phase observations of the KVN+VERA have been carried out every several months. Currently, two observing bands of 22 and 43 GHz are available. We are aiming for early realization of science observations with the 1-Gbps recording system from 2012.

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

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.

We searched a short-term radio variability in an active galactic nucleus PKS 1510-089. A daily flux monitoring for 143 days at 8.4 GHz was performed, and VLBI observations at 8.4, 22, and 43 GHz were carried out 4 times during the flux monitoring period. As a result, variability with time scale of 20 to 30 days was detected. The variation patterns were well alike on three frequencies, moreover those at 22 and 43 GHz were synchronized. These properties support that this short-term variability is an intrinsic one. The Doppler factor estimated from the variability time scale is 47. Since the Doppler factor is not extraordinary large for AGN, such intrinsic variability with time scale less than 30 days would exist in other AGNs.

We have measured the annual parallax of the H2O maser source associated with an infrared dark cloud MSXDC G034.43+00.24 from the observations with VERA (VLBI Exploration of Radio Astrometry). The parallax is 0.643 +/- 0.049 mas, corresponding to the distance of 1.56 +0.12/-0.11 kpc. This value is less than the half of the previous kinematic distance of 3.7 kpc. We revise the core mass estimates of MSXDC G034.43+00.24, based on virial masses, LTE masses and dust masses and show that the core masses decrease from the previous estimations of ~1000 Mo to hundreds of Mo. The spectral type derived from the luminosity also changes from O9.5 to B1 in the case of MM1. This spectral type is still consistent with that of the massive star. The radial velocity derived from the flat rotation model is smaller than the observed velocity, which corresponds to the peculiar motion of ~40 km/s in the line-of-sight direction.

We present multi-frequency simultaneous VLBA observations at 15, 22 and 43 GHz towards the nucleus of the nearby radio galaxy NGC 1052. These three continuum images reveal a double-sided jet structure, whose relative intensity ratios imply that the jet axis is oriented close to the sky plane. The steeply rising spectra at 15-43 GHz at the inner edges of the jets strongly suggest that synchrotron emission is absorbed by foreground thermal plasma. We detected H_2O maser emission in the velocity range of 1550-1850 km/s, which is redshifted by 50-350 km/s with respect to the systemic velocity of NGC 1052. The redshifted maser gas appears projected against both sides of the jet, similar to the HI seen in absorption. The H_2O maser gas is located where the free-free absorption opacity is large. This probably implies that the masers in NGC 1052 are associated with a circumnuclear torus or disk as in the nucleus of NGC 4258. Such circumnuclear structure can be the sense of accreting onto the central engine.

We have used the Submillimeter Array (SMA) to make the first interferometric observations (beam size ~1") of the 12CO J=6-5 line and 435 micron (690 GHz) continuum emission toward the central region of the nearby ULIRG Arp 220. These observations resolve the eastern and western nuclei from each other, in both the molecular line and dust continuum emission. At 435 micron, the peak intensity of the western nucleus is stronger than the eastern nucleus, and the difference in peak intensities is less than at longer wavelengths. Fitting a simple model to the dust emission observed between 1.3 mm and 435 micron suggests that dust emissivity power law index in the western nucleus is near unity and steeper in the eastern nucleus, about 2, and that the dust emission is optically thick at the shorter wavelength. Comparison with single dish measurements indicate that the interferometer observations are missing ~60% of the dust emission, most likely from a spatially extended component to which these observations are not sensitive. The 12CO J=6-5 line observations clearly resolve kinematically the two nuclei. The distribution and kinematics of the 12CO J=6-5 line appear to be very similar to lower J CO lies observed at similar resolution. Analysis of multiple 12CO line intensities indicates that the molecular gas in both nuclei have similar excitation conditions, although the western nucleus is warmer and denser. The excitation conditions are similar to those found in other extreme environments, including M82, Mrk 231, and BR 1202-0725. Simultaneous lower resolution observations of the 12CO, 13CO, and C18O J=2-1 lines show that the 13CO and C18O lines have similar intensities, which suggests that both of these lines are optically thick, or possibly that extreme high mass star formation has produced in an overabundance of C18O.

We have mapped the central region of the Seyfert 1 galaxy NGC 1097 in 12CO(J=2-1) with the Submillieter Array (SMA). The 12CO(J=2-1) map shows a central concentration and a surrounding ring, which coincide respectively with the Seyfert nucleus and a starburst ring. The line intensity peaks at the nucleus, whereas in a previously published 12CO(J=1-0) map the intensity peaks at the starburst ring. The molecular ring has an azimuthally averaged 12CO(J=2-1)/(J=1-0) intensity ratio (R21) of about unity, which is similar to those in nearby active star forming galaxies, suggesting that most of the molecular mass in the ring is involved in fueling the starburst. The molecular gas can last for only about 1.2\times10^8 years without further replenishment assuming a constant star formation rate and a perfect conversion of gas to stars. The velocity map shows that the central molecular gas is rotating with the molecular ring in the same direction, while its velocity gradient is much steeper than that of the ring. This velocity gradient of the central gas is similar to what is usually observed in some Seyfert 2 galaxies. To view the active nucleus directly in the optical, the central molecular gas structure can either be a low-inclined disk or torus but not too low to be less massive than the mass of the host galaxy itself, be a highly-inclined thin disk or clumpy and thick torus, or be an inner part of the galactic disk. The R21 value of ~1.9 of the central molecular gas component, which is significantly higher than the value found at the molecular gas ring, indicates that the activity of the Seyfert nucleus may have a significant influence on the conditions of the molecular gas in the central component.

We present multi-frequency simultaneous VLBA observations at 15, 22 and 43 GHz towards the nucleus of the nearby radio galaxy NGC 1052. These three continuum images reveal a double-sided jet structure, whose relative intensity ratios imply that the jet axis is oriented close to the sky plane. The steeply rising spectra at 15-43 GHz at the inner edges of the jets strongly suggest that synchrotron emission is absorbed by foreground thermal plasma. We detected H2O maser emission in the velocity range of 1550-1850 km/s, which is redshifted by 50-350 km/s with respect to the systemic velocity of NGC 1052. The redshifted maser gas appears projected against both sides of the jet, in the same manner as the HI seen in absorption. The H2O maser gas are located where the free-free absorption opacity is large. This probably imply that the masers in NGC 1052 are associated with a circumnuclear torus or disk as in the nucleus of NGC 4258. Such circumnuclear structure can be the sence of accreting onto the central engine.

We have mapped the central region of NGC 4945 in the $J=2\to1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the continuum at 1.3 mm, at an angular resolution of $5\farc \times 3\farc$ with the Submillimeter Array. The relative proximity of NGC 4945 (distance of only 3.8 Mpc) permits a detailed study of the circumnuclear molecular gas and dust in a galaxy exhibiting both an AGN (classified as a Seyfert 2) and a circumnuclear starburst in an inclined ring with radius $\sim$2\farcs5 ($\sim$50 pc). We find that all three molecular lines trace an inclined rotating disk with major axis aligned with that of the starburst ring and large-scale galactic disk, and which exhibits solid-body rotation within a radius of $\sim$5\farc ($\sim$95 pc). We infer an inclination for the nuclear disk of $62^{\circ} \pm 2^{\circ}$, somewhat smaller than the inclination of the large-scale galactic disk of $\sim$$78^{\circ}$. The continuum emission at 1.3 mm also extends beyond the starburst ring, and is dominated by thermal emission from dust. If it traces the same dust emitting in the far-infrared, then the bulk of this dust must be heated by star-formation activity rather than the AGN. We discover a kinematically-decoupled component at the center of the disk with a radius smaller than $1\farcs4$ (27 pc), but which spans approximately the same range of velocities as the surrounding disk. This component has a higher density than its surroundings, and is a promising candidate for the circumnuclear molecular torus invoked by AGN unification models.

We present 12CO (J=2-1) observations towards the central region of the Seyfert 2 galaxy NGC 4258 with the Submillimeter Array (SMA). Our interferometric maps show two arm-like elongated components along the major axis of the galaxy, with no strong nuclear concentration. The CO (2-1) morphology and kinematics are similar to previous CO (1-0) results. The velocity field of the components agrees with the general galactic rotation, except for the east elongated component, which shows a significant velocity gradient along the east-west direction. In order to account for the velocity field, we propose the kinematical model where the warped rotating disk is also expanding. The line ratio of CO(2-1)/CO(1-0) reveals that the eastern component with the anomalous velocity gradient appears to be warmer and denser. This is consistent with the gas in this component being closer to the center, being heated by the central activities, and possibly interacted by expanding motions from the nuclear region.

We report results of pentachromatic VLBI survey for 18 GHz-peaked spectrum sources. Spectral fitting at every pixel across five frequencies allows us to illustrate distribution of optical depth in terms of free-free absorption or synchrotron self absorption. Quasars and Seyfert 1 sources show one-sided morphology with a core at the end where the optical depth peaks. Radio galaxies and Seyfert 2 show symmetric double-sided jets with a optically thick core at the center.

We report Very-Long-Baseline Array (VLBA) observations at 2.3, 8.4, and 15.4 GHz towards nine GHz-Peaked Spectrum (GPS) sources. One Seyfert 1 galaxy, one Seyfert 2 galaxy, three radio galaxies, and four quasars were included in our survey. We obtained spatial distributions of the Free-Free Absorption (FFA) opacity with milliarcsec resolutions for all sources. It is found that type-1 (Seyfert 1 and quasars) and type-2 (Seyfert 2 and radio galaxies) sources showed different distributions of the FFA opacities. The type-1 sources tend to show more asymmetric opacity distributions towards a double lobe, while those of the type-2 sources are rather symmetric. Our results imply that the different viewing angle of the jet causes the difference of FFA opacity along the external absorber. This idea supports the unified scheme between quasars and radio galaxies, proposed by Barthel (1989).

We report results from nearly simultaneous pentachromatic VLBI observations towa rds a nearby GPS galaxy NGC 1052. The observations at 1.6 and 4.8 GHz with VSOP, and at 2.3, 8.4, and 15.4 GHz wit h VLBA, provide linear resolutions of $\sim 0.1$ pc. Convex spectra of a double-sided jet imply that synchrotron emission is obscured through foreground cold dense plasma, in terms of free--free absorption (FFA). We found a central condensation of the plasma which covers about 0.1 pc and 1 pc of the approaching and receding jets, respectively. A simple model with a geometrically thick plasma torus perpendicular to the jets is established to explain the asymmetric distribution of FFA opacities.

A subparsec-scale dense plasma torus around an active galactic nucleus (AGN) is unveiled. We report on very-long-baseline interferometry (VLBI) observations at 2.3, 8.4, and 15.4 GHz towards the active galaxy NGC 1052. The convex spectra of the double-sided jets and the nucleus imply that synchrotron emission is obscured through free--free absorption (FFA) by the foreground cold dense plasma. A trichromatic image was produced to illustrate the distribution of the FFA opacity. We found a central condensation of the plasma which covers about 0.1 pc and 0.7 pc of the approaching and receding jets, respectively. A simple explanation for the asymmetric distribution is the existence of a thick plasma torus perpendicular to the jets. We also found an ambient FFA absorber, whose density profile can be ascribed to a spherical distribution of the isothermal King model. The coexistence of torus-like and spherical distributions of the plasma suggests a transition from radial accretion to rotational accretion around the nucleus.

In order to investigate the genesis of powerful radio jet, we have mapped the central 10 pc region of the nearby radio galaxy NGC 6251 with a 0.2 pc resolution using Very Long Baseline Interferometer (VLBI) at two radio frequencies, 5 GHz and 15 GHz, we have found the sub-parsec-scale counterjet for the first time in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic beaming model.

Mapping the central 5 pc region of the nearby radio galaxy NGC 6251 with a 0.2 pc resolution using VLBI at two radio frequencies, 5 GHz and 15 GHz, we have found the sub-pc-scale counter jet for the first time in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic beaming model (Ghisellini et al. 1993).