Many studies have shown that there are clear sequences in the period-luminosity relationship (PLR) for Mira variables and semiregular variables (SRVs) in the Large Magellanic Cloud (LMC). To investigate the PLR for SRVs in our galaxy, we examined the annual parallax measurement and conducted K'-band photometric monitoring of an SRV star SV Pegasus (SV Peg). We measured the position change of the associating H$_2$O maser spots by phase-referencing VLBI observations with VERA at 22 GHz, spanning approximately 3 yr, and detected an annual parallax of $\pi = 3.00 \pm 0.06$ mas, corresponding to a distance of $D=333 \pm 7$ pc. This result is in good agreement with the Hipparcos parallax and improves the accuracy of the distance from 35 $\%$ to 2 $\%$. However, the GAIA DR2 catalog gave a parallax of $\pi=1.12\pm0.28$ mas for SV Peg. This indicates that the GAIA result might be blurred by the effect of the stellar size because the estimated stellar radius was $\sim 5$ mas, which is comparable to the parallax. We obtained a K'-band mean magnitude of $m_{K'} = -0.48$ mag and a period of $P=177$ days from our photometric monitoring with a 1-m telescope. Using the trigonometric distance, we derived an absolute magnitude of $M_{K'}=-8.09 \pm 0.05$ mag. This result shows that the position of SV Peg in the PLR falls on the C' sequence found in the PLR in the LMC, which is similar to other SRVs in our galaxy.

H$_2$O maser emission at 22 GHz in the circumstellar envelope is one of the good tracers of detailed physics and kinematics in the mass loss process of asymptotic giant branch stars. Long-term monitoring of an H$_2$O maser spectrum with high time resolution enables us to clarify acceleration processes of the expanding shell in the stellar atmosphere. We monitored the H$_2$O maser emission of the semi-regular variable R Crt with the Kagoshima 6-m telescope, and obtained a large data set of over 180 maser spectra during 1.3 years with an observational span of a few days. Based on an automatic peak detection method based on a least-squares fitting with a Gaussian basis function model, we exhaustively detected peaks as significant velocity components with the radial velocity on a 0.1 km/s scale. This analysis result shows that the radial velocity of red-shifted and blue-shifted components exhibits a change between acceleration and deceleration on the time scale of a few hundred days. These velocity variations are likely to correlate with the intensity variations, in particular during flaring state of H$_2$O masers. It seems reasonable to consider that the velocity variation of the maser source is caused by the shock propagation in the envelope due to stellar pulsation. However, it is difficult to explain the relationship between the velocity variation and the intensity variation only from shock propagation effects. We found that a time delay of the integrated maser intensity with respect to the optical light curve is about 150 days.

We present measurements of proper motion of the sub-parsec scale jet at 22 GHz in the nearby FR I galaxy 3C 66B. Observations were made using the VLBA at six epochs over four years. A phase-referencing technique was used to improve the image quality of the weak and diffuse jet components. We find that the inner knots are almost stationary, though one of them was expected to be detected with the apparent speed of 0.2 mas/yr according to 8 GHz monitoring at the same observation epochs. Clear flux variations are not observed in the core at 22 GHz, in contrast, clear flux enhancement is observed in the core at 8 GHz. We discussed that this can be explained, if the jet has helical structure, that the viewing angle of the jet between 8 and 22 GHz differs by a few degree in case the jet direction is almost along our line of sight. Although these results may imply the existence of a two-zone jet, which has been suggested in certain radio galaxies, it cannot explain the fact that the jet at the higher frequency jet is slower than that at the lower frequency.

We report multifrequency phase-referenced observations of the nearby radio galaxy NGC 4261, which has prominent two-sided jets, using the Very Long Baseline Array at 1.4-43 GHz. We measured radio core positions showing observing frequency dependences (known as "core shift") in both approaching jets and counter jets. The limit of the core position as the frequency approaches infinity, which suggests a jet base, is separated by 82$\pm$16 ${\mu}$as upstream in projection, corresponding to (310$\pm$60)Rs (Rs: Schwarzschild radius) as a deprojected distance, from the 43 GHz core in the approaching jet. In addition, the innermost component at the counter jet side appeared to approach the same position at infinity of the frequency, indicating that cores on both sides are approaching the same position, suggesting a spatial coincidence with the central engine. Applying a phase referencing technique, we also obtained spectral index maps, which indicate that emission from the counter jet is affected by free-free absorption (FFA). The result of the core shift profile on the counter jet also requires FFA because the core positions at 5-15GHz cannot be explained by a simple core shift model based on synchrotron self-absorption (SSA). Our result is apparently consistent with the SSA core shift with an additional disk-like absorber over the counterjet side. Core shift and opacity profiles at the counter jet side suggest a two-component accretion: a radiatively inefficient accretion flow at the inner region and a truncated thin disk in the outer region. We proposed a possible solution about density and temperature profiles in the outer disk on the basis of the radio observation.

