A far-infrared counterpart to the west hot spot of the radio galaxy Pictor A is discovered with the Spectral and Photometric Imaging REceiver (SPIRE) onboard Herschel. The color-corrected flux density of the source is measured as $70.0 \pm 9.9$ mJy at the wavelength of 350 $\mu$m. A close investigation into its radio-to-optical spectrum indicates that the mid-infrared excess over the radio synchrotron component, detected with WISE and Spitzer, significantly contributes to the far-infrared band. Thanks to the SPIRE data, it is revealed that the spectrum of the excess is described by a broken power-law model subjected to a high-energy cutoff. By applying the radiative cooling break under continuous energy injection ($\Delta \alpha = 0.5$), the broken power-law model supports an idea that the excess originates in 10-pc scale substructures within the hot spot. From the break frequency, $\nu_{\rm b} = 1.6_{-1.0}^{+3.0} \times 10^{12}$ Hz, the magnetic field was estimated as $B\simeq1$-$4$ mG. This is higher than the minimum-energy magnetic field of the substructures by a factor of $3$--$10$. Even if the origin of the excess is larger than $\sim 100$ pc, the magnetic field stronger than the minimum-energy field is confirmed. It is proposed that regions with a magnetic field locally boosted via plasma turbulence are observed as the substructures. The derived energy index below the break, $\alpha \sim 0.22$ (conservatively $<0.42$), is difficult to be attributed to the strong-shock acceleration ($\alpha = 0.5$). Stochastic acceleration and magnetic reconnection are considered as a plausible alternative mechanism.

The supermassive black holes (SMBHs) of narrow-line Seyfert 1 galaxies (NLS1s) are at the lowest end of mass function of active galactic nuclei (AGNs) and preferentially reside in late-type host galaxies with pseudobulges, which are thought to be formed by internal secular evolution. On the other hand, the population of radio-loud NLS1s presents a challenge for the relativistic jet paradigm that powerful radio jets are exclusively associated with very high mass SMBHs in elliptical hosts, which are built-up through galaxy mergers. We investigated distorted radio structures associated with the nearest gamma-ray emitting, radio-loud NLS1 1H 0323+342. This provides supporting evidence for the merger hypothesis based on the past optical/near-infrared observations of its host galaxy. The anomalous radio morphology consists of two different structures, the inner curved structure of currently active jet and the outer linear structure of low-brightness relics. Such a coexistence might be indicative of the stage of an established black hole binary with precession before the black holes coalesce in the galaxy merger process. 1H 0323+342 and other radio-loud NLS1s under galaxy interactions may be extreme objects on the evolutionary path from radio-quiet NLS1s to normal Seyfert galaxies with larger SMBHs in classical bulges through mergers and merger-induced jet phases.

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

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.

A new method for measuring the global magnetic field structure of the Galactic plane is presented. We have determined the near-infrared polarization of field stars around 52 Cepheids found in recent surveys toward the Galactic plane. The Cepheids are located at the galactic longitudes $-10^{\circ}\leq \, l\, \leq +10.5^{\circ}$ and latitudes $-0.22^{\circ}\leq \, l\, \leq +0.45^{\circ}$, and their distances are mainly in the range of 10 to 15 kpc from the Sun. Simple classification of the sightlines is made with the polarization behavior vs. $H-K_{\mathrm S}$ color of field stars, and typical examples of three types are presented. Then, division of the field stars in each line of sight into (a) foreground, (b) bulge, and (c) background is made with the $Gaia$ DR2 catalog, the peak of the $H-K_{\mathrm S}$ color histogram, and $H-K_{\mathrm S}$ colors consistent with the distance of the Cepheid in the center, respectively. Differential analysis between them enables us to examine the magnetic field structure more definitely than just relying on the $H-K_{\mathrm S}$ color difference. In one line of sight, the magnetic field is nearly parallel to the Galactic plane and well aligned all the way from the Sun to the Cepheid position on the other side of the Galactic center. Contrary to our preconceived ideas, however, sightlines having such well-aligned magnetic fields in the Galactic plane are rather small in number. At least 36 Cepheid fields indicate random magnetic field components are significant. Two Cepheid fields indicate that the magnetic field orientation changes more than 45 in the line of sight. The polarization increase per color change $P$/ ($H-K_{\mathrm S}$) varies from region to region, reflecting the change in the ratio of the magnetic field strength and the turbulence strength.

We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

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 Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

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

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events.

Imaging the immediate vicinity of supermassive black holes (SMBH) and extracting a BH-spin signature is one of the grand challenges in astrophysics. M87 is known as one of the best targets for imaging the BH shadow and it can be partially thick against synchrotron self-absorption (SSA), particularly in a flaring state with high mass accretion rate. However, little is known about influences of the SSA-thick region on BH shadow images. Here we investigate BH shadow images of M87 at 230 GHz properly taking into account the SSA-thick region. When the BH has a high spin value, the corresponding BH shadow image shows the positional offset between the center of the photon ring and that of the SSA-thick ring at the innermost stable circular orbit (ISCO) due to the frame-dragging effect in the Kerr spacetime. As a result, we find that a dark-crescent structure is generally produced between the photon ring and the SSA-thick ISCO ring in the BH shadow image. The scale size of the dark-crescent increase with BH spin: its width reaches up to $\sim 2$ gravitational radius when the BH spin is 99.8% of its maximum value. The dark crescent is regarded as a new signature of a highly spinning BH. This feature is expected to appear in flaring states with relatively high mass accretion rate rather than the quiescent states. We have simulated the image reconstruction of our theoretical image by assuming the current and future Event Horizon Telescope (EHT) array, and have found that the future EHT including space-very long baseline interferometry in 2020s can detect the dark crescent.

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.

The narrow-line Seyfert 1 (NLS1) galaxy PMN J0948+0022 is an archetype of gamma-ray-emitting NLS1s in active galactic nuclei (AGNs). In this study, we investigate its radio structures using archival data obtained using the Karl G.~Jansky very large array (VLA) and the very long baseline array (VLBA). The new VLA images reveal an outermost radio emission separated by 9.1 arcsec. Its resolved structure and steep spectrum suggest a terminal shock in a radio lobe energized by the jet from the PMN J0948+0022 nucleus. This large-scale radio component is found at almost the same position angle as that of the pc-scale jet, indicating a stable jet direction up to ~1 Mpc. Its apparent one-sidedness implies a moderate advancing speed ($\beta>0.27$). The kinematic age of 1 x 10^7 year of the jet activity is consistent with the expected NLS1 phase of ~10^7 year in the AGN lifetime. The VLBA image reveals the jet structure at distances ranging from r=0.82 milliarcsec to 3.5 milliarcsec, corresponding to approximately 100 pc, where superluminal motions were found. The jet width profile ($\propto r^{1.12}$) and flux-density profile ($\propto r^{-1.44}$) depending on the distance $r$ along the jet suggest that the jet kinetic energy is converted into internal energy in this region. The jet is causally connected in a nearly conical streamline, which is supported by ambient pressure at 100 pc scales in the host galaxy of PMN J0948+0022.

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.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

We study Faraday rotation in the jet of M87 inside the Bondi radius using eight Very Long Baseline Array data sets, one at 8 GHz, four at 5 GHz, and three at 2 GHz. We obtain Faraday rotation measures (RMs) measured across the bandwidth of each data set. We find that the magnitude of RM systematically decreases with increasing distance from the black hole from 5,000 to 200,000 Schwarzschild radii. The data, showing predominantly negative RM sign without significant difference of the RMs on the northern and southern jet edges, suggest that the spatial extent of the Faraday screen is much larger than the jet. We apply models of hot accretion flows, thought to be prevalent in active galactic nuclei having relatively low luminosity such as M87, and find that the decrease of RM is described well by a gas density profile $\rho \propto r^{-1}$. This behavior matches the theoretically expected signature of substantial winds, nonrelativistic un-collimated gas outflows from hot accretion flows, which is consistent with the results of various numerical simulations. The pressure profile inferred from the density profile is flat enough to collimate the jet, which can result in gradual acceleration of the jet in a magneto-hydrodynamical process. This picture is in good agreement with the observed gradual collimation and acceleration of the M87 jet inside the Bondi radius. The dominance of negative RMs suggests that jet and wind axis are misaligned such that the jet emission exposes only one side of the toroidal magnetic fields permeating the winds.

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.

This paper theoretically investigates the relations between the structure of relativistic jets and produced synchrotron images, by using a steady axisymmetric force-free jet model. We especially focus on the limb-brightened jets that are largely symmetric to the jet axes and observed in some active galactic nuclei such as M87, Mrk~501, Cyg~A, and 3C84. We find that symmetrically limb-brightened images can be produced when magnetic field lines of the jet penetrate a fast-spinning black hole as motivated by the Blandford-Znajek mechanism. On the other hand, jets with magnetic field lines that pass through a slowly spinning black hole or the Keplerian accretion disk produce highly asymmetric radio images. In addition, the edge of a counter jet tends to be luminous in the accretion-disk model even for rather small viewing angles, which may be problematic for some observed jets. We also suggest that the site of particle accelerations in relativistic jets can be constrained by fitting the radio images to observations. This kind of study focusing on the jet images far away from the central engine is complimentary to those concentrating directly on the innermost region with upcoming data of the Event Horizon Telescope.

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

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

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 present a kinematic study of the subparsec-scale radio jet of the radio galaxy 3C 84/NGC 1275 with the VLBI Exploration of Radio Astrometry (VERA) array at 22 GHz for 80 epochs from 2007 October to 2013 December. The averaged radial velocity of the bright component "C3" with reference to the radio core is found to be $0.27 pm 0.02c$ between 2007 October and 2013 December. This constant velocity of C3 is naturally explained by the advancing motion of the head of the mini-radio lobe. We also find a non-linear component in the motion of C3 with respect to the radio core. We briefly discuss possible origins of this non-linear motion.

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

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

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

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.

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

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

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

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

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

We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Array) and three additional telescopes in Japan (Takahagi 32m, Yamaguchi 32m and Nobeyama 45m). By comparing the independently obtained gain values based on these two methods, we found that the gain values from each method were consistent within 10% at elevations higher than 10 degree. We also found that the total flux densities of two images produced from the different amplitude calibrations were in agreement within 10% at both 22 and 43 GHz. By using the template spectrum method, furthermore, the additional radio telescopes can participate in the KaVA (i.e. EAVN) so that it can give a notable sensitivity increase. Therefore, our results will constrain the detailed conditions to reliably measure the VLBI amplitude using EAVN and give a potential to extend possible telescopes comprising EAVN.

