We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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.

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 POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

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.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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 report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

The supergiant VX Sagittarii is a strong emitter of both H2O and SiO masers. However, previous VLBI observations have been performed separately, which makes it difficult to spatially trace the outward transfer of the material consecutively. Here we present the astrometrically registered, simultaneous maps of 22.2 GHz H2O and 43.1/42.8/86.2/129.3 GHz SiO masers toward VX Sagittarii. The H2O masers detected above the dust-forming layers have an asymmetric distribution. The multi-transition SiO masers are nearly circular ring, suggesting spherically symmetric wind within a few stellar radii. These results provide the clear evidence that the asymmetry in the outflow is enhanced after the smaller molecular gas clump transform into the inhomogeneous dust layers. The 129.3 GHz maser arises from the outermost region compared to that of 43.1/42.8/86.2 GHz SiO masers. The ring size of the 129.3 GHz maser is maximized around the optical maximum, suggesting that radiative pumping is dominant.

As part of the search for the "dark molecular gas (DMG)," we report on the results of HCO$^+$ J = 1-0 absorption observations toward nine bright extragalactic millimeter wave continuum sources. The extragalactic sources are at high Galactic latitudes (|b| > 10$\deg$) and seen at small extinction (E(B-V) $\lesssim$ 0.1 mag). We have detected the HCO$^+$ absorption lines toward two sources, B0838+133 and B2251+158. The absorption toward B2251+158 was previously reported, while the absorption toward B0838+133 is a new detection. We derive hydrogen column densities or their upper limits toward the nine sources from our observations and compare them to those expected from CO line emission and far-infrared dust continuum emission. Toward the seven sources with no HCO$^+$ detection, CO emission has not been detected, either. Thus the sight lines are likely to be filled with almost pure atomic gas. Toward the two sources with HCO$^+$ detection, CO emission has been also detected. Comparison of the H$_2$ column densities from HCO$^+$ absorption and CO emission suggests a non-negligible amount of DMG toward B0838+133.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We study the linear polarization of the radio cores of eight blazars simultaneously at 22, 43, and 86 GHz with observations obtained by the Korean VLBI Network (KVN) in three epochs between late 2016 and early 2017 in the frame of the Plasma-physics of Active Galactic Nuclei (PAGaN) project. We investigate the Faraday rotation measure (RM) of the cores; the RM is expected to increase with observing frequency if core positions depend on frequency due to synchrotron self-absorption. We find a systematic increase of RMs at higher observing frequencies in our targets. The RM--$\nu$ relations follow power-laws with indices distributed around 2, indicating conically expanding outflows serving as Faraday rotating media. Comparing our KVN data with contemporaneous optical polarization data from the Steward Observatory for a few sources, we find indication that the increase of RM with frequency saturates at frequencies of a few hundreds GHz. This suggests that blazar cores are physical structures rather than simple $\tau=1$ surfaces. A single region, e.g. a recollimation shock, might dominate the jet emission downstream of the jet launching region. We detect a sign change in the observed RMs of CTA 102 on a time scale of $\approx$1 month, which might be related to new superluminal components emerging from its core undergoing acceleration/deceleration and/or bending. We see indication for quasars having higher core RMs than BL Lac objects, which could be due to denser inflows/outflows in quasars.

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 present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

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.