We present the observational results of the Gamma-ray blazar, 3C 66A, at 2.3, 8.4, and 22 GHz at 4 epochs during 2004-05 with the VLBA. The resulting images show an overall core-jet structure extending roughly to the south with two intermediate breaks occurring in the region near the core. By model-fitting to the visibility data, the northmost component, which is also the brightest, is identified as the core according to its relatively flat spectrum and its compactness. As combined with some previous results to investigate the proper motions of the jet components, it is found the kinematics of 3C 66A is quite complicated with components of inward and outward, subluminal and superluminal motions all detected in the radio structure. The superluminal motions indicate strong Doppler boosting exists in the jet. The apparent inward motions of the innermost components last for at least 10 years and could not be caused by new-born components. The possible reason could be non-stationarity of the core due to opacity change.

We present proper motion of the sub-parsec scale jet in a nearby elliptical galaxy 3C 66B. Observations were made using the VLBA and partly Effelsburg 100-m telescope at 2.3 GHz and 8.4 GHz at 10 epochs over 4 years. The 8.4 GHz images showed that a proper motion increases from 0.21 to 0.70 mas/year, corresponding to an apparent speed of 0.30 c to 0.96 c, with a distance from the core on a sub-parsec scale. Our investigation suggests that the apparent increase of the proper motion can be explained by changes in the viewing angle, according to a relativistic beaming model. However, we still cannot eliminate the possibility that acceleration of the jet outflow speed or of changes of emissivity profile in the two-zone jet might be found in 3C 66B.

Recent observational results provide possible evidence that binary black holes (BBHs) exist in the center of giant galaxies and may merge to form a supermassive black hole in the process of their evolution. We first detected a periodic flux variation on a cycle of $93\pm1$ days from the 3-mm monitor observations of a giant elliptical galaxy \object{3C 66B} for which an orbital motion with a period of $1.05\pm0.03$ years had been already observed. The detected signal period being shorter than the orbital period can be explained by taking into consideration the Doppler-shifted modulation due to the orbital motion of a BBH. Assuming that the BBH has a circular orbit and that the jet axis is parallel to the binary angular momentum, our observational results demonstrate the presence of a very close BBH that has the binary orbit with an orbital period of $1.05\pm0.03$ years, an orbital radius of $(3.9\pm1.0) \times 10^{-3}$ pc, an orbital separation of $(6.1^{+1.0}_{-0.9}) \times 10^{-3}$ pc, the larger black hole mass of $(1.2^{+0.5}_{-0.2}) \times 10^9$ $M_{\sun}$, and the smaller black hole mass of $(7.0^{+4.7}_{-6.4}) \times 10^8$ $M_{\sun}$. The BBH decay time of $(5.1^{+60.5}_{-2.5})\times 10^2$ years provides evidence for the occurrence of black hole mergers. This Letter will demonstrate the interesting possibility of black hole collisions to form a supermassive black hole in the process of evolution, one of the most spectacular natural phenomena in the universe.

The next-generation space VLBI mission, VSOP-2, is expected to provide unprecedented spatial resolutions at 8.4, 22, and 43GHz. In this report, phase referencing with VSOP-2 is examined in detail based on a simulation tool called ARIS. The criterion for successful phase referencing was to keep the phase errors below one radian. Simulations with ARIS reveal that phase referencing achieves good performance at 8.4GHz, even under poor tropospheric conditions. At 22 and 43GHz, it is recommended to conduct phase referencing observations under good or typical tropospheric conditions. The satellite is required to have an attitude-switching capability with a one-minute or shorter cycle, and an orbit determination accuracy higher than about 10cm at apogee; the phase referencing calibrators are required to have a signal-to-noise ratio larger than four for a single scan. The probability to find a suitable phase referencing calibrator was estimated by using VLBI surveys. From the viewpoint of calibrator availability, VSOP-2 phase referencing at 8.4GHz is promising. However, the change of finding suitable calibrators at 22 and 43GHz is significantly reduced; it is important to conduct specific investigations for each target at those frequencies.

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

