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 study the properties of the innermost jet of the flat spectrum radio quasar 1633+382 (4C~38.41) based on VLBI data from the radio monitoring observations of the Boston University VLBI program at 43~GHz. Analysis of the components suggests a semi-parabolic jet geometry with jet radius $R$ following the relation $R\propto r^{0.7}$ with distance $r$, with indications of a jet geometry break towards a conical geometry. Brightness temperature falls with distance following $T_B\propto r^{-2.1}$. Combining this information, magnetic field and electron densities are found to fall along the jet as $B\propto r^{-1.5}$ and $n\propto r^{-1.1}$ respectively, suggesting that the magnetic configuration in the jet may be dominated by the poloidal component. Our analysis of the jet structure suggests that the innermost jet regions do not follow a ballistic trajectory and, instead, match a sinusoidal morphology which could be due to jet precession from a helical pattern or Kelvin-Helmholtz instabilities.

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.

In the current paradigm, it is believed that the compact VLBI radio core of radio-loud active galactic nuclei (AGNi) represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with a roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral distribution properties and it has been recently found that, at least in some nearby radio galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNi. Data combined from the sources considered here are not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNi needs to be re-examined by millimeter/sub-millimeter VLBI observations.

PKS 1510--089 is a bright and active $\gamma$-ray source that showed strong and complex $\gamma$-ray flares in mid-2015 during which the Major Atmospheric Gamma Imaging Cherenkov telescopes detected variable very high energy (VHE; photon energies $>$100 GeV) emission. We present long-term multi-frequency radio, optical, and $\gamma$-ray light curves of PKS 1510--089 from 2013 to 2018, and results of an analysis of the jet kinematics and linear polarization using 43 GHz Very Long Baseline Array data observed between late 2015 and mid-2017. We find that a strong radio flare trails the $\gamma$-ray flares in 2015, showing an optically thick spectrum at the beginning and becoming optically thin over time. Two laterally separated knots of emission are observed to emerge from the radio core nearly simultaneously during the $\gamma$-ray flares. We detect an edge-brightened linear polarization near the core in the active jet state in 2016, similar to the quiescent jet state in 2008--2013. These observations indicate that the $\gamma$-ray flares may originate from compression of the knots by a standing shock in the core and the jet might consist of multiple complex layers showing time-dependent behavior, rather than of a simple structure of a fast jet spine and a slow jet sheath.

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.

A long standing problem in the study of Active Galactic Nuclei (AGNs) is that the observed VLBI core is in fact a blending of the actual AGN core (classically defined by the $\tau=1$ surface) and the upstream regions of the jet or optically thin emitting region flows. This blending may cause some biases towards the observables of the core, such as its flux density, size or brightness temperature, which may lead to misleading interpretation of the derived quantities and physics. We study the effects of such blending under the view of the Korean VLBI Network (KVN) for a sample of AGNs at 43 GHz by comparing their observed properties with observations with the Very Large Baseline Array (VLBA). Our results suggest that the observed core sizes are a factor $\sim11$ larger than these of VLBA, which is similar to the factor expected by considering the different resolutions of the two facilities. We suggest the use of this factor to consider blending effects in KVN measurements. Other parameters, such as flux density or brightness temperature, seem to possess a more complicated dependence.

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

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

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

We present the results of very long baseline interferometry (VLBI) observations of gamma-ray bright blazar S5 0716+714 using the Korean VLBI Network (KVN) at the 22, 43, 86, and 129 GHz bands, as part of the Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA) KVN key science program. Observations were conducted in 29 sessions from January 16, 2013 to March 1, 2016, with the source being detected and imaged at all available frequencies. In all epochs, the source was compact on the milliarcsecond (mas) scale, yielding a compact VLBI core dominating the synchrotron emission on these scales. Based on the multi-wavelength data between 15 GHz (Owens Valley Radio Observatory) and 230 GHz (Submillimeter Array), we found that the source shows multiple prominent enhancements of the flux density at the centimeter (cm) and millimeter (mm) wavelengths, with mm enhancements leading cm enhancements by -16$\pm$8 days. The turnover frequency was found to vary between 21 to 69GHz during our observations. By assuming a synchrotron self-absorption model for the relativistic jet emission in S5 0716+714, we found the magnetic field strength in the mas emission region to be $\le$5 mG during the observing period, yielding a weighted mean of 1.0$\pm$0.6 mG for higher turnover frequencies (e.g., >45 GHz).

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.

