We discovered a new growth mode of dust grains to km-sized bodies in protoplanetary disks that evolve by viscous accretion and magnetically driven disk winds (MDWs). We solved an approximate coagulation equation of dust grains with time-evolving disks that consist of both gas and solid components by a one dimensional model. With the collisional growth of dust grains, all solid particles initially drift inward toward the central star by the gas drag force. However, the radial profile of gas pressure, $P$, is modified by the MDW that disperses the gas in an inside-out manner. As a result, a local concentration of solid particles is created by the converging radial flux of drifting dust grains at the location with the convex upward profile of $P$. When the dimensionless stopping time, or the Stokes number, St, there exceeds unity, the solid particles spontaneously reach the growth dominated state because of the positive feedback between the suppressed radial drift and the enhanced accumulation of dust particles that drift from the outer part. These consecutive processes proceed even for spherical dust grains under the Epstein gas drag. Once the solid particles are in the drift limited state, the above mentioned condition of St $\gtrsim 1$ for the dust growth is equivalent with \begin{equation} \Sigma_{\rm d}/\Sigma_{\rm g} \gtrsim \eta, \nonumber \end{equation} where $\Sigma_{\rm d}/\Sigma_{\rm g}$ is the dust-to-gas surface-density ratio and $\eta$ is dimensionless radial pressure-gradient force. As a consequence of the successful growth of dust grains, a ring-like structure containing planetesimal-sized bodies is formed at the inner part of the protoplanetary disks. Such a ring-shaped concentration of planetesimals is expected to play a vital role in the subsequent planet formation.

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disruptions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no flat remnants with the ratio of minor to major axis lengths c/a < 0.5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain significantly flat asteroids. We suggest that these flat asteroids may be interlopers or formed through low-velocity collisions between member asteroids.

Terrestrial planets are believed to be formed via giant impacts of Mars-sized protoplanets. Planets formed via giant impacts have highly eccentric orbits. A swarm of planetesimals around the planets may lead to eccentricity damping for the planets via the equipartition of random energies (dynamical friction). However, dynamical friction increases eccentricities of planetesimals, resulting in high velocity collisions between planetesimals. The collisional cascade grinds planetesimals to dust until dust grains are blown out due to radiation pressure. Therefore, the total mass of planetesimals decreases due to collisional fragmentation, which weakens dynamical friction. We investigate the orbital evolution of protoplanets in a planetesimal disk, taking into account collisional fragmentation of planetesimals. For 100 km-sized or smaller planetesimals, dynamical friction is insignificant for eccentricity damping of planets because of collisional fragmentation. On the other hand, giant impacts eject collisional fragments. Although the total mass of giant impact ejecta is 0.1-0.3 Earth masses, the largest impact ejecta are ~ 1,000 km in size. We also investigate the orbital evolution of single planets with initial eccentricities 0.1 in a swarm of such giant impact ejecta. Although the total mass of giant impact ejecta decreases by a factor of 3 in 30 Myrs, eccentricities of planets are damped down to the Earth level (~0.01) due to interaction with giant impact ejecta. Therefore, giant impact ejecta play an important role for determination of terrestrial planet orbits.

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 present the first sub-arcsecond images of 49 Ceti in the [C I]3P1-3P0 emission and the 614 micron dust continuum emission observed with ALMA, as well as that in the CO(J=3-2) emission prepared by using the ALMA archival data. The spatial distribution of the 614 micron dust continuum emission is found to have a broad-ring structure with a radius of about 100 au around the central star. A substantial amount of gas is also associated with 49 Ceti. The [C I] emission map shows two peaks inside the dust ring, and its overall extent is comparable to that of the dust continuum emission and the CO emission. We find that the [C I]/CO(J=3-2) intensity ratio significantly varies along the major axis. The ratio takes the minimum value of 1.8 around the dust peak position, and increases inwards and outwards. The enhanced ratio around the central star (~3) likely originates from the stellar UV radiation, while that in the outer disk (~10) from the interstellar UV radiation. Such complex distributions of the [C I] and CO(J=3-2) emission will be a key to understand the origin of the gas in 49 Ceti, and will also provide a stringent constraint on physical and chemical models of gaseous debris disks.

We report the first detection of C$_2$ $A^1\Pi_u$--$X^1\Sigma_g^+$ (0,0) and CN $A^2\Pi_u$--$X^2\Sigma^+$ (0,0) absorption bands in the interstellar medium. The detection was made using the near-infrared (0.91--1.35 $\mu$m) high-resolution ($R=20,000$ and 68,000) spectra of Cygnus OB2 No.\,12 collected with the WINERED spectrograph mounted on the 1.3 m Araki telescope. The $A$--$X$ (1,0) bands of C$_2$ and CN were detected simultaneously. These near-infrared bands have larger oscillator strengths, compared with the $A$--$X$ (2,0) bands of C$_2$ and CN in the optical. In the spectrum of the C$_2$ (0,0) band with $R=68,000$, three velocity components in the line of sight could be resolved and the lines were detected up to high rotational levels ($J''\sim20$). By analyzing the rotational distribution of C$_2$, we could estimate the kinetic temperature and gas density of the clouds with high accuracy. Furthermore, we marginally detected weak lines of $^{12}$C$^{13}$C for the first time in the interstellar medium. Assuming that the rotational distribution and the oscillator strengths of the relevant transitions of $^{12}$C$_2$ and $^{12}$C$^{13}$C are the same, the carbon isotope ratio was estimated to be $^{12}\text{C}/^{13}\text{C}=50$--100, which is consistent with the ratio in the local interstellar medium. We also calculated the oscillator strength ratio of the C$_2$ (0,0) and (1,0) bands from the observed band strengths. Unfortunately, our result could not discern theoretical and experimental results because of the uncertainties. High-resolution data to resolve the velocity components will be necessary for both bands in order to put stronger constraints on the oscillator strength ratios.

Context. Although heavily obscured active galactic nuclei (AGNs) have been found by many observational studies, the properties of the surrounding dust are poorly understood. Using AKARI/IRC spectroscopy, we discover a new sample of a heavily obscured AGN in LEDA 1712304 which shows a deep spectral absorption feature due to silicate dust. Aims. We study the infrared (IR) spectral properties of circumnuclear silicate dust in LEDA 1712304. Methods. We perform IR spectral fitting, considering silicate dust properties such as composition, porosity, size and crystallinity. Spectral energy distribution (SED) fitting is also performed to the flux densities in the UV to sub-millimeter range to investigate the global spectral properties. Results. The best-fit model indicates 0.1 $\mu$m-sized porous amorphous olivine (${\rm Mg_{2x}Fe_{2-2x}SiO_4}$; $x=0.4$) with $4\%$ crystalline pyroxene. The optical depth is $\tau_{\rm sil}{\sim}2.3$, while the total IR luminosity and stellar mass are estimated to be $L_{\rm IR}=(5\pm1){\times}10^{10}\,L_{\odot}$ and $M_{\rm star}=(2.7\pm0.8){\times}10^{9}\,M_{\odot}$, respectively. In such low $L_{\rm IR}$ and $M_{\rm star}$ ranges, there are few galaxies which show that large ${\tau}_{\rm sil}$. Conclusions. The silicate dust in the AGN torus of LEDA 1712304 has properties notably similar to those in other AGNs as a whole, but slightly different in the wing shape of the absorption profile. The porosity of the silicate dust suggests dust coagulation or processing in the circumnuclear environments, while the crystallinity suggests that the silicate dust is relatively fresh.

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.

Light curve observations of a recently discovered interstellar object 1I/`Oumuamua suggest that this object has an extremely elongated shape with the axis ratio 0.3 or smaller. Planetesimal collisions can produce irregular shapes including elongated shapes. In this paper, we suggest that the extremely elongated shape of 1I/`Oumuamua may be the result of such an impact. To find detailed impact conditions to form the extremely elongated objects, we conduct numerical simulations of planetesimal collisions using Smoothed Particle Hydrodynamics method for elastic dynamics with self-gravity and interparticle friction. Impacts into strengthless target planetesimals with radius 50 m are conducted with various ratios of impactor mass to target mass q, friction angles phi_d, impact velocities v_imp, and impact angles theta_imp. We find that impacts with q \geq 0.5, phi_d \geq 40 degrees, v_imp \leq 40 degrees, and theta_imp \leq 30 degrees produce remnants with the ratio of intermediate to major axis length less than 0.3. This impact condition suggests that the parent protoplanetary disk in the planetesimal collision stage was weakly turbulent (alpha < 10^{-4} for the inner disk) and composed of planetesimals smaller than ~ 7 km to ensure small impact velocity.

More than a half of asteroids in the main belt have irregular shapes with the ratios of the minor to major axis lengths less than 0.6. One of the mechanisms to create such shapes is collisions between asteroids. The relationship between shapes of collisional outcomes and impact conditions such as impact velocities may provide information on the collisional environments and its evolutionary stages when those asteroids are created. In this study, we perform numerical simulations of collisional destruction of asteroids with radii 50 km and subsequent gravitational reaccumulation using Smoothed Particle Hydrodynamics for elastic dynamics with self-gravity, a model of fracture of rock, and a model of friction of completely damaged rock. We systematically vary the impact velocity from 50 m/s to 400 m/s and the impact angle from 5 degrees to 45 degrees. We investigate shapes of the largest remnants resulting from collisional simulations. As a result, various shapes (bilobed, spherical, flat, elongated, and hemispherical shapes) are formed through equal-mass and low-velocity (50 - 400 m/s) impacts. We clarify a range of the impact angle and velocity to form each shape. Our results indicate that irregular shapes, especially flat shapes, of asteroids with diameters larger than 80 km are likely to be formed through similar-mass and low-velocity impacts, which are likely to occur in the primordial environment prior to the formation of Jupiter.

In the core accretion scenario, gas giant planets are formed form solid cores with several Earth masses via gas accretion. We investigate the formation of such cores via collisional growth from kilometer-sized planetesimals in turbulent disks. The stirring by forming cores induces collisional fragmentation and surrounding planetesimals are ground down until radial drift. The core growth is therefore stalled by the depletion of surrounding planetesiamls due to collisional fragmentation and radial drift. The collisional strength of planetesimals determines the planetesimal-depletion timescale, which is prolonged for large planetesiamls. The size of planetesiamls around growing cores is determined by the planetesimal size distribution at the onset of runaway growth. Strong turbulence delays the onset of runaway growth, resulting in large planetesimals. Therefore, the core mass evolution depends on turbulent parameter $\alpha$; the formation of cores massive enough without significant depletion of surrounding planetesimals needs strong turbulence of $\alpha \gtrsim 10^{-3}$. However, the strong turbulence with $\alpha \gtrsim 10^{-3}$ leads to a significant delay of the onset of runaway growth and prevents the formation of massive cores within the disk lifetime. The formation of cores massive enough within several millions years therefore requires the several times enhancement of the solid surface densities, which is achieved in the inner disk $\lesssim 10$AU due to pile-up of drifting dust aggregates. In addition, the collisional strength $Q_{\rm D}^*$ even for kilometer-sized or smaller bodies affects the growth of cores; $Q_{\rm D}^* \gtrsim 10^7 \,{\rm erg/g}$ for bodies $\lesssim 1$ km is likely for this gas giant formation.

In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which is a grid-based Eulerian hydrocode. We examine the dependence of the critical specific impact energy Q_RD* on impact conditions for a wide range of specific impact energy (Q_R) from disruptive collisions to erosive collisions, and compare our results with previous studies. We find collision outcomes of the iSALE simulation agree well with those of the SPH simulation. Detailed analysis mainly gives three results. (1) The value of Q_RD* depends on numerical resolution, and is close to convergence with increasing numerical resolution. The difference in converged value of Q_RD* between the iSALE code and the SPH code is within 30%. (2) Ejected mass normalized by total mass (M_ej/M_tot) generally depends on various impact conditions. However, when Q_R is normalized by Q_RD* that is calculated for each impact simulation, M_ej/M_tot can be scaled by Q_R/Q_RD*, and is independent of numerical resolution, impact velocity and target size. (3) This similarity law for Q_R/Q_RD* is confirmed for a wide range of specific impact energy. We also derive a semi-analytic formula for Q_RD* based on the similarity law and the crater scaling law. We find that the semi-analytic formula for the case with a non-porous object is consistent with numerical results.

