The nearby bright radio galaxy 3C 84 at the center of Perseus cluster is one of the ideal targets to explore the jet in active galactic nuclei (AGNs) and its parsec-scale environment. A recent research of Fujita & Nagai revealed the existence of the northern counter-jet component (N1) located at $\sim$ 2 mas north from the central core in VLBI images at 15 and 43 GHz and they are explained by the free-free absorption (FFA) due to an ionized plasma foreground. Here we report a new quasi-simultaneous observation of 3C 84 with the Korean VLBI Network (KVN) at 86 GHz and the KVN and VERA (VLBI Exploration of Radio Astrometry) Array (KaVA) at 43 GHz in 2016 February. We succeeded the first detection of N1 at 86 GHz and the data show that N1 still has an inverted spectrum between 43 and 86 GHz with its spectral index $\alpha$ ($S_{\nu} \propto \nu^{\alpha}$) of $1.19 \pm 0.43$, while the approaching lobe component has the steep spectrum with the index of $-0.54 \pm 0.30$. Based on the measured flux asymmetry between the counter and approaching lobes, we constrain the averaged number density of the FFA foreground $n_{\rm e}$ as $1.8 \times 10^{4}~{\rm cm^{-3}} \lesssim n_{\rm e} \lesssim 1.0 \times 10^{6}~{\rm cm^{-3}}$. Those results suggest that the observational properties of the FFA foreground can be explained by the dense ionized gas in the circumnuclear disk and/or assembly of clumpy clouds at the central $\sim 1$ pc region of 3C 84.

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 present ALMA Band 6 observations of the CO(2-1), HCN(3-2), and HCO$^{+}$(3-2) lines in the nearby radio galaxy / brightest cluster galaxy (BCG) of NGC 1275 with the spatial resolution of $\sim20$ pc. In the previous observations, CO(2-1) emission was detected as radial filaments lying in the east-west direction. We resolved the inner filament and found that the filament cannot be represented by a simple infalling stream both morphologically and kinematically. The observed complex nature of the filament resembles the cold gas structure predicted by recent numerical simulations of cold chaotic accretion. A crude estimate suggests that the accretion rate of the cold gas can be higher than that of hot gas. Within the central 100 pc, we detected a rotational disk of the molecular gas whose mass is $\sim10^{8} M_{\sun}$. This is the first evidence of the presence of massive cold gas disk on this spatial scale for BCGs. The disk rotation axis is approximately consistent with the axis of the radio jet on subpc scales. This probably suggests that the cold gas disk is physically connected to the innermost accretion disk which is responsible for jet launching. We also detected absorption features in the HCN(3-2) and HCO$^{+}$(3-2) spectra against the radio continuum emission mostly radiated by $\sim1.2$-pc size jet. The absorption features are blue-shifted from the systemic velocity by $\sim$300-600~km~s$^{-1}$, which suggests the presence of outflowing gas from the active galactic nucleus (AGN). We discuss the relation of the AGN feeding with cold accretion, the origin of blue-shifted absorption, and estimate of black hole mass using the molecular gas dynamics.

The flat spectrum radio quasar 4C 38.41 showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional monitoring observations at various bands across the electromagnetic spectrum. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components and the evolution of the physical properties seem to be in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1 mG, do not seem to dramatically change during the flares, and are much smaller, by a factor 10,000, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting 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.

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.

Radio jets in active galaxies have been expected to interact with circumnuclear environments in their early phase evolutions. By performing the multi-epoch monitoring observation with the KVN and VERA Array (KaVA) at 43~GHz, we investigate the kinematics of the notable newborn bright component C3 located at the tip of the recurrent jet of 3C~84. During 2015 August-September, we discover the flip of C3 and the amount of the flip is about 0.4~milli-arcsecond in angular scale, which corresponds to 0.14 parsec in physical scale. After the flip of C3, it wobbled at the same location for a few months and then it restarted to propagate towards the southern direction. The flux density of C3 coherently showed the monotonic increase during the observation period. The flip is in good agreement with hydrodynamical simulations of jets in clumpy ambient medium. We estimate the number density of the putative clump based on the momentum balance between the jet thrust and the ram pressure from the clump and it is about $10^{3-5}~{\rm cm^{-3}}$. We briefly discuss possible origins of the clump.