We study the centimeter- to millimeter-wavelength synchrotron spectrum of the core of the radio galaxy M87 at $\lesssim0.8\,{\rm mas}~\sim110R_{s}$ spatial scales using four years of fully simultaneous, multi-frequency VLBI data obtained by the Korean VLBI Network (KVN). We find a core spectral index $\alpha$ of $\gtrsim-0.37$ ($S\propto \nu^{+\alpha}$) between 22GHz and 129GHz. By combining resolution-matched flux measurements from the Very Long Baseline Array (VLBA) at 15GHz and taking the Event Horizon Telescope (EHT) 230GHz core flux measurements in epochs 2009 and 2012 as lower limits, we find evidence of a nearly flat core spectrum across 15GHz and 129GHz, which could naturally connect the 230GHz VLBI core flux. The extremely flat spectrum is a strong indication that the jet base does not consist of a simple homogeneous plasma, but of inhomogeneous multi-energy components, with at least one component with the turn-over frequency $\gtrsim100$GHz. The spectral shape can be qualitatively explained if both the strongly (compact, optically thick at $>$100GHz) and the relatively weakly magnetized (more extended, optically thin at $<$100GHz) plasma components are colocated in the footprint of the relativistic jet.

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 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 analysis of near-infrared, high-resolution spectroscopy towards the Flat-spectrum YSO (Young Stellar Object) ESO H$\alpha$ 279a (1.5 solar mass) in the Serpens star forming region, at the distance of 429 pc. Using the Immersion GRating INfrared Spectrometer (IGRINS, R=45,000), we detect emission lines originating from the accretion channel flow, jet, and inner disk. Specifically, we identify hydrogen Brackett series recombination, [Fe II], [Fe III], [Fe IV], Ca I, Na I, H2, H2O and CO overtone emission lines. By modeling five bands of CO overtone emission lines, and the symmetric double-peaked line profile for Na I emission lines, we find that ESO H$\alpha$ 279a has an actively accreting Keplerian disk. From our Keplerian disk model, we find that Na I emission lines originate between 0.04 AU and 1.00 AU, while CO overtone emission lines are from the outer part of disk, in the range between 0.22 AU and 3.00 AU. It reveals that the neutral atomic Na gas is a good tracer of the innermost region of the actively accreting disk. We derive a mass accretion rate of 2-10x10^{-7} M_solar/yr from the measured Br_gamma emission luminosity of 1.78(+-0.31)x10^{31} erg/s.

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.

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 the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

This study aims to search for the existence of intraday variability (IDV) of BL Lac object S5 0716+714 at high radio frequencies for which the interstellar scintillation effect is not significant. Using the 21-meter radio telescope of the Korean VLBI Network (KVN), we present results of multi-epoch simultaneous dual-frequency radio observations. Single-dish observations of S5 0716+714 were simultaneously conducted at 21.7 GHz (K-band) and 42.4 GHz (Q-band), with a high cadence of 30-60 minute intervals.We observed four epochs between December 2009 and June 2010. Over the whole set of observation epochs, S5 0716+714 showed significant inter-month variations in flux density at both the K- and Q-bands, with modulation indices of approximately 19% for the K-band and approximately 36% for the Q-band. In all epochs, no clear intraday variability was detected at either frequency. The source shows monotonic flux density increase in epochs 1 and 3 and monotonic flux density decrease in epochs 2 and 4. In the flux density increasing phases, the flux densities at the Q-band increase more rapidly. In the decreasing phase, no significant flux density difference is seen at the two frequencies. The situation could be different close to flux density peaks that we did not witness in our observations. We find an inverted spectrum with mean spectral indices of -0.57+-0.13 in epoch 1 and -0.15+-0.11 in epoch 3. On the other hand, we find relatively steep indices of +0.24+-0.14 and +0.17+-0.18 in epochs 2 and 4, respectively. We conclude that the frequency dependence of the variability and the change of the spectral index are caused by source-intrinsic effects rather than by any extrinsic scintillation effect.

We present the catalog of the KVN Calibrator Survey (KVNCS). This first part of the KVNCS is a single dish radio survey conducted at 22 (K band) and 43 GHz (Q band) simultaneously using the Korean VLBI Network (KVN) from 2009 to 2011. A total 2045 sources selected from the VLBA Calibrator Survey (VCS) with an extrapolated flux density limit of 100 mJy at K band. The KVNCS contains 1533 sources in the K band with a flux density limit of 70 mJy and 553 sources in the Q band with a flux density limit of 120 mJy; it covers the whole sky down to $-32.^\circ5$ in declination. Five hundred thirteen sources were detected in the K and Q bands, simultaneously; $\sim76\%$ of them are flat-spectrum sources ($-0.5 \leq \alpha \leq 0.5$). From the flux--flux relationship, we anticipated that the most of the radiation of many of the sources comes from the compact components. Therefore, the sources listed in the KVNCS are strong candidates for high frequency VLBI calibrators.