Recent radio observations towards nearby galaxies started to map the whole disk and to identify giant molecular clouds (GMCs) even in the regions between galactic spiral structures. Observed variations of GMC mass functions in different galactic environment indicates that massive GMCs preferentially reside along galactic spiral structures whereas inter-arm regions have many small GMCs. Based on the phase transition dynamics from magnetized warm neutral medium to molecular clouds, Kobayashi et al. 2017 proposes a semi-analytical evolutionary description for GMC mass functions including cloud-cloud collision (CCC) process. Their results show that CCC is less dominant in shaping the mass function of GMCs compared with the accretion of dense HI gas driven by the propagation of supersonic shock waves. However, their formulation does not take into account the possible enhancement of star formation by CCC. Radio observations within the Milky Way indicate the importance of CCC for the formation of star clusters and massive stars. In this article, we reformulate the time evolution equation largely modified from Kobayashi et al. 2017 so that we additionally compute star formation subsequently taking place in CCC clouds. Our results suggest that, although CCC events between smaller clouds outnumber the ones between massive GMCs, CCC-driven star formation is mostly driven by massive GMCs > 10^5.5 Msun (where Msun is the solar mass). The resultant cumulative CCC-driven star formation may amount to a few 10 percent of the total star formation in the Milky Way and nearby galaxies.

We perform global three-dimensional (3D) radiation-hydrodynamic (RHD) simulations of out- flow from supercritical accretion flow around a 10 Msun black hole. We only solve the outflow part, starting from the axisymmetric 2D simulation data in a nearly steady state but with small perturbations in a sinusoidal form being added in the azimuthal direction. The mass accretion rate onto the black hole is ~10^2 L_E/c^2 in the underlying 2D simulation data and the outflow rate is ~10 L_E/c^2 (with LE and c being the Eddington luminosity and speed of light, respectively). We first confirm the emergence of clumpy outflow, which was discovered by the 2D RHD simulations, above the photosphere located at a few hundreds of Schwarzschild radii (r_S) from the central black hole. As prominent 3D features we find that the clumps have the shape of a torn sheet, rather than a cut string, and that they are rotating around the central black hole with a sub-Keplerian velocity at a distance of ~10^3 r_S from the center. The typical clump size is ~30 r_S or less in the radial direction, and is more elongated in the angular directions, ~hundreds of r_S at most. The sheet separation ranges from 50 to 150 r_S. We expect stochastic time variations when clumps pass across the line of the sight of a distant observer. Variation timescales are estimated to be several seconds for a black hole with mass of ten to several tens of Msun, in rough agreement with the observations of some ultra-luminous X-ray sources.

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 investigate the formation of hot and massive circumplanetary disks (CPDs) and the orbital evolution of satellites formed in these disks. Because of the comparatively small size-scale of the sub-disk, quick magnetic diffusion prevents the magnetorotational instability (MRI) from being well-developed at ionization levels that would allow MRI in the parent protoplanetary disk. In the absence of significant angular momentum transport, continuous mass supply from the parental protoplanetary disk leads to the formation of a massive CPD. We have developed an evolutionary model for this scenario and have estimated the orbital evolution of satellites within the disk. We find, in a certain temperature range, that inward migration of a satellite can be stopped by a change in the structure due to the opacity transitions. Moreover, by capturing second and third migrating satellites in mean motion resonances, a compact system in Laplace resonance can be formed in our disk models.

Context. We are creating the AKARI mid-infrared all-sky diffuse maps. Through a foreground removal of the zodiacal emission, we serendipitously detected a bright residual component whose angular size is about 50 x 20 deg. at a wavelength of 9 micron. Aims. We investigate the origin and the physical properties of the residual component. Methods. We measured the surface brightness of the residual component in the AKARI mid-infrared all-sky maps. Results. The residual component was significantly detected only in 2007 January, even though the same region was observed in 2006 July and 2007 July, which shows that it is not due to the Galactic emission. We suggest that this may be a small cloud passing near the Earth. By comparing the observed intensity ratio of I_9um/I_18um with the expected intensity ratio assuming thermal equilibrium of dust grains at 1 AU for various dust compositions and sizes, we find that dust grains in the moving cloud are likely to be much smaller than typical grains that produce the bulk of the zodiacal light. Conclusions. Considering the observed date and position, it is likely that it originates in the solar coronal mass ejection (CME) which took place on 2007 January 25.

We have detected [C I] 3P1-3P0 emissions in the gaseous debris disks of 49 Ceti and Beta Pictoris with the 10 m telescope of the Atacama Submillimeter Telescope Experiment, which is the first detection of such emissions. The line profiles of [C I] are found to resemble those of CO(J=3-2) observed with the same telescope and the Atacama Large Millimeter/submillimeter Array. This result suggests that atomic carbon (C) coexists with CO in the debris disks, and is likely formed by the photodissociation of CO. Assuming an optically thin [C I] emission with the excitation temperature ranging from 30 to 100 K, the column density of C is evaluated to be (2.2+-0.2)x10^17 and (2.5+-0.7)x10^16 cm^-2 for 49 Ceti and Beta Pictoris, respectively. The C/CO column density ratio is thus derived to be 54+-19 and 69+-42 for 49 Ceti and Beta Pictoris, respectively. These ratios are higher than those of molecular clouds and diffuse clouds by an order of magnitude. The unusually high ratios of C to CO are likely attributed to a lack of H2 molecules needed to reproduce CO molecules efficiently from C. This result implies a small number of H2 molecules in the gas disk; i.e., there is an appreciable contribution of secondary gas from dust grains.

Disruptive collisions have been regarded as an important process for planet formation, while non-disruptive, small-scale collisions (hereafter called erosive collisions) have been underestimated or neglected by many studies. However, recent studies have suggested that erosive collisions are also important to the growth of planets, because they are much more frequent than disruptive collisions. Although the thresholds of the specific impact energy for disruptive collisions (Q_RD^*) have been investigated well, there is no reliable model for erosive collisions. In this study, we systematically carried out impact simulations of gravity-dominated planetesimals for a wide range of specific impact energy (Q_R) from disruptive collisions (Q_R ~ Q_RD^*) to erosive ones (Q_R << Q_RD^*) using the smoothed particle hydrodynamics method. We found that the ejected mass normalized by the total mass (M_ej/M_tot) depends on the numerical resolution, the target radius (R_tar) and the impact velocity (v_imp), as well as on Q_R, but that it can be nicely scaled by Q_RD^* for the parameter ranges investigated (R_tar = 30-300 km, v_imp = 2-5 km/s). This means that M_ej/M_tot depends only on Q_R/Q_RD^* in these parameter ranges. We confirmed that the collision outcomes for much less erosive collisions (Q_R < 0.01 Q_RD^*) converge to the results of an impact onto a planar target for various impact angles and that M_ej/M_tot = C * QR/QRD* holds. For disruptive collisions (Q_R ~ Q_RD^*), the curvature of the target has a significant effect on Mej/Mtot. We also examined the angle-averaged value of M_ej/M_tot and found that the numerically obtained relation between angle-averaged M_ej/M_tot and Q_R/Q_RD^* is very similar to the cases for 45-degree impacts. We proposed a new erosion model based on our numerical simulations for future research on planet formation with collisional erosion.

In the Solar System, interplanetary dust particles (IDPs) originating mainly from asteroid collisions and cometary activities drift to the Earth orbit due to the Poynting-Robertson drag. We analyzed the thermal emission from IDPs that was observed by the first Japanese infrared astronomical satellite, AKARI. The observed surface brightness in the trailing direction of the Earth orbit is 3.7% greater than that in the leading direction in the $9{\rm \mu m}$ band and 3.0% in the $18{\rm \mu m}$ band. In order to reveal dust properties causing the leading-trailing surface brightness asymmetry, we numerically integrated orbits of the Sun, the Earth, and a dust particle as a restricted three-body problem including radiation from the Sun. The initial orbits of particles are determined according to the orbits of main-belt asteroids or Jupiter-family comets. The orbital trapping in mean motion resonances results in a significant leading-trailing asymmetry so that intermediate sized dust (~10-100${\rm \mu m}$) produces a greater asymmetry than the zodiacal light has. The leading-trailing surface brightness difference integrated over the size distribution of the asteroidal dust is obtained to be the values of 27.7% and 25.3% in the $9{\rm \mu m}$ and $18{\rm \mu m}$ bands, respectively. In contrast, the brightness difference for cometary dust is calculated as the values of 3.6% and 3.1% in the $9{\rm \mu m}$ and $18{\rm \mu m}$ bands, respectively, if the maximum dust radius is set to be $s_{\rm max} = 3000{\rm \mu m}$. Taking into account these values and their errors, we conclude that the contribution of asteroidal dust to the zodiacal infrared emission is less than ~10%, while cometary dust of the order of 1 mm mainly accounts for the zodiacal light in infrared.

Recent radio observations show that the giant molecular cloud (GMC) mass functions noticeably vary across galactic disks. High-resolution magnetohydrodynamics simulations show that multiple episodes of compression are required for creating a molecular cloud in the magnetized interstellar medium. In this article, we formulate the evolution equation for the GMC mass function to reproduce the observed profiles, for which multiple compression are driven by the network of expanding shells due to HII regions and supernova remnants. We introduce the cloud-cloud collision (CCC) terms in the evolution equation in contrast to the previous work (Inutsuka et al. 2015). The computed time evolution suggests that the GMC mass function slope is governed by the ratio of GMC formation timescale to its dispersal timescale, and that the CCC effect is limited only in the massive-end of the mass function. In addition, we identify a gas resurrection channel that allows the gas dispersed by massive stars to regenerate GMC populations or to accrete onto the pre-existing GMCs. Our results show that almost all of the dispersed gas contribute to the mass growth of pre-existing GMCs in arm regions whereas less than 60 per cent in inter-arm regions. Our results also predict that GMC mass functions have a single power-law exponent in the mass range < 10^5.5 Msun (where Msun represents the solar mass), which is well characterized by GMC self-growth and dispersal timescales. Measurement of the GMC mass function slope provides a powerful method to constrain those GMC timescales and the gas resurrecting factor in various environment across galactic disks.

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

Standard galaxy formation models predict that large-scale double-lobed radio sources, known as DRAGNs, will always be hosted by elliptical galaxies. In spite of this, in recent years a small number of spiral galaxies have also been found to host such sources. These so-called spiral DRAGNs are still extremely rare, with only $\sim 5$ cases being widely accepted. Here we report on the serendipitous discovery of a new spiral DRAGN in data from the Giant Metrewave Radio Telescope (GMRT) at 322 MHz. The host galaxy, MCG+07-47-10, is a face-on late-type Sbc galaxy with distinctive spiral arms and prominent bulge suggesting a high black hole mass. Using WISE infra-red and GALEX UV data we show that this galaxy has a star formation rate of 0.16-0.75 M$_{\odot}$yr$^{-1}$, and that the radio luminosity is dominated by star-formation. We demonstrate that this spiral DRAGN has similar environmental properties to others of this class, but has a comparatively low radio luminosity of $L_{\rm 1.4GHz}$ = 1.12$\times$10$^{22}$ W Hz$^{-1}$, two orders of magnitude smaller than other known spiral DRAGNs. We suggest that this may indicate the existence of a previously unknown low-luminosity population of spiral DRAGNS.