Recent Fermi-Large Area Telescope (LAT) light curves indicate an active $\gamma$-ray state spanning about five months from 2016 June to 2016 October in the BL Lac object 1749+096 (OT 081). During this period, we find two notable $\gamma$-ray events: an exceptionally strong outburst followed by a significant enhancement (local peak). In this study, we analyze multi-waveband light curves (radio, optical, X-ray, and $\gamma$-ray) plus very-long baseline interferometry (VLBI) data to investigate the nature of the $\gamma$-ray events. The $\gamma$-ray outburst coincides with flux maxima at longer wavelengths. We find a spectral hardening of the $\gamma$-ray photon index during the $\gamma$-ray outburst. The photon index shows a transition from a softer-when-brighter to a harder-when-brighter trend at around 1.8 $\times$ $10^{-7}$ ph cm$^{-2}$ s$^{-1}$. We see indication that both the $\gamma$-ray outburst and the subsequent enhancement precede the propagation of a polarized knot in a region near the VLBI core. The highest polarized intensity, 230\,mJy, and an electric vector position angle rotation, by $\sim$32$^{\circ}$, are detected about 12 days after the $\gamma$-ray outburst. We conclude that both $\gamma$-ray events are caused by the propagation of a disturbance in the mm-wave core.

We investigated the detailed radio structure of the jet of 1H 0323+342 using high-resolution multi-frequency Very Long Baseline Array observations. This source is known as the nearest $\gamma$-ray emitting radio-loud narrow-line Seyfert 1 (NLS1) galaxy. We discovered that the morphology of the inner jet is well characterized by a parabolic shape, indicating the jet being continuously collimated near the jet base. On the other hand, we found that the jet expands more rapidly at larger scales, resulting in a conical-like shape. The location of the "collimation break" is coincident with a bright quasi-stationary feature at 7 mas from core (corresponding to a deprojected distance of the order of $\sim$100pc), where the jet width locally contracts together with highly polarized signals, suggesting a recollimation shock. We found that the collimation region is coincident with the region where the jet speed gradually accelerates, suggesting the coexistence of the jet acceleration and collimation zone, ending up with the recollimation shock, which could be a potential site of high-energy $\gamma$-ray flares detected by the Fermi-LAT. Remarkably, these observational features of the 1H 0323+342 jet are overall very similar to those of the nearby radio galaxy M87 and HST-1 as well as some blazars, suggesting that a common jet formation mechanism might be at work. Based on the similarity of the jet profile between the two sources, we also briefly discuss the mass of the central black hole of 1H 0323+342, which is also still highly controversial on this source and NLS1s in general.

The flat spectrum radio quasar 1633+382 (4C~38.41) showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional radio monitoring observations with the OVRO 40m telescope, the Boston University VLBI program, and the Submillimeter Array. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components. Analysis of the spectral energy distribution indicates a higher turnover frequency after the flaring events. The evolution of the flare in the turnover frequency - turnover flux density plane probes the adiabatic losses in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1mG, do not seem to dramatically change during the flares, and are much smaller, by a factor $10^4$, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

We report the discovery of a local convergence of a jet cross section in the quasi-stationary jet feature in the gamma-ray-emitting narrow-line Seyfert 1 galaxy (NLS1) 1H 0323+342. The convergence site is located at approximately 7 mas (corresponding to the order of 100 pc in deprojection) from the central engine. We also found limb-brightened jet structures at both the upstream and downstream of the convergence site. We propose that the quasi-stationary feature showing the jet convergence and limb-brightening occurs as a consequence of recollimation shock in the relativistic jets. The quasi-stationary feature is one of the possible gamma-ray-emitting sites in this NLS1, in analogy with the HST-1 complex in the M87 jet. Monitoring observations have revealed that superluminal components passed through the convergence site and the peak intensity of the quasi-stationary feature, which showed apparent coincidences with the timing of observed gamma-ray activities.