Class I methanol masers are found to be good tracers of the interaction between outflows from massive young stellar objects with their surrounding media. Although polarization observations of Class II methanol masers have been able to provide information about magnetic fields close to the central (proto)stars, polarization observations of Class I methanol masers are rare, especially at 44 and 95 GHz. We present the results of linear polarization observations of 39 Class I methanol maser sources at 44 and 95 GHz. These two lines are observed simultaneously with one of the 21 m Korean VLBI Network telescopes in single-dish mode. Approximately 60% of the observed sources have fractional polarizations of a few percent in at least one transition. This is the first reported detection of linear polarization of the 44 GHz methanol maser. The two maser transitions show similar polarization properties, indicating that they trace similar magnetic environments, although the fraction of the linear polarization is slightly higher at 95 GHz. We discuss the association between the directions of polarization angles and outflows. We also discuss some targets having different polarization properties at both lines, including DR21(OH) and G82.58+0.20, which show the 90{\deg} polarization angle flip at 44 GHz.

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.

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.

We present the results of a broadband radio-to-GeV observing campaign organized to get a better understanding of the radiation processes responsible for the $\gamma$-ray flares observed in 3C 279. The total intensity and polarization observations of the source were carried out between December 28, 2013 and January 03, 2014 using the Fermi-LAT, Swift-XRT, Swift-UVOT, and KVN telescopes. A prominent flare observed in the optical/near-UV passbands was found to be correlated with a concurrent $\gamma$-ray flare at a confidence level $>$95$\%$, which suggests a co-spatial origin of the two. Moreover, the flaring activity in the two regimes was accompanied by no significant spectral variations. A peak in the X-ray light curve coincides with the peaks of the fractional polarization curves at 43 and 86 GHz radio bands. No prominent variation was noticed for the total intensity and the electric vector position angle (EVPA) observations at radio bands during this period. We noticed a possible hint of steepening of the radio spectrum with an increase in percentage polarization, which suggests that the radio polarization variations could be simply due to a spectral change. In a simple scenario, the correlated optical/$\gamma$-ray flares could be caused by the same population of emitting particles. The coincidence of the increase in radio polarization with the X-ray flux supports the picture that X-rays are produced via inverse-Compton scattering of radio photons. The observed fractional variability for the $\gamma$-ray flare $\sim$0.23 does not exceed that in the optical regime, which is inconsistent with what we usually observe for 3C 279; it could be due to different dependencies of the magnetic field and the external radiation field energy density profiles along the jet.

The calibration of Very Long Baseline Interferometry (VLBI) data has long been a time consuming process. The Korean VLBI Network (KVN) is a simple array consisting of three identical antennas. Because four frequencies are observed simultaneously, phase solutions can be transferred from lower frequencies to higher frequencies in order to improve phase coherence and hence sensitivity at higher frequencies. Due to the homogeneous nature of the array, the KVN is also well suited for automatic calibration. In this paper we describe the automatic calibration of single-polarisation KVN data using the KVN Pipeline and comparing the results against VLBI data that has been manually reduced. We find that the pipelined data using phase transfer produces better results than a manually reduced dataset not using the phase transfer. Additionally we compared the pipeline results with a manually reduced phase-transferred dataset and found the results to be identical.