Context: Debris disks are important observational clues for understanding planetary-system formation process. In particular, faint warm debris disks may be related to late planet formation near 1 AU. A systematic search of faint warm debris disks is necessary to reveal terrestrial planet formation. Aims: Faint warm debris disks show excess emission that peaks at mid-IR wavelengths. Thus we explore debris disks using the AKARI mid-IR all-sky point source catalog (PSC), a product of the second generation unbiased IR all-sky survey. Methods : We investigate IR excess emission for 678 isolated main-sequence stars for which there are 18 micron detections in the AKARI mid-IR all-sky catalog by comparing their fluxes with the predicted fluxes of the photospheres based on optical to near-IR fluxes and model spectra. The near-IR fluxes are first taken from the 2MASS PSC. However, 286 stars with Ks<4.5 in our sample have large flux errors in the 2MASS photometry due to saturation. Thus we have measured accurate J, H, and Ks band fluxes, applying neutral density (ND) filters for Simultaneous InfraRed Imager for Unbiased Survey (SIRIUS) on IRSF, the \phi 1.4 m near-IR telescope in South Africa, and improved the flux accuracy from 14% to 1.8% on average. Results: We identified 53 debris-disk candidates including eight new detections from our sample of 678 main-sequence stars. The detection rate of debris disks for this work is ~8%, which is comparable with those in previous works by Spitzer and Herschel. Conclusion: The importance of this study is the detection of faint warm debris disks around nearby field stars. At least nine objects have a large amount of dust for their ages, which cannot be explained by the conventional steady-state collisional cascade model.

Adhesion forces between submicrometer-sized silicate grains play a crucial role in the formation of silicate dust agglomerates, rocky planetesimals, and terrestrial planets. The surface energy of silicate dust particles is the key to their adhesion and rolling forces in a theoretical model based on the contact mechanics. Here we revisit the cohesion of amorphous silica spheres by compiling available data on the surface energy for hydrophilic amorphous silica in various circumstances. It turned out that the surface energy for hydrophilic amorphous silica in a vacuum is a factor of 10 higher than previously assumed. Therefore, the previous theoretical models underestimated the critical velocity for the sticking of amorphous silica spheres, as well as the rolling friction forces between them. With the most plausible value of the surface energy for amorphous silica spheres, theoretical models based on the contact mechanics are in harmony with laboratory experiments. Consequently, we conclude that silicate grains with a radius of $0.1~\mu$m could grow to planetesimals via coagulation in a protoplanetary disk. We argue that the coagulation growth of silicate grains in a molecular cloud is advanced either by organic mantles rather than icy mantles or, if there are no mantles, by nanometer-sized grain radius.

The latest observation of HL Tau by ALMA revealed spectacular concentric dust rings in its circumstellar disk. We attempt to explain the multiple ring structure as a consequence of aggregate sintering. Sintering is known to reduce the sticking efficiency of dust aggregates and occurs at temperatures slightly below the sublimation point of their constituent material. We here present a dust growth model incorporating sintering and use it to simulate global dust evolution due to sintering, coagulation, fragmentation, and radial inward drift in a modeled HL Tau disk. We show that aggregates consisting of multiple species of volatile ices experience sintering, collisionally disrupt, and pile up at multiple locations slightly outside the snow lines of the volatiles. At wavelengths of 0.87--1.3 mm, these sintering zones appear as bright, optically thick rings with a spectral slope of $\approx 2$, whereas the non-sintering zones as darker, optically thinner rings of a spectral slope of $\approx$ 2.3--2.5. The observational features of the sintering and non-sintering zones are consistent with those of the major bright and dark rings found in the HL Tau disk, respectively. Radial pileup and vertical settling occur simultaneously if disk turbulence is weak and if monomers constituting the aggregates are $\sim 1~{\rm \mu m}$ in radius. For the radial gas temperature profile of $T = 310(r/1~{\rm AU})^{-0.57}~{\rm K}$, our model perfectly reproduces the brightness temperatures of the optically thick bright rings, and reproduces their orbital distances to an accuracy of $\lesssim$ 30%.

The zodiacal emission, which is the thermal infrared (IR) emission from the interplanetary dust (IPD) in our Solar System, has been studied for a long time. Nevertheless, accurate modeling of the zodiacal emission has not been successful to reproduce the all-sky spatial distribution of the zodiacal emission, especially in the mid-IR where the zodiacal emission peaks. We therefore aim to improve the IPD cloud model based on Kelsall et al. 1998, using the AKARI 9 and 18 micron all-sky diffuse maps. By adopting a new fitting method based on the total brightness, we have succeeded in reducing the residual levels after subtraction of the zodiacal emission from the AKARI data and thus in improving the modeling of the zodiacal emission. Comparing the AKARI and the COBE data, we confirm that the changes from the previous model to our new model are mostly due to model improvements, but not temporal variations between the AKARI and the COBE epoch, except for the position of the Earth-trailing blob. Our results suggest that the size of the smooth cloud, a dominant component in the model, is by about 10% more compact than previously thought, and that the dust sizes are not large enough to emit blackbody radiation in the mid-IR. Furthermore we significantly detect an isotropically-distributed IPD component, owing to accurate baseline measurement with AKARI.

When planetesimals grow via collisions in a turbulent disk, stirring through density fluctuation caused by turbulence effectively increases the relative velocities between planetesimals, which suppresses the onset of runaway growth. We investigate the onset of runaway growth in a turbulent disk through simulations that calculate the mass and velocity evolution of planetesimals. When planetesimals are small, the average relative velocity between planetesimals, $v_{\rm r}$, is much greater than their surface escape velocity, $v_{\rm esc}$, so that runaway growth does not occur. As planetesimals become large via collisional growth, $v_{\rm r}$ approaches $v_{\rm esc}$. When $v_{\rm r} \approx 1.5 v_{\rm esc}$, runaway growth of the planetesimals occurs. During the oligarchic growth subsequent to runaway growth, a small number of planetary embryos produced via runaway growth become massive through collisions with planetesimals with radii of that at the onset of runaway growth, $r_{\rm p,run}$. We analytically derive $r_{\rm p,run}$ as a function of the turbulent strength. Growing $\sim 10\,M_\oplus$ embryos that are suitable to become the cores of Jupiter and Saturn requires $r_{\rm p,run} \sim 100$\,km, which is similar to the proposed fossil feature in the size distribution of main belt asteroids. In contrast, the formation of Mars as quickly as suggested from Hf-W isotope studies requires small planetesimals at the onset of runaway growth. Thus, the conditions required to form Mars, Jupiter, and Saturn and the size distribution of the main-belt asteroids indicate that the turbulence increased in amplitude relative to the sound speed with increasing distance from the young Sun.

The East-Asian VLBI Network (EAVN) is the international VLBI facility in East Asia and is conducted in collaboration with China, Japan, and Korea. The EAVN consists of VLBI arrays operated in each East Asian country, containing 21 radio telescopes and three correlators. The EAVN will be mainly operated at 6.7 (C-band), 8 (X-band), 22 (K-band), and 43 GHz (Q-band), although the EAVN has an ability to conduct observations at 1.6 - 129 GHz. We have conducted fringe test observations eight times to date at 8 and 22 GHz and fringes have been successfully detected at both frequencies. We have also conducted science commissioning observations of 6.7 GHz methanol masers in massive star-forming regions. The EAVN will be operational from the second half of 2017, providing complementary results with the FAST on AGNs, massive star-forming regions, and evolved stars with high angular resolution at cm- to mm-wavelengths.

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

Collisions are a fundamental process in planet formation. If colliding objects simply merge, a planetary object can grow. However, if the collision is disruptive, planetary growth is prevented. Therefore, the impact conditions under which collisions are destructive are important in understanding planet formation. So far, the critical specific impact energy for a disruptive collision Q_D^* has been investigated for various types of collisions between objects ranging in scale from centimeters to thousands of kilometers. Although the values of Q_D^* have been calculated numerically while taking into consideration various physical properties such as self-gravity, material strength, and porosity, the dependence of Q_D^* on numerical resolution has not been sufficiently investigated. In this paper, using the smoothed particle hydrodynamics (SPH) method, we performed numerical simulations of collisions between planetesimals at various numerical resolutions (from 2 x 10^4 to 5 x 10^6 SPH particles) and investigated the resulting variation in Q_D^*. The value of Q_D^* is shown to decrease as the number of SPH particles increases, and the difference between the Q_D^* values for the lowest and highest investigated resolutions is approximately a factor of two. Although the results for 5 x 10^6 SPH particles do not fully converge, higher-resolution simulations near the impact site show that the value of Q_D^* for the case with 5 x 10^6 SPH particles is close to the expected converged value. Although Q_D^* depends on impact parameters and material parameters, our results indicate that at least 5 x 10^6 SPH particles are required for numerical simulations in disruptive collisions to obtain the value of Q_D^* within 20% error.

In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-type stars have been observed. Here, we quantitatively estimate the total mass of ejected materials during the giant impact stages. We found that $\sim$0.4 times the Earth's mass is ejected in total throughout the giant impact stage. Ejected materials are ground down by collisional cascade until micron-sized grains are blown out by radiation pressure. The depletion timescale of these ejected materials is determined primarily by the mass of the largest body among them. We conducted high-resolution simulations of giant impacts to accurately obtain the mass of the largest ejected body. We then calculated the evolution of the debris disks produced by a series of giant impacts and depleted by collisional cascades to obtain the infrared excess evolution of the debris disks. We found that the infrared excess is almost always higher than the stellar infrared flux throughout the giant impact stage ($\sim$100 Myr) and is sometimes $\sim$10 times higher immediately after a giant impact. Therefore, giant impact stages would explain the infrared excess from most observed warm debris disks. The observed fraction of stars with warm debris disks indicates that the formation probability of our solar system-like terrestrial planets is approximately 10%.

The bulk chemical compositions of planets are uncertain, even for major elements such as Mg and Si. This is due to the fact that the samples available for study all originate from relatively shallow depths. Comparison of the stable isotope compositions of planets and meteorites can help overcome this limitation. Specifically, the non-chondritic Si isotope composition of the Earth's mantle was interpreted to reflect the presence of Si in the core, which can also explain its low density relative to pure Fe-Ni alloy. However, we have found that angrite meteorites display a heavy Si isotope composition similar to the lunar and terrestrial mantles. Because core formation in the angrite parent-body (APB) occurred under oxidizing conditions at relatively low pressure and temperature, significant incorporation of Si in the core is ruled out as an explanation for this heavy Si isotope signature. Instead, we show that equilibrium isotopic fractionation between gaseous SiO and solid forsterite at 1370 K in the solar nebula could have produced the observed Si isotope variations. Nebular fractionation of forsterite should be accompanied by correlated variations between the Si isotopic composition and Mg/Si ratio following a slope of 1, which is observed in meteorites. Consideration of this nebular process leads to a revised Si concentration in the Earth's core of 3.6 (+6.0/-3.6) wt% and provides estimates of Mg/Si ratios of bulk planetary bodies.

Planets are formed from collisional growth of small bodies in a protoplanetary disk. Bodies much larger than approximately $1$\,m are mainly controlled by the gravity of the host star and experience weak gas drag; their orbits are mainly expressed by orbital elements: semimajor axes $a$, eccentricities $e$, and inclinations $i$, which are modulated by gas drag. In a previous study, $\dot a$, $\dot e$, and $\dot i$ were analytically derived for $e \ll 1$ and $i \ll H/a$, where $H$ is the scale height of the disk. Their formulae are valid in the early stage of planet formation. However, once massive planets are formed, $e$ and $i$ increase greatly. Indeed, some small bodies in the solar system have very large $e$ and $i$. Therefore, in this paper, I analytically derive formulae for $\dot a$, $\dot e$, and $\dot i$ for $1-e^2 \ll 1$ and $i \ll H/a$ and for $i \gg H/a$. The formulae combined from these limited equations will represent the results of orbital integration unless $e \geq 1$ or $i > \pi - H/a$. Since the derived formulae are applicable for bodies not only in a protoplanetary disk but also in a circumplanetary disk, I discuss the possibility of the capture of satellites in a circumplanetary disk using the formulae.