3C 84 (NGC 1275) is a well-studied mis-aligned Active Galactic Nucleus (AGN), which has been active in Gamma rays since at least 2008. We have monitored the source at four wavelengths (14 mm, 7 mm, 3 mm and 2 mm) using the Korean VLBI network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C 84 exhibits bright radio emission both near the central supermassive black hole (SMBH) feature known as C1 and from a moving feature located to the south known as C3. Other facilities have also detected these short-term variations above a slowly rising trend at shorter wavelengths, such as in Gamma ray and 1 mm total intensity light-curves. We find that the variations in the $\gamma$ rays and 1 mm total intensity light-curves are correlated, with the $\gamma$ rays leading and lagging the radio emission. Analysis of the 2 mm KVN data shows that both the Gamma rays and 1 mm total intensity short-term variations are better correlated with the SMBH region than C3, likely placing the short-term variations in C1. We interpret the emission as being due to the random alignment of spatially separated emission regions. We place the slowly rising trend in C3, consistent with previous results. Additionally, we report that since mid-2015, a large mm-wave radio flare has been occurring in C3, with a large Gamma ray flare coincident with the onset of this flare at all radio wavelengths.

Atmospheric propagation effects at millimeter wavelengths can significantly alter the phases of radio signals and reduce the coherence time, putting tight constraints on high frequency Very Long Baseline Interferometry (VLBI) observations. In previous works, it has been shown that non-dispersive (e.g. tropospheric) effects can be calibrated with the frequency phase transfer (FPT) technique. The coherence time can thus be significantly extended. Ionospheric effects, which can still be significant, remain however uncalibrated after FPT as well as the instrumental effects. In this work, we implement a further phase transfer between two FPT residuals (i.e. so-called FPT-square) to calibrate the ionospheric effects based on their frequency dependence. We show that after FPT-square, the coherence time at 3 mm can be further extended beyond 8~hours, and the residual phase errors can be sufficiently canceled by applying the calibration of another source, which can have a large angular separation from the target (>20 deg) and significant temporal gaps. Calibrations for all-sky distributed sources with a few calibrators are also possible after FPT-square. One of the strengths and uniqueness of this calibration strategy is the suitability for high-frequency all-sky survey observations including very weak sources. We discuss the introduction of a pulse calibration system in the future to calibrate the remaining instrumental effects and allowing the possibility of imaging the source structure at high frequencies with FPT-square, where all phases are fully calibrated without involving any additional sources.

We present multi-frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C~38.41) as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program combined with additional observations in radio, optical, X-rays and $\gamma-$rays carried out between the period 2012 March - 2015 August. The monitoring of this source reveals a significant long-lived increase in its activity since approximately two years in the radio bands, which correlates with a similar increase on all other bands from sub-millimeter to $\gamma-$rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function (DCF) indicates time lags smaller than the $\sim40$ days uncertainties among both radio bands and also high-energy bands, and a time lag of $\sim$70 days, with $\gamma-$rays leading radio. We interpret that the high-energy and radio fluxes are arising from different emitting regions, located at $1\pm12$ and $40\pm13$ pc from the central engine respectively.