Mid-infrared properties are reported of the west hot spot of the radio galaxy Pictor A with the Wide-field Infrared Survey Explorer (WISE). The mid-infrared counterpart to the hot spot, WISE J051926.26-454554.1, is listed in the AllWISE source catalog. The source was detected in all the four WISE photometric bands. A comparison between the WISE and radio images reinforces the physical association of the \wise\ source to the hot spot. The WISE flux density of the source was carefully evaluated. A close investigation of the multi-wavelength synchrotron spectral energy distribution from the object reveals a mid-infrared excess at the wavelength of $\lambda=22$ $\mu$m with a statistical significance of $4.8 \sigma$ over the simple power-law extrapolation from the synchrotron radio spectrum. The excess is reinforced by single and double cutoff power-law modeling of the radio-to-optical spectral energy distribution. The synchrotron cutoff frequency of the main and excess components was evaluated as $7.1 \times 10^{14}$ Hz and $5.5 \times 10^{13}$ Hz, respectively. From the cutoff frequency, the magnetic field of the emission region was constrained as a function of the region size. In order to interpret the excess component, an electron population different from the main one dominating the observed radio spectrum is necessary. The excess emission is proposed to originate in a sub structure within the hot spot, in which the magnetic field is by a factor of a few stronger than that in the minimum energy condition. The relation of the mid-infrared excess to the X-ray emission is briefly discussed.

We present Very Long Baseline Array polarimetric observations of the innermost jet of 3C$\sim$84 (NGC$\sim$1275) at 43$\sim$GHz. A significant polarized emission is detected at the hotspot of the innermost re-started jet, which is located $\sim$1 pc south from the radio core. While the previous report presented a hotspot at the southern end of the western limb, the hotspot location has been moved to the southern end of the eastern limb. Faraday rotation is detected within an entire bandwidth of the 43-GHz band. The measured rotation measure (RM) is at most (6.3$\pm$1.9)$\times10^{5}$$\sim$rad$\sim$m$^{-2}$ and might be slightly time variable on the timescale of a month by a factor of a few. Our measured RM and the RM previously reported by the CARMA and SMA observations cannot be consistently explained by the spherical accretion flow with a power-law profile. We propose that a clumpy/inhomogeneous ambient medium is responsible for the observed rotation measure. Using equipartition magnetic field, we derive the electron density of $2\times10^{4}$$\sim$cm$^{-3}$. Such an electron density is consistent with the cloud of narrow line emission region around the central engine. We also discuss the magnetic field configuration from black hole scale to pc scale and the origin of low polarization.

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

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

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

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

High-resolution millimeter wave very-long-baseline interferometry (mm-VLBI) is an ideal tool for probing the structure at the base of extragalactic jets in detail. The TeV blazar Mrk 501 is one of the best targets among BL Lac objects for studying the nature of off-axis jet structures because it shows different jet position angles at different scales. The aim of this study is to investigate the properties of the off-axis jet structure through high-resolution mm-VLBI images at the jet base and physical parameters such as kinematics, flux densities, and spectral indices. We performed Very Long Baseline Array (VLBA) observations over six epochs from 2012 February to 2013 February at 43 GHz. Quasi-simultaneous Global Millimeter VLBI Array (GMVA) observations at 86 GHz were performed in May 2012. We discover a new jet component at the northeast direction from the core in all the images at 43 and 86 GHz. The new component shows the off-axis location from the persistent jet extending to the southeast. The 43 GHz images reveal that the scattering of the positions of the NE component is within ~0.2 mas. The 86 GHz data reveals a jet component located 0.75 mas southeast of the radio core. We also discuss the spectral indices between 43 and 86 GHz, where the northeast component has steeper spectral index and the southeast component has comparable or flatter index than the radio core does.

We examine the plasma composition of relativistic jets in four FRII radio galaxies by analyzing the total cocoon pressure in terms of partial pressures of thermal and non-thermal electrons/positrons and protons. The total cocoon pressure is determined by cocoon dynamics via comparison of theoretical model with the observed cocoon shape. By inserting the observed number density of non-thermal electrons/positrons and the upper limit of thermal electron/positron number density into the equation of state, the number density of protons is constrained. We apply this method to four FRII radio galaxies (Cygnus A, 3C 219, 3C 223 and 3C 284), for which the total cocoon pressures have been already evaluated. We find that the positron-free plasma comprising of protons and electrons is ruled out, when we consider plausible particle distribution functions. In other words, the mixture of positrons is required for all four FRII radio galaxies; the number density ratio of electrons/positrons to protons is larger than two. Thus, we find that the plasma composition is independent of the jet power and the size of cocoons. We also investigate the additional contribution of thermal electrons/positrons and protons on the cocoon dynamics. When thermal electrons/positrons are absent, the cocoon is supported by the electron/ proton plasma pressure, while both electron/positron pressure supported and electron/proton plasma pressure supported cocoons are allowed if the number density of thermal electrons/positrons is about 10 times larger than that of non-thermal ones.

We report on results from new high-sensitivity, high-resolution 86GHz (3.5 millimeter) observations of the jet base in the nearby radio galaxy M87, obtained by the Very Long Baseline Array in conjunction with the Green Bank Telescope. The resulting image has a dynamic range exceeding 1500 to 1, the highest ever achieved for this jet at this frequency, resolving and imaging a detailed jet formation/collimation structure down to ~10 Schwarzschild radii (Rs). The obtained 86GHz image clearly confirms some important jet features known at lower frequencies, i.e., a wide-opening angle jet base, a limb-brightened intensity profile, a parabola-shape collimation profile and a counter jet. The limb-brightened structure is already well developed at < 0.2mas (< 28Rs, projected) from the core, where the corresponding apparent opening angle becomes as wide as ~100 degrees. The subsequent jet collimation near the black hole evolves in a complicated manner; there is a "constricted" structure at tens Rs from the core, where the jet cross section is locally shrinking. We suggest that an external pressure support from the inner part of radiatively-inefficient accretion flow may be dynamically important in shaping/confining the footprint of the magnetized jet. We also present the first VLBI 86GHz polarimetric experiment for this source, where a highly polarized (~20%) feature is detected near the jet base, indicating the presence of a well-ordered magnetic field. As a by-product, we additionally report a 43/86 GHz polarimetric result for our calibrator 3C 273 suggesting an extreme rotation measure near the core.

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

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.

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 investigate the location of the radio jet bases ("radio cores") of blazars in radio images, and their stationarity by means of dense very long baseline interferometry (VLBI) observations. In order to measure the position of a radio core, we conducted 12 epoch astrometric observation of the blazar Markarian 421 with the VLBI Exploration of Radio Astrometry at 22 GHz immediately after a large X-ray flare, which occurred in the middle of 2011 September. For the first time, we find that the radio core is not stationary but rather changes its location toward 0.5 mas downstream. This angular scale corresponds to the de-projected length of a scale of $10^5$ Schwarzschild radii (Rs) at the distance of Markarian~421. This radio-core wandering may be a new type of manifestation associated with the phenomena of large X-ray flares.

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 examine the fate of a dead radio source in which jet injection from the central engine has stopped at an early stage of its evolution ($t = t_j \lesssim 10^5$ yr). To this aim, we theoretically evaluate the evolution of the emission from both the lobe and the shell, which are composed of shocked jet matter and a shocked ambient medium, respectively. Based on a simple dynamical model of expanding lobe and shell, we clarify how the broadband spectrum of each component evolves before and after the cessation of the jet activity. It is shown that the spectrum is strongly dominated by the lobe emission while the jet is active ($t \leq t_j$). On the other hand, once the jet activity has ceased ($t > t_j$), the lobe emission fades out rapidly, since fresh electrons are no longer supplied from the jet. Meanwhile, shell emission only shows a gradual decrease, since accelerated electrons are continuously supplied from the bow shock that is propagating into the ambient medium. As a result, overall emission from the shell overwhelms that from the lobe at wide range of frequencies from radio up to gamma-ray soon after the jet activity has ceased. Our result predicts a new class of dead radio sources that are dominated by shell emission. We suggest that the emission from the shell can be probed in particular at a radio wavelengths with the Square Kilometer Array (SKA) phase 1.

We present results from 43 GHz (VLBA, six epochs from 2012.2 to 2013.2) and 86 GHz (GMVA, one epoch in 2012.4) observations toward the basis of the jet in the TeV Blazar Mrk 501. The 43-GHz data analysis reveals a new feature located northeast of the radio core, with a flux density of several tens of mJy, perpendicularly to the jet axis. The 86-GHz image shows the jet feature located 0.75 mas southeast of the radio core, which is consistent with the previous result. The location of Gaussian model for 0.75 mas feature does not coincide with those for the jets in the 43-GHz image, however, a distribution of emission is found. We also discuss the spectral indices of the core, the northeast feature, and the jet feature between 43 GHz and 86 GHz, which show flat-to-steep, steep, and flat-to-invert, respectively.

We report on VLBI observations of NGC 1275 with VLBI Exploration of Radio Astrometry (VERA) and Very Long Baseline Array (VLBA) during 2008 and 2013 at 15, 22 and 43 GHz. Our observations provide long-term variations in the radio flux and structure of the inner jet in NGC 1275 on sub-parsec scales, particularly at C3 which is a new radio component accountable for the recent active state in radio regime since early-2000s. We found the apparent velocity of C3 was $\beta_{\rm{app}} = 0.267 \pm 0.007 $, which was consistent with typical velocities of the hot spot and/or radio lobes in the young radio sources. Furthermore, the radio flux and size of C3 had simultaneously increased during September 2007 and May 2012 with keeping optically-thin spectra, suggesting that the particle acceleration occurred inside C3 as generally seen in the hot spots of the radio lobes. Our new observations suggest that C3 is a new-born hot spot and/or radio lobe at a very early stage of its evolution. Our results imply the radio lobes in the radio galaxies might be already formed on sub-parsec scales from the central SMBH and might play an important role to produce $\rm \gamma$-ray emission in the young radio sources.

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

In this work we present a polarization analysis at radio frequencies of Markarian 421, one of the closest (z=0.03) TeV blazars. The observations were obtained, both in total and in polarized intensity, with the Very Long Baseline Array (VLBA) at 15, 24, and 43 GHz throughout 2011, with one observation per month (for a total of twelve epochs). We investigate the magnetic field topology and the polarization structure on parsec scale and their evolution with time. We detect polarized emission both in the core and in the jet region, and it varies with frequency, location and time. In the core region we measure a mean fractional polarization of about 1-2%, with a peak of about 4% in March at 43 GHz; the polarization angle is almost stable at 43 GHz, but it shows significant variability in the range 114-173 deg at 15 GHz. In the jet region the polarization properties show a more stable behavior; the fractional polarization is about 16% and the polarization angle is nearly perpendicular to the jet axis. The higher EVPA variability observed at 15 GHz is due both to a variable Faraday rotation effect and to opacity. The residual variability observed in the intrinsic polarization angle, together with the low degree of polarization in the core region, could be explained with the presence of a blend of variable cross-polarized subcomponents within the beam.

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.