Recent three-dimensional magnetohydrodynamical simulations have identified a disk wind by which gas materials are lost from the surface of a protoplanetary disk, which can significantly alter the evolution of the inner disk and the formation of terrestrial planets. A simultaneous description of the realistic evolution of the gaseous and solid components in a disk may provide a clue for solving the problem of the mass concentration of the terrestrial planets in the solar system. We simulate the formation of terrestrial planets from planetary embryos in a disk that evolves via magnetorotational instability and a disk wind. The aim is to examine the effects of a disk wind on the orbital evolution and final configuration of planetary systems. We perform N-body simulations of sixty 0.1 Earth-mass embryos in an evolving disk. The evolution of the gas surface density of the disk is tracked by solving a one-dimensional diffusion equation with a sink term that accounts for the disk wind. We find that even in the case of a weak disk wind, the radial slope of the gas surface density of the inner disk becomes shallower, which slows or halts the type I migration of embryos. If the effect of the disk wind is strong, the disk profile is significantly altered (e.g., positive surface density gradient, inside-out evacuation), leading to outward migration of embryos inside ~ 1 AU. Disk winds play an essential role in terrestrial planet formation inside a few AU by changing the disk profile. In addition, embryos can undergo convergent migration to ~ 1 AU in certainly probable conditions. In such a case, the characteristic features of the solar system's terrestrial planets (e.g., mass concentration around 1 AU, late giant impact) may be reproduced.

We present a comprehensive survey of diffuse interstellar bands (DIBs) in $0.91-1.32\mu$m with a newly-developed near-infrared (NIR) spectrograph, WINERED, mounted on the Araki 1.3 m Telescope in Japan. We obtained high-resolution ($R=28,300$) spectra of 25 early-type stars with color excesses of $0.07<E(B-V)<3.4$. In addition to the five DIBs previously detected in this wavelength range, we identified 15 new DIBs, 7 of which were reported as DIB "candidates" by Cox. We analyze the correlations among NIR DIBs, strong optical DIBs, and the reddening of the stars. Consequently, we found that all NIR DIBs show weaker correlations with the reddening rather than the optical DIBs, suggesting that the equivalent widths of NIR DIBs depend on some physical conditions of the interstellar clouds, such as UV flux. Three NIR DIBs, $\lambda\lambda$10780, 10792, and 11797, are found to be classifiable as a "family," in which the DIBs are well correlated with each other, suggesting that the carriers of these DIBs are connected with some chemical reactions and/or have similar physical properties such as ionization potential. We also found that three strongest NIR DIBs $\lambda\lambda10780, 11797, 13175$ are well correlated with the optical DIB $\lambda 5780.5$, whose carrier is proposed to be a cation molecule with high ionization potential, indicating that the carriers of the NIR DIBs could be cation molecules.

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

Several hundred stars older than 10 million years have been observed to have infrared excesses. These observations are explained by dust grains formed by the collisional fragmentation of hidden planetesimals. Such dusty planetesimal discs are known as debris discs. In a dynamically cold planetesimal disc, collisional coagulation of planetesimals produces planetary embryos which then stir the surrounding leftover planetesimals. Thus, the collisional fragmentation of planetesimals that results from planet formation forms a debris disc. We aim to determine the properties of the underlying planetesimals in debris discs by numerically modelling the coagulation and fragmentation of planetesimal populations. The brightness and temporal evolution of debris discs depend on the radial distribution of planetesimal discs, the location of their inner and outer edges, their total mass, and the size of planetesimals in the disc. We find that a radially narrow planetesimal disc is most likely to result in a debris disc that can explain the trend of observed infrared excesses of debris discs around G-type stars, for which planet formation occurs only before 100 million years. Early debris disc formation is induced by planet formation, while the later evolution is explained by the collisional decay of leftover planetesimals around planets that have already formed. Planetesimal discs with underlying planetesimals of radii $\sim 100\,$km at $\approx 30$ AU most readily explain the Spitzer Space Telescope 24 and 70$ \mu$m fluxes from debris discs around G-type stars.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We report the results of the measurement of the trigonometric parallax of an H$_2$O maser source in \r22 with the VLBI Exploration of Radio Astrometry (VERA). The annual parallax was determined to be 0.314$\pm$0.070~mas, corresponding to a distance of 3.18$^{+0.90}_{-0.66}$~kpc. Our results confirm \r22 to be located in the Perseus arm. We computed the peculiar motion of \r22 to be $(U_{src}, V_{src}, W_{src}) = (4\pm4, -32\pm6, 8\pm6)$~\kms, where $U_{src}$, $V_{src}$, and $W_{src}$ are directed toward the Galactic center, in the direction of Galactic rotation and toward the Galactic north pole, respectively. \r22, NGC\,7538 and Cepheus A lie along the same line of sight, and are within 2$^{\circ}$ on the sky. Their parallax distances with which we derived their absolute position in the Milky Way show that \r22 and \ngc~are associated with the Perseus Arm, while \cepa~is located in the Local Arm. We compared the kinematic distances of \iras, \ga, \gb, \gc, \gd, \cepa, \ngc, \ge, \gf~and \r22~derived with flat and non-flat rotation curve with its parallax distance and found the kinematic distance derived from the non-flat rotation assumption ($-$5 to $-$39 \kms~lag) to be consistent with the parallax distance.

We investigate the formation of multiple-planet systems in the presence of a hot Jupiter using extended N-body simulations that are performed simultaneously with semi-analytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of a hot Jupiter. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal-mass; the largest planet constitutes more than 50 percent of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper (Ogihara et al. 2013), we have found a new physical mechanism of induced migration of the hot Jupiter, which is called a crowding-out. If the hot Jupiter opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the hot Jupiter remains. We also discuss angular momentum transfer between the planets and disk.

We conducted phase referencing VLBI observations of the Mira variable T~Lepus (T~Lep) using VERA, from 2003 to 2006. The distance to the source was determined from its annual parallax which was measured to be 3.06$\pm$0.04 mas, corresponding to a distance of 327$\pm$4\,pc. Our observations revealed the distribution and internal kinematics of H$_2$O masers in T~Lep, and we derived a source systemic motions of 14.60$\pm$0.50 mas yr$^{-1}$ and $-$35.43$\pm$0.79 mas yr$^{-1}$ in right ascension and declination, respectively. We also determined a LSR velocity of $V_\mathrm{LSR}^{\ast} = -$27.63 km s$^{-1}$. Comparison of our result with an image recently obtained from the VLTI infrared interferometer reveals a linear scale picture of the circumstellar structure of T~Lep. Analysis of the source systemic motion in the Galacto-centric coordinate frame indicates a large peculiar motion, which is consistent with the general characteristics of AGB stars. This source makes a contribution to the calibration of the period-luminosity relation of Galactic Mira variables. From the compilation of data for nearby Mira variables found in the literature, whose distances were derived from astrometric VLBI observations, we have calibrated the Galactic Mira period-luminosity relation to a high degree of accuracy.

We carried out the photometric observations of the SU UMa-type dwarf nova ER UMa during 2011 and 2012, which showed the existence of persistent negative superhumps even during the superoutburst. We performed two-dimensional period analysis of its light curves by using a method called "least absolute shrinkage and selection operator" (Lasso) and "phase dispersion minimization" (PDM) analysis, and we found that the period of negative superhumps systematically changed between a superoutburst and the next superoutburst. The trend of the period change can beinterpreted as reflecting the change of the disk radius. This change of the disk radius is in good agreement with the predicted change of the disk radius by the thermal-tidal instability (TTI) model. The normal outbursts within a supercycle showed a general trend that the rising rate to maximum becomes slower as the next superoutburst approaches. The change can be interpreted as the consequence of the increased gas-stream flow onto the inner region of the disk as the result of the tilted disk. Some of the superoutbursts were found to be triggered by a precursor normal outburst when the positive superhumps appeared to develop. The positive and negative superhumps co-existed during the superoutburst. The positive superhumps were prominent only during four or five days after the supermaximum, while the signal of the negative superhumps became strong after the middle phase of the superoutburst plateau. A simple combination of the positive and negative superhumps was found to be insufficient in reproducing the complex profile variation. We were able to detect the developing phase of positive superhumps (stage A superhumps) for the first time in ER UMa-type dwarf novae. Using the period of stage A superhumps, we obtained a mass ratio of 0.100(15), which indicates that ER UMa is on the ordinary evolutional track of CVs.

We have measured the internal proper motions of the 6.7 GHz methanol masers associated with Cepheus A (Cep A) HW2 using Very Long Baseline Interferometery (VLBI) observations. We conducted three epochs of VLBI monitoring observations of the 6.7 GHz methanol masers in Cep A-HW2 with the Japanese VLBI Network (JVN) over the period between 2006-2008. In 2006, we were able to use phase-referencing to measure the absolute coordinates of the maser emission with an accuracy of a few milliarcseconds. We compared the maser distribution with other molecular line observations that trace the rotating disk. We measured the internal proper motions for 29 methanol maser spots, of which 19 were identified at all three epochs and the remaining ten at only two epochs. The magnitude of proper motions ranged from 0.2 to 7.4 km/s, with an average of 3.1 km/s. Although there are large uncertainties in the observed internal proper motions of the methanol maser spots in Cep A, they are well fitted by a disk that includes both rotation and infall velocity components. The derived rotation and infall velocities at the disk radius of 680 au are 0.5 +- 0.7 and 1.8 +- 0.7 km/s, respectively. Assuming that the modeled disk motion accurately represents the accretion disk around the Cep A-HW2 high-mass YSO, we estimated the mass infall rate to be 3 x 10^{-4} n_8 Msun/yr (n_8 is the gas volume density in units of 10^{8} cm^{-3}). The combination of the estimated mass infall rate and the magnitude of the fitted infall velocity suggests that Cep A-HW2 is at an evolutionary phase of active gas accretion from the disk onto the central high-mass YSO. The infall momentum rate is estimated to be 5 x 10^{-4} n_8 Msun/yr km/s, which is larger than the estimated stellar radiation pressure of the HW2 object, supporting the hypothesis that this object is in an active gas accretion phase.

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion.

Continuing the project described by Kato et al. (2009a, arXiv:0905.1757), we collected times of superhump maxima for SU UMa-type dwarf novae mainly observed during the 2012-2013 season. We found three objects (V444 Peg, CSS J203937 and MASTER J212624) having strongly positive period derivatives despite the long orbital period (Porb). By using the period of growing stage (stage A) superhumps, we obtained mass ratios for six objects. We characterized nine new WZ Sge-type dwarf novae. We made a pilot survey of the decline rate of slowly fading part of SU UMa-type and WZ Sge-type outbursts. The decline time scale was found to generally follow the expected Porb^(1/4) dependence and WZ Sge-type outbursts also generally follow this trend. There are some objects which show slower decline rates, and we consider these objects good candidates for period bouncers. We also studied unusual behavior in some objects, including BK Lyn which made a transition from an ER UMa-type state to the novalike (standstill) state in 2013 and unusually frequent occurrence of superoutbursts in NY Ser and CR Boo. We applied least absolute shrinkage and selection operator (Lasso) power spectral analysis, which has been proven to be very effective in analyzing the Kepler data, to ground-based photometry of BK Lyn and detected the dramatic disappearance of the signal of negative superhumps in 2013. We suggested that the mass-transfer rates did not vary strongly between the ER UMa-type state and novalike state in BK Lyn, and this transition was less likely caused by a systematic variation of the mass-transfer rate.