We present a study of the inexplicit connection between radio jet activity and gamma-ray emission of BL Lacertae (BL Lac; 2200+420). We analyze the long-term millimeter activity of BL Lac via interferometric observations with the Korean VLBI Network (KVN) obtained at 22, 43, 86, and 129 GHz simultaneously over three years (from January 2013 to March 2016); during this time, two gamma-ray outbursts (in November 2013 and March 2015) can be seen in gamma-ray light curves obtained from Fermi observations. The KVN radio core is optically thick at least up to 86 GHz; there is indication that it might be optically thin at higher frequencies. To first order, the radio light curves decay exponentially over the time span covered by our observations, with decay timescales of 411+/-85 days, 352+/-79 days, 310+/-57 days, and 283+/-55 days at 22, 43, 86, and 129 GHz, respectively. Assuming synchrotron cooling, a cooling time of around one year is consistent with magnetic field strengths B~2microT and electron Lorentz factors gamma~10,000. Taking into account that our formal measurement errors include intrinsic variability and thus over-estimate the statistical uncertainties, we find that the decay timescale tau scales with frequency nu like tau~nu^{-0.2}. This relation is much shallower than the one expected from opacity effects (core shift), but in agreement with the (sub-)mm radio core being a standing recollimation shock. We do not find convincing radio flux counterparts to the gamma$ray outbursts. The spectral evolution is consistent with the `generalized shock model' of Valtaoja et al. (1992). A temporary increase in the core opacity and the emergence of a knot around the time of the second gamma-ray event indicate that this gamma-ray outburst might be an `orphan' flare powered by the `ring of fire' mechanism.

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 explore physical properties of the shocked external medium (i.e., a shell) in 3C 84 associated with the recurrent radio lobe born around 1960. In the previous work of Ito et al., we investigated a dynamical and radiative evolution of such a shell after the central engine stops the jet launching and we found that a fossil shell emission overwhelms that of the rapidly fading radio lobe. We apply this model to 3C 84 and find the followings: (i) The fossil shell made of shocked diffuse ambient matter with the number density of 0.3 cm$^{-3}$ radiates bright Inverse-Compton (IC) emission with the seed photons of the radio emission from the central compact region and the IC emission is above the sensitivity threshold of the Cherenkov Telescope Array (CTA). (ii) When the fossil shell is produced in a geometrically thick ionized plasma with the number density of $10^{3}$ cm$^{-3}$ and the field strength in the shell may reach about 17 mG in the presence of magnetic fields amplification and the radio emission becomes comparable to the sensitivity of deep imaging VLBI observations. A possible production of ultra high energy cosmic-rays (UHECRs) in the dense shocked plasma is also argued.

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

We present results of single-epoch very long baseline interferometry (VLBI) observations of gamma-ray bright active galactic nuclei (AGNs) using the Korean VLBI Network (KVN) at 22, 43, 86, and 129~GHz bands, which are part of a KVN key science program, Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA). We selected a total of 34 radio-loud AGNs of which 30 sources are gamma-ray bright AGNs with flux densities of $>6\times10^{-10}$~ph~cm$^{-2}$~s$^{-1}$. Single-epoch multi-frequency VLBI observations of the target sources were conducted during a 24-hr session on 2013 November 19 and 20. All observed sources were detected and imaged at all frequency bands with or without a frequency phase transfer technique which enabled the imaging of 12 faint sources at 129~GHz, except for one source. Many of the target sources are resolved on milliarcsecond scales, yielding a core-jet structure with the VLBI core dominating the synchrotron emission on the milliarcsecond scale. CLEAN flux densities of the target sources are 0.43-28~Jy, 0.32-21~Jy, 0.18-11~Jy, and 0.35-8.0~Jy in the 22, 43, 86, and 129~GHz bands, respectively.

The East-Asian VLBI Network (EAVN) is the international VLBI facility in East Asia and is conducted in collaboration with China, Japan, and Korea. The EAVN consists of VLBI arrays operated in each East Asian country, containing 21 radio telescopes and three correlators. The EAVN will be mainly operated at 6.7 (C-band), 8 (X-band), 22 (K-band), and 43 GHz (Q-band), although the EAVN has an ability to conduct observations at 1.6 - 129 GHz. We have conducted fringe test observations eight times to date at 8 and 22 GHz and fringes have been successfully detected at both frequencies. We have also conducted science commissioning observations of 6.7 GHz methanol masers in massive star-forming regions. The EAVN will be operational from the second half of 2017, providing complementary results with the FAST on AGNs, massive star-forming regions, and evolved stars with high angular resolution at cm- to mm-wavelengths.