We present results of long-term multi-wavelength polarization observations of the powerful blazar 3C~279 after its $\gamma$-ray flare on 2013~December 20. We followed up this flare with single-dish polarization observations using two 21-m telescopes of the Korean VLBI Network. Observations carried out weekly from 2013~December~25 to 2015~January~11, at 22~GHz, 43~GHz, 86~GHz simultaneously, as part of the Monitoring Of GAmma-ray Bright AGN (MOGABA) program. We measured 3C~279 total flux densities of 22--34~Jy at 22~GHz, 15--28~Jy (43~GHz), and 10--21~Jy (86~GHz), showing mild variability of $\leq 50\,\%$ over the period of our observations. The spectral index between 22~GHz and 86~GHz ranged from $-0.13$ to $-0.36$. Linear polarization angles were 27$^{\circ}$--38$^{\circ}$, 30$^{\circ}$--42$^{\circ}$, and 33$^{\circ}$--50$^{\circ}$ at 22~GHz, 43~GHz, and 86~GHz, respectively. The degree of linear polarization was in the range of 6--12\,\%, and slightly decreased with time at all frequencies. We investigated Faraday rotation and depolarization of the polarized emission at 22--86~GHz, and found Faraday rotation measures (RM) of $-300$ to $-1200$~rad~m$^{-2}$ between 22~GHz and 43~GHz, and $-800$ to $-5100$~rad~m$^{-2}$ between 43~GHz and 86~GHz. The RM values follow a power law with a mean power law index $a$ of $2.2$, implying that the polarized emission at these frequencies travels through a Faraday screen in or near the jet. We conclude that the regions emitting polarized radio emission may be different from the region responsible for the 2013 December $\gamma$-ray flare and are maintained by the dominant magnetic field perpendicular to the direction of the radio jet at milliarcsecond scales.

Active galactic nuclei (AGN) are known for irregular variability on all time scales, down to intra-day variability with relative variations of a few percent within minutes to hours. On such short timescales, unexplored territory, such as the possible existence of a shortest characteristic time scale of activity and the shape of the high frequency end of AGN power spectra, still exists. We present the results of AGN single-dish fast photometry performed with the Korean VLBI Network (KVN). Observations were done in a "anti-correlated" mode using two antennas, with always at least one antenna pointing at the target. This results in an effective time resolution of less than three minutes. We used all four KVN frequencies, 22, 43, 86, and 129 GHz, in order to trace spectral variability, if any. We were able to derive high-quality light curves for 3C 111, 3C 454.3, and BL Lacertae at 22 and 43 GHz, and for 3C 279 at 86 GHz, between May 2012 and April 2013. We performed a detailed statistical analysis in order to assess the levels of variability and the corresponding upper limits. We found upper limits on flux variability ranging from $\sim$1.6% to $\sim$7.6%. The upper limits on the derived brightness temperatures exceed the inverse Compton limit by three to six orders of magnitude. From our results, plus comparison with data obtained by the University of Michigan Radio Astronomy Observatory, we conclude that we have not detected source-intrinsic variability which would have to occur at sub-per cent levels.

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.

CIT 6 is a carbon star in the transitional phase from the asymptotic giant branch (AGB) to the protoplanetary nebulae (pPN). Observational evidences of two point sources in the optical, circumstellar arc segments in an HC$_3$N line emission, and a bipolar nebula in near-infrared provide strong support for the presence of a binary companion. Hence, CIT 6 is very attractive for studying the role of companions in the AGB-pPN transition. We have carried out high resolution $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ observations of CIT 6 with the Submillimeter Array combined with the Submillimeter Telescope (single-dish) data. The $^{12}$CO channel maps reveal a spiral-shell pattern connecting the HC$_3$N segments in a continuous form, and an asymmetric outflow corresponding to the near-infrared bipolar nebula. Rotation of the $^{12}$CO channel peak position may be related to the inner spiral winding and/or the bipolar outflow. An eccentric orbit binary is suggested for the presences of an anisotropic mass loss to the west and a double spiral pattern. The lack of interarm emission to the west may indicate a feature corresponding to the periastron passage of a highly eccentric orbit of the binary. Spatially-averaged radial and spectral profiles of $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ are compared with simple spherical radiative transfer models, suggesting a change of $^{12}$CO/$^{13}$CO abundance ratio from $\sim30$ to $\sim50$ inward in the CSE of CIT 6. The millimeter continuum emission is decomposed into extended dust thermal emission (spectral index $\sim-2.4$) and compact emission from radio photosphere (spectral index $\sim-2.0$).