We investigate formation of close-in terrestrial planets from planetary embryos under the influence of a hot Jupiter (HJ) using gravitational N-body simulations that include gravitational interactions between the gas disk and the terrestrial planet (e.g., type I migration). Our simulations show that several terrestrial planets efficiently form outside the orbit of the HJ, making a chain of planets, and all of them gravitationally interact directly or indirectly with the HJ through resonance, which leads to inward migration of the HJ. We call this mechanism of induced migration of the HJ as "crowding out." The HJ is eventually lost by collision with the central star, and only several terrestrial planets remain. We also find that the efficiency of the crowding-out effect depends on model parameters; for example, the heavier the disk is, the more efficient the crowding out is. When planet formation occurs in a massive disk, the HJ can be lost to the central star and is never observed. On the other hand, for a less massive disk, the HJ and terrestrial planets can coexist; however, the companion planets can be below the detection limit of current observations. In both cases, systems with the HJ and terrestrial planets have little chance for detection. Therefore, our model naturally explains the lack of companion planets in HJ systems regardless of the disk mass. In effect, our model provide a theoretical prediction for future observations; additional planets can be discovered just outside the HJ, and their masses should generally be small.

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

We report the results of multi-epoch observations of Sgr A* with VLBI Exploration of Radio Astrometry (VERA) at 43 GHz, carried out from 2004 to 2008. We detected a time variation of flux at 11 % level and intrinsic size at 19 %. In addition, comparisons with previous Very Long Baseline Array (VLBA) results shows that Sgr A* underwent the flaring event at least longer than 10 days in May 2007. The intrinsic size of Sgr A* remained unchanged within 1 $\sigma$ level from the size before/after the flaring event, indicating that the brightness temperature of Sgr A* was increased. The flaring event occurred within 31 d, which is shorter than the refractive time scale. Moreover it is difficult to explain the increase in the spectral index at the flaring event by the simple interstellar scattering model. Hence, the flaring event is most likely associated with the changes in intrinsic properties of Sgr A*. We considered the origin of the brightness temperature variation, and concluded that the flaring event of Sgr A* could be explained by the continuous heating of electrons, such as a standing shock in accretion flow.

A maintenance mechanism of an approximately linear velocity profile of the Venus zonal flow or superrotation is explored, with the aid of a Reynolds-averaged turbulence modeling approach. The basic framework is similar to that of Gierasch (1975) in the sense that the mechanism is examined under a given meridional circulation. The profile mimicking the observations of the flow is initially assumed, and its maintenance mechanism in the presence of turbulence effects is investigated from a viewpoint of the suppression of energy cascade. In the present work, the turbulent viscosity is regarded as an indicator of the intensity of the cascade. A novelty of this formalism is the use of the isotropic turbulent viscosity based on a nonlocal time scale linked to a large-scale flow structure. The mechanism is first discussed qualitatively. On the basis of these discussions, the two-dimensional numerical simulation of the proposed model is performed, with an initially assumed superrotation, and the fast zonal flow is shown to be maintained, compared with the turbulent viscosity lacking the nonlocal time scale. The relationship of the present model with the current general-circulation-model simulation is discussed in light of a crucial role of the vertical viscosity.

We investigate the condition for capture into first-order mean motion resonances using numerical simulations with a wide range of various parameters. In particular, we focus on deriving the critical migration timescale for capture into the 2:1 resonance; additional numerical experiments for closely spaced resonances (e.g., 3:2) are also performed. We find that the critical migration timescale is determined by the planet-to-stellar mass ratio, and its dependence exhibits power-law behavior with index -4/3. This dependence is also supported by simple analytic arguments. We also find that the critical migration timescale for systems with equal-mass bodies is shorter than that in the restricted problem; for instance, for the 2:1 resonance between two equal-mass bodies, the critical timescale decreases by a factor of 10. In addition, using the obtained formula, the origin of observed systems that include first-order commensurabilities is constrained. Assuming that pairs of planets originally form well separated from each other and then undergo convergent migration and are captured in resonances, it is possible that a number of exoplanets experienced rapid orbital migration. For systems in closely spaced resonances, the differential migration timescale between the resonant pair can be constrained well; it is further suggested that several exoplanets underwent migration that can equal or even exceed the type I migration rate predicted by the linear theory. This implies that some of them may have formed in situ. Future observations and the use of our model will allow us to statistically determine the typical migration speed in a protoplanetary disk.

Mars is likely to be a planetary embryo formed through collisions with planetesimals, which can explain its small mass and rapid formation timescale obtained from 182Hf-182$W chronometry. In the classical theory of planet formation, the final embryo mass is determined only by the solid surface density. However, embryos can stir surrounding planetesimals, leading to fragmentation through erosive (cratering) collisions. We find that radial drift of small fragments can drastically reduce the solid surface density. On the other hand, embryo growth is accelerated by fragment accretion. Since collisional fragmentation efficiency depends on the initial size of planetesimals, the final embryo mass and its growth time are determined by the initial planetesimal size and disk surface density. We have investigated the effect of these two parameters on the mass of Mars and the predicted radiogenic excess of 182W in the martian mantle. Two scenarios can explain the rapid formation of small Mars: (i) it formed by accretion of small planetesimals in a massive disk or (ii) it formed from large planetesimals but its growth was arrested by the inward then outward migration of Jupiter. Taking into account all constraints, we conclude that Mars is likely to have formed in a massive disk of about ~ 0.1 solar mass from planetesimals smaller than ~ 10 km in radius. Such small planetesimal size cannot explain core accretion of Jupiter, suggesting that there may have been a heliocentric gradient in planetesimal size in the solar nebula.

Shattering of dust grains in the interstellar medium is a viable mechanism of small grain production in galaxies. We examine the robustness or uncertainty in the theoretical predictions of shattering. We identify $P_1$ (the critical pressure above which the deformation destroys the original lattice structures) as the most important quantity in determining the timescale of small grain production, and confirm that the same $P_1/t$ ($t$ is the duration of shattering) gives the same grain size distribution [$n(a)$, where $a$ is the grain radius] after shattering within a factor of 3. The uncertainty in the fraction of shocked material that is eventually ejected as fragments causes uncertainties in $n(a)$ by a factor of 1.3 and 1.6 for silicate and carbonaceous dust, respectively. The size distribution of shattered fragments have minor effects as long as the power index of the fragment size distribution is less than ~ 3.5, since the slope of grain size distribution $n(a)$ continuously change by shattering and becomes consistent with $n(a)\propto a^{-3.5}$. The grain velocities as a function of grain radius can have an imprint in the grain size distribution especially for carbonaceous dust. We also show that the formulation of shattering can be simplified without losing sufficient precision.

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we studied 86 SU UMa-type dwarf novae. We confirmed the general trends such as the relation between period derivatives and orbital periods. There are some systems showing positive period derivatives despite the long orbital periods. We observed the 2011 outburst of the WZ Sge-type dwarf nova BW Scl, and recorded an O-C diagram similar to those of previously known WZ Sge-type dwarf novae. The WZ Sge-type dwarf nova OT J184228.1+483742 showed an unusual pattern of double outbursts composed of an outburst with early superhumps and one with ordinary superhumps. We propose an interpretation that a very small growth rate of the 3:1 resonance due to an extremely low mass-ratio led to a quenching of the superoutburst before the ordinary superhumps appeared. We studied ER UMa-type dwarf novae and found that V1159 Ori showed positive superhumps similar to ER UMa in the 1990s. The recently recognized ER UMa-type object BK Lyn dominantly showed negative superhumps, and its behavior was very similar to the present-day state of ER UMa. The pattern of period variations in AM CVn-type objects were very similar to short-period hydrogen-rich SU UMa-type dwarf novae, making them helium analogue of hydrogen-rich SU UMa-type dwarf novae. SBS 1108+574, a peculiar hydrogen-rich dwarf nova below the period minimum, showed a very similar pattern of period variations to those of short-period SU UMa-type dwarf novae. The mass-ratio derived from the detected orbital period suggests that this secondary is a somewhat evolved star whose hydrogen envelope was mostly stripped during the mass-exchange. CC Scl, MASTER OT J072948.66+593824.4 and OT J173516.9+154708 showed only low-amplitude superhumps with complex profiles. These superhumps are likely a combination of closely separated two periods.

This paper reports very-long-baseline interferometry observations of the radio-loud broad absorption line (BAL) quasar J1020+4320 at 1.7, 2.3, 6.7, and 8.4 GHz using the Japanese VLBI network (JVN) and European VLBI network (EVN). The radio morphology is compact with a size of ~10 pc. The convex radio spectrum is stable over the last decade; an observed peak frequency of 3.2 GHz is equivalent to 9.5 GHz in the rest frame, suggesting an age of the order of ~100 years as a radio source, according to an observed correlation between linear size and peak frequency of compact steep spectrum (CSS) and giga-hertz peaked spectrum (GPS) radio sources. A low-frequency radio excess suggests relic of past jet activity. J1020+4320 may be one of the quasars with recurrent and short-lived jet activity during a BAL-outflowing phase.

The trajectories of dust particles ejected from a comet are affected by solar radiation pressure as a function of their ratios of radiation pressure cross section to mass. Therefore, a study on the orbital evolution of the particles caused by the radiation pressure reveals the physical properties of dust on the surface of the comet nucleus. In the course of NASA's Deep Impact mission, the ejecta plume evolved under the influence of the radiation pressure. From the evolution and shape of the plume, we have succeeded in obtaining $\beta \approx 0.4$, where $\beta$ is the ratio of the radiation pressure to the solar gravity. Taking into account $\beta \approx 0.4$ as well as the observational constraints of a high color temperature and a small silicate-feature strength, dust particles ejected from the surface of comet 9P/Tempel 1 are likely compact dust aggregates of sizes $\approx 20\,\mu$m (mass $\sim 10^{-8}$\,g). This is comparable to the major dust on the surface of comet 1P/Halley ($\sim 10\mu$m) inferred from in-situ measurements and theoretical considerations. Since such dust aggregates with $\beta \approx 0.4$ must have survived on the surface against jets due to ice sublimation, the temperature of ice in the nucleus must be kept below 145\,K, which is much lower than equilibrium temperature determined by solar irradiation and thermal emission. These facts indicate that 9P/Tempel 1 has a dust mantle composed of $20\,\mu$m-sized dust aggregates with low thermal conductivities $\sim 1 \, {\rm erg\, cm}^{-1} \, {\rm K}^{-1}\,{\rm s}^{-1}$.

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 present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

We present measurements of the linear diameter of the emission region of the Vela pulsar at observing wavelength lambda=18 cm. We infer the diameter as a function of pulse phase from the distribution of visibility on the Mopra-Tidbinbilla baseline. As we demonstrate, in the presence of strong scintillation, finite size of the emission region produces a characteristic W-shaped signature in the projection of the visibility distribution onto the real axis. This modification involves heightened probability density near the mean amplitude, decreased probability to either side, and a return to the zero-size distribution beyond. We observe this signature with high statistical significance, as compared with the best-fitting zero-size model, in many regions of pulse phase. We find that the equivalent full width at half maximum of the pulsar's emission region decreases from more than 400 km early in the pulse to near zero at the peak of the pulse, and then increases again to approximately 800 km near the trailing edge. We discuss possible systematic effects, and compare our work with previous results.

We compare the noise in interferometric measurements of the Vela pulsar from ground- and space-based antennas with theoretical predictions. The noise depends on both the flux density and the interferometric phase of the source. Because the Vela pulsar is bright and scintillating, these comparisons extend into both the low and high signal-to-noise regimes. Furthermore, our diversity of baselines explores the full range of variation in interferometric phase. We find excellent agreement between theoretical expectations and our estimates of noise among samples within the characteristic scintillation scales. Namely, the noise is drawn from an elliptical Gaussian distribution in the complex plane, centered on the signal. The major axis, aligned with the signal phase, varies quadratically with the signal, while the minor axis, at quadrature, varies with the same linear coefficients. For weak signal, the noise approaches a circular Gaussian distribution. Both the variance and covariance of the noise are also affected by artifacts of digitization and correlation. In particular, we show that gating introduces correlations between nearby spectral channels.