We have made multi-epoch VLBI observations of H2O maser emission in the massive star forming region IRAS 06061+2151 with the Japanese VLBI network (JVN) from 2005 May to 2007 October. The detected maser features are distributed within an 1\arcsec$\times$1\arcsec (2000 au$\times$2000 au at the source position) around the ultra-compact H {\small\bf II} region seen in radio continuum emission. Their bipolar morphology and expanding motion traced through their relative proper motions indicate that they are excited by an energetic bipolar outflow. Our three-dimensional model fitting has shown that the maser kinematical structure in IRAS 06061+2151 is able to be explained by a biconical outflow with a large opening angle ($>$ 50\degr). The position angle of the flow major axis coincides very well with that of the large scale jet seen in 2.1$\:\mu\rmn{m}$ hydrogen emission. This maser geometry indicates the existence of dual structures composed of a collimated jet and a less collimated massive molecular flow. We have also detected a large velocity gradient in the southern maser group. This can be explained by a very small (on a scale of several tens of au) and clumpy (the density contrast by an order of magnitude or more) structure of the parental cloud. Such a structure may be formed by strong instability of shock front or splitting of high density core.

We present the results of the multi-frequency observations of radio outburst of the microquasar Cyg X-3 in February and March 2006 with the Nobeyama 45-m telescope, the Nobeyama Millimeter Array, and the Yamaguchi 32-m telescope. Since the prediction of a flare by RATAN-600, the source has been monitored from Jan 27 (UT) with these radio telescopes. At the eighteenth day after the quench of the activity, successive flares exceeding 1 Jy were observed successfully. The time scale of the variability in the active phase is presumably shorter in higher frequency bands. We also present the result of a follow-up VLBI observation at 8.4 GHz with the Japanese VLBI Network (JVN) 2.6 days after the first rise. The VLBI image exhibits a single core with a size of <8 mas (80 AU). The observed image was almost stable, although the core showed rapid variation in flux density. No jet structure was seen at a sensitivity of $T_b = 7.5\times 10^5$ K.

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

We present the results of six-epoch Very Long Baseline Array (VLBA) observations of 3C~66A. The high-resolution Very Long Baseline Interferometer (VLBI) maps obtained at multi-frequency (2.3, 8.4, and 22.2 GHz) simultaneously enabled us to identify the brightest compact component with the core. We find that the spectrum of the core can be reasonably fitted by the synchrotron self-absorption model. Our VLBA maps show that the jet of 3C~66A has two bendings at about 1.2 and 4 mas from the core. We also give possible identifications of our jet components with the components in previous VLBA observations by analysing their proper motions. We find consistent differences of the position from the core in one component between different frequencies at six epochs.

We present the basic features and the activities of Japanese VLBI network (JVN), a newly-established VLBI network with baselines ranging from 50 to 2560 km spreading across the Japanese islands, and capable of observing at 6.7, 8.4, and 22 GHz. We show a number of results of JVN observations: 8.4-GHz continuum images of a Giga-hertz Peaked Spectrum (GPS) source and radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s), the spatial and velocity structures of water masers in NML Cygni as well as methanol masers in Cep A, and demonstrative observations with the bigradient phase referencing.

We study the accretion flows from the circumbinary disks onto the supermassive binary black holes in a subparsec scale of the galactic center, using a smoothed particles hydrodynamics (SPH) code. Simulation models are presented in four cases of a circular binary with equal and unequal masses, and of an eccentric binary with equal and unequal masses. We find that the circumblack-hole disks are formed around each black holes regardless of simulation parameters. There are two-step mechanisms to cause an accretion flow from the circumbinary disk onto supermassive binary black holes: First, the tidally induced elongation of the circumbinary disk triggers mass inflow towards two closest points on the circumbinary disk from the black holes. Then, the gas is increasingly accumulated on these two points owing to the gravitational attraction of black holes. Second, when the gas can pass across the maximum loci of the effective binary potential, it starts to overflow via their two points and freely infalls to each black hole. In circular binaries, the gas continues to be supplied from the circumbinary disk (i.e. the gap between the circumbinary disk and the binary black hole is always closed.) In eccentric binaries, the mass supply undergoes the periodic on/off transitions during one orbital period because of the variation of periodic potential. The gap starts to close after the apastron and to open again after the next periastron passage. Due to this gap closing/opening cycles, the mass-capture rates are eventually strongly phase dependent. This could provide observable diagnosis for the presence of supermassive binary black holes in merged galactic nuclei.