We present a simultaneous single-dish survey of 22 GHz water maser and 44 GHz and 95 GHz class I methanol masers toward 77 6.7 GHz class II methanol maser sources, which were selected from the Arecibo methanol maser Galactic plane survey (AMGPS) catalog.Water maser emission is detected in 39 (51%) sources, of which 15 are new detections. Methanol maser emission at 44 GHz and 95 GHz is found in 25 (32%) and 19 (25%) sources, of which 21 and 13 sources are newly detected, respectively. We find 4 high-velocity (> 30 km/s) water maser sources, including 3 dominant blue- or redshifted outflows.The 95 GHz masers always appear with the 44 GHz maser emission. They are strongly correlated with 44 GHz masers in velocity, flux density, and luminosity, while they are not correlated with either water or 6.7 GHz class II methanol masers. The average peak flux density ratio of 95 GHz to 44 GHz masers is close to unity, which is two times higher than previous estimates. The flux densities of class I methanol masers are more closely correlated with the associated BGPS core mass than those of water or class II methanol masers. Using the large velocity gradient (LVG) model and assuming unsaturated class I methanol maser emission, we derive the fractional abundance of methanol to be in a range of 4.2*10^-8 to 2.3*10^-6, with a median value of 3.3\pm2.7*10^-7.

We report results of investigation of amplitude calibration for very long baseline interferometry (VLBI) observations with Korean VLBI Network (KVN). Amplitude correction factors are estimated based on comparison of KVN observations at 22~GHz correlated by Daejeon hardware correlator and DiFX software correlator in Korea Astronomy and Space Science Institute (KASI) with Very Long Baseline Array (VLBA) observations at 22~GHz by DiFX software correlator in National Radio Astronomy Observatory (NRAO). We used the observations for compact radio sources, 3C 454.3, NRAO 512, OJ 287, BL Lac, 3C 279, 1633+382, and 1510-089, which are almost unresolved for baselines in a range of 350-477 km. Visibility data of the sources obtained with similar baselines at KVN and VLBA are selected, fringe-fitted, calibrated, and compared for their amplitudes. We found that visibility amplitudes of KVN observations should be corrected by factors of 1.10 and 1.35 when correlated by DiFX and Daejeon correlators, respectively. These correction factors are attributed to the combination of two steps of 2-bit quantization in KVN observing systems and characteristics of Daejeon correlator.

We report the first detection of linearly polarized emission at an observing wavelength of 350 $\mu$m from the radio-loud Active Galactic Nucleus 3C 279. We conducted polarization observations for 3C 279 using the SHARP polarimeter in the Caltech Submillimeter Observatory on 2014 March 13 and 14. For the first time, we detected the linear polarization with the degree of polarization of 13.3\%$\pm$3.4\% ($3.9\sigma$) and the Electric Vector Position Angle (EVPA) of 34.7$^\circ\pm5.6^\circ$. We also observed 3C~279 simultaneously at 13, 7, and 3.5 mm in dual polarization with the Korean very long baseline interferometry (VLBI) Network on 2014 March 6 (single dish) and imaged in milliarcsecond (mas) scales at 13, 7, 3.5, and 2.3 mm on March 22 (VLBI). We found that the degree of linear polarization increases from 10\% to 13\% at 13 mm to 350 $\mu$m and the EVPAs at all observing frequencies are parallel within $<10^\circ$ to the direction of the jet at mas scale, implying that the integrated magnetic fields are perpendicular to the jet in the innermost regions. We also found that the Faraday rotation measures RM are in a range of $-6.5\times10^2 \sim -2.7\times10^3$ rad m$^{-2}$ between 13 and 3.5 mm, and are scaled as a function of wavelength: $|{\rm RM}|\propto\lambda^{-2.2}$. These results indicate that the mm and sub-mm polarization emission are generated in the compact jet within 1~mas scale and affected by a Faraday screen in or in the close proximity of the jet.

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.