The relativistic jet in M87 offers a unique opportunity for understanding the detailed jet structure and emission processes due to its proximity. In particular, the peculiar jet region HST-1 at ~1 arcsecond (or 80 pc, projected) from the nucleus has attracted a great deal of interest in the last decade because of its superluminal motion and broadband radio-to-X-ray outbursts, which may be further connected to the gamma-ray productions up to TeV energies. Over the last five years, we have been doing an intensive monitoring of HST-1 with EVN at 5GHz in order to examine the detailed structural evolution and its possible connection to high-energy activities. While this program already yielded interesting results in terms of the detailed mas-scale structure, proper motion measurements and structural variations, the recent HST-1 brightness is continuously decreasing at this frequency. To counter this, we have shifted our monitoring frequency to 1.7GHz from October 2013. This strategy successfully recovered the fainter emission that was missed in the last 5GHz session. Moreover, we again discovered the sudden emergence of a new component at the upstream edge of HST-1, demonstrating that the use of EVN 1.7GHz is indeed powerful to probe the current weak nature of HST-1. Here we report early results from the 1.7GHz monitoring as well as further progress on the long-term kinematic study.

The gas cloud G2 falling toward Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, is supposed to provide valuable information on the physics of accretion flows and the environment of the black hole. We observed Sgr A* with four European stations of the Global Millimeter Very Long Baseline Interferometry Array (GMVA) at 86 GHz on 1 October 2013 when parts of G2 had already passed the pericenter. We searched for possible transient asymmetric structure -- such as jets or winds from hot accretion flows -- around Sgr A* caused by accretion of material from G2. The interferometric closure phases remained zero within errors during the observation time. We thus conclude that Sgr A* did not show significant asymmetric (in the observer frame) outflows in late 2013. Using simulations, we constrain the size of the outflows that we could have missed to ~2.5 mas along the major axis, ~0.4 mas along the minor axis of the beam, corresponding to approximately 232 and 35 Schwarzschild radii, respectively; we thus probe spatial scales on which the jets of radio galaxies are suspected to convert magnetic into kinetic energy. As probably less than 0.2 Jy of the flux from Sgr A* can be attributed to accretion from G2, one finds an effective accretion rate eta*Mdot < 1.5*10^9 kg/s ~ 7.7*10^-9 Mearth/yr for material from G2. Exploiting the kinetic jet power--accretion power relation of radio galaxies, one finds that the rate of accretion of matter that ends up in jets is limited to Mdot < 10^17 kg/s ~ 0.5 Mearth/yr, less than about 20% of the mass of G2. Accordingly, G2 appears to be largely stable against loss of angular momentum and subsequent (partial) accretion at least on time scales < 1 year.

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

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 investigate the position of the radio core in a blazar by multi-epoch astrometric observations at 43 GHz. Using the VLBI Exploration of Radio Astrometry (VERA), we have conducted four adjacent observations in February 2011 and another four in October 2011, and succeeded in measuring the position of the radio core in the TeV blazar Mrk 501 relative to a distant compact quasar NRAO 512. During our observations, we find that (1) there is no positional change within ~0.2 mas or ~2.0 pc de-projected with $\pm1\sigma$ error for the weighted-mean phase-referenced positions of Mrk 501 core relative to NRAO 512 over four adjacent days, and (2) there is an indication of position change for 3C 345 core relative to NRAO 512. By applying our results to the standard internal shock model for blazars, we constrain the bulk Lorenz factors of the ejecta.

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

We analyzed data in polarized intensity obtained with the Very Long Baseline Array (VLBA) at twelve epochs (one observation per month from January to December 2011) at 15, 24, and 43 GHz. For the absolute orientation of the electric vector position angles (EVPA) we used the D-terms method. We also used gamma-ray data from the Fermi Large Area Telescope on weekly time bins throughout 2011. The source shows polarized emission, and its properties vary with time, frequency, and location along the jet. The core mean polarization fraction is generally between 1% and 2%, with a 4% peak at 43 GHz in March; the polarization angle is variable, mainly at 15 GHz, where it changes frequently, and less so at 43 GHz, where it oscillates in the range 114 - 173 deg. The jet polarization properties are more stable, with a fractional polarization of around 16% and a polarization angle nearly perpendicular to the jet axis. The average flux and photon index at gamma-ray energies are (17.7+-0.5)x10^8 ph cm^-2 s^-1 and the photon index is 1.77+-0.02. The gamma-ray light curve shows variability, with a main peak that appears to be associated with the peak in the core polarized emission at 43 GHz, as well as with the total intensity light curve. A discrete correlation function analysis yields a correlation coefficient of 0.54 at zero delay, with a significance level above 99.7%. We accurately determine the polarization properties of Mrk 421, both in the core and in the jet region. The radio and gamma-ray light curves are correlated. The observed EVPA variability at 15 GHz is partly due to opacity and partly to a variable Faraday rotation effect. To explain the residual variability of the intrinsic polarization angle and the low degree of polarization in the core region, we invoke a blend of variable cross-polarized subcomponents with different polarization properties within the beam.

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

We explore energy densities of magnetic field and relativistic electrons in the M87 jet. Since the radio core at the jet base is identical to the optically thick surface against synchrotron self absorption (SSA), the observing frequency is identical to the SSA turnover frequency. As a first step, we assume the radio core as a simple uniform sphere geometry. Using the observed angular size of the radio core measured by the Very Long Baseline Array at 43 GHz, we estimate the energy densities of magnetic field ($U_{B}$) and relativistic electrons ($U_{e}$) based on the standard SSA formula. Imposing the condition that the Poynting power and relativistic electron one should be smaller than the total power of the jet, we find that (i) the allowed range of the magnetic field strength ($B_{tot}$) is from 1 G to 15 G, and that (ii) $1 times 10^{-5} < U_{e}/U_{B} < 6 times 10^{2}$ holds. The uncertainty of $U_{e}/U_{B}$ comes from the strong dependence on the angular size of the radio core and the minimum Lorentz factor of non-thermal electrons ($gamma_{e,min}$) in the core. It is still open that the resultant energetics is consistent with either the magnetohydrodynamic jet or with kinetic power dominated jet even on ~10 Schwarzschild radii scale.

The gamma-ray-detected blazar 3C 454.3 exhibits dramatic flux and polarization variations in the optical and near-infrared bands. In December 2010, the object emitted a very bright outburst. We monitored it for approximately four years (including the 2010 outburst) by optical and near-infrared photopolarimetry. During the 2010 outburst, the object emitted two rapid, redder brightenings, at which the polarization degrees (PDs) in both bands increased significantly and the bands exhibited a frequency-dependent polarization. The observed frequency-dependent polarization leads us to propose that the polarization vector is composed of two vectors. Therefore, we separate the observed polarization vectors into short and long-term components that we attribute to the emissions of the rapid brightenings and the outburst that varied the timescale of days and months, respectively. The estimated PD of the short-term component is greater than the maximum observed PD and is close to the theoretical maximum PD. We constrain the bulk Lorentz factors and inclination angles between the jet axis and the line of sight from the estimated PDs. In this case, the inclination angle of the emitting region of short-term component from the first rapid brightening should be equal to 90$^{\circ}$, because the estimated PD of the short-term component was approximately equal to the theoretical maximum PD. Thus, the Doppler factor at the emitting region of the first rapid brightening should be equal to the bulk Lorentz factor.

We present a study of sub-pc scale radio structure of the radio galaxy 3C 84/NGC 1275 based on the Very Long Baseline Array (VLBA) data at 43 GHz. We discover a limb-brightening in the "restarted" jet associated with the 2005 radio outburst. In the 1990s, the jet structure was ridge-brightening rather than limb-brightening, despite the observations being done with similar angular resolution. This indicates that the transverse jet structure has changed recently. This change in the morphology shows an interesting agreement with the $\gamma$-ray flux increase, i.e., the $\gamma$-ray flux in 1990s was at least seven times lower than the current one. One plausible explanation for the limb-brightening is the velocity structure of the jet in the context of the stratified jet, which is a successful scenario to explain the $\gamma$-ray emission in some active galactic nuclei (AGNs). If this is the case, the change in apparent transverse structure might be caused by the change in the transverse velocity structure. We argue the possibility that the transition from ridge-brightening to limb-brightening is related to the $\gamma$-ray time variability on the timescale of decades. We also discuss the collimation profile of the jet.

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 report the results obtained for the AGN Markarian 421 by model-fitting the data in the visibility plane, studing the proper motion of jet components, the light curve, and the spectral index of the jet features. We compare the radio data with optical light curves obtained at the Steward Observatory, considering also the optical polarization information. Mrk 421 has a bright nucleus and a one-sided jet extending towards the north-west for a few parsecs. The model-fits show that brightness distribution is well described using 6-7 circular Gaussian components, four of which are reliably identified at all epochs; all components are effectively stationary except for component D, at ~0.4 mas from the core, whose motion is however subluminal. Analysis of the light curve shows two different states, with the source being brighter and more variable in the first half of 2011 than in the second half. The highest flux density is reached in February. A comparison with the optical data reveals an increase of the V magnitude and of the fractional polarization simultaneous with the enhancement of the radio activity.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We investigated the detailed inner jet structure of M87 using the Very Long Baseline Array data at 2, 5, 8.4, 15, 23.8, 43, and 86 GHz, especially focusing on the multi-frequency properties of the radio core at the jet base. First, we measured a size of the core region transverse to the jet axis, defined as $W_{\rm c}$, at each frequency \nu, and found a relation between $W_{\rm c}$ and $\nu$ as $W_{\rm c}(\nu) \propto \nu^{-0.71\pm0.05}$. Then, by combining $W_{\rm c}(\nu)$ and the frequency dependence of the core position $r_{\rm c}(\nu)$, which was obtained by our previous study, we have constructed a collimation profile of the innermost jet $W_{\rm c}(r)$ down to ~10 Schwarzschild radii ($R_{\rm s}$) from the central black hole. We found that $W_{\rm c}(r)$ smoothly connects with the width profile of the outer edge-brightened, parabolic jet, and then follows a similar radial dependence down to several tens of $R_{\rm s}$. Closer to the black hole, the measured radial profile suggests a possible change of the jet collimation shape from the outer parabolic one, in which the jet shape tends to become more radially-oriented. This could be related to a magnetic collimation process or/and interaction with surrounding materials at the jet base. The present results shed light on the importance of higher-sensitivity/resolution imaging studies for M87 at 86, 43 and also 22 GHz, and should be examined more rigorously.

We explore non-thermal emission from a shocked interstellar medium, which is identified as an expanding shell, driven by a relativistic jet in active galactic nuclei (AGNs). In this work, we particularly focus on parsec-scale size mini-shells surrounding mini radio lobes. From radio to X-ray band, the mini radio lobe emission dominates the faint emission from the mini-shell. On the other hand, we find that inverse-Compton (IC) emission from the shell can overwhelm the associated lobe emission at very high energy (VHE; E>100 GeV) gamma-ray range, because energy densities of synchrotron photons from the lobe and/or soft photons from the AGN nucleus are large and IC scattering effectively works. The predicted IC emission from nearby mini-shells can be detected with the Cherencov Telescope Array (CTA) and they are potentially a new class of VHE gamma-ray emitters.

