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.

Observations of gamma-ray bursts (GRBs) at Very High Energy (VHE) offer a unique opportunity to investigate particle acceleration processes, magnetic fields and radiation fields in these events. Very Long Baseline Interferometry (VLBI) observations have been proven to be a powerful tool providing unique information on the source size of the GRBs at mas scales, as well as their accurate positions and possible expansion speeds. This paper reports on the follow-up observations of GRB 190114C, the first ever GRB detected with high significance at TeV photon energies by the MAGIC telescope, conducted with the East Asia VLBI Network (EAVN) at 22 GHz on three epochs, corresponding to 6, 15 and 32 days after the burst. The derived maps do not show any significant source above 5 sigma. The inferred upper limits on the GRB 190114C flux density at 22 GHz are used here to constrain the allowable two-dimensional parameter space for the afterglow emission. We find that our limits are consistent with most afterglow parameter combinations proposed so far in the literature. This is the first effort for the EAVN to search and monitor a radio transient in the Target of Opportunity mode. In addition to the useful constraints on GRB 190114C radio emission, experience gained from these observations is very helpful for future routine operation of EAVN transient program.

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.

This paper presents a survey of microwave front-end receivers installed at radio telescopes throughout the World. This unprecedented analysis was conducted as part of a review of front-end developments for Italian radio telescopes, initiated by the Italian National Institute for Astrophysics in 2016. Fifteen international radio telescopes have been selected to be representative of the instrumentation used for radio astronomical observations in the frequency domain from 300 MHz to 116 GHz. A comprehensive description of the existing receivers is presented and their characteristics are compared and discussed. The observing performances of the complete receiving chains are also presented. An overview of on-going developments illustrates and anticipates future trends in front-end projects to meet the most ambitious scientific research goals.

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.

Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to shed light on the crucial physical processes giving rise to the jet formation, as well as to its extraordinary radiation output up to gamma-ray energies. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm-VLBI polarimetric observations, reaching an unprecedented resolution (~50 microarcsec). We also investigate the variability and physical process occurring in the source during the observing period which coincides with a very active state of the source till high-energies. The Faraday rotation analysis between 3 and 7mm shows a gradient in rotation measure with a maximum value of ~6X10^4 rad/m^2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7mm EVPAs orientation is different when the component is exiting the core or crossing a stationary feature at ~0.1 mas. The interaction between the superluminal component and a recollimation shock at ~0.1 mas could have triggered the multi-wavelengths flares. The variability Doppler factor associated with such interaction is large enough to explain the high energy emission, as we infer from the analysis of gamma-ray and X-ray data, and it is in agreement with the Doppler factor obtained to explain the extraordinary optical flare by Raiteri et al.(2017).

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.

Understanding the launching, acceleration, and collimation of jets powered by active galactic nuclei remains an outstanding problem in relativistic astrophysics. This is partly because observational tests of jet formation models suffer from the limited angular resolution of ground-based very long baseline interferometry that has thus far been able to probe the transverse jet structure in the acceleration and collimation zone of only two sources. Here we report radio interferometric observations of 3C 84 (NGC 1275), the central galaxy of the Perseus cluster, made with an array including the orbiting radio telescope of the RadioAstron mission. The obtained image transversely resolves the edge-brightened jet in 3C 84 only 30 microarcseconds from the core, which is ten times closer to the central engine than what has been possible in previous ground-based observations, and it allows us to measure the jet collimation profile from ~ 100 to ~10000 gravitational radii from the black hole. The previously found, almost cylindrical jet profile on scales larger than a few thousand r_g is now seen to continue at least down to a few hundred r_g from the black hole and we find a broad jet with a transverse radius larger than about 250 r_g at only 350 r_g from the core. If the bright outer jet layer is launched by the black hole ergosphere, it has to rapidly expand laterally on scales smaller than 100 r_g. If this is not the case, then this jet sheath is likely launched from the accretion disk.