We investigate the accretion flows onto the supermassive binary black holes (SMBBHs) from the circumbinary disk with the equal mass, eccentric binary on the subparsec scale, using Smoothed Particle Hydrodynamics (SPH) code. We find that the material can be supplied from circumbinary disk, which leads to the formation of two accretion disks around the SMBBHs. The mass accretion rates significantly modulate with the binary orbital motion. These could provide the observable diagnosis of the existence of the supermassive binary black holes (e.g. OJ287) on the subparsec scale in merged galactic nuclei.

We propose bigradient phase referencing (BPR), a new radio-observation technique, and report its performance using the Japanese very-long-baseline-interferometry network (JVN). In this method, a weak source is detected by phase-referencing using a primary calibrator, in order to play a role as a secondary calibrator for phase-referencing to a weak target. We will be given the opportunity to select a calibrator from lots of milli-Jansky sources, one of which may be located at the position closer to the target. With such a smaller separation, high-quality phase-referencing can be achieved. Furthermore, a subsequent more-sophisticated calibration can relocate array's focus to a hypothetical point much closer to the target; a higher quality of phase referencing is available. Our demonstrative observations with strong radio sources have proved the capabilities of BPR in terms of image dynamic ranges and astrometric reproducibility. The image dynamic range on a target has been improved with a factor of about six compared to that of normal phase-referencing; the resultant position difference of target's emission between two epochs was only 62+-50 micro-arcsecond, even with less than 2300-km baselines at 8.4 GHz and fast-switching of a target-calibrator pair of a 2.1-degree separation.

Supermassive black hole binaries may exist in the centers of active galactic nuclei like quasars and radio galaxies and mergers between galaxies may result in the formation of supermassive binaries during the course of galactic evolution. Using the very-long-baseline interferometer, we imaged the radio galaxy 3C 66B at radio frequencies and found that the unresolved radio core of 3C 66B shows well-defined elliptical motions with a period of 1.05 +/- 0.03 years, which provides a direct detection of a supermassive black hole binary.

Water masers around an AGB star, IRC+60169, were observed at four epochs using the Japanese VLBI networks. The distribution of the maser features is limited in a thick-shell region, which has inner and outer expansion velocities of 7 km/s and 14 km/s at radii of 25 mas and 120 mas, respectively. The distribution of the red-shifted features exhibits a ring-like structure, the diameter of which is 30 mas, and corresponds to the inner radius of the maser shell. This implies that dense gas around the star obscures red-shifted emission. Although a position--radial velocity diagram for the maser features is consistent with a spherical shell model, the relative proper motions do not indicate an expansion motion of the shell. A remarkable property has been found that is a possible periodic change of the alignment pattern of water maser spots.

We present new optical spectroscopy of the high-redshift powerful radio galaxy MRC 0406$-$244 at redshift of 2.429. We find that the two extensions toward NW and SE probed in the rest-frame ultraviolet image are heated mainly by the nonthermal continuum of the active galactic nucleus. However, each extension shows a shell-like morphology, suggesting that they are a pair of superbubbles induced by the superwind activity rather than by the interaction between the radio jet and the ambient gas clouds. If this is the case, the intense starburst responsible for the formation of superbubbles could occur $\sim 1 \times 10^9$ yr ago. On the other hand, the age of the radio jets may be of the order of $\sim 10^6$ yr, being much shorter than the starburst age. Therefore, the two events, i.e., the starburst and the radio-jet activities, are independent phenomena. However, their directions of the expanding motions could be governed by the rotational motion of the gaseous component in the host galaxy. This idea appears to explain the alignment effect of MRC 0406$-$244.

In order to investigate the genesis of powerful radio jet, we have mapped the central 10 pc region of the nearby radio galaxy NGC 6251 with a 0.2 pc resolution using Very Long Baseline Interferometer (VLBI) at two radio frequencies, 5 GHz and 15 GHz, we have found the sub-parsec-scale counterjet for the first time in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic beaming model.

We report on large-scale, regular morphological patterns found in the radio jet of the nearby radio galaxy NGC 6251. Investigating morphological properties of this radio jet from the nucleus to a radial distance of $\sim$ 300 arcsec ($\approx$ 140 kpc) mapped at 1662 MHz and 4885 MHz by Perley, Bridle, & Willis, we find three chains, each of which consists of five radio knots. We also find that eight radio knots in the first two chains consist of three small sub-knots (the triple-knotty substructures). We discuss the observational properties of these regular morphological patterns.

Mapping the central 5 pc region of the nearby radio galaxy NGC 6251 with a 0.2 pc resolution using VLBI at two radio frequencies, 5 GHz and 15 GHz, we have found the sub-pc-scale counter jet for the first time in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic beaming model (Ghisellini et al. 1993).