We investigate shell emission associated with dying radio loud AGNs. First, based on our recent work by Ito et al. (2015), we describe the dynamical and spectral evolutions of shells after stopping the jet energy injection. We find that the shell emission overwhelms that of the radio lobes soon after stopping the jet energy injection because fresh electrons are continuously supplied into the shell via the forward shock while the radio lobes rapidly fade out without jet energy injection. We find that such fossil shells can be a new class of target sources for SKA telescope. Next, we apply the model to the nearby radio source 3C84. Then, we find that the fossil shell emission in 3C84 is less luminous in radio band while it is bright in TeV gamma-ray band and it can be detectable by CTA.

The Interferometric Monitoring of Gamma-ray Bright Active galactic nuclei (iMOGABA) program provides not only simultaneous multifrequency observations of bright gamma-ray detected active galactic nuclei (AGN), but also covers the highest Very Large Baseline Interferometry (VLBI) frequencies ever being systematically monitored, up to 129 GHz. However, observation and imaging of weak sources at the highest observed frequencies is very challenging. In the second paper in this series, we evaluate the viability of the frequency phase transfer technique to iMOGABA in order to obtain larger coherence time at the higher frequencies of this program (86 and 129 GHz) and image additional sources that were not detected using standard techniques. We find that this method is applicable to the iMOGABA program even under non-optimal weather conditions.

We present the detection of compact radio structures of fourteen radio-loud narrow line Seyfert 1 (NLS1) galaxies from Very Long Baseline Array observations at 5 GHz, which were performed in 2013. While 50\% of the sources of our sample show a compact core only, the remaining 50\% exhibit a core-jet structure. The measured brightness temperatures of the cores range from $10^{8.4}$ to $10^{11.4}$ K with a median value of $10^{10.1}$ K, indicating that the radio emission is from non-thermal jets, and that, likely, most sources are not strongly beamed, then implying a lower jet speed in these radio-loud NLS1 galaxies. In combination with archival data taken at multiple frequencies, we find that seven sources show flat or even inverted radio spectra, while steep spectra are revealed in the remaining seven objects. Although all these sources are very radio-loud with $R > 100$, their jet properties are diverse, in terms of their milli-arcsecond (mas) scale (pc scale) morphology and their overall radio spectral shape. The evidence for slow jet speeds (i.e., less relativistic jets), in combination with the low kinetic/radio power, may offer an explanation for the compact VLBA radio structure in most sources. The mildly relativistic jets in these high accretion rate systems are consistent with a scenario, where jets are accelerated from the hot corona above the disk by the magnetic field and the radiation force of the accretion disk. Alternatively, a low jet bulk velocity can be explained by low spin in the Blandford-Znajek mechanism.

We report results of the performance evaluation of a new hardware correlator in Korea, the Daejeon correlator, developed by the Korea Astronomy and Space Science Institute (KASI) and the National Astronomical Observatory of Japan (NAOJ). We conducted Very Long Baseline Interferometry (VLBI) observations at 22~GHz with the Korean VLBI Network (KVN) in Korea and the VLBI Exploration of Radio Astrometry (VERA) in Japan, and correlated the aquired data with the Daejeon correlator. For evaluating the performance of the new hardware correlator, we compared the correlation outputs from the Daejeon correlator for KVN observations with those from a software correlator, the Distributed FX (DiFX). We investigated the correlated flux densities and brightness distributions of extragalactic compact radio sources. The comparison of the two correlator outputs show that they are consistent with each other within $<8\%$, which is comparable with the amplitude calibration uncertainties of KVN observations at 22~GHz. We also found that the 8\% difference in flux density is caused mainly by (a) the difference in the way of fringe phase tracking between the DiFX software correlator and the Daejeon hardware correlator, and (b) an unusual pattern (a double-layer pattern) of the amplitude correlation output from the Daejeon correlator. The visibility amplitude loss by the double-layer pattern is as small as 3\%. We conclude that the new hardware correlator produces reasonable correlation outputs for continuum observations, which are consistent with the outputs from the DiFX software correlator.