Oxygen-rich Asymptotic Giant Branch (AGB) stars can be intense emitters of SiO ($v$=1 and 2, J=1$\rightarrow$0) and H$_2$O maser lines at 43 and 22 GHz, respectively. VLBI observations of the maser emission provide a unique tool to probe the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned \water\ and $v$=1 and $v$=2 SiO maser maps have traditionally limited the physical constraints that can be placed on the SiO maser pumping mechanism. We present phase referenced simultaneous spectral-line VLBI images for the SiO $v$=1 and $v$=2, J=1$\rightarrow$0, and H$_2$O maser emission around the AGB star R\,LMi, obtained from the Korean VLBI Network (KVN). The simultaneous multi-channel receivers of the KVN offer great possibilities for astrometry in the frequency domain. With this facility we have produced images with bona-fide absolute astrometric registration between high frequency maser transitions of different species to provide the positions of the \water\ maser emission, and the centre of the SiO maser emission, and hence reducing the uncertainty in the proper motion for R\,LMi by an order of magnitude over that from Hipparcos. This is the first successful demonstration of source frequency phase referencing for mm-VLBI spectral-line observations and also where the ratio between the frequencies is not an integer.

We present the results of unbiased 22 GHz H2O water and 44 GHz class I CH3OH methanol maser surveys in the central 7x10 arcmin area of NGC 1333 and two additional mapping observations of a 22 GHz water maser in a ~3x3arcmin area of the IRAS4A region. In the 22 GHz water maser survey of NGC 1333 with sensitivity of sigma~0.3Jy, we confirmed masers toward H2O(B) in the region of HH 7-11 and IRAS4B. We also detected new water masers at ~20arcsec away in the western direction of IRAS4B or ~25arcsec away in the southern direction of IRAS4A. We could not however find young stellar objects or molecular outflows associated with them. They showed two different velocity components of ~0 and ~16 km/s, which are blue- and red-shifted relative to the adopted systemic velocity of ~7 km/s for NGC 1333. They also showed time variabilities in both intensity and velocity from multi-epoch observations and an anti-correlation between the intensities of the blue- and the red-shifted velocity components. We suggest that the unidentified powering source of these masers might be in the earliest evolutionary stage of star formation before the onset of molecular outflows. Finding this kind of water masers is only possible by an unbiased blind survey. In the 44 GHz methanol maser survey with sensitivity of sigma~0.5 Jy, we confirmed masers toward the IRAS4A2 and the eastern shock region of the IRAS2A. Both sources are also detected in 95 and 132 GHz methanol maser lines. In addition, we had new detections of methanol masers at 95 and 132 GHz toward IRAS4B. In terms of the isotropic luminosity, we detected the methanol maser sources brighter than ~5x1025 erg/s from our unbiased survey.

The Korean VLBI Network (KVN) is a new mm-VLBI dedicated array with capability for simultaneous observations at multiple frequencies, up to 129 GHz. The innovative multi-channel receivers present significant benefits for astrometric measurements in the frequency domain. The aim of this work is to verify the astrometric performance of the KVN using a comparative study with the VLBA, a well established instrument. For that purpose, we carried out nearly contemporaneous observations with the KVN and the VLBA, at 14/7 mm, in April 2013. The KVN observations consisted of simultaneous dual frequency observations, while the VLBA used fast frequency switching observations. We used the Source Frequency Phase Referencing technique for the observational and analysis strategy. We find that having simultaneous observations results in a superior performance for compensation of all atmospheric terms in the observables, in addition to offering other significant benefits for astrometric analysis. We have compared the KVN astrometry measurements to those from the VLBA. We find that the structure blending effects introduce dominant systematic astrometric shifts and these need to be taken into account. We have tested multiple analytical routes to characterize the impact of the low resolution effects for extended sources in the astrometric measurements. The results from the analysis of KVN and full VLBA datasets agree within 2-$\sigma$ of the thermal error estimate. We interpret the discrepancy as arising from the different resolutions. We find that the KVN provides astrometric results with excellent agreement, within 1-$\sigma$, when compared to a VLBA configuration which has a similar resolution. Therefore this comparative study verifies the astrometric performance of KVN using SFPR at 14/7 mm, and validates the KVN as an astrometric instrument.

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 started two observing programs with the Korean VLBI Network (KVN) monitoring changes in the flux density and polarization of relativistic jets in gamma-ray bright AGNs simultaneously at 22, 43, 86, 129 GHz. One is a single-dish weekly-observing program in dual polarization with KVN 21-m diameter radio telescopes beginning in 2011 May. The other is a VLBI monthly-observing program with the three-element VLBI network at an angular resolution range of 1.0--9.2 mas beginning in 2012 December. The monitoring observations aim to study correlation of variability in gamma-ray with that in radio flux density and polarization of relativistic jets when they flare up. These observations enable us to study the origin of the gamma-ray flares of AGNs.