The very long baseline interferometry (VLBI) technique offers angular resolutions superior to any other instruments at other wavelengths, enabling unique science applications of high-resolution imaging of radio sources and high-precision astrometry. The East Asia VLBI Network (EAVN) is a collaborative effort in the East Asian region. The EAVN currently consists of 21 telescopes with diverse equipment configurations and frequency setups, allowing flexible subarrays for specific science projects. The EAVN provides the highest resolution of 0.5 mas at 22 GHz, allowing the fine imaging of jets in active galactic nuclei, high-accuracy astrometry of masers and pulsars, and precise spacecraft positioning. The soon-to-be-operational Five-hundredmeter Aperture Spherical radio Telescope (FAST) will open a new era for the EAVN. This state-of-the-art VLBI array also provides easy access to and crucial training for the burgeoning Asian astronomical community. This Perspective summarizes the status, capabilities and prospects of the EAVN.

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 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 results from our Giant Metrewave Radio Telescope (GMRT) HI observations of the Arp 305 system. The system consists of two interacting spiral galaxies NGC 4016 and NGC 4017, a large amount of resultant tidal debris and a prominent tidal dwarf galaxy (TDG) candidate projected within the tidal bridge between the two principal galaxies. Our higher resolution GMRT HI mapping, compared to previous observations, allowed detailed study of smaller scale features. Our HI analysis supports the conclusion in Hancock et al. (2009) that the most recent encounter between the pair occurred $\sim$ 4 $\times$ 10$^8$ yrs ago. The GMRT observations also show HI features near NGC 4017 which may be remnants of an earlier encounter between the two galaxies. The HI properties of the Bridge TDG candidate include: M(HI) $\sim$ 6.6 $\times$ 10$^8$ msolar and V(HI) = 3500$\pm$ 7 km/s, which is in good agreement with the velocities of the parent galaxies. Additionally the TDG's HI linewidth of 30 km/s and a modest velocity gradient together with its SFR of 0.2 msolar/yr add to the evidence favouring the bridge candidate being a genuine TDG. The Bridge TDG's \textit{Spitzer} 3.6 $\mu$m and 4.5 $\mu$m counterparts with a [3.6]--[4.5] colour $\sim$ -0.2 mag suggests stellar debris may have seeded its formation. Future spectroscopic observations could confirm this formation scenario and provide the metallicity which is a key criteria for the validation for TDG candidates.

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.

Radio-loud high-redshift quasars (HRQs), although only a few of them are known to date, are crucial for the studies of the growth of supermassive black holes (SMBHs) and the evolution of active galactic nuclei (AGN) at early cosmological epochs. Radio jets offer direct evidence of SMBHs, and their radio structures can be studied with the highest angular resolution using Very Long Baseline Interferometry (VLBI). Here we report on the observations of three HRQs (J0131-0321, J0906+6930, J1026+2542) at z>5 using the Korean VLBI Network and VLBI Exploration of Radio Astrometry Arrays (together known as KaVA) with the purpose of studying their pc-scale jet properties. The observations were carried out at 22 and 43 GHz in 2016 January among the first-batch open-use experiments of KaVA. The quasar J0906+6930 was detected at 22 GHz but not at 43 GHz. The other two sources were not detected and upper limits to their compact radio emission are given. Archival VLBI imaging data and single-dish 15-GHz monitoring light curve of J0906+6930 were also acquired as complementary information. J0906+6930 shows a moderate-level variability at 15 GHz. The radio image is characterized by a core-jet structure with a total detectable size of ~5 pc in projection. The brightness temperature, 1.9x10^{11} K, indicates relativistic beaming of the jet. The radio properties of J0906+6930 are consistent with a blazar. Follow-up VLBI observations will be helpful for determining its structural variation.

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.

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