We report finding kiloparsec-scale radio emissions aligned with parsec-scale jet structures in the narrow-line Seyfert 1 (NLS1) galaxy Mrk 1239 using the Very Large Array and the Very Long Baseline Array. Thus, this radio-quiet NLS1 has a jet-producing central engine driven by essentially the same mechanism as that of other radio-loud active galactic nuclei (AGNs). Most of the radio luminosity is concentrated within 100 parsecs and overall radio morphology looks edge-darkened; the estimated jet kinetic power is comparable to Fanaroff--Riley Type I radio galaxies. The conversion from accretion to jet power appears to be highly inefficient in this highly accreting low-mass black hole system compared with that in a low-luminosity AGN with similar radio power driven by a sub-Eddington, high-mass black hole. Thus, Mrk 1239 is a crucial probe to the unexplored parameter spaces of central engines for a jet formation.

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 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.

We made simultaneous single-dish and very long baseline interferometer (VLBI) observations of a narrow-line Seyfert 1 galaxy (NLS1) 1H 0323+342, showing gamma-ray activity revealed by Fermi/LAT observations. We found significant variation of the total flux density at 8 GHz on the time scale of one month by the single-dish monitoring. The total flux density varied by 5.5% in 32 days, which is comparable to the gamma-ray variability time scale, corresponding to the variability brightness temperature of $7.0 \times 10^{11}$ K. The source consists of central and southeastern components on the parsec (pc) scale. The flux of only the central component decreased in the same way as the total flux density, indicating that the short-term radio variability, and probably the gamma-ray emitting region, is associated with this component. From the VLBI observations we obtained the brightness temperatures of greater than $(5.2 \pm 0.3) \times 10^{10}$ K, and derived the equipartition Doppler factor of greater than 1.7, the variability Doppler factor of 2.2, and the 8 GHz radio power of $10^{24.6}$ W Hz$^{-1}$. Combining them we conclude that acceleration of radio jets and creation of high-energy particles are ongoing in the central engine, and that the apparent very radio-loud feature of the source is due to the Doppler-boosting effect, resulting in the intrinsic radio loudness to be an order of magnitude smaller than the observed values. We also conclude that the pc-scale jet represents recurrent activity from the the spectral fitting and the estimated kinematic age of pc- and kpc-scale extended components with different position angle.

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 conducted Very Long Baseline Array (VLBA) observations of seven nearby narrow-line Seyfert 1 (NLS1) galaxies at 1.7 GHz (18cm) with milli-arcsecond resolution. This is the first systematic very long baseline interferometry (VLBI) study focusing on the central parsec-scale regions of radio-quiet NLS1s. Five of the seven were detected at a brightness temperature of >~5x10^6 K and contain radio cores with high brightness temperatures of >6x10^7 K, indicating a nonthermal process driven by jet-producing central engines as is observed in radio-loud NLS1s and other active galactic nucleus (AGN) classes. VLBA images of MRK 1239, MRK 705, and MRK 766 exhibit parsec-scale jets with clear linear structures. A large portion of the radio power comes from diffuse emission components that are distributed within the nuclear regions (<~300 pc), which is a common characteristic throughout the observed NLS1s. Jet kinetic powers limited by the Eddington limit may be insufficient to allow the jets escape to kiloparsec scales for these radio-quiet NLS1s with low-mass black holes of <~10^7 M_sun.

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.