We have investigated Saturn's core formation at a radial pressure maximum in a protoplanetary disk, which is created by gap opening by Jupiter. A core formed via planetesimal accretion induces the fragmentation of surrounding planetesimals, which generally inhibits further growth of the core by removal of the resulting fragments due to radial drift caused by gas drag. However, the emergence of the pressure maximum halts the drift of the fragments, while their orbital eccentricities and inclinations are efficiently damped by gas drag. As a result, the core of Saturn rapidly grows via accretion of the fragments near the pressure maximum. We have found that in the minimum-mass solar nebula, kilometer sized planetesimals can produce a core exceeding 10 Earth masses within two million years. Since Jupiter may not have undergone significant type II inward migration, it is likely that Jupiter's formation was completed when the local disk mass has already decayed to a value comparable to or less than Jovian mass. The expected rapid growth of Saturn's core on a timescale comparable to or shorter than observationally inferred disk lifetime enables Saturn to acquire the current amount of envelope gas before the disk gas is completely depleted. The high heat energy release rate onto the core surface due to the rapid accretion of the fragments delays onset of runaway gas accretion until the core mass becomes somewhat larger than that of Jupiter, which is consistent with the estimate based on interior modeling. Therefore, the rapid formation of Saturn induced by gap opening of Jupiter can account for the formation of multiple gas giants (Jupiter and Saturn) without significant inward migration and larger core mass of Saturn than that of Jupiter.

Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles' aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.

(Abridged) We access the expected EKB dust disk properties by modeling. We treat the debiased population of the known transneptunian objects (TNOs) as parent bodies and generate the dust with our collisional code. The resulting dust distributions are modified to take into account the influence of gravitational scattering and resonance trapping by planets on migrating dust grains as well as the effect of sublimation. A difficulty is that the amount and distribution of dust are largely determined by sub-kilometer-sized bodies. These are directly unobservable, and their properties cannot be accessed by collisional modeling, because objects larger than 10...60m in the present-day EKB are not in a collisional equilibrium. To place additional constraints, we use in-situ measurements of the New Horizons spacecraft within 20AU. We show that the TNO population has to have a break in the size distribution at s<70km. However, even this still leaves us with several models that all correctly reproduce a nearly constant dust impact rates in the region of giant planet orbits and do not violate the constraints from the non-detection of the EKB dust thermal emission by the COBE spacecraft. The modeled EKB dust disks, which conform to the observational constraints, can either be transport-dominated or intermediate between the transport-dominated and collision-dominated regime. The in-plane optical depth of such disks is tau(r>10AU)~10^-6 and their fractional luminosity is f_d~10^-7. Planets and sublimation are found to have little effect on dust impact fluxes and dust thermal emission. The spectral energy distribution of an EKB analog, as would be seen from 10pc distance, peaks at wavelengths of 40...50\mum at F~0.5mJy, which is less than 1% of the photospheric flux at those wavelengths. Therefore, exact EKB analogs cannot be detected with present-day instruments such as Herschel/PACS.

We present distance measurement of the semiregular variable RX Bootis (RX Boo) with its annual parallax. Using the unique dual-beam system of the VLBI Exploration of Radio Astrometry (VERA) telescope, we conducted astrometric observations of a water maser spot accompanying RX Boo referred to the quasar J1419+2706 separated by 1.69 degrees from RX Boo. We have measured the annual parallax of RX Boo to be 7.31 +/- 0.50 mas, corresponding to a distance of 136 +10/-9 pc, from the one-year monitoring observation data of one maser spot at VLSR = 3.2 km/s. The distance itself is consistent with the one obtained with Hipparcos. The distance uncertainty is reduced by a factor of two, allowing us to determine the stellar properties more accurately. Using our distance, we discuss the location of RX Boo in various sequences of Period-Luminosity (PL) relations. We found RX Boo is located in the Mira sequence of PL relation. In addition, we calculated the radius of photosphere and the mass limitation of RX Boo and discussed its evolutionary status.