Recently, both stellar mass-segregation and binary-fractions were uniformly measured on relatively large samples of Galactic Globular Clusters (GCs). Simulations show that both sizeable binary-star populations and Intermediate-Mass Black Holes (IMBHs) quench mass-segregation in relaxed GCs. Thus mass-segregation in GCs with a reliable binary-fraction measurement is a valuable probe to constrain IMBHs. In this paper we combine mass-segregation and binary-fraction measurements from the literature to build a sample of 33 GCs (with measured core-binary fractions), and a sample of 43 GCs (with a binary fraction measurement in the area between the core radius and the half-mass radius). Within both samples we try to identify IMBH-host candidates. These should have relatively low mass-segregation, a low binary fraction (< 5%), and short (< 1 Gyr) relaxation time. Considering the core binary fraction sample, no suitable candidates emerge. If the binary fraction between the core and the half-mass radius is considered, two candidates are found, but this is likely due to statistical fluctuations. We also consider a larger sample of 54 GCs where we obtained an estimate of the core binary fraction using a predictive relation based on metallicity and integrated absolute magnitude. Also in this case no suitable candidates are found. Finally, we consider the GC core- to half-mass radius ratio, that is expected to be larger for GCs containing either an IMBH or binaries. We find that GCs with large core- to half-mass radius ratios are less mass-segregated (and show a larger binary fraction), confirming the theoretical expectation that the energy sources responsible for the large core are also quenching mass-segregation

We present the results of Very Long Baseline Interferometry (VLBI) observations using the phase reference technique to detect weak Active Galactic Nuclei (AGN) cores in the Virgo cluster. Our observations were carried out using the Korean VLBI Network (KVN). We have selected eight representative radio galaxies, seven Virgo cluster members and one galaxy (NGC 4261) that is likely to be in the background. The selected galaxies are located in a range of density regions showing various morphology in 1.4 GHz continuum. Since half of our targets are too weak to be detected at K-band we applied a phase referencing technique to extend the source integration time by calibrating atmospheric phase fluctuations. We discuss the results of the phase referencing method at high frequency observations and we compare them with self-calibration on the relatively bright AGNs, such as M87, M84 and NGC 4261. In this manuscript we present the radio intensity maps at 22 GHz of the Virgo cluster sample while we demonstrate for first time the capability of KVN phase referencing technique.

Simultaneous observations at multiple frequency bands have the potential to overcome the fundamental limitation imposed by the atmospheric propagation in mm-VLBI observations. The propagation effects place a severe limit in the sensitivity achievable in mm-VLBI, reducing the time over which the signals can be coherently combined, and preventing the use of phase referencing and astrometric measurements. We carried out simultaneous observations at 22, 43, 87 and 130 GHz of a group of five AGNs, the weakest of which is ca. 200 mJy at 130 GHz, with angular separations ranging from 3.6 to 11 degrees, using the KVN. We analysed this data using the Frequency Phase Transfer (FPT) and the Source Frequency Phase Referencing (SFPR) techniques, which use the observations at a lower frequency to correct those at a higher frequency. The results of the analysis provide an empirical demonstration of the increase in the coherence times at 130 GHz from a few tens of seconds to about twenty minutes, with FPT, and up to many hours with SFPR. Moreover the astrometric analysis provides high precision relative position measurements between two frequencies, including, for the first time, astrometry at 130 GHz. Finally we demonstrate a method for the generalised decomposition of the relative position measurements into absolute position shifts for bona fide astrometric registration of the maps of the individual sources at multiple frequencies, up to 130 GHz.

Active Galactic Nuclei (AGN) with bright radio jets offer the opportunity to study the structure of and physical conditions in relativistic outflows. For such studies, multi-frequency polarimetric very long baseline interferometric (VLBI) observations are important as they directly probe particle densities, magnetic field geometries, and several other parameters. We present results from first-epoch data obtained by the Korean VLBI Network (KVN) within the frame of the Plasma Physics of Active Galactic Nuclei (PAGaN) project. We observed seven radio-bright nearby AGN at frequencies of 22, 43, 86, and 129 GHz in dual polarization mode. Our observations constrain apparent brightness temperatures of jet components and radio cores in our sample to $>10^{8.01}$ K and $>10^{9.86}$ K, respectively. Degrees of linear polarization $m_{L}$ are relatively low overall: less than 10%. This indicates suppression of polarization by strong turbulence in the jets. We found an exceptionally high degree of polarization in a jet component of BL Lac at 43 GHz, with $m_{L} \sim$ 40%. Assuming a transverse shock front propagating downstream along the jet, the shock front being almost parallel to the line of sight can explain the high degree of polarization.