We present the results of CO emission line observations toward the semicircular Galactic supernova remnant (SNR) 3C434.1 (G94.0+1.0). We mapped an area covering the whole SNR in the 12CO J=1-0 emission line using the Seoul Radio Astronomy Observatory (SRAO) 6-m telescope, and found a large molecular cloud superposed on the faint western part of the SNR. The cloud was elongated along the north-south direction and showed a very good spatial correlation with the radio features of the SNR. We carried out 12CO J=2-1 line observations of this cloud using the K\"{o}lner Observatorium f\"{u}r Sub-Millimeter Astronomie (KOSMA) 3-m telescope and found a region in which the 12CO J=2-1 to J=1-0 ratio was high (~ 1.6). This higher excitation, together with the morphological relation, strongly suggested that the molecular cloud was interacting with the SNR. The systemic velocity of the molecular cloud (-13 km/s) gave a kinematic distance of 3.0 kpc to the SNR-molecular cloud system. We derived the physical parameters of the SNR based on this new distance. We examined the variation of the radio spectral index over the remnant and found that it was flatter in the western part, wherein the SNR was interacting with the molecular cloud. We therefore propose that 3C434.1 is the remnant of a supernova explosion that occurred just outside the boundary of a relatively thin, sheetlike molecular cloud. We present a hydrodynamic model showing that its asymmetric radio morphology can result from its interaction with this blocking molecular cloud.

This paper presents the catalog of correlated flux densities in three ranges of baseline projection lengths of 637 sources from a 43 GHz (Q-band) survey observed with the Korean VLBI Network. Of them, 623 sources have not been observed before at Q-band with VLBI. The goal of this work in the early science phase of the new VLBI array is twofold: to evaluate the performance of the new instrument that operates in a frequency range of 22-129 GHz and to build a list of objects that can be used as targets and as calibrators. We have observed the list of 799 target sources with declinations down to -40 degrees. Among them, 724 were observed before with VLBI at 22 GHz and had correlated flux densities greater than 200 mJy. The overall detection rate is 78%. The detection limit, defined as the minimum flux density for a source to be detected with 90% probability in a single observation, was in a range of 115-180 mJy depending on declination. However, some sources as weak as 70 mJy have been detected. Of 623 detected sources, 33 objects are detected for the first time in VLBI mode. We determined their coordinates with the median formal uncertainty 20 mas. The results of this work set the basis for future efforts to build the complete flux-limited sample of extragalactic sources at frequencies 22 GHz and higher at 3/4 of the celestial sphere.

The NGC 2024 FIR 6 region was observed in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. The water maser spectra displayed several velocity components and month-scale time variabilities. Most of the velocity components may be associated with FIR 6n, while one component was associated with FIR 4. A typical life time of the water-maser velocity-components is about 8 months. The components showed velocity fluctuations with a typical drift rate of about 0.01 km/s/day. The methanol class I masers were detected toward FIR 6. The methanol emission is confined within a narrow range around the systemic velocity of the FIR 6 cloud core. The methanol masers suggest the existence of shocks driven by either the expanding H II region of FIR 6c or the outflow of FIR 6n.