3C 66A is an intermediate-frequency-peaked BL Lac object detected by the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope. We present a study of the long-term variations of this blazar seen over 2 years at GeV energies with Fermi and in the optical (flux and polarization) and near infrared with the Kanata telescope. In 2008, the first year of the study, we find a correlation between the gamma-ray flux and the measurements taken with the Kanata telescope. This is in contrast to the later measurements performed during 2009--2010 which show only a weak correlation along with a gradual increase of the optical flux. We calculate an external seed photon energy density assuming that the gamma-ray emission is due to external Compton scattering. The energy density of the external photons is found to be higher by a factor of two in 2008 compared to 2009--2010. We conclude that the different behaviors observed between the first year and the later years might be explained by postulating two different emission components.

We investigate observational properties of a kilo-parsec scale knot in radio-loud quasar 3C 380 by using two epoch archival data obtained by Very Long Baseline Interferometry (VLBI) at 5 GHz on 1998 July and 2001 April. We succeed in obtaining the highest spatial resolution image of the bright knot K1 located at 732 milliarcseconds, or more than 20 kpc de-projected, downstream from the nucleus three times better than previously obtained highest resolution image by Papageorgiou et al. (2006). Our images reveal, with new clarity, "inverted bow-shock" structure in K1 facing the nucleus and its morphology resembles a conical shock wave. By comparing the two epoch images directly, we explore the kinematics of K1 and obtain the upper limit of apparent velocity, 0.25 mas/yr or 9.8 c of K1 for the first time. The upper limit of apparent velocity is marginally smaller than superluminal motions seen in the core region. Further new epoch VLBI observations are necessary to measure the proper motion at K1.

We report on the detailed radio status of the M87 jet during the Very-High-Energy (VHE) gamma-ray flaring event in April 2010, obtained from high-resolution, multi-frequency, phase-referencing VLBA observations. We especially focus on the properties for the jet base (the radio core) and the peculiar knot HST-1, which are currently favored as the gamma-ray emitting sites. During the VHE flaring event, the HST-1 region remains stable in terms of its structure and flux density in the optically thin regime above 2GHz, being consistent with no signs of enhanced activities reported at X-ray for this feature. The radio core shows an inverted spectrum at least up to 43GHz during this event. Astrometry of the core position, which is specified as ~20Rs from the central engine in our previous study, shows that the core position is stable on a level of 4Rs. The core at 43 and 22GHz tends to show slightly (~10%) higher flux level near the date of the VHE flux peak compared with the epochs before/after the event. The size of the 43-GHz core is estimated to be ~17Rs, which is close to the size of the emitting region suggested from the observed time scale of rapid variability at VHE. These results tend to favor the scenario that the VHE gamma-ray flare in 2010 April is associated with the radio core.

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

We introduce the GENJI program (Gamma-ray Emitting Notable AGN Monitoring by Japanese VLBI), which is a monitoring program of gamma-ray bright AGNs with the VERA array (VLBI Exploration of Radio Astrometry). The GENJI programme aims a dense monitoring at 22 GHz towards the $\gamma$-ray emitting active galactic nuclei (AGNs) to investigate the radio time variation of the core and possible ejection of new radio component, motion of jets, and their relation with the emission at other wavelengths especially in $\gamma$-rays. Currently we are monitoring 8 $\gamma$-ray-emitting notable AGNs (DA 55, 3C 84, M 87, PKS 1510-089, DA 406, NRAO 530, BL Lac, 3C 454.3) about once every two weeks. This programme is promising to trace the trend of radio time variation on shorter timescale than conventional VLBI monitoring programme and to provide complimentary data with them (e.g., MOJAVE, Boston University Blazar Project). In particular, we successfully coordinated quick follow-up observations after the GeV $\gamma$-ray flare in NRAO 530 and 3C 454.3 reported by the Fermi Gamma-ray Space Telescope. Here we present the initial results of morphology and light curves for the first 7-month operation.

We report the finding of kiloparsec (kpc)-scale radio structures in three radio-loud narrow-line Seyfert 1 (NLS1) galaxies from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) of the Very Large Array (VLA), which increases the number of known radio-loud NLS1s with kpc-scale structures to six, including two gamma-ray emitting NLS1s (PMN J0948+0022 and 1H 0323+342) detected by the Fermi Gamma-ray Space Telescope. The detection rate of extended radio emissions in NLS1s is lower than that in broad-line active galactic nuclei (AGNs) with a statistical significance. We found both core-dominated (blazar-like) and lobe-dominated (radio-galaxy-like) radio structures in these six NLS1s, which can be understood in the framework of the unified scheme of radio-loud AGNs that considers radio galaxies as non-beamed parent populations of blazars. Five of the six NLS1s have (i) extended radio luminosities suggesting jet kinetic powers of >~10^44 erg/s, which is sufficient to make jets escape from hosts' dense environments, (ii) black holes of >~10^7 solar mass, which can generate the necessary jet powers from near-Eddington mass accretion, and (iii) two-sided radio structures at kpc scales, requiring expansion rates of ~0.01c--0.3c and kinematic ages of >~10^7 years. On the other hand, most typical NLS1s would be driven by black holes of <~10^7 solar mass in a limited lifetime of ~10^7 years. Hence the kpc-scale radio structures may originate in a small window of opportunity during the final stage of the NLS1 phase just before growing into broad-line AGNs.

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

We report on the VLBI follow-up observations using the Japanese VLBI Network (JVN) array at 22 GHz for the largest X-ray flare of TeV blazar Mrk 421 that occurred in mid-February, 2010. The total of five epochs of observations were performed at intervals of about 20 days between March 7 and May 31, 2010. No new-born component associated with the flare was seen directly in the total intensity images obtained by our multi-epoch VLBI observations. However, one jet component located at ~1 mas north-west from the core was able to be identified, and its proper motion can be measured as -1.66+/-0.46 mas yr^-1, which corresponds to an apparent velocity of -3.48+/-0.97 c. Here, this negative velocity indicates that the jet component was apparently moving toward the core. As the most plausible explanation, we discuss that the apparent negative velocity was possibly caused by the ejection of a new component, which could not be resolved with our observations. In this case, the obtained Doppler factor of the new component is around 10 to 20, which is consistent with the ones typically estimated by model fittings of spectral energy distribution for this source.

High-resolution radio observations are ideal for constraining the value of physical parameters in the inner regions of active-galactic-nucleus jets and complement results on multiwavelength (MWL) observations. This study is part of a wider multifrequency campaign targeting the nearby TeV blazar Markarian 421 (z=0.031), with observations in the sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and gamma rays (Fermi-LAT, MAGIC, VERITAS). We investigate the jet's morphology and any proper motions, and the time evolution of physical parameters such as flux densities and spectral index. The aim of our wider multifrequency campaign is to try to shed light on questions such as the nature of the radiating particles, the connection between the radio and gamma-ray emission, the location of the emitting regions and the origin of the flux variability. We consider data obtained with the Very Long Baseline Array (VLBA) over twelve epochs (one observation per month from January to December 2011) at 15 GHz and 24 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch. The structure of Mrk 421 is dominated by a compact (~0.13 mas) and bright component, with a one-sided jet detected out to ~10 mas. We identify 5-6 components in the jet that are consistent with being stationary during the 12-month period studied here. Measurements of the spectral index agree with those of other works: they are fairly flat in the core region and steepen along the jet length. Significant flux-density variations are detected for the core component. From our results, we draw an overall scenario in which we estimate a viewing angle 2{\deg} < theta < 5{\deg} and a different jet velocity for the radio and the high-energy emission regions, such that the respective Doppler factors are {\delta}r ~3 and {\delta}h.e. ~14.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We present a preliminary analysis of new high resolution radio observations of the nearby TeV blazar Markarian 421 (z=0.031). This study is part of an ambitious multifrequency campaign, with observations in sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and gamma rays (Fermi-LAT, MAGIC, VERITAS). In this manuscript we consider only data obtained with the Very Long Baseline Array (VLBA) at seven epochs (one observation per month from January to July 2011) at 15 and 23.8 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch. We identified 5-6 components which are consistent with being stationary during the 6-month period reported here. The aim is to try to shed light on questions such as the nature of radiating particles, the connection between radio and gamma-ray emission, the location of emitting regions and the origin of the flux variability.

We examine plasma composition of jets in active galactic nuclei through the comparison of the total pressure ($P$) with partial pressures of electrons and protons in a cocoon. The total pressure is estimated from the analysis of an expanding cocoon dynamics. We determine the average kinetic energy per particle for several representative cases of particle energy distribution such as one- and two-temperature thermal plasmas and non-thermal electrons by evaluating the dissipation of total kinetic energy of the jet into the internal energy of cocoon plasma. The number density of the total electrons/positrons ($n_{\pm}$) in the cocoon is constrained by using the particle supply from hot spots and the absence of thermal bremsstrahlung emission from radio lobes. By inserting $P$, $n_{\pm}$ and the particle energy of each population into the equation of state, the number density ($n_{p}$) and pressure ($P_{p}$) of protons in the cocoon can be constrained. Applying this method to Cygnus A, we find that (i) electron/positron ($e^{\pm}$) pairs always dominate in terms of number density, but that (ii) either an "$e^{\pm}$-supported cocoon (i.e., $P_{\pm} >P_{p}$)" or "proton-supported one (i.e, $P_{\pm} <P_{p}$)" is possible.

The radio galaxy 3C 84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGNs) and one of the best laboratories to study the radio properties of the sub-pc jet in connection with the gamma-ray emission. In order to identify possible radio counterparts of the gamma-ray emissions in 3C 84, we study the change in structure within the central 1 pc and the light curve of sub-pc-size components C1, C2, and C3. We search for any correlation between changes in the radio components and the gamma-ray flares by making use of VLBI and single dish data. Throughout the radio monitoring spanning over two GeV gamma-ray flares detected by the {\it Fermi}-LAT and the MAGIC Cherenkov Telescope in the periods of 2009 April to May and 2010 June to August, total flux density in radio band increases on average. This flux increase mostly originates in C3. Although the gamma-ray flares span on the timescale of days to weeks, no clear correlation with the radio light curve on this timescale is found. Any new prominent components and change in morphology associated with the gamma-ray flares are not found on the VLBI images.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

Aims: We aim to constrain the structural variations within the HST-1 region downstream of the radio jet of M87, in general as well as in connection to the episodes of activity at very high energy (VHE). Methods: We analyzed and compared 26 VLBI observations of the M87 jet, obtained between 2006 and 2011 with the Very Long Baseline Array (VLBA) at 1.7 GHz and the European VLBI Network (EVN) at 5 GHz. Results: HST-1 is detected at all epochs; we model-fitted its complex structure with two or more components, the two outermost of which display a significant proper motion with a superluminal velocity around ~4c. The motion of a third feature that is detected upstream is more difficult to characterize. The overall position angle of HST-1 has changed during the time of our observations from -65deg to -90deg, while the structure has moved by over 80 mas downstream. Our results on the component evolution suggest that structural changes at the upstream edge of HST-1 can be related to the VHE events.