It seems to be a widely accepted opinion that the types of accretion disks (or flows) generally realized in the nuclei of radio galaxies and in further lower mass-accretion rate nuclei are inner, hot, optically thin, radiatively inefficient accretion flows (RIAFs) surrounded by outer, cool, optically thick, standard type accretion disks. However, observational evidence for the existence of such outer cool disks in these nuclei is rather poor. Instead, recent observations sometimes suggest the existence of inner cool disks of non-standard type, which develop in the region very close to their central black holes. Taking NGC 4261 as a typical example of such light eating nuclei, for which both flux data ranging from radio to X-ray and data for the counterjet occultation are available, we examine the plausibility of such a picture for the accretion states as mentioned above, based on model predictions. It is shown that the explanation of the gap seen in the counterjet emission in terms of the free-free absorption by an outer standard disk is unrealistic, and moreover, the existence itself of such an outer standard disk seems very implausible. Instead, the model of RIAF in an ordered magnetic field (so called resistive RIAF model) can well serve to explain the emission gap in terms of the synchrotron absorption, as well as to reproduce the observed features of the overall spectral energy distribution (SED). This model also predicts that the RIAF state starts directly from an interstellar hot gas phase at around the Bondi radius and terminates at the inner edge whose radius is about 100 times the Schwartzschild radii. Therefore, there is a good possibility for a cool disk to develop within this innermost region.

We conducted radio detection observations at 8.4 GHz for 22 radio-loud broad absorption line (BAL) quasars, selected from the Sloan Digital Sky Survey (SDSS) Third Data Release, by a very-long-baseline interferometry (VLBI) technique. The VLBI instrument we used was developed by the Optically ConnecTed Array for VLBI Exploration project (OCTAVE), which is operated as a subarray of the Japanese VLBI Network (JVN). We aimed at selecting BAL quasars with nonthermal jets suitable for measuring their orientation angles and ages by subsequent detailed VLBI imaging studies to evaluate two controversial issues of whether BAL quasars are viewed nearly edge-on, and of whether BAL quasars are in a short-lived evolutionary phase of quasar population. We detected 20 out of 22 sources using the OCTAVE baselines, implying brightness temperatures greater than 10^5 K, which presumably come from nonthermal jets. Hence, BAL outflows and nonthermal jets can be generated simultaneously in these central engines. We also found four inverted-spectrum sources, which are interpreted as Doppler-beamed, pole-on-viewed relativistic jet sources or young radio sources: single edge-on geometry cannot describe all BAL quasars. We discuss the implications of the OCTAVE observations for investigations for the orientation and evolutionary stage of BAL quasars.

We performed two types of radiation testing on high-speed LSI chips to test their suitability for use in wideband observations by the Japanese next space VLBI mission, VSOP-2. In the total ionization dose experiment we monitored autocorrelation spectra which were taken with irradiated LSI chips and the source current at intervals up to 1,000 hours from the ionization dose, but we could not see any change of these features for the chips irradiated with dose rates expected in the VSOP-2 mission. In the single event effect experiment, we monitored the cross correlation phase and power spectra between the data from radiated and non-radiated devices, and the source current during the irradiation of heavy-ions. We observed a few tens of single event upsets as discrete delay jumps for each LSI. We estimated the occurrence rate of single events in space as between once a few days to once a month. No single event latch-up was seen in any of the LSIs. These results show that the tested LSIs have sufficient tolerance to the environment for space VLBI observations.