We report first results from KVN and VERA Array (KaVA) VLBI observations obtained in the frame of our Plasma-physics of Active Galactic Nuclei (PAGaN) project. We observed eight selected AGN at 22 and 43 GHz in single polarization (LCP) between March 2014 and April 2015. Each source was observed for 6 to 8 hours per observing run to maximize the $uv$ coverage. We obtained a total of 15 deep high-resolution images permitting the identification of individual circular Gaussian jet components and three spectral index maps of BL Lac, 3C 111 and 3C 345 from simultaneous dual-frequency observations. The spectral index maps show trends in agreement with general expectations -- flat core and steep jets -- while the actual value of the spectral index for jets shows indications for a dependence on AGN type. We analyzed the kinematics of jet components of BL Lac and 3C 111, detecting superluminal proper motions with maximum apparent speeds of about $5c$. This constrains the lower limits of the intrinsic component velocities to $\sim0.98c$ and the upper limits of the angle between jet and line of sight to $\sim$20$\deg$. In agreement with global jet expansion, jet components show systematically larger diameters $d$ at larger core distances $r$, following the global relation $d\approx0.2r$, albeit within substantial scatter.

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.

Context: Pre-merger interactions between galaxies can induce significant changes in the morphologies and kinematics of the stellar and ISM components. Large amounts of gas and stars are often found to be disturbed or displaced as tidal debris. This debris then evolves, sometimes forming stars and occasionally tidal dwarf galaxies. Here we present results from our HI study of Arp 65, an interacting pair hosting extended HI tidal debris. Aims: In an effort to understand the evolution of tidal debris produced by interacting pairs of galaxies, including in situ star and tidal dwarf galaxy formation, we are mapping HI in a sample of interacting galaxy pairs. The Arp 65 pair is one of them. Methods: Our resolved HI 21 cm line survey is being carried out using the Giant Metrewave Radio Telescope (GMRT). We used our HI survey data as well as available SDSS optical, Spitzer infra-red and GALEX UV data to study the evolution of the tidal debris and the correlation of HI with the star-forming regions within it. Results: In Arp 65 we see a high impact pre-merger interaction involving a pair of massive galaxies (NGC 90 and NGC 93) that have a stellar mass ratio of ~ 1:3. The interaction, which probably occurred ~ 1.0 -- 2.5 $\times$ 10$^8$ yr ago, appears to have displaced a large fraction of the HI in NGC 90 (including the highest column density HI) beyond its optical disk. We also find extended ongoing star formation in the outer disk of NGC 90. In the major star-forming regions, we find the HI column densities to be ~ 4.7 $\times$ 10$^{20}$ cm$^{-2}$ or lower. But no signature of star formation was found in the highest column density HI debris, SE of NGC 90. This indicates conditions within the highest column density HI debris remain hostile to star formation and it reaffirms that high HI column densities may be a necessary but not sufficient criterion for star formation.

KaVA (KVN and VERA Array) is a new combined VLBI array with KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). First, we briefly review the imaging capabilities of KaVA array which actually achieves more than three times better dynamic range than that achieved by VERA alone. The KaVA images clearly show detailed structures of extended radio jets in AGNs. Next, we represent the key science program to be led by KaVA AGN sub working group. We will conduct the monitoring observations of Sgr A* and M87 because of the largeness of their central super-massive black hole angular sizes. The main science goals of the program are (i) testing magnetically-driven-jet paradigm by mapping velocity fields of the M87 jet, and (ii) obtaining tight constraints on physical properties of radio emitting region in Sgr A*.

The radio galaxy 3C84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGN) and one of the best laboratories to study the radio properties of subparsec scale jets. We discuss here the past and present activity of the nuclear region within the central 1pc and the properties of subparsec-sized components C1, C2 and C3. We compare these results with the high resolution space-VLBI image at 5GHz obtained with the RadioAstron satellite and we shortly discuss the possible correlation of radio emission with the gamma-ray emission.