In December 2009, the bright blazar, 3C 454.3 exhibited a strong outburst in the optical, X-ray and gamma-ray regions. We performed photometric and polarimetric monitoring of this outburst in the optical and near-infrared bands with TRISPEC and HOWPol attached to the Kanata telescope. We also observed this outburst in the infrared band with AKARI, and the radio band with the 32-m radio telescope of Yamaguchi University. The object was in an active state from JD 2455055 to 2455159. It was 1.3 mag brighter than its quiescent state before JD 2455055 in the optical band. After the end of the active state in JD 2455159, a prominent outburst was observed in all wavelengths. The outburst continued for two months. Our optical and nearinfrared polarimetric observations revealed that the position angle of the polarization (PA) apparently rotated clockwise by 240 degrees within 11 d in the active state (JD 2455063-2455074), and after this rotation, PA remained almost constant during our monitoring. In the outburst state, PA smoothly rotated counterclockwise by 350 degrees within 35 d (JD 2455157-2455192). Thus, we detected two distinct rotation events of polarization vector in opposite directions. We discuss these two events compared with the past rotation events observed in 2005, 2007 and 2008.

Following the discovery of a new radio component right before the GeV \gamma-ray detection since 2008 August by Fermi Gamma-ray Space Telescope, we present a detailed study of the kinematics and lightcurve on the central sub-pc scale of 3C 84 using the archival VLBA 43-GHz data covering the period between 2002 January to 2008 November. We find that the new component "C3", previously reported by the observations with the VLBI Exploration of Radio Astrometry (VERA), was already formed in 2003. The flux density of C3 increases moderately until 2008, and then it becomes brighter rapidly after 2008. The radio core, C1, also shows a similar trend. The apparent speed of C3 with reference to the core C1 shows moderate acceleration from 0.10c to 0.47c between 2003 November to 2008 November, but is still sub-relativistic. We further try to fit the observed broadband spectrum by the one-zone synchrotron self-Compton (SSC) model using the measured apparent speed of C3. The fit can reproduce the observed \gamma-ray emission, but does not agree with the observed radio spectral index between 22 and 43 GHz.

We report on the correlation between the flux, color and polarization variations on time scales of days--months in blazars, and discuss their universal aspects. We performed monitoring of 42 blazars in the optical and near-infrared bands from 2008 to 2010 using TRISPEC attached to the "Kanata" 1.5-m telescope. We found that 28 blazars exhibited "bluer-when-brighter" trends in their whole or a part of time-series data sets. This corresponds to 88% of objects that were observed for >10 days. Thus, our observation unambiguously confirmed that the "bluer-when-brighter" trend is common in the emission from blazar jets. This trend was apparently generated by a variation component with a constant and relatively blue color and an underlying red component. Prominent short-term flares on time scales of days--weeks tended to exhibit a spectral hysteresis; their rising phases were bluer than their decay phases around the flare maxima. In contrast to the strong flux--color correlation, the correlation of the flux and polarization degree was relatively weak; only 10 objects showed significant positive correlations. Rotations of polarization were detected only in three objects: PKS 1510-089, 3C 454.3, and PKS 1749+096, and possibly in S5 0716+714. We also investigated the dependence of the degree of variability on the luminosity and the synchrotron peak frequency, \nu_peak. As a result, we found that lower luminosity and higher \nu_peak objects had smaller variations in their amplitudes both in the flux, color, and polarization degree. Our observation suggests the presence of several distinct emitting sources, which have different variation time-scales, colors, and polarizations. We propose that the energy injection by, for example, internal shocks in relativistic shells is a major factor for blazar variations on time scales of both days and months.

There is growing evidence of relativistic jets in radio-loud narrow-line Seyfert 1 (RL-NLS1) galaxies. We constrain the observational properties of the radio emission in the first RL-NLS1 galaxy ever detected in gamma-rays, PMN J0948+0022, i.e., its flux density and structure in total intensity and in polarization, its compactness, and variability. We performed three real-time e-VLBI observations of PMN J0948+0022 at 22 GHz, using a global array including telescopes in Europe, East Asia, and Australia. These are the first e-VLBI science observations ever carried out with a global array, reaching a maximum baseline length of 12458 km. The observations were part of a large multiwavelength campaign in 2009. The source is detected at all three epochs. The structure is dominated by a bright component, more compact than 55 microarcsec, with a fainter component at a position angle theta~ 35deg. Relativistic beaming is required by the observed brightness temperature of 3.4x10^11 K. Polarization is detected at a level of about 1%. The parameters derived by the VLBI observations, in addition to the broad-band properties, confirm that PMN J0948+0022 is similar to flat spectrum radio quasars. Global e-VLBI is a reliable and promising technique for future studies.

We report on multi-band photometric and polarimetric observations of the blazars AO 0235+164 and PKS 1510-089. These two blazars were active in 2008 and 2009, respectively. In these active states, prominent short flares were observed in both objects, having amplitudes of >1 mag within 10 d. The $V-J$ color became bluer when the objects were brighter in these flares. On the other hand, the color of PKS 1510-089 exhibited a trend that it became redder when it was brighter, except for its prominent flare. This redder-when-brighter trend can be explained by the strong contribution of thermal emission from an accretion disk. The polarization degree increased at the flares, and reached >25 % at the maxima. We compare these flares in AO 0235+164 and PKS 1510-089 with other short flares which were detected by our monitoring of 41 blazars. Those two flares had one of the largest variation amplitudes in both flux and polarization degree. Furthermore, we found a significant positive correlation between the amplitudes of the flux and polarization degree in the short flares. It indicates that the short flares originate from the region where the magnetic field is aligned.

We explore the evolution of the emissions by accelerated electrons in shocked shells driven by jets in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, we evaluated the broadband emission spectra by properly taking into account adiabatic and radiative cooling effects on the electron distribution. The synchrotron radiation and inverse Compton (IC) scattering of various photons that are mainly produced in the accretion disc and dusty torus are considered as radiation processes. We show that the resultant radiation is dominated by the IC emission for compact sources (< 10kpc), whereas the synchrotron radiation is more important for larger sources. We also compare the shell emissions with those expected from the lobe under the assumption that a fractions of the energy deposited in the shell and lobe carried by the non-thermal electrons are $\epsilon_e \sim 0.01$ and $\epsilon_{e, lobe} \sim 1$, respectively. Then, we find that the shell emissions are brighter than the lobe ones at infra-red and optical bands when the source size is > 10kpc, and the IC emissions from the shell at > 10 GeV can be observed with the absence of contamination from the lobe irrespective of the source size. In particular, it is predicted that, for most powerful nearby sources ($L_j \sim 10^{47} ergs s^{-1}$), TeV gamma-rays produced via the IC emissions can be detected by the modern Cherenkov telescopes such as MAGIC, HESS and VERITAS.

Recent radio observations reveal the existence of mini radio lobes in active galaxies with their scales of $\sim 10 {\rm pc}$. The lobes are expected to be filled with shock accelerated electrons and protons. In this work, we examine the photon spectra from the mini lobes, properly taking the hadronic processes into account. We find that the resultant broadband spectra contain the two distinct hadronic bumps in $\gamma$-ray bands, i.e., the proton synchrotron bump at $\sim$ MeV and the synchrotron bump at $\sim$ GeV due to the secondary electrons/positrons produced via photo-pion cascade. Especially when the duration of particle injection is shorter than the lobe age, radio-dark $\gamma$-ray lobes are predicted. The existence of the $\gamma$-ray lobes could be testable with the future TeV-$\gamma$ telescope {\it CTA}.

The detection of gamma-ray emission by Fermi-LAT from the radio loud Narrow Line Seyfert 1 PMN J0948+0022 (Abdo et al. 2009, ApJ 699, 976) triggered a multi-wavelength campaign between March and July 2009. Given its high compactness (Doi et al. 2006, PASJ 58, 829), inverted spectrum, and 0deg declination, the source was an ideal target to observe at 22 GHz with a Global VLBI array extending from Europe to East Asia and Australia. In order to deliver prompt results to be analysed in combination with the other instruments participating in the campaign, the observations were carried out with real time VLBI, for the first time on a Global scale. Indeed, the main results have been published just a few months after the campaign (Abdo et al. 2009, ApJ 707, 727). Here we present additional details about the e-VLBI observations.

Aims. We report on a detailed study of the optical afterglow of GRB 061121 with our original time-series photometric data. In conjunction with X-ray observations, we discuss the origin of its optical and X-ray afterglows. Methods. We observed the optical afterglow of Swift burst GRB 061121 with the Kanata 1.5-m telescope at Higashi-Hiroshima Observatory. Our observation covers a period just after an X-ray plateau phase. We also performed deep imaging with the Subaru telescope in 2010 in order to estimate the contamination of the host galaxy. Results. In the light curve, we find that the optical afterglow also exhibited a break as in the X-ray afterglow. However, our observation suggests a possible hump structure or a flattening period before the optical break in the light curve. There is no sign of such a hump in the X-ray light curve. Conclusions. This implies that the emitting region of optical was distinct from that of X-rays. The hump in the optical light curve was possibly caused by the passage of the typical frequency of synchrotron emission from another forward shock distinct from the early afterglow. The observed decay and spectral indices are inconsistent with the standard synchrotron-shock model. Hence, the observation requires a change in microphysical parameters in the shock region or a prior activity of the central engine. Alternatively, the emission during the shallow decay phase may be a composition of two forward shock emissions, as indicated by the hump structure in the light curve.

We present optical and near-infrared (NIR) photometry of a classical nova, V2362 Cyg (= Nova Cygni 2006). V2362 Cyg experienced a peculiar rebrightening with a long duration from 100 to 240 d after the maximum of the nova. Our multicolor observation indicates an emergence of a pseudophotosphere with an effective temperature of 9000 K at the rebrightening maximum. After the rebrightening maximum, the object showed a slow fading homogeneously in all of the used bands for one week. This implies that the fading just after the rebrightening maximum ( less or equal 1 week ) was caused by a slowly shrinking pseudophotosphere. Then, the NIR flux drastically increased, while the optical flux steeply declined. The optical and NIR flux was consistent with blackbody radiation with a temperature of 1500 K during this NIR rising phase. These facts are likely to be explained by dust formation in the nova ejecta. Assuming an optically thin case, we estimate the dust mass of 10^(-8) -- 10^(-10) M_solar, which is less than those in typical dust-forming novae. These results support the senario that a second, long-lasting outflow, which caused the rebrightening, interacted with a fraction of the initial outflow and formed dust grains.

Hadronic emission from parsec size radio lobes in active galactic nuclei (AGN) is discussed. The lobes are composed of shocked jet plasmaand expected to be filled with high energy particles. By using the Monte Carlo simulation, we calculate the photon spectra from the lobes including photo-meson interaction processes. When the synchrotron emission from primary electrons is bright, synchrotron-self-Compton component is dominant in gamma-ray bands. The hadronic emission from the lobes can be dominated in gamma-ray bands when the primary emission is not very bright. Proton synchrotron component arises at sub MeV band. The synchrotron emission radiated from secondary electron/positron pairs produced via photo-meson cascade emerges in GeV-TeV energy ranges. These high energy emission signatures provide a test for proton accelerations in young AGN jets.