The structure of the NGC 6251 jet on the milliarcsecond scale is investigated using images taken with the European VLBI Network and the Very Long Baseline Array. We detect a structural transition of the jet from a parabolic to a conical shape at a distance of (1-2) x 10^5 times the Schwarzschild radius from the central engine, which is close to the sphere of gravitational influence of the supermassive black hole (SMBH). We also examine the jet pressure profiles with the synchrotron minimum energy assumption to discuss the physical origin of the structural transition. The NGC 6251 jet, together with the M87 jet, suggests a fundamental process of structural transition in the jets of active galactic nuclei (AGNs). Collimated AGN jets are characterized by their external galactic medium, showing that AGN jets interplay with the SMBH and its host galaxy.

The mis-aligned active galactic nuclei (AGN), 3C\,84, has been observed approximately monthly at four frequencies simultaneously (22\,GHz -- 129\,GHz) on the Korean VLBI Network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C\,84 is known to have to sites of bright radio emission, with the first thought to be near the central super-massive black hole and the second in a slowly moving emission feature, several milliarcseconds south of the core. Analysis of this data suggests that the highly variable Gamma-ray emission originates near the SMBH and the slowly rising $\gamma$-ray emission originates in the slow moving feature.

In the current paradigm, it is believed that the compact VLBI radio core of radio-loud AGNs represents the innermost upstream regions of relativistic outflows. These regions of AGN jets have generally been modeled by a conical outflow with roughly constant opening angle and flow speed. Nonetheless, some works suggest that a parabolic geometry would be more appropriate to fit the high energy spectral energy distribution properties and it has been recently found that, at least in some nearby radio-galaxies, the geometry of the innermost regions of the jet is parabolic. We compile here multi-frequency core sizes of archival data to investigate the typically unresolved upstream regions of the jet geometry of a sample of 56 radio-loud AGNs. Data combined from the sources considered here is not consistent with the classic picture of a conical jet starting in the vicinity of the super-massive black hole (SMBH), and may exclude a pure parabolic outflow solution, but rather suggest an intermediate solution with quasi-parabolic streams, which are frequently seen in numerical simulations. Inspection of the large opening angles near the SMBH and the range of the Lorentz factors derived from our results support our analyses. Our result suggests that the conical jet paradigm in AGNs needs to be re-examined by mm/sub-mm VLBI observations.

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 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 investigate the magnetic field configuration of the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. By stacking images over three years in order to enhance the sensitivity, we reveal, for the first time, systematic transverse gradients of the Faraday rotation measure in several knots along the jet. Combining this result with polarization properties and the dynamics of the jet, we suggest the magnetic structure in several knots at kiloparsec scales consists of a systematically wrapped, tightly wound helical configuration. Our analysis brings us a new paradigm where the M87 jet is a fundamentally current carrying system produced in the vicinity of the supermassive black hole, transferring a huge amount of the electromagnetic energy over the host galaxy scale.