We explore the degree of magnetization at the jet base of M87 by using the observational data of the event horizon telescope (EHT) at 230~GHz obtained by Doeleman et al. By utilizing the method in Kino et al., we derive the energy densities of magnetic fields ($U_{B}$) and electrons and positrons ($U_{\pm}$) in the compact region detected by EHT (EHT-region) with its full-width-half-maximum size $40~{\rm \mu as}$. First, we assume that an optically-thick region for synchrotron self absorption (SSA) exists in the EHT-region. Then, we find that the SSA-thick region should not be too large not to overproduce the Poynting power at the EHT-region. The allowed ranges of the angular size and the magnetic field strength of the SSA-thick region are $21~{\rm \mu as} \le \theta_{\rm thick}\le 26.3~{\rm \mu as}$ and $50~{\rm G}\le B_{\rm tot}\le 124~{\rm G}$, respectively. Correspondingly $U_{B}\gg U_{\pm} $ is realized in this case. We further examine the composition of plasma and energy density of protons by utilizing the Faraday rotation measurement ($RM$) at 230~GHz obtained by Kuo et al. Then, we find that $U_{B}\gg U_{\pm}+U_{p} $ still holds in the SSA-thick region. Second, we examine the case when EHT-region is fully SSA-thin. Then we find that $U_{B}\gg U_{\pm}$ still holds unless protons are relativistic. Thus, we conclude that magnetically driven jet scenario in M87 is viable in terms of energetics close to ISCO scale unless the EHT-region is fully SSA-thin and relativistic protons dominated.

We study the warping and tearing of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on an eccentric orbit. The circumbinary disk is significantly misaligned with the binary orbital plane, and is subject to the time-dependent tidal torques. In principle, such a disk is warped and precesses, and is torn into mutually misaligned rings in the region, where the tidal precession torques are stronger than the local viscous torques. We derive the tidal-warp and tearing radii of the misaligned circumbinary disks around eccentric SMBH binaries. We find that in disks with the viscosity parameter, alpha, larger than a critical value depending on the disk aspect ratio, the disk warping appears outside the tearing radius. This condition is expressed as alpha > sqrt{H/3r} for H/r ~<0.1, where H is the disk scale height. If alpha < sqrt{H/3r}, only the disk tearing occurs because the tidal warp radius is inside the tearing radius, where most of disk material is likely to rapidly accrete onto SMBHs. In warped and torn disks, both the tidal-warp and the tearing radii most strongly depend on the binary semi-major axis, although they also mildly depend on the other orbital and disk parameters. This strong dependence enables us to estimate the semi-major axis, once the tidal warp or tearing radius is determined observationally: For the tidal warp radius of 0.1 pc, the semi-major axis is estimated to be ~10^{-2} pc for 10^7 Msun black hole with typical orbital and disk parameters. We also briefly discuss the possibility that central objects of observed warped maser disks in active galactic nuclei are supermassive black hole binaries.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

We study a warping instability of a geometrically thin, non-self-gravitating, circumbinary disk around young binary stars on an eccentric orbit. Such a disk is subject to both the tidal torques due to a time-dependent binary potential and the radiative torques due to radiation emitted from each star. The tilt angle between the circumbinary disk plane and the binary orbital plane is assumed to be very small. We find that there is a radius within/beyond which the circumbinary disk is unstable to radiation-driven warping, depending on the disk density and temperature gradient indices. This marginally stable warping radius is very sensitive to viscosity parameters, a fiducial disk radius and the temperature measured there, the stellar luminosity, and the disk surface density at a radius where the disk changes from the optically thick to thin for the irradiation from the central stars. On the other hand, it is insensitive to the orbital eccentricity and binary irradiation parameter, which is a function of the binary mass ratio and luminosity of each star. Since the tidal torques can suppress the warping in the inner part of the circumbinary disk, the disk starts to be warped in the outer part. While the circumbinary disks are most likely to be subject to the radiation-driven warping on a AU to kilo-AU scale for binaries with young massive stars more luminous than 10^4Lsun, the radiation driven warping does not work for those around young binaries with the luminosity comparable to the solar luminosity.

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 results from our Giant Metrewave Radio Telescope (GMRT) HI observations of the interacting pair Arp 202 (NGC 2719 and NGC 2719A). Earlier deep UV(GALEX) observations of this system revealed a tidal tail like extension with a diffuse object towards its end, proposed as a tidal dwarf galaxy (TDG) candidate. We detect HI emission from the Arp 202 system, including HI counterparts for the tidal tail and the TDG candidate. Our GMRT HI morphological and kinematic results clearly link the HI tidal tail and the HI TDG counterparts to the interaction between NGC 2719 and NGC 2719A, thus strengthening the case for the TDG. The Arp 202 TDG candidate belongs to a small group of TDG candidates with extremely blue colours. In order to gain a better understanding of this group we carried out a comparative study of their properties from the available data. We find that HI (and probably stellar) masses of this extremely blue group are similar to the lowest HI mass TDGs in the literature. However the number of such blue TDG candidates examined so far is too small to conclude whether or not their properties justify them to be considered as a subgroup of TDGs.

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.