Continuing the project described by Kato et al. (2009, PASJ 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 51 SU UMa-type dwarf novae mainly observed during the 2010-2011 season. Although most of the new data for systems with short superhump periods basically confirmed the findings by Kato et al. (2009) and Kato et al. (2010, PASJ 62, 1525, arXiv:1009.5444), the long-period system GX Cas showed an exceptionally large positive period derivative. An analysis of public Kepler data of V344 Lyr and V1504 Cyg yielded less striking stage transitions. In V344 Lyr, there was prominent secondary component growing during the late stage of superoutbursts, and the component persisted at least for two more cycles of successive normal outbursts. We also investigated the superoutbursts of two conspicuous eclipsing objects: HT Cas and the WZ Sge-type object SDSS J080434.20+510349.2. Strong beat phenomena were detected in both objects, and late-stage superhumps in the latter object had an almost constant luminosity during the repeated rebrightenings. The WZ Sge-type object SDSS J133941.11+484727.5 showed a phase reversal around the rapid fading from the superoutburst. The object showed a prominent beat phenomenon even after the end of the superoutburst. A pilot study of superhump amplitudes indicated that the amplitudes of superhumps are strongly correlated with orbital periods, and the dependence on the inclination is weak in systems with inclinations smaller than 80 deg.

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.

The formation of giant planets requires accumulation of ~10 Earth mass in solids; but how do protoplanets acquire their mass? There are many, often competing processes that regulate the accretion rate of protoplanets. To assess their effects we present a new, publicly-available toy model. The rationale behind the toy model is that it encompasses as many physically-relevant processes as possible, but at the same time does not compromise its simplicity, speed, and physical insight. The toy model follows a modular structure, where key features -- e.g. planetesimal fragmentation, radial orbital decay, nebula turbulence -- can be switched on or off. Our model assumes three discrete components (fragments, planetesimals, and embryos) and is zero dimensional in space. We have tested the outcomes of the toy model against literature results and generally find satisfactory agreement. We include, for the first time, model features that capture the three-way interactions among small particles, gas, and protoplanets. Collisions among planetesimals will result in fragmentation, transferring a substantial amount of the solid mass to small particles, which couple strongly to the gas. Our results indicate that the efficiency of the accretion process then becomes very sensitive to the gas properties -- especially to the turbulent state and the magnitude of the disk headwind (the decrease of the orbital velocity of the gas with respect to Keplerian) -- as well as to the characteristic fragment size.

We report the detection of bursts of 6.7 GHz methanol maser emission in a high-mass star-forming region, G33.64-0.21. One of the spectral components of the maser in this source changed its flux density by 7 times that of the previous day, and it decayed with a timescale of 5 days. The burst occurred repeatedly in the spectral component, and no other components showed such variability. A VLBI observation with the Japanese VLBI Network (JVN) showed that the burst location was at the southwest edge of a spatial distribution, and the bursting phenomenon occurred in a region much smaller than 70 AU. We suggest an impulsive energy release like a stellar flare as a possible mechanism for the burst. These results imply that 6.7 GHz methanol masers could be a useful new probe for studying bursting activity in the process of star formation of high-mass YSOs with a high-resolution of AU scale.

In February 2011, a burst event of the H$_{2}$O maser in Orion KL (Kleinmann-Low object) has started after 13-year silence. This is the third time to detect such phenomena in Orion KL, followed by those in 1979-1985 and 1998. We have carried out astrometric observations of the bursting H$_{2}$O maser features in Orion KL with VERA (VLBI Exploration of Radio Astrometry), a Japanese VLBI network dedicated for astrometry. The total flux of the bursting feature at the LSR velocity of 7.58 km s$^{-1}$ reaches 4.4$\times10^{4}$ Jy in March 2011. The intensity of the bursting feature is three orders of magnitudes larger than that of the same velocity feature in the quiescent phase in 2006. Two months later, another new feature appears at the LSR velocity of 6.95 km s$^{-1}$ in May 2011, separated by 12 mas north of the 7.58 km s$^{-1}$ feature. Thus, the current burst occurs at two spatially different features. The bursting masers are elongated along the northwest-southeast direction as reported in the previous burst in 1998. We determine the absolute positions of the bursting features for the first time ever with a submilli-arcsecond (mas) accuracy. Their positions are coincident with the shocked molecular gas called the Orion Compact Ridge. We tentatively detect the absolute proper motions of the bursting features toward southwest direction. It is most likely that the outflow from the radio source I or another young stellar object interacting with the Compact Ridge is a possible origin of the H$_{2}$O maser burst.

We surveyed the (1,1), (2,2), and (3,3) lines of NH3 and the H2O maser toward the molecular cloud L935 in the extended HII region NGC 7000 with an angular resolution of 1.6' using the Kashima 34-m telescope. We found five clumps in the NH3 emission with a size of 0.2--1 pc and mass of 9--452 M_sun. The molecular gas in these clumps has a similar gas kinetic temperature of 11--15 K and a line width of 1--2 km/s. However, they have different star formation activities such as the concentration of T-Tauri type stars and the association of H2O maser sources. We found that these star formation activities are related to the geometry of the HII region. The clump associated with the T-Tauri type star cluster has a high star formation efficiency of 36--62%. This clump is located near the boundary of the HII region and molecular cloud. Therefore, we suggest that the star formation efficiency increases because of the triggered star formation.

Massive planetary cores ($\sim 10$ Earth masses) trigger rapid gas accretion to form gas giant planets \rev{such as} Jupiter and Saturn. We investigate the core growth and the possibilities for cores to reach such a critical core mass. At the late stage, planetary cores grow through collisions with small planetesimals. Collisional fragmentation of planetesimals, which is induced by gravitational interaction with planetary cores, reduces the amount of planetesimals surrounding them, and thus the final core masses. Starting from small planetesimals that the fragmentation rapidly removes, less massive cores are formed. However, planetary cores acquire atmospheres that enlarge their collisional cross section before rapid gas accretion. Once planetary cores exceed about Mars mass, atmospheres significantly accelerate the growth of cores. We show that, taking into account the effects of fragmentation and atmosphere, initially large planetesimals enable formation of sufficiently massive cores. On the other hand, because the growth of cores is slow for large planetesimals, a massive disk is necessary for cores to grow enough within a disk lifetime. If the disk with 100\,km-sized initial planetesimals is 10 times as massive as the minimum mass solar nebula, planetary cores can exceed 10 Earth masses in the Jovian planet region ($>5\,$AU).

We have performed VLBI and single-dish monitoring of 22 GHz H$_{2}$O maser emission from the high mass young stellar object G353.273+0.641 with VERA (VLBI Exploration of Radio Astrometry) and Tomakamai 11-m radio telescope. Two maser flares have been detected, separated almost two years. Frequent VLBI monitoring has revealed that these flare activities have been accompanied by structural change of the prominent shock front traced by H2O maser alignments. We have detected only blue-shifted emissions and all maser features have been distributed within very small area of 200 $\times$ 200 au$^{2}$ in spite of wide velocity range (> 100 km s$^{-1}$). The light curve shows notably intermittent variation and suggests that the H$_{2}$O masers in G353.273+0.641 are excited by episodic radio jet. The time-scale of \sim2 yr and characteristic velocity of \sim500 km s$^{-1}$ also support this interpretation. Two isolated velocity components of C50 (-53 \pm 7 km s$^{-1}$) and C70 (-73 \pm 7 km s$^{-1}$) have shown synchronised linear acceleration of the flux weighted $V_{\rmn{LSR}}$ values (\sim-5 km s$^{-1}$ yr$^{-1}$) during the flare phase. This can be converted to the lower-limit momentum rate of 1.1 \times 10$^{-3}$ M$_{\sun}$ km s$^{-1}$ yr$^{-1}$. Maser properties are quite similar to that of IRAS 20126+4104 especially. This corroborates the previous suggestion that G353.273+0.641 is a candidate of high mass protostellar object. The possible pole-on geometry of disc-jet system can be suitable for direct imaging of the accretion disc in this case.

Dust particles in orbit around a star drift toward the central star by the Poynting-Robertson effect and pile up by sublimation. We analytically derive the pile-up magnitude, adopting a simple model for optical cross sections. As a result, we find that the sublimation temperature of drifting dust particles plays the most important role in the pile-up rather than their optical property does. Dust particles with high sublimation temperature form a significant dust ring, which could be found in the vicinity of the sun through in-situ spacecraft measurements. While the existence of such a ring in a debris disk could not be identified in the spectral energy distribution (SED), the size of a dust-free zone shapes the SED. Since we analytically obtain the location and temperature of sublimation, these analytical formulae are useful to find such sublimation evidences.

We performed the astrometry of H2O masers in the Galactic star-forming region G48.61+0.02 with the VLBI Exploration of Radio Astrometry (VERA). We derived a trigonometric parallax of 199+/-7 micro as, which corresponds to a distance of 5.03+/-0.19 kpc. The distance to G48.61+0.02 is about a half of its far kinematic distance, which was often assumed previously. This distance places G48.61+0.02 in the Sagittarius-Carina arm and near the active star forming region and the supernova remnant W51. We also obtained the three dimensional motion of G48.61+0.02, and found that it has a large peculiar motion of 40+/-5 km/s. This peculiar motion would be originated with the multiple supernovae explosions in W51, or the streaming motion across the Sagittarius-Carina arm.

We present the results of astrometric observations with VERA toward the H2O maser sources in IRAS 05137+3919, which is thought to be located in the far outer Galaxy. We have derived the parallax of \pi = 0.086 +/- 0.027 mas, which corresponds to the source distance of D=11.6+5.3-2.8 kpc. Although the parallax measurement is only 3-sigma level and thus the distance uncertainty is considerably large, we can strongly constrain the minimum distance to this source, locating the source at the distance from the Sun greater than 8.3 kpc (or 16.7 kpc from the Galaxy's center) at 90% confidence level. Our results provide an astrometric confirmation that this source is located in the far outer Galaxy beyond 15 kpc from the Galaxy center, indicating that IRAS 05137+3919 is one of the most distant star-forming regions from the Galaxy center.

We conducted the astrometry of H2O masers in the Galactic star-forming region Onsala 1 (ON1) with VLBI Exploration of Radio Astrometry (VERA). We measured a trigonometric parallax of 0.404+/-0.017 mas, corresponding to a distance of 2.47+/-0.11 kpc. ON1 is appeared to be located near the tangent point at the Galactic longitude of 69.54 deg. We estimate the angular velocity of the Galactic rotation at Sun, the ratio of the distance from Sun to the Galactic center and the Galactic rotation velocity at Sun, to be 28.7+/-1.3 km/s/kpc using the measured distance and proper motion of ON1. This value is larger than the IAU recommended value of 25.9 km/s/kpc, but consistent with other results recently obtained with the VLBI technique.

We have performed the astrometry of H2O masers in the Galactic star-forming region Onsala 2 North (ON2N) with the VLBI Exploration of Radio Astrometry.We obtained the trigonometric parallax of 0.261+/-0.009 mas, corresponding to the heliocentric distance of 3.83+/-0.13 kpc. ON2N is expected to be on the solar circle, because its radial velocity with respect to the Local Standard of Rest (LSR) is nearly zero. Using present parallax and proper motions of the masers, the Galactocentric distance of the Sun and the Galactic rotation velocity at the Sun are found to be R_0 = 7.80+/-0.26 kpc and Theta_0 = 213+/-5 km/s, respectively. The ratio of Galactic constants, namely the angular rotation velocity of the LSR can be determined more precisely, and is found to be Omega_0=Theta_0/R_0 = 27.3+/-0.8 km/s/kpc, which is consistent with the recent estimations but different from 25.9 km/s/kpc derived from the recommended values of Theta_0 and R_0 by the International Astronomical Union (1985).

We report on the astrometric observations of the 22 GHz H2O masers in the high mass star-forming region G5.89-0.39 with VERA (VLBI Exploration of Radio Astrometry). Newly derived distance of 1.28^{+0.09}_{-0.08} kpc is the most precise and significantly nearer than previous values. We revised physical parameters and reconsidered nature of G5.89-0.39 based on the new distance as follows. (1) The ionizing star of the ultra compact (UC) HII region is a late O-type (O8 - 8.5) zero age main sequence (ZAMS) star, consistent with previously established limits based on its infrared spectral line emission. (2) Crescent-like maser alignment at the position of the O type ZAMS star may trace accretion disk (or its remnant), which suggests that the star is still young and before complete evaporation of circumstellar materials. (3) Although the revised mass for the east-west outflow has been reduced, it still quite large (100 Msun) which indicates that a significant fraction of the mass is entrained material and that the dynamical age significantly underestimates the actual outflow age. Our newly-derived distance emphasizes that G5.89-0.39 is one of the nearest targets to investigate ongoing high-mass star formation and evolution in a compact cluster containing a young O-type star.

As an extension of the project in Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 61 SU UMa-type dwarf novae mainly observed during the 2009-2010 season. The newly obtained data confirmed the basic findings reported in Kato et al. (2009): the presence of stages A-C, as well as the predominance of positive period derivatives during stage B in systems with superhump periods shorter than 0.07 d. There was a systematic difference in period derivatives for systems with superhump periods longer than 0.075 d between this study and Kato et al. (2009). We suggest that this difference is possibly caused by the relative lack of frequently outbursting SU UMa-type dwarf novae in this period regime in the present study. We recorded a strong beat phenomenon during the 2009 superoutburst of IY UMa. The close correlation between the beat period and superhump period suggests that the changing angular velocity of the apsidal motion of the elliptical disk is responsible for the variation of superhump periods. We also described three new WZ Sge-type objects with established early superhumps and one with likely early superhumps. We also suggest that two systems, VX For and EL UMa, are WZ Sge-type dwarf novae with multiple rebrightenings. The O-C variation in OT J213806.6+261957 suggests that the frequent absence of rebrightenings in very short-Porb objects can be a result of sustained superoutburst plateau at the epoch when usual SU UMa-type dwarf novae return to quiescence preceding a rebrightening. We also present a formulation for a variety of Bayesian extension to traditional period analyses.

We measured the trigonometric annual parallax of H$_2$O maser source associated with the massive star-forming regions IRAS 06061+2151 with VERA. The annual parallax of $0.496\pm0.031$ mas corresponding to a distance of $2.02^{+0.13}_{-0.12}$ kpc was obtained by 10 epochs of observations from 2007 October to 2009 September. This distance was obtained with a higher accuracy than the photometric distance previously measured, and places IRAS 06061+2151 in the Perseus spiral arm. We found that IRAS 06061+2151 also has a peculiar motion of larger than 15 km s$^{-1}$ counter to the Galactic rotation. That is similar to five sources in the Perseus spiral arm, whose parallaxes and proper motions have already been measured with higher accuracy. Moreover, these sources move at on average 27 km s$^{-1}$ toward the Galactic center and counter to the Galactic rotation.

We have carried out multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22~GHz H$_{2}$O masers associated with a Class 0 protostar L1448C in the Perseus molecular cloud. The maser features trace the base of collimated bipolar jet driven by one of the infrared counter parts of L1448C named as L1448C(N) or L1448-mm A. We detected possible evidences for apparent acceleration and precession of the jet according to the three-dimensional velocity structure. Based on the phase-referencing VLBI astrometry, we have successfully detected an annual parallax of the H$_{2}$O maser in L1448C to be 4.31$\pm$0.33~milliarcseconds (mas) which corresponds to a distance of 232$\pm$18~pc from the Sun. The present result is in good agreement with that of another H$_{2}$O maser source NGC~1333 SVS13A in the Perseus molecular cloud, 235~pc. It is also consistent with the photometric distance, 220~pc. Thus, the distance to the western part of the Perseus molecular cloud complex would be constrained to be about 235~pc rather than the larger value, 300~pc, previously reported.

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct $N$-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about $\sim 0.1$ Earth mass within $10^7$ years. This brings another difficulty in formation of gas giant planets, which requires cores with $\sim 10$ Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos' collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the solar system.

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 on trigonometric parallax measurements for the Galactic star forming region G14.33-0.64 toward the Sagittarius spiral arm. We conducted multi-epoch phase-referencing observations of an H2O maser source in G14.33-0.64 with the Japanese VLBI array VERA. We successfully detected a parallax of 0.893+/-0.101 mas, corresponding to a source distance of 1.12+/-0.13 kpc, which is less than half of the kinematic distance for G14.33-0.64. Our new distance measurement demonstrates that the Sagittarius arm lies at a closer distance of ~1 kpc, instead of previously assumed ~2-3 kpc from kinematic distances. The previously suggested deviation of the Sagittarius arm toward the Galactic center from the symmetrically fitted model (Taylor & Cordes 1993) is likely due to large errors of kinematic distances at low galactic longitudes. G14.33-0.64 most likely traces the near side of the Sagittarius arm. We attempted fitting the pitch angle of the arm with other parallax measurements along the arm, which yielded two possible pitch angles of i=34.7+/-2.7 degrees and i=11.2+/-10.5 degrees. Our proper motion measurements suggest G14.33-0.64 has no significant peculiar motion relative to the differential rotation of the Galaxy (assumed to be in a circular orbit), indicating that the source motion is in good agreement with the Galactic rotation.

It is an open question why the temperature of molecular gas in the Galactic center region is higher than that of dust. To address this problem, we made simultaneous observations in the NH_3 (J,K) = (1,1), (2,2), and (3,3) lines of the central molecular zone (CMZ) using the Kagoshima 6 m telescope. The ortho-to-para ratio of NH_3 molecules in the CMZ is 1.5--3.5 at most observed area. This ratio is higher than the statistical equilibrium value, and suggests that the formation temperature of NH_3 is 11--20 K. This temperature is similar to the dust temperature estimated from the submillimeter and infrared continuum. This result suggests that the NH_3 molecules in the CMZ were produced on dust grains with the currently observed temperature (11--20 K), and they were released into the gas phase by supernova shocks or collisions of dust particles. The discrepancy between warm molecular gas and cold dust can be explained by the transient heating of the interstellar media in the CMZ approximately 10^5 years ago when the NH_3 molecules were released from the dust.

Joint analysis of Cosmic Microwave Background, Baryon Acoustic Oscillation, and supernova data has enabled precision estimation of cosmological parameters. New programs will push to 1% uncertainty in the dark energy equation of state and tightened constraint on curvature, requiring close attention to systematics. Direct 1% measurement of the Hubble constant (H0) would provide a new constraint. It can be obtained without overlapping systematics directly from recessional velocities and geometric distance estimates for galaxies via the mapping of water maser emission that traces the accretion disks of nuclear black holes. We identify redshifts 0.02<z<0.06 as best for small samples, e.g., 10 widely distributed galaxies, each with 3% distance uncertainty. Knowledge of peculiar radial motion is also required. Mapping requires very long baseline interferometry (VLBI) with the finest angular resolution, sensitivity to individual lines of a few mJy-km/s, and baselines that can detect a complex of ~10 mJy lines (peak) in < 1 min. For 2010-2020, large ground apertures (50-100m diameter) augmenting the VLBA are critical, such as EVLA, GBT, Effelsberg, and the Large Millimeter Telescope, for which we propose a 22 GHz receiver and VLBI instrumentation. A space-VLBI aperture may be required, thus motivating US participation in the Japanese VSOP-2 mission (launch c.2013). This will provide 3-4x longer baselines and ~5x improvement in distance uncertainty. There are now 5 good targets at z>0.02, out of ~100 known masers. A single-dish discovery survey of >10,000 nuclei (>2500 hours on the GBT) would build a sample of tens of potential distance anchors. Beyond 2020, a high-frequency SKA could provide larger maser samples, enabling estimation of H0 from individually less accurate distances, and possibly without the need for peculiar motion corrections.

Oxygen-rich Asymptotic Giant Branch (AGB) stars are intense emitters of SiO and H$_2$O maser lines at 43 (J=1-0, v=1 and 2) and 22 GHz, respectively. VLBI observations of the maser emission provides a unique tool to sample the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned v=1 and v=2 SiO maser maps have traditionally prevented a unique interpretation of the observations in terms of physical underlying conditions, which depend on the nature of the SiO pumping mechanism. We have carried out observations of the SiO and H$_2$O maser emission towards RLMi, using the astrometric capabilities of VERA. Due to the too-weak emission of the reference calibrator we had to develop a special method to accurately relate the coordinates for both transitions. We present relative astrometrically aligned v=1 and v=2 J=1-0 SiO maser maps, at multiple epochs, and discuss the astrophysical results. The incorporation of astrometric information into the maps of SiO masers challenges the weak points in the current theoretical models, which will need further refinements to address the observations results.

According to our common understandings, the original surface of a short-period comet nucleus has been lost by sublimation processes during its close approaches to the Sun. Sublimation results in the formation of a dust mantle on the retreated surface and in chemical differentiation of ices over tens or hundreds of meters below the mantle. In the course of NASA's Deep Impact mission, optical and infrared imaging observations of the ejecta plume were conducted by several researchers, but their interpretations of the data came as a big surprise: (1) The nucleus of comet 9P/Tempel 1 is free of a dust mantle, but maintains its pristine crust of submicron-sized carbonaceous grains; (2) Primordial materials are accessible already at a depth of several tens of cm with abundant silicate grains of submicrometer sizes. In this study, we demonstrate that a standard model of cometary nuclei explains well available observational data: (1) A dust mantle with a thickness of ~1-2 m builds up on the surface, where compact aggregates larger than tens of micrometers dominate; (2) Large fluffy aggregates are embedded in chemically differentiated layers as well as in the deepest part of the nucleus with primordial materials. We conclude that the Deep Impact results do not need any peculiar view of a comet nucleus.

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 a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

We present the results of daily monitoring of 6.7 GHz methanol maser in Cepheus A (Cep A) using Yamaguchi 32-m radio telescope as well as the results of imaging observations conducted with the JVN (Japanese VLBI Network). We indentified five spectral features, which are grouped into red-shifted (-1.9 and -2.7 km/s) and blue-shifted (-3.8, -4.2, and -4.9 km/s), and we detected rapid variabilities of these maser features within a monitoring period of 81 days. The red-shifted features decreased in flux density to 50% of its initial value, while the flux density of the blue-shifted features rapidly increased within a 30 days. The time variation of these maser features showed two remarkable properties; synchronization and anti-correlation between the red-shifted and the blue-shifted. The spatial distribution of the maser spots obtained by the JVN observation showed an arclike structure with a scale of $\sim$1400 AU, and separations of the five maser features were found to be larger than 100 AU. The absolute position of the methanol maser was also obtained based on the phase-referencing observations, and the arclike structure were found to be associated with the Cep A-HW2 object, with the elongation of the arclike structure nearly perpendicularly to the radio continuum jet from the Cep A-HW2 object. These properties of the masers, namely, the synchronization of flux variation, and the spectral and spatial isolation of features, suggest that the collisional excitation by shock wave from a common exciting source is unlikely. Instead, the synchronized time variation of the masers can be explained if all the maser features are excited by infrared radiation from dust which is heated by a common exciting source with a rapid variability.

Comet 73P-B/Schwassmann-Wachmann 3 was observed with IRCS/Subaru at geocentric distance of 0.074 AU on UT 10 May 2006. Multiple H2O emission lines were detected in non-resonant fluorescence near 2.9 microns. No significant variation in total H2O production rate was found during the (3 hour) duration of our observations. H2O rotational temperatures and ortho-to-para abundance ratios were measured for several positions in the coma. The temperatures extracted from two different time intervals show very similar spatial distributions. For both, the rotational temperature decreased from ~110 K to ~90 K as the projected distance from the nucleus increased from ~5 to ~30 km. We see no evidence for OPR change in the coma. The H2O ortho-para ratio is consistent with the statistical equilibrium value (3.0) for all spatially resolved measurements. This implies a nuclear spin temperature higher than ~ 45 K.

We have observed 13 methanol maser sources associated with massive star-forming regions; W3(OH), Mon R2, S 255, W 33A, IRAS 18151-1208, G 24.78+0.08, G 29.95-0.02, IRAS 18556+0136, W 48, OH 43.8-0.1, ON 1, Cep A and NGC 7538 at 6.7 GHz using the Japanese VLBI Network (JVN). Twelve of the thirteen sources were detected at our longest baseline of $\sim$50 M$\lambda$, and their images are presented. Seven of them are the first VLBI images at 6.7 GHz. This high detection rate and the small fringe spacing of $\sim$4 milli-arcsecond suggest that most of the methanol maser sources have compact structure. Given this compactness as well as the known properties of long-life and small internal-motion, this methanol maser line is suitable for astrometry with VLBI.

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 report on the results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H2O masers associated with the young stellar object SVS 13 in the NGC 1333 region. We have carried out phase-referencing VLBI astrometry and measured an annual parallax of the maser features in SVS 13 of 4.25+/-0.32 mas, corresponding to the distance of 235+/-18 pc from the Sun. Our result is consistent with a photometric distance of 220 pc previously reported. Even though the maser features were detectable only for 6 months, the present results provide the distance to NGC 1333 with much higher accuracy than photometric methods. The absolute positions and proper motions have been derived, revealing that the H2O masers with the LSR (local standard of rest) velocities of 7-8 km s-1 are most likely associated with VLA4A, which is a radio counterpart of SVS 13. The origin of the observed proper motions of the maser features are currently difficult to attribute to either the jet or the rotating circumstellar disk associated with VLA4A, which should be investigated through future high-resolution astrometric observations of VLA4A and other radio sources in NGC 1333.

We have performed high-precision astrometry of H2O maser sources in Galactic star forming region Sharpless 269 (S269) with VERA. We have successfully detected a trigonometric parallax of 189+/-8 micro-arcsec, corresponding to the source distance of 5.28 +0.24/-0.22 kpc. This is the smallest parallax ever measured, and the first one detected beyond 5 kpc. The source distance as well as proper motions are used to constrain the outer rotation curve of the Galaxy, demonstrating that the difference of rotation velocities at the Sun and at S269 (which is 13.1 kpc away from the Galaxy's center) is less than 3%. This gives the strongest constraint on the flatness of the outer rotation curve and provides a direct confirmation on the existence of large amount of dark matter in the Galaxy's outer disk.

The Far-Infrared Surveyor (FIS) is one of two focal plane instruments on the AKARI satellite. FIS has four photometric bands at 65, 90, 140, and 160 um, and uses two kinds of array detectors. The FIS arrays and optics are designed to sweep the sky with high spatial resolution and redundancy. The actual scan width is more than eight arcmin, and the pixel pitch is matches the diffraction limit of the telescope. Derived point spread functions (PSFs) from observations of asteroids are similar to the optical model. Significant excesses, however, are clearly seen around tails of the PSFs, whose contributions are about 30% of the total power. All FIS functions are operating well in orbit, and its performance meets the laboratory characterizations, except for the two longer wavelength bands, which are not performing as well as characterized. Furthermore, the FIS has a spectroscopic capability using a Fourier transform spectrometer (FTS). Because the FTS takes advantage of the optics and detectors of the photometer, it can simultaneously make a spectral map. This paper summarizes the in-flight technical and operational performance of the FIS.

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from the mid- to far-infrared. The instruments also have the capability for imaging and spectroscopy in the wavelength range 2 - 180 micron in the pointed observation mode, occasionally inserted into the continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90 percent of the whole sky with higher spatial resolution and wider wavelength coverage than that of the previous IRAS all-sky survey. Point source catalogues of the All-Sky Survey will be released to the astronomical community. The pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

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 report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.

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.

This paper presents results of a survey search for bright compact radio sources at 22 GHz with the VERA radio-interferometer. Each source from a list of 2494 objects was observed in one scan for 2 minutes. The purpose of this survey was to find compact extragalactic sources bright enough at 22 GHz to be useful as phase calibrators. Observed sources were either a) within 6 degrees of the Galactic plane, or b) within 11 degrees from the Galactic center; or c) within 2 degrees from known water masers. Among the observed sources, 549 were detected, including 180 extragalactic objects which were not previously observed with the very long baseline interferometry technique. Estimates of the correlated flux densities of the detected sources are presented.

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.