We report multi-epoch Very Long Baseline Interferometry (VLBI) study of the innermost jet of 3C 84. We carried out 14-epoch VLBI observations during 2006-2009 with the Japanese VLBI Network (JVN) and the VLBI Exploration of Radio Astrometry (VERA), immediately following the radio outburst that began in 2005. Comparison between VLBI lightcurve and single-dish lightcurve indicates that the outburst was associated with the central ~1 pc core. We found that this outburst accompanied the emergence of a new component, and the projected speed of this new component was 0:23c from 2007/297 (2007 October 24) to 2009/114 (2009 April 24). We argue the site of gamma-ray emission detected by Fermi/LAT and jet kinematics in connection with gamma-ray emission mechanism.

We investigate the origin of short X-ray flares which are occasionally observed in early stages of afterglows of gamma-ray bursts (GRBs). We observed two events, GRB 071112C and GRB 080506, before the start of X-ray flares in the optical and near-infrared (NIR) bands with the 1.5-m Kanata telescope. In conjunction with published X-ray and optical data, we analyzed densely sampled light curves of the early afterglows and spectral energy distributions (SEDs) in the NIR-X-ray ranges. We found that the SEDs had a break between the optical and X-ray bands in the normal decay phases of both GRBs regardless of the model for the correction of the interstellar extinction in host galaxies of GRBs. In the X-ray flares, X-ray flux increased by 3 and 15 times in the case of GRB 071112C and 080506, respectively, and the X-ray spectra became harder than those in the normal decay phases. No significant variation in the optical-NIR range was detected together with the X-ray flares. These results suggest that the X-ray flares were associated with either late internal shocks or external shocks from two-component jets.

We explore the emissions by accelerated electrons in shocked shells driven by jets in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, the synchrotron radiation and inverse Compton (IC) scattering of various photons that are mainly produced in the core are considered as radiation processes. We show that the radiative output is dominated by the IC emission for compact sources (< 30kpc), whereas the synchrotron radiation is more important for larger sources. It is predicted that, for powerful sources ($L_j \sim 10^{47} ergs s^{-1}$), GeV-TeV gamma-rays produced via the IC emissions can be detected by the Fermi satellite and modern Cherenkov telescopes such as MAGIC, HESS and VERITAS if the source is compact.

Using deep Chandra ACIS observation data for Cygnus A, we report evidence of non-thermal X-ray emission from radio lobes surrounded by a rich intra-cluster medium (ICM). The diffuse X-ray emission, which are associated with the eastern and western radio lobes, were observed in a 0.7--7 keV Chandra$ ACIS image. The lobe spectra are reproduced with not only a single-temperature Mekal model, such as that of the surrounding ICM component, but also an additional power-law (PL) model. The X-ray flux densities of PL components for the eastern and western lobes at 1 keV are derived as 77.7^{+28.9}_{-31.9} nJy and 52.4^{+42.9}_{-42.4} nJy, respectively, and the photon indices are 1.69^{+0.07}_{-0.13} and 1.84^{+2.90}_{-0.12}, respectively. The non-thermal component is considered to be produced via the inverse Compton (IC) process, as is often seen in the X-ray emission from radio lobes. From a re-analysis of radio observation data, the multiwavelength spectra strongly suggest that the seed photon source of the IC X-rays includes both cosmic microwave background radiation and synchrotron radiation from the lobes. The derived parameters indicate significant dominance of the electron energy density over the magnetic field energy density in the Cygnus A lobes under the rich ICM environment.

We report on optical-near-infrared photopolarimetric observations of a blazar 3C 454.3 over 200 d. The object experienced an optical outburst in July 2007. This outburst was followed by a short state fainter than $V=15.2$ mag lasting $\sim 25$ d. The object, then, entered an active state during which we observed short flares having a timescale of 3-10 d. The object showed two types of features in the color-magnitude relationship. One is a "bluer-when-brighter" trend in the outburst state, and the other is a "redder-when-brighter" trend in the faint state. These two types of features suggest a contribution of a thermal emission to the observed flux, as suspected in previous studies. Our polarimetric observation detected two episodes of the rotation of the polarization vector. The first one was a counterclockwise rotation in the $QU$ plane during the outburst state. After this rotation event of the polarization vector, the object entered a rapidly fading stage. The second one was seen in a series of flares during the active state. Each flare had a specific position angle of polarization, and it apparently rotated clockwise from the first to the last flares. Thus, the object exhibited rotations of the polarization vector in opposite directions. We estimated a decay timescale of the short flares during the active state, and then calculated an upper limit of the strength of the magnetic field, $B$=0.2 G, assuming a typical beaming factor of blazars, $\delta=20$. This upper limit of $B$ is smaller than those previously estimated from spectral analysis.

Multi-epoch Very Long Baseline Interferometry (VLBI) study of the sub-pc scale jet of 3C 84 is presented. We carried out 14-epoch VLBI observations during 2006-2009 with the Japanese VLBI Network (JVN) and the VLBI Exploration of Radio Astrometry (VERA), immediately following the radio outburst that began in 2005. We confirmed that the outburst was associated with the central ~1 pc core, accompanying the emergence of a new component. This is striking evidence of the recurrence of jet activity. The new component became brighter during 2008, in contrast to the constant gamma-ray emission that was observed with the Fermi Gamma-ray Space Telescope during the same time. We found that the projected speed of the new component is 0.23c from 2007/297 (2007 October 24) to 2009/114 (2009 April 24). The direction of movement of this component differs from that of the pre-existing component by ~40 degree. This is the first measurement of kinematics of a sub-pc jet in a gamma-ray active phase. Possible detection of jet deceleration and the jet kinematics in connection with the gamma-ray emission is discussed.

We report the observation of a new dwarf nova, OT J055717+683226, during its first-ever recorded superoutburst in December 2006. Our observation shows that this object is an SU UMa-type dwarf nova having a very short superhump period of 76.67+/- 0.03 min (0.05324+/-0.00002 d). The next superoutburst was observed in March 2008. The recurrence time of superoutbursts (supercycle) is, hence, estimated to be ~480 d. The supercycle is much shorter than those of WZ Sge-type dwarf novae having supercycles of >~ 10 yr, which are a major population of dwarf novae in the shortest orbital period regime (<~85 min). Using a hierarchical cluster analysis, we identified seven groups of dwarf novae in the shortest orbital period regime. We identified a small group of objects that have short supercycles, small outburst amplitudes, and large superhump period excesses, compared with those of WZ Sge stars. OT J055717+683226 probably belongs to this group.

We report two different types of backflow from jets by performing 2D special relativistic hydrodynamical simulations. One is anti-parallel and quasi-straight to the main jet (quasi-straight backflow), and the other is bent path of the backflow (bent backflow). We find that the former appears when the head advance speed is comparable to or higher than the local sound speed at the hotspot while the latter appears when the head advance speed is slower than the sound speed bat the hotspot. Bent backflow collides with the unshocked jet and laterally squeezes the jet. At the same time, a pair of new oblique shocks are formed at the tip of the jet and new bent fast backflows are generated via these oblique shocks. The hysteresis of backflow collisions is thus imprinted in the jet as a node and anti-node structure. This process also promotes broadening of the jet cross sectional area and it also causes a decrease in the head advance velocity. This hydrodynamic process may be tested by observations of compact young jets.

We developed a method to separate a long-term trend from observed temporal variations of polarization in blazars using a Bayesian approach. The temporal variation of the polarization vector is apparently erratic in most blazars, while several objects occasionally exhibited systematic variations, for example, an increase of the polarization degree associated with a flare of the total flux. We assume that the observed polarization vector is a superposition of distinct two components, a long-term trend and a short-term variation component responsible for short flares. Our Bayesian model estimates the long-term trend which satisfies the condition that the total flux correlates with the polarized flux of the short-term component. We demonstrate that assumed long-term polarization components are successfully separated by the Bayesian model for artificial data. We applied this method to photopolarimetric data of OJ 287, S5 0716+714, and S2 0109+224. Simple and systematic long-term trends were obtained in OJ 287 and S2 0109+224, while no such a trend was identified in S5 0716+714. We propose that the apparently erratic variations of polarization in OJ 287 and S2 0109+224 are due to the presence of the long-term polarization component. The behavior of polarization in S5 0716+714 during our observation period implies the presence of a number of polarization components having a quite short time-scale of variations.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

There has not been a comprehensive framework for comparing spectral data from different planets.Such a framework is needed for the study of extrasolar planets and objects within the solar system. We have undertaken observations to compile a library of planet spectra for all planets, some moons, and some dwarf planets in the solar system to study their general spectroscopic and photometric natures. During May and November of 2008, we acquired spectra for the planets using TRISPEC, which is capable of simultaneous three-band spectroscopy in a wide wavelength range of 0.45 - 2.5 microns with low resolving power (lambda-over-Delta-lambda is 140 - 360). Patterns emerge from comparing the spectra. Analyzing their general spectroscopic and photometric natures, we show that it is possible to distinguish between gas planets, soil planets and ice planets. These methods can be applied to extrasolar observations using low resolution spectrography or broad-band filters. The present planet spectral library is the first library to contain observational spectra for all of the solar system planets, based on simultaneous observations in visible and near infrared wavelengths. This library will be a useful reference for analyzing extrasolar planet spectra, and for calibrating planetary data sets.