We present our work towards using the Korean VLBI (Very Long Baseline Interferometer) Network (KVN) and VLBI Exploration of Radio Astronomy (VERA) arrays combined into the KVN and VERA Array (KaVA) for observations of radio pulsars at high frequencies ($\simeq$22-GHz). Pulsar astronomy is generally focused at frequencies approximately 0.3 to several GHz and pulsars are usually discovered and monitored with large, single-dish, radio telescopes. For most pulsars, reduced radio flux is expected at high frequencies due to their steep spectrum, but there are exceptions where high frequency observations can be useful. Moreover, some pulsars are observable at high frequencies only, such as those close to the Galactic Center. The discoveries of a radio-bright magnetar and a few dozen extended Chandra sources within 15 arc-minute of the Galactic Center provide strong motivations to make use of the KaVA frequency band for searching pulsars in this region. Here, we describe the science targets and report progresses made from the KVN test observations for known pulsars. We then discuss why KaVA pulsar observations are compelling.

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 study a warping instability of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on a circular orbit. Such a circumbinary disk is subject to not only tidal torques due to the binary gravitational potential but also radiative torques due to radiation emitted from an accretion disk around each black hole. We find that a circumbinary disk initially aligned with the binary orbital plane is unstable to radiation-driven warping beyond the marginally stable warping radius, which is sensitive to both the ratio of vertical to horizontal shear viscosities and the mass-to-energy conversion efficiency. As expected, the tidal torques give no contribution to the growth of warping modes but tend to align the circumbinary disk with the orbital plane. Since the tidal torques can suppress the warping modes in the inner part of circumbinary disk, the circumbinary disk starts to be warped at radii larger than the marginally stable warping radius. If the warping radius is of the order of 0.1 pc, a resultant semi-major axis is estimated to be of the order of 10^-2 pc to 10^-4 pc for 10^7 Msun black hole. We also discuss the possibility that the central objects of observed warped maser disks in active galactic nuclei are binary supermassive black holes with a triple disk: two accretion disks around the individual black holes and one circumbinary disk surrounding them.

We present the results of radio monitoring observations of Sgr A* at 7 mm (i.e. 43 GHz) with VLBI Exploration of Radio Astrometry (VERA), which is a VLBI array in Japan. VERA provides angular resolutions on millisecond scales, resolving structure within ~100 Schwarzschild radii of Sgr A* similar to Very Large Baseline Array (VLBA). We performed multi-epoch observations of Sgr A* in 2005 - 2008, and started monitoring it again with VERA from January 2013 for tracing the current G2 encounter event. Our preliminary results in 2013 show that Sgr A* on mas scales has been in ordinary state as of August 2013, although some fraction of the G2 cloud already passed pericenter of Sgr A* in April 2013. We will continue on monitoring Sgr A* with VERA and newly developed KaVA (KVN and VERA Array).