We have made a high-resolution VLBI observation of the gamma-ray loud quasar PKS 1622-297 with the HALCA spacecraft and ground radio telescopes at 5 GHz in 1998 February, almost three years after the source exhibited a spectacular GeV gamma-ray flare. The source shows an elongated structure toward the west on the parsec scale. The visibility data are well modeled by three distinct components; a bright core and two weaker jet components. Comparison with previous observations confirms that the jet components have an apparent superluminal motion up to 12.1 h^{-1}c, with the inner jet components having lower superluminal speeds. We apply the inverse Compton catastrophe model and derive a Doppler factor, \delta, of 2.45, which is somewhat lower than that of other gamma-ray loud active galactic nuclei (AGNs), suggesting the source was in a more quiescent phase at the epoch of our observation. As an alternative probe of the sub-parsec scale structure, we also present the results from multi-epoch ATCA total flux monitoring, which indicate the presence of persistent intraday variability consistent with refractive interstellar scintillation. We examine the gamma-ray emission mechanism in the light of these observations.

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

We report on a strong outburst of the W49N H2O maser observed with VERA. Single-dish monitoring with VERA 20 m telescopes detected a strong outburst of the maser feature at V_LSR = -30.7 km/s in 2003 October. The outburst had a duration of ~100 days and a peak intensity of 7.9 x 10^4 Jy, being one of the strongest outbursts in W49N observed so far. VLBI observations with the VERA array were also carried out near to the maximum phase of the outburst, and the outburst spot was identified in the VLBI map. While the map was in good agreement with previous studies, showing three major concentrations of maser spots, we found a newly formed arc-like structure in the central maser concentration, which may be a shock front powered by a forming star or a star cluster. The outburst spot was found to be located on the arc-like structure, indicating a possible connection of the present outburst to a shock phenomenon.

We present the results of first dual-beam observations with VERA (VLBI Exploration of Radio Astrometry). The first dual-beam observations of a pair of H2O maser sources W49N and OH43.8-0.1 have been carried out on 2002 May 29 and July 23, and fringes of H2O maser lines at 22 GHz have been successfully detected. While the residual fringe phases of both sources show rapid variations over 360 degree due to the atmospheric fluctuation, the differential phase between the two sources remains constant for 1 hour with r.m.s. of 8 degree, demonstrating that the atmospheric phase fluctuation is removed effectively by the dual-beam phase referencing. The analysis based on Allan standard deviation reveals that the differential phase is mostly dominated by white phase noise, and the coherence function calculated from the differential phase shows that after phase referencing the fringe visibility can be integrated for arbitrarily long time. These results demonstrate VERA's high capability of phase referencing, indicating that VERA is a promising tool for phase referencing VLBI astrometry at 10 microarcsec level accuracy.

We investigate the effects of a passing stellar encounter on a planetesimal disk through analytical calculations and numerical simulations, and derive the boundary radius ($a_{\rm planet}$) outside which planet formation is inhibited by disruptive collisions with high relative velocities.

We describe measurements of the size of the Vela pulsar via scintillation, using both fits to the distribution of intensity and measurements of the modulation index. We briefly discuss systematic effects other than source size that can affect the distribution, including gain variations, self-noise, scintillation shot noise, and correlator saturation. Modulation index, a single number, can be biased by all of these, whereas the distribution of intensity is affected in different ways by different effects, providing means of distinguishing among them. Self-noise and gain variations are likely more important at long observing wavelengths, and correlator saturation and scintillation shot noise at short wavelengths. We find a size of about 500 km at decimeter wavelengths. Interestingly, this agrees with measurements of modulation index by Roberts & Ables at the same wavelength. Their results (and more recently that reported by Macquart et al.) suggest that size decreases with increasing wavelength. Although consistent with the observations, this conclusion is perhaps surprising from the standpoint of the traditional radius-to-frequency mapping. However, these measurements are of size rather than height; and of course systematic effects may play a role.