We present the results of 12CO J = 1-0 line observations of eleven Galactic supernova remnants (SNRs) obtained using the Seoul Radio Astronomy Observatory (SRAO) 6-m radio telescope. The observation was made as a part of the SRAO CO survey of SNRs between l = 70 and 190 deg, which is intended to identify SNRs interacting with molecular clouds. The mapping areas for the individual SNRs are determined to cover their full extent in the radio continuum. We used halfbeam grid spacing (60") for 9 SNRs and full-beam grid spacing (120") for the rest. We detected CO emission towards most of the remnants. In six SNRs, molecular clouds showed a good spatial relation with their radio morphology, although no direct evidence for the interaction was detected. Two SNRs are particularly interesting: G85.4+0.7, where there is a filamentary molecular cloud along the radio shell, and 3C434.1, where a large molecular cloud appears to block the western half of the remnant. We briefly summarize the results obtained for individual SNRs.

We report simultaneous multi-frequency observing performance at 22 and 43 GHz of the 21-m shaped-Cassegrain radio telescopes of the Korean VLBI Network (KVN). KVN is the first millimeter-dedicated VLBI network in Korea having a maximum baseline length of 480 km. It currently operates at 22 and 43 GHz and planed to operate in four frequency bands, 22, 43, 86, and 129 GHz. The unique quasioptics of KVN enable simultaneous multi-frequency observations based on efficient beam filtering and accuarate antenna-beam alignment at 22 and 43 GHz. We found that the offset of the beams is within <5 arcseconds over all pointing directions of antenna. The dual polarization, cooled HEMT receivers at 22 and 43 GHz result in receiver noise temperatures less than 40 K at 21.25-23.25 GHz and 80 K at 42.11-44.11 GHz. The pointing accuracies have been measured to be 3 arcseconds in azimuth and elevation for all antennas. The measured aperture efficiencies are 65%(K)/67%(Q), 62%(K)/59%(Q), and 66%(K)/60%(Q) for the three KVN antennas, KVNYS, KVNUS, and KVNTN, respectively. The main-beam efficiencies are measured to be 50%(K)/52%(Q), 48%(K)/50%(Q), and 50%(K)/47%(Q) for KVNYS, KVNUS, and KVNTN, respectively. The estimated Moon efficiencies are 77%(K)/90%(Q), 74%(K)/79%(Q), and 80%(K)/86%(Q) for KVNYS, KVNUS, KVNTN, respectively. The elevation dependence of the aperture efficiencies is quite flat for elevations > 20 degrees.

We report a multi-epoch, simultaneous 22 GHz H2O and 44 GHz class I CH3OH maser line survey towards 180 intermediate-mass young stellar objects, including 14 Class 0, 19 Class I objects, and 147 Herbig Ae/Be stars. We detected H2O and CH3OH maser emission towards 16 (9%) and 10 (6%) sources with one new H2O and six new CH3OH maser sources. The detection rates of both masers rapidly decrease as the central (proto)stars evolve, which is contrary to the trends in high-mass star-forming regions. This suggests that the excitations of the two masers are closely related to the evolutionary stage of the central (proto)stars and the circumstellar environments. H2O maser velocities deviate on average 9 km s^-1 from the ambient gas velocities whereas CH3OH maser velocities match quite well with the ambient gas velocities. For both maser emissions, large velocity differences (|v_{H2O} - v_{sys} | > 10 km s^-1 and |v_{CH3OH} - v_{sys}| > 1 km s^-1) are mostly confined to Class 0 objects. The formation and disappearance of H2O masers is frequent and their integrated intensities change by up to two orders of magnitude. In contrast, CH3OH maser lines usually show no significant change in intensity, shape, or velocity. This is consistent with the previous suggestion that H2O maser emission originates from the base of an outflow while 44 GHz class I CH3OH maser emission arises from the interaction region of the outflow with the ambient gas. The isotropic maser luminosities are well correlated with the bolometric luminosities of the central objects. The fitted relations are L_{H2O} = 1.71 * 10^{-9} (L_{bol})^{0.97} and L_{CH3OH} = 1.71 * 10^{-10} (L_{bol})^{1.22}.