We investigate the relationship between radio power and properties related to active galactic nuclei (AGNs). Radio power at 1.4 or 5 GHz, which has been used in many studies, can be affected by synchrotron self-absorption and free-free absorption in a dense region. On the other hand, these absorption effects get smaller at higher frequencies. Thus, we performed simultaneous observations at 22 and 43 GHz using the Korean VLBI Network (KVN) radio telescope based on a sample of 305 AGN candidates residing in elliptical galaxies from the overlap between the Sloan Digital Sky Survey (SDSS) Data Release 7 and Faint Images of the Radio Sky at Twenty-Centimeters (FIRST). About 37% and 22% of the galaxies are detected at 22 and 43 GHz, respectively. Assuming no flux variability between the FIRST and KVN observation, spectral indices were derived from FIRST and KVN data and we found that over 70% of the detected galaxies have flat or inverted spectra, implying the presence of optically thick compact regions near the centres of the galaxies. Core radio power does not show a clear dependence on black hole mass at either low (1.4 GHz) or high (22 and 43 GHz) frequencies. However, we found that the luminosity of the [OIII] $\lambda$5007 emission line and the Eddington ratio correlate with radio power more closely at high frequencies than at low frequencies. This suggests that radio observation at high frequencies can be an appropriate tool for unveiling the innermost region. In addition, the luminosity of the [OIII] $\lambda$5007 emission line and the Eddington ratio can be used as a tracer of AGN activity. Our study suggests a causal connection between high frequency radio power and optical properties of AGNs.

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.

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.

The BL Lac object S5 0716+71 was observed in a global multi-frequency campaign to search for rapid and correlated flux density variability and signatures of an inverse-Compton (IC) catastrophe during the states of extreme apparent brightness temperatures. The observing campaign involved simultaneous monitoring at radio to IR/optical wavelengths centered around a 500-ks INTEGRAL pointing (November 10-17, 2003). We present the combined analysis and results of the cm- to sub-mm observations including a detailed study of the inter- to intra-day variability and spectral characteristics of 0716+714. We further constrain the variability brightness temperatures (T_B) and Doppler factors (delta) comparing the radio-bands with the hard X-ray emission (3-200 keV). 0716+714 was in an exceptionally high state (outburst) and different (slower) phase of short-term variability. The flux density variability in the cm- to mm-bands is dominated by a correlated, ~4 day time scale amplitude increase of up to ~35% systematically more pronounced towards shorter wavelengths. This contradicts expectations from standard interstellar scintillation (ISS) and suggests a source-intrinsic origin of the variability. The derived lower limits to T_B exceed the 10^12 K IC-limit by up to 3-4 orders of magnitude. Assuming relativistic boosting, we obtain robust and self-consistent lower limits of delta >= 5-33, in good agreement with delta_VLBI obtained from VLBI studies and the IC-Doppler factors delta_IC > 14-16 obtained from the INTEGRAL data. Since a strong contribution from ISS can be excluded and a simultaneous IC catastrophe was not observed, we conclude that relativistic Doppler boosting naturally explains the apparent violation of the theoretical limits within standard synchrotron-self-Compton (SSC) jet models of AGN.

The capabilities of the Global mm-VLBI Array are summarized and demonstrated through actual images from our monitoring of extragalactic radio jets. This sensitive 3mm-VLBI interferometer is able to provide images of up to 50 microarcseconds resolution. For the near future, ALMA, the GBT, the LMT, CARMA, SRT, Yebes, Nobeyama and Noto are some of the most sensitive stations suitable to participate in mm-VLBI. This future array, together with the present Global mm-VLBI Array, would achieve 10 times better sensitivities than nowadays. Image fidelity would also largely increased. T he addition of ALMA would improve the (u,v)-coverage for sources with low declination (<20 deg.) and facilitate the VLBI imaging of the Galactic Centre source SgrA*.