We report on optical and infrared photometric observations of a WZ Sge-type dwarf nova, V455 And during a superoutburst in 2007. These observations were performed with the KANATA (V, J, and K_s bands) and MITSuME (g', Rc, and Ic bands) telescopes. Our 6-band simultaneous observations allowed us to investigate the temporal variation of the temperature and the size of the emitting region associated with the superoutburst and short-term modulations, such as early and ordinary superhumps. A hot (>11000 K) accretion disk suddenly disappeared when the superoutburst finished, while blackbody emission, probably from the disk, still remained dominant in the optical region with a moderately high temperature (~8000 K). This indicates that a substantial amount of gas was stored in the disk even after the outburst. This remnant matter may be a sign of an expected mass-reservoir which can trigger echo outbursts observed in several WZ Sge stars. The color variation associated with superhumps indicates that viscous heating in a superhump source stopped on the way to the superhump maximum, and a subsequent expansion of a low-temperature region made the maximum. The color variation of early superhumps was totally different from that of superhumps: the object was bluest at the early superhump minimum. The temperature of the early superhump light source was lower than that of an underlying component, indicating that the early superhump light source was a vertically expanded low-temperature region at the outermost part of the disk.

We argue that the origin of "FRI/FRI{-.1em}I dichotomy" -- the division between Fanaroff-Riley class I (FRI) with subsonic lobes and class I{-.1em}I (FRI{-.1em}I) radio sources with supersonic lobes is sharp in the radio-optical luminosity plane (Owen-White diagram) -- can be explained by the deceleration of advancing radio lobes. The deceleration is caused by the growth of the effective cross-sectional area of radio lobes. We derive the condition in which an initially supersonic lobe turns into a subsonic lobe, combining the ram-pressure equilibrium between the hot spots and the ambient medium with the relation between "the hot spot radius" and "the linear size of radio sources" obtained from the radio observations. We find that the dividing line between the supersonic lobes and subsonic ones is determined by the ratio of the jet power $L_{\rm j}$ to the number density of the ambient matter at the core radius of the host galaxy $\bar{n}_{\rm a}$. It is also found that there exists the maximal ratio of $(L_{\rm j}/\bar{n}_{\rm a})$ and its value resides in $(L_{\rm j}/\bar{n}_{\rm a})_{\rm max}\approx 10^{44-47} {\rm erg} {\rm s}^{-1} {\rm cm}^{3}$, taking account of considerable uncertainties. This suggests that the maximal value $(L_{\rm j}/\bar{n}_{\rm a})_{\rm max}$ separates between FRIs and FRI{-.1em}Is.

We report detailed, long term near-infrared (NIR) light curves of GRS 1915+105 in 2007-2008, covering its long "soft state" for the first time. From our NIR monitoring and the X-ray data of the All Sky Monitor (ASM) onboard Rossi X-ray Timing Explorer (RXTE), we discovered that the NIR flux dropped by > 1 mag during short X-ray flares with a time-scale of days. With the termination of the soft state, the H-Ks color reddened and the anti-correlation pattern was broken. The observed H-Ks color variation suggests that the dominant NIR source was an accretion disk during the soft state. The short X-ray flares during the soft state were associated with spectral hardening in X-rays and increasing radio emission indicating jet ejection. The temporal NIR fading during the X-ray flares, hence, implies a sudden decrease of the contribution of the accretion disk when the jet is ejected.

We present a new prediction of GeV $\gamma$-ray emission from radio lobes of young AGN jets. In the previous work of Kino et al. (2007), MeV $\gamma$-ray bremsstrahlung emission was predicted from young cocoons/radio-lobes in the regime of no coolings. In this study, we include cooling effects of bremsstrahlung emission and adiabatic loss. With the initial conditions determined by observed young radio lobes, we solve a set of equations describing the expanding lobe evolution. Then we find that the lobes initially have electron temperature of $\sim$GeV, and they cool down to $\sim$MeV by the adiabatic loss. Correspondingly, the lobes initially yield bright bremsstrahlung luminosity in $\sim$GeV range and they fade out. We estimate these $\gamma$-ray emissions and show that nearby young radio lobes could be detected with Fermi Gamma-ray Space Telescope.

Gamma-ray emission from cocoons of young radio galaxies is predicted. Considering the process of adiabatic injection of the shock dissipation energy and mass of the relativistic jet into the cocoon, we find that the thermal electron temperature of the cocoon is typically predicted to be of the order of $\sim$ MeV, and is determined only by the bulk Lorentz factor of the jet. Together with the time-dependent dynamics of the cocoon expansion, we find that young cocoons can yield thermal Bremsstrahlung emissions at energies $\sim$MeV. Hotter cocoons (i.e., GeV) for younger sources are also discussed.

In order to reveal the long-term evolution of relativistic jets in active galactic nuclei (AGNs), we examine the dynamical evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects (MSOs), Fanaroff-Riley type II radio galaxies (FRIIs)]. By comparing the observed relation between the hot spot size and the linear size of radio source with a coevolution model of hot spot and cocoon, we find that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRII phase). The deceleration is caused by the growth of the cross-sectional area of the cocoon head. Moreover, by comparing the hot spot speed with the sound speed of the ambient medium, we predict that only CSOs whose initial advance speed is higher than 0.3-0.5c can evolve into FRIIs.

We present the result of near-infrared and optical observations of the BL Lac object S5 0716$ + $714 carried out by the KANATA telescope. S5 0716$ + $714 has both a long term high-amplitude variability and a short-term variability within a night. The shortest variability (microvariability) time-scale is important for understanding the geometry of jets and magnetic field, because it provides a possible minimum size of variation sources. Here, we report the detection of 15-min variability in S5 0716$ + $714, which is one of the shortest time-scales in optical and near-infrared variations observed in blazars. The detected microvariation had an amplitude of $0.061{\pm}0.005$ mag in $V$ band and a blue color of $\Delta(V-J)=-0.025{\pm}0.011$. Furthermore, we successfully detected an unprecedented, short time-scale polarimetric variation which correlated with the brightness change. We revealed that the microvariation had a specific polarization component. The polarization degree of the variation component was higher than that of the total flux. These results suggest that the microvariability originated from a small and local region where the magnetic field is aligned.

We investigate the total kinetic powers (L_j) and ages (t_age) of powerful jets in FR II radio galaxies by comparison of the dynamical model of expanding cocoons with observations. We select four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219), for which the mass-density profiles of intracluster medium (ICM) are known in the literature. It is found that large fractions > 0.02 - 0.7 of the Eddington luminosity (L_Edd) are carried away as a kinetic power of jet. The upper limit of estimated 2 L_j / L_Edd are larger than unity (< 10) for some sources, suggesting a possibility of super-Eddington mass accretions. As a consequence of the large powers, we also find that the total energy stored in the cocoon (E_c) exceeds the energy derived from the minimum energy condition for the energy of radiating non-thermal electrons and magnetic fields (E_min): 4< E_c /E_min <310. This implies that most of the energy in cocoon is carried by invisible components such as thermal leptons (electron and positron) and/or protons.

We examine the evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects (MSOs), Fanaroff-Riley type II radio galaxies (FRIIs)], by comparing the relation between the hot spot size and the projected linear size with a coevolution model of hot spots and a cocoon. We take account of the deceleration effect by the cocoon head growth. We find that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRII phase). This is ascribed to the change of the power-law index of ambient density profile in the MSO phase ($\sim $1 kpc). It is also found that the cocoon shape becomes nearly spherical or disrupted for MSOs, while an elongated morphology is predicted for CSOs and FRIIs. This seems to be consistent with the higher fraction of distorted morphology of MSOs than that of CSOs and FRI. Finally, we predict that only CSOs whose initial advance speed is higher than about 0.1c can evolve into FRIIs, comparing the hot spot speed with the sound speed of the ambient medium.

We report time-resolved CCD photometry of the cataclysmic variable EG Aquarii during the 2006 November outburst During the outburst, superhumps were unambiguously detected with a mean period of 0.078828(6) days, firstly classifying the object as an SU UMa-type dwarf nova. It also turned out that the outburst contained a precursor. At the end of the precursor, immature profiles of humps were observed. By a phase analysis of these humps, we interpreted the features as superhumps. This is the second example that the superhumps were shown during a precursor. Near the maximum stage of the outburst, we discovered an abrupt shift of the superhump period by ${\sim}$ 0.002 days. After the supermaximum, the superhump period decreased at the rate of $\dot{P}/P$=$-8.2{\times}10^{-5}$, which is typical for SU UMa-type dwarf novae. Although the outburst light curve was characteristic of SU UMa-type dwarf novae, long-term monitoring of the variable shows no outbursts over the past decade. We note on the basic properties of long period and inactive SU UMa-type dwarf novae.

In this paper, we initiate a new study of steady funnel-flow accretion onto strongly magnetized neutron stars, including a full treatment of shock generation. As a first step, we adopt a simplified model considering the flow within Newtonian theory and neglecting radiative pressure and cooling. The flow is taken to start from an accretion disc and then to follow magnetic field lines, forming a transonic funnel flow onto the magnetic poles. A standing shock occurs at a certain point in the flow and beyond this material accretes subsonically onto the star with high pressure and density. We calculate the location of the standing shock and all other features of the flow within the assumptions of our model. Applications to observed X-ray pulsars are discussed.

Invisible plasma content in blazar jets such as protons and/or thermal electron-positron ($e^{\pm}$) pairs is explored through combined arguments of dynamical and radiative processes. By comparing physical quantities required by the internal shock model with those obtained through the observed radio-to-gamma-ray spectra for Mrk 421, we find the existence of a copious amount of invisible plasma in the jet. We speculate that the blazar sequence could arise from variations of total amount and/or blending ratio of $e^{\pm}$ pair and electron-proton plasma.

We report photometric observations of AL Com during its rare outburst in 2007. The light curve is reminiscent of its past superoutbursts in 1995 and 2001, except for the rebrightening phase after the main superoutburst. During the rebrightening phase in 2007, we found clear modulations between V=16.2-15.2. In conjunction with the lack of prominent superhumps in our time-series observations, the modulations can most naturally be interpreted as repetitive short rebrightenings with a cycle of 1-2 days. The rebrightening characteristics in 2007 are different from those in 1995 and 2001. This indicates that the type of rebrightenings in WZ Sge stars depends not on binary parameters of objects, but on the mass-accretion process for each outburst.

Several SU UMa-type dwarf novae, in particular, WZ Sge-type stars tend to exhibit rebrightenings after superoutbursts. The rebrightening phenomenon is problematic for the disk instability theory of dwarf novae since it requires a large amount of remnant matter in the disk even after superoutbursts. Here, we report our optical and infrared observations during the first-ever outburst of a new dwarf nova, SDSS J102146.44+234926.3. During the outburst, we detected superhumps with a period of 0.056281 +/- 0.000015 d, which is typical for superhump periods in WZ Sge stars. In conjunction with the appearance of a long-lived rebrightening, we conclude that the object is a new member of WZ Sge stars. Our observations, furthermore, revealed infrared behaviors for the first time in the rebrightening phase of WZ Sge stars. We discovered prominent infrared superhumps. We calculate the color temperature of the infrared superhump source to be 4600-6400 K. These temperatures are too low to be explained with a fully-ionized disk appearing during dwarf nova outbursts. We also found a Ks-band excess over the hot disk component. These unprecedented infrared activities provide evidence for the presence of mass reservoir at the outermost part of the accretion disk. We propose that a moderately high mass-accretion rate at this infrared active region leads to the long-lived rebrightening observed in SDSS J102146.44+234926.3.

Invisible plasma content in blazar jets such as protons and/or thermal electron-positron ($e^{\pm}$) pairs is explored through combined arguments of dynamical and radiative processes. By comparing physical quantities required by the internal shock model with those obtained through the observed broadband spectra for Mrk 421, we obtain that the ratio of the Lorentz factors of a pair of cold shells resides in about $2\sim 20$, which implies that the shocks are at most mildly relativistic. Using the obtained Lorentz factors, the total mass density $\rho$ in the shocked shells is investigated. The upper limit of $\rho$ is obtained from the condition that thermal bremsstrahlung emission should not exceed the observed $\gamma$-ray luminosity, whilst the lower limit is constrained from the condition that the energy density of non-thermal electrons is smaller than that of the total plasma. Then we find $\rho$ is $10^2$-$10^3$ times heavier than that of non-thermal electrons for pure $e^{\pm}$ pairs, while $10^2$-$10^6$ times heavier for pure electron-proton ($e/p$) content, implying the existence of a large amount of invisible plasma. The origin of the continuous blazar sequence is shortly discussed and we speculate that the total mass density and/or the blending ratio of $e^{\pm}$ pairs and $e/p$ plasma could be new key quantities for the origin of the sequence.