We performed phase-reference very long baseline interferometry (VLBI) observations on five radio-loud narrow-line Seyfert 1 galaxies (NLS1s) at 8.4 GHz with the Japanese VLBI Network (JVN). Each of the five targets (RXS J08066+7248, RXS J16290+4007, RXS J16333+4718, RXS J16446+2619, and B3 1702+457) in milli-Jansky levels were detected and unresolved in milli-arcsecond resolutions, i.e., with brightness temperatures higher than 10^7 K. The nonthermal processes of active galactic nuclei (AGN) activity, rather than starbursts, are predominantly responsible for the radio emissions from these NLS1s. Out of the nine known radio-loud NLS1s, including the ones chosen for this study, we found that the four most radio-loud objects exclusively have inverted spectra. This suggests a possibility that these NLS1s are radio-loud due to Doppler beaming, which can apparently enhance both the radio power and the spectral frequency.

We observed the narrow-line quasar SDSS J094857.3+002225, which has the highest known radio loudness for a narrow-line Seyfert~1 galaxy (NLS1), at 1.7--15.4 GHz with the Very Long Baseline Array (VLBA). This is the first very-long-baseline interferometry (VLBI) investigation for a radio-loud NLS1. We independently found very high brightness temperatures from (1) its compactness in a VLBA image and (2) flux variation among the VLBA observation, our other observations with the VLBA, and the Very Large Array (VLA). A Doppler factor larger than 2.7--5.5 was required to meet an intrinsic limit of brightness temperature in the rest frame. This is evidence for highly relativistic nonthermal jets in an NLS1. We suggest that the Doppler factor is one of the most crucial parameters determining the radio loudness of NLS1s. The accretion disk of SDSS J094857.3+002225 is probably in the very high state, rather than the high/soft state, by analogy with X-ray binaries with strong radio outbursts and superluminal jets such as GRS 1915+105.

Korean VLBI Network (KVN) is the first dedicated mm-wavelength VLBI Network in East Asia and will be available from the middle of 2008. KVN consists of three stations and has the maximum observation frequency of 129 GHz with the maximum baseline length of 480 km. KVN has unique characteristics in the multifrequency, simultaneous observing system. By taking advantage of this we are considering various science topics, including not only maser emitting regions and young stellar objects in our galaxy, but also extragalactic objects. Construction of the first site is in progress. We are concurrently developing components, including receivers, data acquisition systems, and a correlator, and also arranging the international collaboration.

We have made a high-resolution VLBI observation of the gamma-ray loud quasar PKS 1622-297 with the HALCA spacecraft and ground radio telescopes at 5 GHz in 1998 February, almost three years after the source exhibited a spectacular GeV gamma-ray flare. The source shows an elongated structure toward the west on the parsec scale. The visibility data are well modeled by three distinct components; a bright core and two weaker jet components. Comparison with previous observations confirms that the jet components have an apparent superluminal motion up to 12.1 h^{-1}c, with the inner jet components having lower superluminal speeds. We apply the inverse Compton catastrophe model and derive a Doppler factor, \delta, of 2.45, which is somewhat lower than that of other gamma-ray loud active galactic nuclei (AGNs), suggesting the source was in a more quiescent phase at the epoch of our observation. As an alternative probe of the sub-parsec scale structure, we also present the results from multi-epoch ATCA total flux monitoring, which indicate the presence of persistent intraday variability consistent with refractive interstellar scintillation. We examine the gamma-ray emission mechanism in the light of these observations.

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.

We investigate whether a pc-scale jet of 3C345 is dominated by a normal plasma or an electron-positron plasma. We present a general condition that a jet component becomes optically thick for synchrotron self-absorption, by extending the method originally developed by Reynolds et al. The general condition gives a lower limit of the electron number density, with the aid of the surface brightness condition, which enables us to compute the magnetic field density. Comparing the lower limit with another independent constraint for the electron density that is deduced from the kinetic luminosity, we can distinguish the matter content. We apply the procedure to the five components of 3C345 (C2, C3, C4, C5, and C7) of which angular diameters and radio fluxes at the peak frequencies were obtainable from literature. Evaluating the representative values of Doppler beaming factors by their equipartition values, we find that all the five components are likely dominated by an electron-positron plasma. The conclusion does not depend on the lower cutoff energy of the power-law distribution of radiating particles.