The 18 mag QSO 1633+382 (4C38.41, z=1.807) showed a very pronounced outburst in 2001/2002. With a peak amplitude of more than 9 Jy at 90GHz, this flare was brighter than any known previous flare in this source (data available since 1980).During onset, the mm-flare was particulary fast, with an increase of more than 2 Jy at 230 GHz in less than 8 days. Since January 2002, the mm-flux of 1633+382 is decaying. During this decline, however, local flux variations with amplitudes of 1-3 Jy were seen, indicative of underlying and more rapid source activity on time scales of 1-2 months. After the main peak occurring in 2001.99, the 90 GHz flux showed secondary maxima at approximately half year intervals in 2002.3, 2002.7 and 2003.13. This kind of periodicity might be explained via the lighthouse model (Camenzind and Krockenberger 1992), which is based on the magnetic accelerator of Blanford and Payne (1982). At present the millimeter flux is nearly back to its quiescent level of 2-2.5 Jy, which the source had before the flare began. Our VLBA Polarimetry monitoring started June 2002 during the onset of the flare. At cm wavelength, the flare is only marginally detected which implies very high opacity of the source.

Korea's new VLBI project to construct the Korean VLBI Network (KVN) started in 2001, as a 7-year project that is fully funded by our government. We plan to build 3 new high-precision radio telescopes of 21-m diameter in 3 places in Korea, which will be exclusively used for VLBI observations. We will install the 2/8, 22 and 43 GHz HEMT receivers within 2007 as a first target, and later we will expand the receiving frequency up to 86 and 129 GHz for astronomical, geodetic, and earth science VLBI research. The millimeter-wave VLBI will be the ultimate goal of KVN. For the front-ends, we are going to install a multi-channel receiver system that employs low-pass filters within a quasi-optical beam transportation system. This receiver system will give reliable phase calibrations for millimeter-wave VLBI as well as enable simultaneous multi-frequency band observations. The hard-disk type new Mark 5 will be used as the main recorder of KVN. We have completed the design of the KVN DAS system of 2 Gsps sampling rate, which will use 4 data streams to meet the multi-channel requirement. A VERA type DAS modified for Mark 5 recorder is also under consideration. A new correlator project for KVN was recently approved from Korean government, and will start in the second half of 2004.

We summarize the present status of VLBI experiments at 3 mm (86 GHz), 2 mm (129-150 GHz) and 1.3 mm (215-230 GHz). We present and discuss a new 3 mm VLBI map of M87 (Virgo A), which has a spatial resolution of only approx. 20 Schwarzschild radii. We discuss recent VLBI results for SgrA* and argue in favor of new observations within an extended European mm-VLBI network, in order to search for variability. We discuss the possibilities to image the `event horizon' of a super-massive black hole at wavelengths < 2mm, and conclude that the addition of large and sensitive millimetre telescopes such as CARMA, the SMA, the LMT and ALMA will be crucial for this.

The so-called Spectral Curvature Parameter(SCP), when plotted versus the high-frequency spectral index ($alpha_{high}$) of synchrotron sources, provides crucial parameters on the continuum spectrum of synchrotron radiation without the more complex modeling of spectral ageing scenarios. An important merit of the SCP-$alpha$ diagram is the enhanced reliability of extracting multiple injection spectra, $alpha_{inj}$. Different from the colour-colour diagram, tracks of different $alpha_{inj}$s, especially when the synchrotron particles are young, exhibit less overlap and less smearing in the SCP-$alpha$ diagram. Three giant radio galaxies(GRGs) and a sample of Compact steep spectrum(CSS) souces are presented. GRGs exhibit asymmetries of their injection spectral indices $alpha_{inj}$ in the SCP-$alpha_{high}$ diagram. The obtained $alpha_{inj}$s and the trends in the sources are cross-checked with the literature and show remarkable confidence. Besides the spectral steepening, spectral flattening is prominent in the radio lobes. The spectral flattening is a clue to efficient re-acceleration processes in the lobes. It implies interaction with the surrounding intergalactic or intra-cluster medium is an important characteristic of GRGs. In the SW lobe of DA240, there is a clear sign of CI and KP/JP bifurcation at the source extremity. This indicates a highly relativistic energy transportation from the core or in situ acceleration in this typical FR I lobe. Our analysis proves, if exists, KP spectra imply the existence of strong $B_{sync}$ field with $B_{sync} > B_{CMB}$. In the CSS sources, our result confirms the CI model and $B_{sync} >> B_{CMB}$. The synchrotron self-absorption is significant in the CSS sample.