We investigate the total kinetic powers (L_{j}) and ages (t_{age}) of powerful jets of four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219) by the detail comparison of the dynamical model of expanding cocoons with observed ones. It is found that these sources have quite large kinetic powers with the ratio of L_{j} to the Eddington luminosity (L_{Edd}) resides in $0.02 <L_{j}/L_{Edd} <10$. Reflecting the large kinetic powers, we also find that the total energy stored in the cocoon (E_{c}) exceed the energy derived from the minimum energy condition (E_{min}): $2< E_{c}/E_{min} <160$. This implies that a large amount of kinetic power is carried by invisible components such as thermal leptons (electron and positron) and/or protons.

Recently, strong AGN outbursts at the centers of galaxy clusters have been found. Using a simple model, we study particle acceleration around a shock excited by an outburst and estimate nonthermal emission from the accelerated particles. We show that emission from secondary electrons is consistent with the radio observations of the minihalo in the Perseus cluster, if there was a strong AGN outburst >~10^8 yrs ago with an energy of ~1.8x10^62 erg. The validity of our model depends on the frequency of the large outbursts. We also estimate gamma-ray emission from the accelerated particles and show that it could be detected with GLAST.

Strong $\gamma$-ray emission from cocoons of young radio galaxies is newly predicted. Considering the process of adiabatic injection of the shock dissipation energy and mass of the relativistic jet in active nuclei (AGNs) into the cocoon, while assuming thermalizing electron plasma interactions, we find that the thermal electron temperature of the cocoon is typically predicted in $\sim$MeV, which is determined only by the bulk Lorentz factor of the relativistic jet. Together with the time-dependent dynamics of the cocoon expansion, we find that young cocoons can yield thermal bremsstrahlung emissions at energies $\sim$ MeV.

We describe the dynamical evolution of hot spots velocity, pressure and mass density in radio loud active galactic nuclei (AGNs), taking proper account of (1) the conservations of the mass, momentum, and kinetic energy flux of the unshocked jet, (2) the deceleration process of the jet by shocks, and (3) the cocoon expansion without assuming the constant aspect ratio of the cocoon. By the detailed comparison with two dimensional relativistic hydrodynamic simulations, we show that our model well reproduces the whole evolution of relativistic jets. Our model can explain also the observational trends of the velocity, the pressure, the size, and mass density of hot spots in compact symmetric objects (CSOs) and FR II radio galaxies.

We examine the constraints imposed on the total kinetic power ($L_{\rm j}$) and the age ($t_{\rm age}$) of relativistic jets in FR II radio sources in a new way. We solve the dynamical expansion of its cocoon embedded in the intra-cluster medium (ICM) and obtain the analytic solution of its physical quantities in terms of $L_{\rm j}$ and $t_{\rm age}$. The estimate of $L_{\rm j}$ and $t_{\rm age}$ is done by the comparison of the model and the observed shape of the cocoon. The analysis is focused on Cygnus A and we find that (i) the source age is $3 {\rm Myr} <t_{\rm age}<30 {\rm Myr}$ and (ii) the total kinetic power of the jet is estimated as $0.2\times 10^{46}{\rm erg s^{-1}} <L_{\rm j}<1\times 10^{48} {\rm erg s^{-1}}$ which is larger than 1% of the Eddington luminosity of Cygnus A.

The matter content of jets in active galactic nuclei is examined in a new way. We model the dynamical expansion of its cocoon embedded in the intra-cluster medium (ICM). By comparing the observed shape of the cocoon with that expected from the theoretical model, we estimate the total pressure ($P_{\rm c}$) and electron temperature ($T_{e}$) of the cocoon. The number density of the total electrons ($n_{e^{-}}$) is constrained by using the non-thermal spectrum of the hot spot and the analysis of the momentum balance between the jet thrust and the rum pressure of ICM. Together with the obtained $P_{\rm c}$, $T_{e}$ and $n_{e^{-}}$, we constrain the matter content in the jets. We find that, in the case of Cygnus A, the ratio of number density of protons to that of electrons is of order of $10^{-3}$ . This implies the existence of a large number of positron in the jet.

We study the hydrodynamical effects of two colliding shells, adopted to model internal shocks in various relativistic outflows such as gamma-ray bursts and blazars. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of the two shells. The shell spreading with a few ten percent of the speed of light is also shown to be a notable aspect. The conversion efficiency of the bulk kinetic energy to internal one shows deviations from the widely-used inelastic two-point-mass-collision model. Observational implications are also shortly discussed.

We study both analytically and numerically hydrodynamical effects of two colliding shells, the simplified models of the internal shock in various relativistic outflows such as gamma-ray bursts and blazars. We pay particular attention to three interesting cases: a pair of shells with the same rest mass density (``{\it equal rest mass density}''), a pair of shells with the same rest mass (``{\it equal mass}''), and a pair of shells with the same bulk kinetic energy (``{\it equal energy}'') measured in the intersteller medium (ISM) frame. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of two shells for the ``equal mass'' case, while no significant split appears for the ``equal energy'' and ``equal rest mass density'' cases. The shell spreading with a few ten percent of the speed of light is also shown as a notable aspect caused by rarefaction waves. The conversion efficiency of the bulk kinetic energy to internal one is numerically evaluated. The time evolutions of the efficiency show deviations from the widely-used inellastic two-point-mass-collision model.

We explore the energetics and plasma composition in FR II sources using a new simple method of combining shock dynamics and radiation spectrum. The hot spots are identified with the reverse shocked region of jets. With the one-dimensional shock jump conditions taking account of the finite pressure of hot ICM, we estimate the rest mass and energy densities of the sum of thermal and non-thermal particles in hot spots. Independently, based on the Synchrotron Self-Compton (SSC) model, we estimate the number and energy densities of {\it non-thermal} electrons using the multi-frequency radiation spectrum of hot spots. We impose the condition that the obtained rest mass, internal energy, and number densities of non-thermal electrons should be lower than those of the total particles determined by shock dynamics. We apply this method to Cygnus A. We examine three extreme cases of pure electron-positron pair plasma (Case I), pure electron-proton plasma with separate thermalization (Case II), and pure electron-proton plasma in thermal-equilibrium (Case III). By detailed SSC analysis for Cygnus A and 3C123, we find that the energy density of non-thermal electrons is about 10 times larger than that of magnetic field. We find that the Case III is not acceptable because predicted photon spectra do not give a good fit to the observed one. We find that Case II can also be ruled out since the number density of non-thermal electrons exceeds that of the total number density. Hence, we find that only pure $e^{\pm}$ plasma (Case I) is acceptable among the three cases. Total kinetic power of jet and electron acceleration efficiency are also constrained by internal energy densities of non-thermal and total particles.

We report the Chandra detection of an X-ray jet in 3C 15. The peak of the X-ray emission in the jet is 4.1'' (a projected distance of 5.1 kpc) from the nucleus, and coincident with a component previously identified in the radio and optical jets. We examine four models for the X-ray jet emission: (I) weak synchrotron cooling in equip., (II) moderate synchrotron cooling in equip., (III) weak synchrotron plus SSC cooling, and (IV) moderate synchrotron plus SSC cooling. We argue that case (II) can most reasonably explain the overall emission from knot C. Case (III) is also possible, but requires a large departure from equipartition and for the jet power to be comparable to that of the brightest quasars. Diffuse X-ray emission has also been detected, distributed widely over the full extent (63kpc x 25kpc) of the radio lobes. We compare the total energy contained in the lobes with the jet power estimated from knot C, and discuss the energetic link between the jet and the lobes. We argue that the fueling time (t_fuel) and the source age (t_src) are comparable for case (II), whereas t_fuel << t_src is likely for case (III). The latter may imply that the jet has a very small filling factor, ~10^{-3}. We consider the pressure balance between the thermal galaxy halo and non-thermal relativistic electrons in the radio lobes. Finally, we show that the X-ray emission from the nucleus is not adequately fitted by a simple absorbed power-law model, but needs an additional power-law with heavy absorption intrinsic to the source. Such a high column density is consistent with the presence of a dense, dusty torus which obscures the quasar nucleus.

The effects of a wind on the emerging spectrum from an inefficiently-radiating accretion flow in a global magnetic field are examined, based on the analytic solution obtained recently by one of the present authors. The results exhibit the steepening of the negative slope appearing in the intermediate frequency range of bremsstrahlung spectrum and the decrease in the luminosity ratio of thermal synchrotron to bremsstrahlung, in accordance with the increasing wind strength. Both effects are due to a suppressed mass accretion rate in the inner disk, caused by a mass loss in terms of wind. In order to demonstrate the reliability of this model, Sagittarius A^* (Sgr A^*) and the nucleus of M 31, both of which have been resolved in an X-ray band by Chandra, are taken up as the best candidates for the broadband spectral fittings. Although the observed X-ray data are reproduced for these objects by both of the inverse-Compton and the bremsstrahlung fittings, some evidence of preference for the latter are recognized. The wind effects are clearly seen in the latter fitting case, in which we can conclude that a widely extending accretion disk is present in each nucleus, with no or only weak wind in Sgr A^* and with a considerably strong wind in the nuclear region of M 31. Especially in Sgr A^*, the inferred mass accretion rates are much smaller than the Bondi rate whose estimate has become reliable due to Chandra. This fact strongly suggests that the accretion in this object does not proceed like Bondi's prediction, though its extent almost reaches the Bondi radius.

Using multi-frequency spectra from TeV blazars in quiescent states, we obtain the physical parameters of the emission region of blazars within the framework of the one-zone synchrotron self-Compton (SSC) model. We numerically calculate the steady-state energy spectra of electrons by self-consistently taking into account the effects of radiative cooling with a proper account of the Klein-Nishina effects. Here electrons are assumed to be injected with a power-law spectrum and to escape on a finite time scale, which naturally leads to the existence of a break energy scale. Although we do not use time variabilities but utilize a model of electron escape to constrain the size of the emission region, the resultant size turns out to be similar to that obtained based on time variabilities. Through detailed comparison of the predicted emission spectra with observations, we find that for Mrk 421, Mrk 501, and PKS 2155--304, the energy density of relativistic electrons is about an order of magnitude larger than that of magnetic fields with an uncertainty within a factor of a few.

We calculate the radiation spectra from advection-dominated accretion flows (ADAFs), taking into account the effects of a global magnetic field. Calculation is based on the analytic model for magnetized ADAFs proposed by Kaburaki, where a large-scale magnetic field controls the accretion process. Adjusting a few parameters, we find that our model can well reproduce the observed spectrum of Sagittarius A$^{*}$. The result is discussed in comparison with those of well-known ADAF models, where the turbulent viscosity controls the accretion process.