We investigate the distribution of Faraday rotation measure (RM) in the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. We resolve the structure of the RM in several knots along the jet for the first time. We derive the power spectrum in the arcsecond scale jet and find indications that the RM cannot be associated with a turbulent magnetic field with 3D Kolmogorov spectrum. Our analysis indicates that the RM probed on jet scales has a significant contribution of a Faraday screen associated with the vicinity of the jet, in contrast with that on kiloparsec scales, typically assumed to be disconnected from the jet. Comparison with previous RM analyses suggests that the magnetic fields giving rise to the RMs observed in jet scales have different properties and are well less turbulent than these observed in the lobes.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

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 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 on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($\chi$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $\chi$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the M 87 core to be between $-$7.5$\times$10$^{5}$ and 3.4$\times$10$^{5}$ rad/m$^{2}$. Assuming a density profile of the accretion flow that follows a power-law distribution and a magnetic field that is ordered, radial, and has equipartition strength, the limit on the rotation measure constrains the mass accretion rate $\dot{M}$ to be below 9.2$\times$10$^{-4}$ M$_{\odot}$~yr$^{-1}$ at a distance of 21 Schwarzchild radii from the central black hole. This value is at least two orders of magnitude smaller than the Bondi accretion rate, suggesting significant suppression of the accretion rate in the inner region of the accretion flow. Consequently, our result disfavors the classical \emph{advection dominated accretion flow} (ADAF) and prefers the \emph{adiabatic inflow-outflow solution} (ADIOS) or \emph{convection-dominated accretion flow} (CDAF) for the hot accretion flow in M 87.

The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $\mu$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $\mu$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016.

We report results from 5-day VLBI observations of the well-known quasar 3C 279 at 1.3 mm (230 GHz) in 2011. The measured nonzero closure phases on triangles including stations in Arizona, California and Hawaii indicate that the source structure is spatially resolved. We find an unusual inner jet direction at scales of $\sim$1 parsec extending along the northwest-southeast direction (PA = $127^{\circ}\pm3^{\circ}$), as opposed to other (previously) reported measurements on scales of a few parsecs showing inner jet direction extending to the southwest. The 1.3 mm structure corresponds closely with that observed in the central region of quasi-simultaneous super-resolution VLBA images at 7 mm. The closure phase changed significantly on the last day when compared with the rest of observations, indicating that the inner jet structure may be variable on daily timescales. The observed new direction of the inner jet shows inconsistency with the prediction of a class of jet precession models. Our observations indicate a brightness temperature of $\sim 8\times10^{10}$ K in the 1.3 mm core, much lower than that at centimeter wavelengths. Observations with better uv coverage and sensitivity in the coming years will allow the discrimination between different structure models and will provide direct images of the inner regions of the jet with 20--30 $\mu$as (5--7 light months) resolution.

We have studied Very Long Baseline Array (VLBA) polarimetric observations of 8 sources including quasars and BL Lacs at 12, 15, 22, 24 and 43 GHz and high frequency rotation measure ($RM$) maps are presented. We find typical values for the $RM$ in the VLBI core of several thousand rad/m$^2$, which are higher than values in the literature at lower frequencies. Assuming a dependence of the form $RM\propto \nu^a$, we obtain an average value of $a=3.6\pm1.3$, which is larger than that expected by theoretical considerations. Rotation measures are detected in the jet of only two sources and we find that only 0906+430 (and possibly 1633+382) show indications of a robust gradient. We discuss the Faraday--corrected polarization properties of the sources. Our interpretation supports the presence of helical magnetic fields with new, unresolved, components affecting the intrinsic direction of polarization close to the base of the jet of some objects.

One of the tools that can provide evidence about the existence of helical magnetic fields in AGN is the observation of rotation measure gradients across the jet. Such observations have been previously made successfully, proving that such gradients are far from being rare, but common and typically persistent over several years, although some of them may show a reversal in the direction along the jet. Further studies of rotation measure gradients can help us in our understanding of the magnetic field properties and structure in the base of the jets. We studied Very Long Baseline Array (VLBA) polarimetric observations of 8 sources consistent of some quasars and BL Lacs at 12, 15, 22, 24 and 43 GHz and we find that all but two sources show indications of rotation measure gradients, either parallel or perpendicular to the jet. We interpret gradients perpendicular to the jet as indications of the change of the line of sight of the magnetic field due to its helicity, and gradients parallel to the jet as the decrease of magnetic field strength and/or electron density as we move along the jet. When comparing our results with the literature, we find tentative evidence of a rotation measure gradient flip, which can be explained as a change of the pitch angle or jet bending.

Previous research showed most BL Lac objects and some quasars have aligned VLBI-core and optical polarizations, although some of the AGNs also showed no obvious relationship between VLBI-core and optical polarization angles. This may indicate that some AGNs have co-spatial regions of optical and radio emission, while others do not. Another possibility is that some of the VLBI cores had Faraday rotations of several tens of thousand of rad/m^2, which were not properly fit using the three-frequency data due to n\pi ambiguities in the observed polarization angles, leading to incorrect subtraction of the effects of the core Faraday rotation, and so incorrect intrinsic radio polarization angles \chi_0. We obtained additional 12+15+22+24+43GHz plus optical observations for 8 of 40 AGNs previously considered. Our results indicate that, although some VLBI radio cores have comparatively high rotation measures, this alone cannot explain the misalignments found between the radio core and optical VLBI polarization angles, \Delta\chi = |\chi_opt - \chi_0|. Comparison between \Delta\chi and (i) the orientation of \chi_0 relative to the jet direction, (ii) the degree of polarization of the core, (iii) a depolarization factor, (iv) the core rotation measures and (v) the core magnetic fields 1pc from the jet base do not yield evidence for any correlations between these properties. There is some evidence that the maximum observed \Delta\chi tends to decrease as the core-region magnetic field increases, suggesting that misalignments in \Delta\chi could be associated in part with relatively low core magnetic fields. Thus, although the overall distribution of \Delta\chi for all 40 sources in our sample does show a significant peak at \Delta\chi\sim0, it remains unclear what distinguishes these AGN cores from those showing appreciable misalignment between optical and VLBI-core polarization position angles.

Although the continua of radio-loud Active Galactic Nuclei (AGN) are typically dominated by synchrotron radiation over virtually the entire spectrum, it is not clear whether the radio and higher-frequency emission originate in the same or different parts of the jet. Several different radio--optical correlations based on polarization data have been found recently, suggesting that the optical and radio polarization may be closely related, and that the corresponding emission regions may be cospatial (Gabuzda et. al2006, Jorstad et al. 2007, D'Arcangelo et al. 2007) Our joint analysis of optical and VLBA polarization data for a sample of about 40 AGNs shows that, after correction for the inferred VLBA core Faraday rotations, most BL Lac objects and some quasars have aligned VLBA-core and optical polarizations, although many quasars also show no obvious relationship between their VLBA-core and optical polarization angles. This may indicate that not all AGNs have cospatial regions of optical and radio emission in their jets. However, another possibility is that some of the 7mm-2cm VLBA cores have Faraday rotations of the order of several tens of thousand of rad/m^2, which were not properly fit by our three-frequency data due to n*pi ambiguities in the observed polarization angles, leading to incorrect subtraction of the effects of the core Faraday rotation, and so incorrect "zero-wavelength" radio polarization angles. The possibility of such high core Faraday rotations is supported by the results of the parsec-scale Faraday-rotation studies of Zavala & Taylor (2004) and Jorstad et al. (2007).

We present an astrometric analysis of the binary systems ABDorA /ABDorC and ABDorBa / ABDorBb. These two systems of well-known late-type stars are gravitationally associated and they constitute the quadruple ABDoradus system. From the astrometric data available at different wavelengths, we report: (i) a determination of the orbit of ABDorC, the very low mass companion to ABDorA, which confirms the mass estimate of 0.090Msun reported in previous works; (ii) a measurement of the parallax of ABDorBa, which unambiguously confirms the long-suspected physical association between this star and ABDorA; and (iii) evidence of orbital motion of ABDorBa around ABDorA, which places an upper bound of 0.4Msun on the mass of the pair ABDorBa / ABDorBb (50% probability). Further astrometric monitoring of the system at all possible wavelengths would determine with extraordinary precision the dynamical mass of its four components.