We investigate properties of the interstellar medium (ISM) interacting with shocks around the Galactic supernova remnant IC443, using the results of near-infrared [FeII] and H2 line mapping with the IRSF/SIRIUS. In the present study, we newly performed H2 1-0 S(1) and 2-1 S(1) line mapping with the narrow-band filters tuned for these lines, covering the entire remnant (30'x35'). Combined with an [FeII] line map, our result shows that the H2 line emission is significantly detected in the southern region, while the [FeII] line emission is detected all over the remnant, suggesting that slow and fast shocks propagate through the southern region and the entire remnant, respectively. In particular, the H2 line emission is relatively strong compared to the [FeII] line emission in the southwestern region, where TeV gamma-ray emission is detected. As the strong H2 line emission indicates the dominance of the dense ISM, this result supports the scenario that the gamma-ray emission is likely to be produced through a heavy interaction between cosmic-ray protons and the dense ISM in the southwestern shell. We also find that the H2 and [FeII] line emissions show an anti-correlated spatial distribution in the same region, suggesting the presence of the clumpy ISM. Such a clumpy morphology of the ISM around IC443 may assist cosmic-ray protons to efficiently interact with large amounts of the ISM protons.

We report results of astrometric VLBI observations toward a Mira variable star BX~Cam using the VLBI array "VERA". The observations were performed from February 2012 to November 2014. Obtained parallax is 1.73$\pm$0.03 mas corresponding to a distance of 0.58$\pm$0.01 kpc. Parallax of this source was also reported in Gaia DR2 as 4.13$\pm$0.25 mas, and there is a 240 \% difference between these two measurements. Astrometric results from our VLBI observations show that we exactly traced angular motions of the seven maser spots in BX~Cam. We calculated stellar luminosities using both parallaxes, and obtained luminosities of $L_{\ast}^{\mathrm{VERA}} = 4950\pm170 L_{\odot}$ and $L_{\ast}^{\mathrm{Gaia}} = 870\pm110 L_{\odot}$. Deduced luminosities also support a validity of the parallax that we determined with VERA. Evaluating the two parallaxes, we concluded that the parallax of 1.73$\pm$0.03 mas from the VERA observations is correct for BX~Cam. We obtained a systemic motion of BX~Cam as (${\mu_{\alpha}}\cos{\delta}^{\mathrm{sys}}$, ${\mu}_{\delta}^{\mathrm{sys}}$) $=$ (13.48$\pm$0.14, $-$34.30$\pm$0.18) mas\,yr$^{-1}$. A total of 73 H$_2$O maser spots detected from our VLBI observations show a spatial distribution of 30 au $\times$ 80 au with a strong elongation along north-south direction. They show outflows with a three-dimensional velocity of 14.79$\pm$1.40 km\,s$^{-1}$. From a comparison between time variations of $V$-band magnitudes and H$_2$O maser, we found that variation of the H$_2$O maser is relevant to that seen in $V$-band even though the H$_2$O maser does not recover its maximum flux in each cycle.

The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $\theta_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $\mu$G, indicating that the BHR 71 core is magnetically supercritical with $\lambda = 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core.

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.

In this study, the detailed magnetic field structure of the dense protostellar core Barnard 335 (B335) was revealed based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. Magnetic fields pervading B335 were mapped using 24 stars after subtracting unrelated ambient polarization components, for the first time revealing that they have an axisymmetrically distorted hourglass-shaped structure toward the protostellar core. On the basis of simple two- and three-dimensional magnetic field modeling, magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $90^{\circ} \pm 7^{\circ}$ and $50^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B335 was determined to be $30.2 \pm 17.7$ $\mu {\rm G}$. The critical mass of B335, evaluated using both magnetic and thermal/turbulent support against collapse, was determined to be $M_{\rm cr} = 3.37 \pm 0.94$ ${\rm M}_{\odot}$, which is identical to the observed core mass of $M_{\rm core}=3.67$ M$_{\odot}$. We thus concluded that B335 started its contraction from a condition near equilibrium. We found a linear relationship in the polarization versus extinction diagram, up to $A_V \sim 15$ mag toward the stars with the greatest obscuration, which verified that our observations and analysis provide an accurate depiction of the core.

The detailed magnetic field structure of the starless dense core CB81 (L1774, Pipe 42) in the Pipe Nebula was determined based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. The magnetic fields pervading CB81 were mapped using 147 stars and axisymmetrically distorted hourglass-like fields were identified. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $4^{\circ} \pm 8^{\circ}$ and $20^{\circ} \pm 20^{\circ}$, respectively. The total magnetic field strength of CB81 was found to be $7.2 \pm 2.3$ $\mu{\rm G}$. Taking into account the effects of thermal/turbulent pressure and magnetic fields, the critical mass of CB81 was calculated to be $M_{\rm cr}=4.03 \pm 0.40$ M$_{\odot}$, which is close to the observed core mass of $M_{\rm core}=3.37 \pm 0.51$ M$_{\odot}$. We thus conclude that CB81 is in a condition close to the critical state. In addition, a spatial offset of $92''$ was found between the center of magnetic field geometry and the dust extinction distribution; this offset structure could not have been produced by self-gravity. The data also indicate a linear relationship between polarization and extinction up to $A_V \sim 30$ mag going toward the core center. This result confirms that near-infrared polarization can accurately trace the overall magnetic field structure of the core.

Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for the flux freezing model were a line-of-sight magnetic inclination angle of $\gamma_{\rm mag} = 35^{\circ} \pm 15^{\circ}$ and a core center to ambient (background) density contrast of $\rho_{\rm c} / \rho_{\rm bkg} = 75$. The initial density for core formation ($\rho_0$) was estimated to be $\rho_{\rm c} / 75 = 4670$ cm$^{-3}$, which is about one order of magnitude higher than the expected density ($\sim 300$ cm$^{-3}$) for the inter-clump medium of the Pipe Nebula. FeSt 1-457 is likely to have been formed from the accumulation of relatively dense gas, and the relatively dense background column density of $A_V \simeq 5$ mag supports this scenario. The initial radius (core formation radius) $R_0$ and the initial magnetic field strength $B_0$ were obtained to be 0.15 pc ($1.64 R$) and $10.8-14.6$ $\mu$G, respectively. We found that the initial density $\rho_0$ is consistent with the mean density of the nearly critical magnetized filament with magnetic field strength $B_0$ and radius $R_0$. The relatively dense initial condition for core formation can be naturally understood if the origin of the core is the fragmentation of magnetized filaments.

On 15 February 2013, the asteroid 367943 Duende (2012 DA14) experienced a near-Earth encounter at an altitude of 27,700 km or 4.2 Earth radii. We present here the results of an extensive, multi-observatory campaign designed to probe for spectral and/or rotational changes to Duende due to gravitational interactions with the Earth during the flyby. Our spectral data reveal no changes within the systematic uncertainties of the data. Post-flyby lightcurve photometry places strong constraints on the rotation state of Duende, showing that it is in non-principal axis rotation with fundamental periods of P_1 = 8.71 +/- 0.03 and P_2 = 23.7 +/- 0.2 hours. Multiple lightcurve analysis techniques, coupled with theoretical considerations and delay-doppler radar imaging, allows us to assign these periods to specific rotational axes of the body. In particular we suggest that Duende is now in a non-principal, short axis mode rotation state with a precessional period equal to P_1 and oscillation about the symmetry axis at a rate equal to P_2. Temporal and signal-to-noise limitations inherent to the pre-flyby photometric dataset make it difficult to definitively diagnose whether these periods represent a change imparted due to gravitational torques during the flyby. However, based on multiple analysis techniques and a number of plausibility arguments, we suggest that Duende experienced a rotational change during the planetary encounter with an increase in its precessional rotation period. Our preferred interpretation of the available data is that the precession rate increased from 8.4 hours prior to the flyby to 8.7 hours afterwards. A companion paper by Benson et al. (2019) provides a more detailed dynamical analysis of this event and compares the data to synthetic lightcurves computed from a simple shape model of Duende. (abbreviated abstract)

The magnetic field structure, kinematical stability, and evolutionary status of the starless dense core Barnard 68 (B68) are revealed based on the near-infrared polarimetric observations of background stars, measuring the dichroically polarized light produced by aligned dust grains in the core. After subtracting unrelated ambient polarization components, the magnetic fields pervading B68 are mapped using 38 stars and axisymmetrically distorted hourglass-like magnetic fields are obtained, although the evidence for the hourglass field is not very strong. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles on the plane-of-sky and in the line-of-sight direction are determined to be $47^{\circ} \pm 5^{\circ}$ and $20^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B68 is obtained to be $26.1 \pm 8.7$ $\mu {\rm G}$. The critical mass of B68, evaluated using both magnetic and thermal/turbulent support, is $M_{\rm cr} = 2.30 \pm 0.20$ ${\rm M}_{\odot}$, which is consistent with the observed core mass of $M_{\rm core}=2.1$ M$_{\odot}$, suggesting nearly critical state. We found a relatively linear relationship between polarization and extinction up to $A_V \sim 30$ mag toward the stars with deepest obscuration. Further theoretical and observational studies are required to explain the dust alignment in cold and dense regions in the core.

We report measurements of parallax and proper motion for four 22 GHz water maser sources as part of VERA Outer Rotation Curve project. All sources show Galactic latitudes of $>$ 2$^{\circ}$ and Galactocentric distances of $>$ 11 kpc at the Galactic longitude range of 95$^{\circ}$ $< l <$ 126$^{\circ}$. The sources trace the Galactic warp reaching to 200$\sim$400 pc, and indicate the signature of the warp to 600 pc toward the north Galactic pole. The new results along with previous results in the literature show the maximum height of the Galactic warp is increased with Galactocentric distance. Also, we examined velocities perpendicular to the disk for the sample, and found an oscillatory behavior between the vertical velocities and Galactic heights. This behavior suggests the existence of the bending (vertical density) waves, possibly induced by a perturbing satellite (e.g. passage of the Sagittarius dwarf galaxy).

As part our investigation into the Galactic rotation curve, we carried out Very Long Baseline Interferometry (VLBI) observations towards the star-forming region IRAS 01123+6430 using VLBI Exploration of Radio Astrometry (VERA) to measure its annual parallax and proper motion. The annual parallax was measured to be 0.151+/-0.042 mas, which corresponds to a distance of D=6.61^{+2.55}_{-1.44} kpc, and the obtained proper motion components were ({\mu}_{\alpha}cos{\delta}, {\mu}_{\delta})=(-1.44+/-0.15, -0.27+/-0.16) mas yr^{-1} in equatorial coordinates. Assuming Galactic constants of (R_0, {\Theta}_0)=(8.05+/-0.45 kpc, 238+/-14 km s^{-1}), the Galactocentric distance and rotation velocity were measured to be (R, {\Theta})=(13.04+/-2.24 kpc, 239+/-22 km s^{-1}), which are consistent with a flat Galactic rotation curve. The newly estimated distance provides a more accurate bolometric luminosity of the central young stellar object, L_Bol=(3.11+/-2.86)\times 10^3 L_solar, which corresponds to a spectral type of B1-B2. The analysis of 12CO(J=1-0) survey data obtained with the Five College Radio Astronomical Observatory (FCRAO) 14 m telescope shows that the molecular cloud associated with IRAS 01123+6430 consists of arc-like and linear components, which well matches a structure predicted by numerical simulation of the cloud-cloud collision (CCC) phenomenon. The coexistence of arc-like and linear components implies that the relative velocity of initial two clouds was as slow as 3-5 km s^{-1}, which meets the expected criteria of massive star formation where the core mass is effectively increased in the presence of low relative velocity (~3-5 km s^{-1}), as suggested by Takahira et al.(2014).

Star formation is thought to be driven by two groups of mechanisms; spontaneous collapse and triggered collapse. Triggered star formation mechanisms further diverge into cloud-cloud collision (CCC), "collect and collapse" (C\&C) and shock induced collapse of pre-existing, gravitationally stable cores, or 'radiation driven implosion' (RDI). To evaluate the contributions of these mechanisms and establish whether these processes can occur together within the same star forming region we performed mapping observations of radio frequency ammonia, and water maser emission lines in the S235 massive star forming region. Via spectral analyses of main, hyperfine and multi-transitional ammonia lines we explored the distribution of temperature and column density in the dense gas in the S235 and S235AB star forming region. The most remarkable result of the mapping observations is the discovery of high density gas in inter-core bridges which physically link dense molecular cores that house young proto-stellar clusters. The presence of dense gas implies the potential for future star formation within the system of cores and gas bridges. Cluster formation implies collapse and the continuous physical links, also seen in re-imaged archival CS and $^{13}$CO maps, suggests a common origin to the molecular cores housing these clusters, i.e the structure condensed from a single, larger parent cloud, brought about by the influence of a local expanding H${\rm II}$ region. An ammonia absorption feature co-locating with the center of the extended H${\rm II}$ region may be attributed to an older gas component left over from the period prior to formation of the H${\rm II}$ region. Our observations also detail known and new sites of water maser emission, highlighting regions of active ongoing star formation.

Aims. We aim to constrain the size and porosity of ejected dust particles from comet 252P/LINEAR and their evolution near the perihelion via near-infrared multiband polarimetry. A close approach of the comet to the Earth in March 2016 (~0.036 au) provided a rare opportunity for the sampling of the comet with a high spatial resolution. Methods. We made NIR JHKS bands polarimetric observations of the comet for 12 days near perihelion, interspersed between broadband optical imaging observations over four months. In addition, dynamical simulation of the comet was performed 1000 yr backward in time. Results. We detected two discontinuous brightness enhancements. Before the first enhancement, the NIR polarization degrees were far lower than those of ordinary comets at a given phase angle. Soon after the activation, however, they increased by ~13 % at most, showing unusual blue polarimetric color over the J and H bands (-2.55 % / um on average) and bluing of both J-H and H-Ks dust color. Throughout the event, the polarization vector was marginally aligned perpendicular to the scattering plane. The subsequent postperihelion reactivation of the comet lasted for approximately 1.5 months, with a factor of ~30 times pre-activation dust mass-loss rates in the Rc band. Conclusions. The marked increase in the polarization degree with blue NIR polarimetric color is reminiscent of the behaviors of a fragmenting comet D/1999 S4 (LINEAR). The most plausible scenario for the observed polarimetric properties of 252P/LINEAR would be an ejection of predominantly large, compact dust particles from the desiccated surface layer. We conjecture that the more intense solar heating that the comet has received in the near-Earth orbit would cause the paucity of small, fluffy dust particles around the nucleus of the comet.

We conducted near-infrared (JHKs) imaging polarimetry toward the infrared dark cloud (IRDC) M17 SWex, including almost all of the IRDC filaments as well as its outskirts, with the polarimeter SIRPOL on the IRSF 1.4 m telescope. We revealed the magnetic fields of M17 SWex with our polarization-detected sources that were selected by some criteria based on their near-IR colors and the column densities toward them, which were derived from the Herschel data. The selected sources indicate not only that the ordered magnetic field is perpendicular to the cloud elongation as a whole, but also that at both ends of the elongated cloud the magnetic field appears to bent toward its central part, i.e., large-scale hourglass-shaped magnetic field perpendicular to the cloud elongation. In addition to this general trend, the elongations of the filamentary subregions within the dense parts of the cloud appear to be mostly perpendicular to their local magnetic fields, while the magnetic fields of the outskirts appear to follow the thin filaments that protrude from the dense parts. The magnetic strengths were estimated to be ~70-300 microG in the subregions, of which lengths and average number densities are ~3-9 pc and ~2-7x10^3 cm^{-3}, respectively, by the Davis-Chandrasekhar-Fermi method with the angular dispersion of our polarization data and the velocity dispersion derived from the C^{18}O (J=1-0) data obtained by the Nobeyama 45 m telescope. These field configurations and our magnetic stability analysis of the subregions imply that the magnetic field have controlled the formation/evolution of the M17 SWex cloud.

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

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

We report X-ray, optical, and near-infrared monitoring of the new X-ray transient MAXI J1820$+$070 discovered with MAXI on 2018 March 11. Its X-ray intensity reached $\sim 2$ Crab in 2--20 keV at the end of March, and then gradually decreased until the middle of June. In this period, the X-ray spectrum was described by Comptonization of the disk emission, with a photon index of $\sim$1.5 and an electron temperature of $\sim$50 keV, which is consistent with a black hole X-ray binary in the low/hard state. The electron temperature and the photon index were slightly decreased and increased with increasing flux, respectively. The source showed significant X-ray flux variation on timescales of seconds. This short-term variation was found to be associated with changes in the spectral shape, and the photon index became slightly harder at higher fluxes. This suggests that the variation was produced by a change in the properties of the hot electron cloud responsible for the strong Comptonization. Modeling a multi-wavelength SED around the X-ray flux peak at the end of March, covering the near-infrared to X-ray bands, we found that the optical and near-infrared fluxes were likely contributed substantially by the jet emission. Before this outburst, the source was never detected in the X-ray band with MAXI (with a 3$\sigma$ upper limit of $\sim$0.2 mCrab in 4--10 keV, obtained from the 7-year data in 2009--2016), whereas weak optical and infrared activity was found at their flux levels $\sim$3 orders of magnitude lower than the peak fluxes in the outburst.

The relationship between submillimeter (submm) dust emission polarization and near-infrared (NIR) $H$-band polarization produced by dust dichroic extinction was studied for the cold starless dense core FeSt 1-457. The distribution of polarization angles ($90^{\circ}$-rotated for submm) and degrees were found to be very different between at submm and NIR wavelengths. The mean polarization angles for FeSt 1-457 at submm and NIR wavelengths are $132.1^{\circ} \pm 22.0^{\circ}$ and $2.7^{\circ} \pm 16.2^{\circ}$, respectively. The correlation between $P_H$ and $A_V$ was found to be linear from outermost regions to relatively dense line of sight of $A_V \approx 25$ mag, indicating that NIR polarization reflects overall polarization (magnetic field) structure of the core at least in this density range. The flat $P_H/A_V$ versus $A_V$ correlations were confirmed, and the polarization efficiency was found to be comparable to the observational upper limit (Jones 1989). On the other hand, as reported by Alves et al., submm polarization degrees show clear linearly decreasing trend against $A_V$ from $A_V \approx 20$ mag to the densest center ($A_V \approx 41$ mag), appearing as "polarization hole" structure. The power law index for the $P_{\rm submm}$ versus $A_V$ relationship was obtained to be $\approx -1$, indicating that the alignment for the submm sensitive dust is lost. These very different polarization distributions at submm and NIR wavelengths suggest that (1) there is different radiation environment at these wavelengths or (2) submm-sensitive dust is localized or the combination of them.

Spectropolarimetry is a powerful technique for investigating the physical properties of gas and solid materials in cometary comae without mutual contamination, but there have been few spectropolarimetric studies to extract each component. We attempt to derive the continuum polarization degree of comet 2P/Encke, free from influence of molecular emissions. The target is unique in that it has an orbit dynamically decoupled from Jupiter like main-belt asteroids, while ejecting gas and dust like ordinary comets. We observed the comet using the Higashi-Hiroshima Optical and Near-Infrared Camera attached to the Cassegrain focus of the 150-cm Kanata telescope on UT 2017 February 21 when the comet was at the solar phase angle of 75.7 deg. We find that the continuum polarization degree with respect to the scattering plane is 33.8+/-2.7 % at the effective wavelength of 0.815 um, which is significantly higher than those of cometary dust in a high-Pmax group at similar phase angles. Assuming that an ensemble polarimetric response of 2P/Encke's dust as a function of phase angle is morphologically similar with those of other comets, its maximum polarization degree is estimated to > 40 % at the phase angle of ~100 deg. In addition, we obtain the polarization degrees of the C2 swan bands (0.51-0.56 um), the NH2 alpha bands (0.62-0.69 um) and the CN-red system (0.78-0.94 um) in a range of 3-19 %, which depend on the molecular species and rotational quantum numbers of each branch. The polarization vector aligns nearly perpendicularly to the scattering plane with the average of 0.4 deg over a wavelength range of 0.50-0.97 um. From the observational evidence, we conjecture that the large polarization degree of 2P/Encke would be attributable to a dominance of large dust particles around the nucleus, which have remained after frequent perihelion passages near the Sun.

A very long term near-infrared variable star survey towards the Large and Small Magellanic Clouds was carried out using the 1.4m InfraRed Survey Facility at the South African Astronomical Observatory. This project was initiated in December 2000 in the LMC, and in July 2001 in the SMC. Since then an area of 3 square degrees along the bar in the LMC and an area of 1 square degree in the central part of the SMC have been repeatedly observed. This survey is ongoing, but results obtained with data taken until December 2017 are reported in this paper. Over more than 15 years we have observed the two survey areas more than one hundred times. This is the first survey that provides near-infrared time-series data with such a long time baseline and on such a large scale. This paper describes the observations in the SMC and publishes a point source photometric catalogue, a variable source catalogue, and time-series data.

In the present study, the magnetic field scaling on density, $|B| \propto \rho^{\kappa}$, was revealed in a single starless core for the first time. The $\kappa$ index of $0.78 \pm 0.10$ was obtained toward the starless dense core FeSt 1-457 based on the analysis of the radial distribution of the polarization angle dispersion of background stars measured at the near-infrared wavelengths. The result prefers $\kappa = 2/3$ for the case of isotropic contraction, and the difference of the observed value from $\kappa = 1/2$ is 2.8 sigma. The distribution of the ratio of mass to magnetic flux was evaluated. FeSt 1-457 was found to be magnetically supercritical near the center ($\lambda \approx 2$), whereas nearly critical or slightly subcritical at the core boundary ($\lambda \approx 0.98$). Ambipolar-diffusion-regulated star formation models for the case of moderate magnetic field strength may explain the physical status of FeSt 1-457. The mass-to-flux ratio distribution for typical dense cores (critical Bonnor--Ebert sphere with central $\lambda=2$ and $\kappa=1/2$--$2/3$) was calculated and found to be magnetically critical/subcritical at the core edge, which indicates that typical dense cores are embedded in and evolve from magnetically critical/subcritical diffuse surrounding medium.

We present optical and near-infrared observations of a low-luminosity Type IIP supernova (SN) 2016bkv from the initial rising phase to the plateau phase. Our observations show that the end of the plateau is extended to $\gtrsim 140$ days since the explosion, indicating that this SN takes one of the longest time to finish the plateau phase. among Type IIP SNe (SNe IIP), including low-luminosity (LL) SNe IIP. The line velocities of various ions at the middle of the plateau phase are as low as 1,000--1,500 km s$^{-1}$, which is the lowest even among LL SNe IIP. These measurements imply that the ejecta mass in SN 2016bkv is larger than that of the well-studied LL IIP SN 2003Z. In the early phase, SN 2016bkv shows a strong bump in the light curve. In addition, the optical spectra in this bump phase exhibit a blue continuum accompanied with a narrow H$\alpha$ emission line. These features indicate an interaction between the SN ejecta and the circumstellar matter (CSM) as in SNe IIn. Assuming the ejecta-CSM interaction scenario, the mass loss rate is estimated to be $\sim 1.7 \times 10^{-2} M_{\odot}$ yr$^{-1}$ within a few years before the SN explosion. This is comparable to or even larger than the largest mass loss rate observed for the Galactic red supergiants ($\sim 10^{-3} M_{\odot}$ yr$^{-1}$ for VY CMa). We suggest that the progenitor star of SN 2016bkv experienced a violent mass loss just before the SN explosion.

The relationship between dust polarization and extinction was determined for the cold dense starless molecular cloud core FeSt 1-457 based on the background star polarimetry of dichroic extinction at near-infrared wavelengths. Owing to the known (three-dimensional) magnetic field structure, the observed polarizations from the core were corrected by considering (a) the subtraction of the ambient polarization component, (b) the depolarization effect of inclined distorted magnetic fields, and (c) the magnetic inclination angle of the core. After these corrections, a linear relationship between polarization and extinction was obtained for the core in the range up to $A_V \approx 20$ mag. The initial polarization vs. extinction diagram changed dramatically after the corrections of (a) to (c), with the correlation coefficient being refined from 0.71 to 0.79. These corrections should affect the theoretical interpretation of the observational data. The slope of the finally obtained polarization--extinction relationship is $P_H / E_{H-K_s} = 11.00 \pm 0.72$ $\%$ ${\rm mag}^{-1}$, which is close to the statistically estimated upper limit of the interstellar polarization efficiency (Jones 1989). This consistency suggests that the upper limit of interstellar polarization efficiency might be determined by the observational viewing angle toward polarized astronomical objects.

We investigate the optical properties of Ophiuchus to obtain clues on the formation epoch of this cluster, and compare them to those of the Coma cluster, which is comparable in mass to Ophiuchus but much more disturbed dynamically. Based on a deep image of the Ophiuchus cluster in the r' band obtained at the Canada France Hawaii Telescope with the MegaCam camera, we have applied an iterative process to subtract the contribution of the numerous stars that pollute the image, due to the low Galactic latitude of the cluster, and obtained a photometric catalogue of 2818 galaxies fully complete at r'=20.5 mag and still 91% complete at r'=21.5 mag. We use this catalogue to derive the cluster Galaxy Luminosity Function (GLF) for the overall image and for a region (hereafter the "rectangle" region) covering exactly the same physical size as the region in which the GLF of the Coma cluster was studied by Adami et al. (2007). We then compute density maps based on an adaptive kernel technique, for different magnitude limits, and define three circular regions covering 0.08, 0.08 and 0.06 deg^2 respectively centered on the cluster (C), northwest (NW) and southeast (SE) of the cluster, in which we compute the GLFs. The GLF fits are much better when a Gaussian is added to the usual Schechter function, to account for the excess of very bright galaxies. Compared to Coma, Ophiuchus shows a strong excess of bright galaxies. The properties of the two nearby very massive clusters Ophiuchus and Coma are quite comparable, though they seem embedded in different large scale environments. Our interpretation is that Ophiuchus has built up long ago, as confirmed by its relaxed state (see paper I) while Coma is still in the process of forming.

ThunderKAT is the image-plane transients programme for MeerKAT. The goal as outlined in 2010, and still today, is to find, identify and understand high-energy astrophysical processes via their radio emission (often in concert with observations at other wavelengths). Through a comprehensive and complementary programme of surveying and monitoring Galactic synchrotron transients (across a range of compact accretors and a range of other explosive phenomena) and exploring distinct populations of extragalactic synchrotron transients (microquasars, supernovae and possibly yet unknown transient phenomena) - both from direct surveys and commensal observations - we will revolutionise our understanding of the dynamic and explosive transient radio sky. As well as performing targeted programmes of our own, we have made agreements with the other MeerKAT large survey projects (LSPs) that we will also search their data for transients. This commensal use of the other surveys, which remains one of our key programme goals in 2016, means that the combined MeerKAT LSPs will produce by far the largest GHz-frequency radio transient programme to date.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

Recent detection of gravitational waves from a neutron star (NS) merger event GW170817 and identification of an electromagnetic counterpart provide a unique opportunity to study the physical processes in NS mergers. To derive properties of ejected material from the NS merger, we perform radiative transfer simulations of kilonova, optical and near-infrared emissions powered by radioactive decays of r-process nuclei synthesized in the merger. We find that the observed near-infrared emission lasting for > 10 days is explained by 0.03 Msun of ejecta containing lanthanide elements. However, the blue optical component observed at the initial phases requires an ejecta component with a relatively high electron fraction (Ye). We show that both optical and near-infrared emissions are simultaneously reproduced by the ejecta with a medium Ye of ~ 0.25. We suggest that a dominant component powering the emission is post-merger ejecta, which exhibits that mass ejection after the first dynamical ejection is quite efficient. Our results indicate that NS mergers synthesize a wide range of r-process elements and strengthen the hypothesis that NS mergers are the origin of r-process elements in the Universe.

Three dimensional (3D) magnetic field information on molecular clouds and cores is important for revealing their kinematical stability (magnetic support) against gravity which is fundamental for studying the initial conditions of star formation. In the present study, the 3D magnetic field structure of the dense starless core FeSt 1-457 is determined based on the near-infrared polarimetric observations of the dichroic polarization of background stars and simple 3D modeling. With an obtained angle of line-of-sight magnetic inclination axis $\theta_{\rm inc}$ of $45^{\circ}\pm10^{\circ}$ and previously determined plane-of-sky magnetic field strength $B_{\rm pol}$ of $23.8\pm12.1$ $\mu{\rm G}$, the total magnetic field strength for FeSt 1-457 is derived to be $33.7\pm18.0$ $\mu{\rm G}$. The critical mass of FeSt 1-457, evaluated using both magnetic and thermal/turbulent support is ${M}_{\rm cr} = 3.70\pm0.92$ ${\rm M}_{\odot}$, which is identical to the observed core mass, $M_{\rm core}=3.55\pm0.75$ ${\rm M}_{\odot}$. We thus conclude that the stability of FeSt 1-457 is in a condition close to the critical state. Without infalling gas motion and no associated young stars, the core is regarded to be in the earliest stage of star formation, i.e., the stage just before the onset of dynamical collapse following the attainment of a supercritical condition. These properties would make FeSt 1-457 one of the best starless cores for future studies of the initial conditions of star formation.

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations.

Magnetic fields are believed to play an important role in controlling the stability and contraction of dense condensations of gas and dust leading to the formation of stars and planetary systems. In the present study, the magnetic field of FeSt 1-457, a cold starless molecular cloud core, was mapped on the basis of the polarized near-infrared light from 185 background stars after being dichroically absorbed by dust aligned with the magnetic field in the core. A distinct "hourglass-shaped" magnetic field was identified in the region of the core, which was interpreted as the first evidence of a magnetic field structure distorted by mass condensation in a starless core. The steep curvature of the magnetic field lines obtained in the present study indicates that the distortion was mainly created during the formation phase of the dense core. The derived mass-to-magnetic flux ratio indicates that the core is in a magnetically supercritical state. However, the stability of the core can be considered to be in a nearly critical state if the additional contributions from the thermal and turbulent support are included. Further diffusion of the magnetic field and/or turbulent dissipation would cause the onset of dynamical collapse of the core. The geometrical relationship between the direction of the magnetic field lines and the elongation of the core was found to be in good agreement with the theoretical predictions for the formation of Sun-like stars under the influence of a magnetic field.

From measuring the annual parallax of water masers over one and a half years with VERA, we present the trigonometric parallax and corresponding distance of another newly identified water maser source in the region of IRAS 20231+3440 as $\pi=0.611\pm0.022$ mas and $D=1.64\pm0.06$ kpc respectively. We measured the absolute proper motions of all the newly detected maser spots (30 spots) and presented two pictures describing the possible spatial distribution of the water maser as the morphology marks out an arc of masers whose average proper motion velocity in the jet direction was 14.26 km s$^{-1}$. As revealed by the ALLWISE composite image, and by applying the colour-colour method of YSO identification and classification on photometric archived data, we identified the driving source of the north maser group to be a class I, young stellar object. To further probe the nature of the progenitor, we used the momentum rate maximum value (1.2$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$ km s$^{-1}$) of the outflow to satisfy that the progenitor under investigation is a low mass young stellar object concurrently forming alongside an intermediate-mass YSO $\sim 60,000$ au ($\sim 37$ arcsecs) away from it.

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths.

We present the results of multi-epoch VLBI observations of water masers in the AGFL 5142 massive star forming region. We measure an annual parallax of $\pi=0.467 \pm 0.010$ mas, corresponding to a source distance of $D=2.14^{+0.051}_{-0.049}$ kpc. Proper motion and line of sight velocities reveal the 3D kinematics of masers in this region, most of which associate with millimeter sources from the literature. In particular we find remarkable bipolar bowshocks expanding from the most massive member, AFGL 5142 MM1, which are used to investigate the physical properties of its protostellar jet. We attempt to link the known outflows in this region to possible progenitors by considering a precessing jet scenario and we discuss the episodic nature of ejections in AFGL 5142.

We investigated the star formation activities in the AFGL333 region, which is in the vicinity of the W4 expanding bubble, by conducting NH3 (1,1), (2,2), and (3,3) mapping observations with the 45 m Nobeyama Radio Telescope at an angular resolution of 75". The morphology of the NH3 (1,1) map shows a bow-shape structure with the size of 2.0 x 0.6 pc as seen in the dust continuum. At the interface between the W4 bubble and the dense NH3 cloud, the compact HII region G134.2+0.8, associated with IRAS02245+6115, is located. Interestingly, just north and south of G134.2+0.8 we found NH3 emission exhibiting large velocity widths of ~ 2.8 km/s, compared to 1.8 km/s at the other positions. As the possibility of mechanical energy injection through the activity of YSO(s) is low, we considered the origin of the large turbulent gas motion as indication of interaction between the compact HII region and the periphery of the dense molecular cloud. We also found expanding motion of the CO emission associated with G134.2+0.8. The overall structure of the AFGL333-Ridge might have been formed by the expanding bubble of W4. However, the small velocity widths observed west of IRAS02245+6115, around the center of the dense molecular cloud, suggest that interaction with the compact HII region is limited. Therefore the YSOs (dominantly Class 0/I) in the core of the AFGL333-Ridge dense molecular cloud most likely formed in quiescent mode. As has been previously suggested for the large scale star formation in the W3 giant molecular cloud, our results show an apparent coexistence of induced and quiescent star formation in this region. It appears that star formation in the AFGL333 region has proceeded without significant external triggers, but accompanying stellar feedback environment.

We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226.

We present $Spitzer$/IRS mid-infrared spectral maps of the Galactic star-forming region M17 as well as IRSF/SIRIUS Br$\gamma$ and Nobeyama 45-m/FOREST $^{13}$CO ($J$=1--0) maps. The spectra show prominent features due to polycyclic aromatic hydrocarbons (PAHs) at wavelengths of 6.2, 7.7, 8.6, 11.3, 12.0, 12.7, 13.5, and 14.2 $\mu$m. We find that the PAH emission features are bright in the region between the HII region traced by Br$\gamma$ and the molecular cloud traced by $^{13}$CO, supporting that the PAH emission originates mostly from photo-dissociation regions. Based on the spatially-resolved $Spitzer$/IRS maps, we examine spatial variations of the PAH properties in detail. As a result, we find that the interband ratio of PAH 7.7 $\mu$m/PAH 11.3 $\mu$m varies locally near M17SW, but rather independently of the distance from the OB stars in M17, suggesting that the degree of PAH ionization is mainly controlled by local conditions rather than the global UV environments determined by the OB stars in M17. We also find that the interband ratios of the PAH 12.0 $\mu$m, 12.7 $\mu$m, 13.5 $\mu$m, and 14.2 $\mu$m features to the PAH 11.3 $\mu$m feature are high near the M17 center, which suggests structural changes of PAHs through processing due to intense UV radiation, producing abundant edgy irregular PAHs near the M17 center.

Accurately measuring the masses of white dwarf stars is crucial in many astrophysical contexts (e.g., asteroseismology and cosmochronology). These masses are most commonly determined by fitting a model atmosphere to an observed spectrum; this is known as the spectroscopic method. However, for cases in which more than one method may be employed, there are well known discrepancies between masses determined by the spectroscopic method and those determined by astrometric, dynamical, and/or gravitational-redshift methods. In an effort to resolve these discrepancies, we are developing a new model of hydrogen in a dense plasma that is a significant departure from previous models. Experiments at Sandia National Laboratories are currently underway to validate these new models, and we have begun modifications to incorporate these models into stellar-atmosphere codes.

As part of our laboratory investigation of the theoretical line profiles used in white dwarf atmosphere models, we extend the electron-density ($n_{\rm e}$) range measured by our experiments to higher densities (up to $n_{e}\sim80\times10^{16}$ cm$^{-3}$). Whereas inferred parameters using the hydrogen-$\beta$ spectral line agree among different line-shape models for $n_{\rm e}\lesssim30\times10^{16}$ cm$^{-3}$, we now see divergence between models. These are densities beyond the range previously benchmarked in the laboratory, meaning theoretical profiles in this regime have not been fully validated. Experimentally exploring these higher densities enables us to test and constrain different line-profile models, as the differences in their relative H-Balmer line shapes are more pronounced at such conditions. These experiments also aid in our study of occupation probabilities because we can measure these from relative line strengths.

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.

Recent large-scale infrared surveys have been revealing stellar populations in the inner Galaxy seen through strong interstellar extinction in the disk. In particular, classical Cepheids with their period-luminosity and period-age relations are useful tracers of Galactic structure and evolution. Interesting groups of Cepheids reported recently include four Cepheids in the Nuclear Stellar Disk (NSD), about 200 pc around the Galactic Centre, found by Matsunaga et al. and those spread across the inner part of the disk reported by Dekany and collaborators. We here report our discovery of nearly thirty classical Cepheids towards the bulge region, some of which are common with Dekany et al., and discuss the large impact of the reddening correction on distance estimates for these objects. Assuming that the four Cepheids in the NSD are located at the distance of the Galactic Centre and that the near-infrared extinction law, i.e. wavelength dependency of the interstellar extinction, is not systematically different between the NSD and other bulge lines-of-sight, most of the other Cepheids presented here are located significantly further than the Galactic Centre. This suggests a lack of Cepheids in the inner 2.5 kpc region of the Galactic disk except the NSD. Recent radio observations show a similar distribution of star-forming regions.

We report spectroscopic and photometric follow-up of the peculiar nova V5852~Sgr (discovered as OGLE-2015-NOVA-01), which exhibits a combination of features from different nova classes. The photometry shows a flat-topped light curve with quasi-periodic oscillations, then a smooth decline followed by two fainter recoveries in brightness. Spectroscopy with the Southern African Large Telescope shows first a classical nova with an Fe II or Fe IIb spectral type. In the later spectrum, broad emissions from helium, nitrogen and oxygen are prominent and the iron has faded which could be an indication to the start of the nebular phase. The line widths suggest ejection velocities around $1000\,{\rm km\,s^{-1}}$. The nova is in the direction of the Galactic bulge and is heavily reddened by an uncertain amount. The $V$ magnitude 16 days after maximum enables a distance to be estimated and this suggests that the nova may be in the extreme trailing stream of the Sagittarius dwarf spheroidal galaxy. If so it is the first nova to be detected from that, or from any dwarf spheroidal galaxy. Given the uncertainty of the method and the unusual light curve we cannot rule out the possibility that it is in the bulge or even the Galactic disk behind the bulge.

We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers.

We report the results of VERA multi-epoch VLBI 22 GHz water maser observations of S255IR-SMA1, a massive young stellar object located in the S255 star forming region. By annual parallax the source distance was measured as D = 1.78 +-0.12 kpc and the source systemic motion was (u alpha cos d, u d) = (-0.13 +- 0.20, -0.06 +- 0.27) mas yr-1. Masers appear to trace a U-shaped bow shock whose morphology and proper motions are well reproduced by a jet-driven outflow model with a jet radius of about 6 AU. The maser data, in the context of other works in the literature, reveal ejections from S255IR-SMA1 to be episodic, operating on timescales of ~1000 years.

We present extensively dense observations of the super-Chandrasekhar supernova (SC SN) candidate SN 2012dn from $-11$ to $+140$ days after the date of its $B$-band maximum in the optical and near-infrared (NIR) wavelengths conducted through the OISTER ToO program. The NIR light curves and color evolutions up to 35 days after the $B$-band maximum provided an excellent match with those of another SC SN 2009dc, providing a further support to the nature of SN 2012dn as a SC SN. We found that SN 2012dn exhibited strong excesses in the NIR wavelengths from $30$ days after the $B$-band maximum. The $H$ and $K_{s}$-band light curves exhibited much later maximum dates at $40$ and $70$ days after the $B$-band maximum, respectively, compared with those of normal SNe Ia. The $H$ and $K_{s}$-band light curves subtracted by those of SN 2009dc displayed plateaued evolutions, indicating a NIR echo from the surrounding dust. The distance to the inner boundary of the dust shell is limited to be $4.8 - 6.4\times10^{-2}$ pc. No emission lines were found in its early phase spectrum, suggesting that the ejecta-CSM interaction could not occur. On the other hand, we found no signature that strongly supports the scenario of dust formation. The mass loss rate of the pre-explosion system is estimated to be $10^{-6}-10^{-5}$ M$_{\odot}$ yr$^{-1}$, assuming that the wind velocity of the system is $10-100$ km~s$^{-1}$, which suggests that the progenitor of SN 2012dn could be a recurrent nova system. We conclude that the progenitor of this SC SN could be explained by the single-degenerate scenario.

We present simultaneous optical and near-infrared (NIR) polarimetric results for the black hole binary V404 Cygni spanning the duration of its 7-day long optically-brightest phase of its 2015 June outburst. The simultaneous R and Ks-band light curves showed almost the same temporal variation except for the isolated (~30 min duration) orphan Ks-band flare observed at MJD 57193.54. We did not find any significant temporal variation of polarization degree (PD) and position angle (PA) in both R and Ks bands throughout our observations, including the duration of the orphan NIR flare. We show that the observed PD and PA are predominantly interstellar in origin by comparing the V404 Cyg polarimetric results with those of the surrounding sources within the 7'x7' field-of-view. The low intrinsic PD (less than a few percent) implies that the optical and NIR emissions are dominated by either disk or optically-thick synchrotron emission, or both. We also present the broadband spectra of V404 Cyg during the orphan NIR flare and a relatively faint and steady state by including quasi-simultaneous Swift/XRT and INTEGRAL fluxes. By adopting a single-zone synchrotron plus inverse-Compton model as widely used in modeling of blazars, we constrained the parameters of a putative jet. Because the jet synchrotron component cannot exceed the Swift/XRT disk/corona flux, the cutoff Lorentz factor in the electron energy distribution is constrained to be <10^2, suggesting particle acceleration is less efficient in this microquasar jet outburst compared to AGN jets. We also suggest that the loading of the baryon component inside the jet is inevitable based on energetic arguments.

We present infrared multi-epoch observations of the dust forming nova V1280 Sco over $\sim$2000 days from the outburst. The temporal evolution of the infrared spectral energy distributions at 1272, 1616 and 1947 days can be explained by the emissions produced by amorphous carbon dust of mass (6.6--8.7)$\times$10$^{-8}$M$_{\odot}$ with a representative grain size of 0.01$~\mu$m and astronomical silicate dust of mass (3.4--4.3)$\times$10$^{-7}$M$_{\odot}$ with a representative grain size of 0.3--0.5$~\mu$m. Both of these dust species travel farther away from the white dwarf without an apparent mass evolution throughout those later epochs. The dust formation scenario around V1280 Sco suggested from our analyses is that the amorphous carbon dust is formed in the nova ejecta followed by the formation of silicate dust in the expanding nova ejecta or as a result of the interaction between the nova wind and the circumstellar medium.

We conducted an optical and near-infrared polarimetric observation of the highly dormant Jupiter-Family Comet, 209P/LINEAR. Because of its low activity, we were able to determine the linear polarization degrees of the coma dust particles and nucleus independently, that is $P_n$=30.3$^{+1.3}_{-0.9}$% at $\alpha$=92.2$^\circ$ and $P_n$=31.0$^{+1.0}_{-0.7}$% at $\alpha$=99.5$^\circ$ for the nucleus, and $P_c$=28.8$^{+0.4}_{-0.4}$% at $\alpha$=92.2$^\circ$ and 29.6$^{+0.3}_{-0.3}$% at $\alpha$=99.5$^\circ$ for the coma. We detected no significant variation in $P$ at the phase angle coverage of 92.2$^\circ$-99.5$^\circ$, which may imply that the obtained polarization degrees are nearly at maximum in the phase-polarization curves. By fitting with an empirical function, we obtained the maximum values of linear polarization degrees $P_\mathrm{max}$=30.8% for the nucleus and $P_\mathrm{max}$=29.6% for the dust coma. The $P_\mathrm{max}$ of the dust coma is consistent with those of dust-rich comets. The low geometric albedo of $P_v$=0.05 was derived from the slope-albedo relationship and was associated with high $P_\mathrm{max}$. We examined $P_\mathrm{max}$-albedo relations between asteroids and 209P, and found that the so-called Umov law seems to be applicable on this cometary surface.

We report on annual parallax and proper motion observations of H2O masers in S235AB-MIR, which is a massive young stellar object in the Perseus Arm. Using multi-epoch VLBI astrometry we measured a parallax of pi = 0.63 +- 0.03 mas, corresponding to a trigonometric distance of D = 1.56+-0.09 kpc, and source proper motion of ( u alpha cos d , u d) = (0.79 +- 0.12, -2.41 +- 0.14) mas/yr. Water masers trace a jet of diameter 15 au which exhibits a definite radial velocity gradient perpendicular to its axis. 3D maser kinematics were well modelled by a rotating cylinder with physical parameters: v_out = 45+-2 km/s, v_rot = 22+-3 km/s, i = 12+-2 degrees, which are the outflow velocity, tangential rotation velocity and line-of-sight inclination, respectively. One maser feature exhibited steady acceleration which may be related to the jet rotation. During our 15 month VLBI programme there were three `maser burst' events caught `in the act' which were caused by the overlapping of masers along the line of sight.

The Ophiuchus cluster, at a redshift z=0.0296, is known from X-rays to be one of the most massive nearby clusters, but due to its very low Galactic latitude its optical properties have not been investigated in detail. We discuss the optical properties of the galaxies in the Ophiuchus cluster, in particular with the aim of understanding better its dynamical properties. We have obtained deep optical imaging in several bands with various telescopes, and applied a sophisticated method to model and subtract the contributions of stars in order to measure galaxy magnitudes as accurately as possible. The colour-magnitude relations obtained show that there are hardly any blue galaxies in Ophiuchus (at least brighter than r'<=19.5), and this is confirmed by the fact that we only detect two galaxies in Halpha. We also obtained a number of spectra with ESO-FORS2, that we combined with previously available redshifts. Altogether, we have 152 galaxies with spectroscopic redshifts in the 0.02<=z<=0.04 range, and 89 galaxies with both a redshift within the cluster redshift range and a measured r' band magnitude (limited to the Megacam 1x1 deg^2 field). A complete dynamical analysis based on the galaxy redshifts available shows that the overall cluster is relaxed and has a mass of 1.1x10^15 solar masses. The Sernal-Gerbal method detects a main structure and a much smaller substructure that are not separated in projection. From its dynamical properties derived from optical data, the Ophiuchus cluster seems to be overall a relaxed structure, or at most a minor merger, though in X-rays the central region (radius ~ 150 kpc) may show evidence for merging effects.

We spectroscopically measure multiple hydrogen Balmer line profiles from laboratory plasmas to investigate the theoretical line profiles used in white dwarf atmosphere models. X-ray radiation produced at the Z Pulsed Power Facility at Sandia National Laboratories initiates plasma formation in a hydrogen-filled gas cell, replicating white dwarf photospheric conditions. Here we present time-resolved measurements of H$\beta$ and fit this line using different theoretical line profiles to diagnose electron density, $n_{\rm e}$, and $n=2$ level population, $n_2$. Aided by synthetic tests, we characterize the validity of our diagnostic method for this experimental platform. During a single experiment, we infer a continuous range of electron densities increasing from $n_{\rm e}\sim4$ to $\sim30\times10^{16}\,$cm$^{-3}$ throughout a 120-ns evolution of our plasma. Also, we observe $n_2$ to be initially elevated with respect to local thermodynamic equilibrium (LTE); it then equilibrates within $\sim55\,$ns to become consistent with LTE. This supports our electron-temperature determination of $T_{\rm e}\sim1.3\,$eV ($\sim15,000\,$K) after this time. At $n_{\rm e}\gtrsim10^{17}\,$cm$^{-3}$, we find that computer-simulation-based line-profile calculations provide better fits (lower reduced $\chi^2$) than the line profiles currently used in the white dwarf astronomy community. The inferred conditions, however, are in good quantitative agreement. This work establishes an experimental foundation for the future investigation of relative shapes and strengths between different hydrogen Balmer lines.

We report observations of the Type Iax supernova (SN Iax) 2012Z at optical and near-infrared wavelengths from immediately after the explosion until $\sim$ $260$ days after the maximum luminosity using the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) Target-of-Opportunity (ToO) program and the Subaru telescope. We found that the near-infrared (NIR) light curve evolutions and color evolutions are similar to those of SNe Iax 2005hk and 2008ha. The NIR absolute magnitudes ($M_{J}\sim-18.1$ mag and $M_{H}\sim-18.3$ mag) and the rate of decline of the light curve ($\Delta$ $m_{15}$($B$)$=1.6 \pm 0.1$ mag) are very similar to those of SN 2005hk ($M_{J}\sim-17.7$ mag, $M_{H}\sim$$-18.0$ mag, and $\Delta$ $m_{15}$($B$)$\sim1.6$ mag), yet differ significantly from SNe 2008ha and 2010ae ($M_{J}\sim-14 - -15$ mag and $\Delta$ $m_{15}$($B$)$\sim2.4-2.7$ mag). The estimated rise time is $12.0 \pm 3.0$ days, which is significantly shorter than that of SN 2005hk or any other Ia SNe. The rapid rise indicates that the $^{56}$Ni distribution may extend into the outer layer or that the effective opacity may be lower than that in normal SNe Ia. The late-phase spectrum exhibits broader emission lines than those of SN 2005hk by a factor of 6--8. Such high velocities of the emission lines indicate that the density profile of the inner ejecta extends more than that of SN 2005hk. We argue that the most favored explosion scenario is a `failed deflagration' model, although the pulsational delayed detonations is not excluded.

We present the discovery of two new X-ray transients in archival Chandra data. The first transient, XRT 110103, occurred in January 2011 and shows a sharp rise of at least three orders of magnitude in count rate in less than 10 s, a flat peak for about 20 s and decays by two orders of magnitude in the next 60 s. We find no optical or infrared counterpart to this event in preexisting survey data or in an observation taken by the SIRIUS instrument at the Infrared Survey Facility 2.1 yr after the transient, providing limiting magnitudes of J>18.1, H>17.6 and Ks>16.3. This event shows similarities to the transient previously reported in Jonker et al. which was interpreted as the possible tidal disruption of a white dwarf by an intermediate mass black hole. We discuss the possibility that these transients originate from the same type of event. If we assume these events are related a rough estimate of the rates gives 1.4*10^5 per year over the whole sky with a peak 0.3-7 keV X-ray flux greater than 2*10^-10 erg cm^-2 s^-1 . The second transient, XRT 120830, occurred in August 2012 and shows a rise of at least three orders of magnitude in count rate and a subsequent decay of around one order of magnitude all within 10 s, followed by a slower quasi-exponential decay over the remaining 30 ks of the observation. We detect a likely infrared counterpart with magnitudes J=16.70+/-0.06, H=15.92+/-0.04 and Ks=15.37+/-0.06 which shows an average proper motion of 74+/-19 milliarcsec per year compared to archival 2MASS observations. The JHKs magnitudes, proper motion and X-ray flux of XRT 120830 are consistent with a bright flare from a nearby late M or early L dwarf.

We present new near-infrared (NIR) light-curve templates for fundamental (FU, JHK) and first overtone (FO, J) Cepheids. The new templates together with PL and PW relations provide Cepheid distances from single-epoch observations with a precision only limited by the intrinsic accuracy of the method adopted. The templates rely on a very large set of Galactic and Magellanic Clouds (MCs) Cepheids (FU,~600; FO,~200) with well sampled NIR (IRSF data) and optical (V,I; OGLE data) light curves. To properly trace the change in the shape of the light curve as a function of period, we split the sample of calibrating Cepheids into 10 different period bins. The templates for the first time cover FO Cepheids and the FU short-period Cepheids (P<5 days). Moreover, the zero-point phase is anchored to the phase of the mean magnitude along the rising branch. The new approach has several advantages in sampling the light curve of bump Cepheids when compared with the phase of maximum light. We also provide new estimates of the NIR-to-optical amplitude ratios for FU and FO Cepheids. We perform detailed analytical fits using both 7th-order Fourier series and multi-Gaussian periodic functions. The latter are characterized by a smaller number of free parameters (9 vs 15). Mean NIR magnitudes based on the new templates are up to 80% more accurate than single-epoch measurements and up to 50% more accurate than mean magnitudes based on previous templates, with typical associated uncertainties ranging from 0.015 mag (J) to 0.019 mag (K). Moreover, the errors on individual distances of Small MC Cepheids derived from NIR PW relations, are essentially reduced to the intrinsic scatter of the adopted relations. Thus, the new templates are the ultimate tool to estimate precise Cepheid distances from NIR single-epoch observations, which can be adopted to derive the 3D structure of the MCs.

Monitor of All sky X-ray Image (MAXI) discovered a new outburst of an X-ray transient source named MAXI J1421-613. Because of the detection of three X-ray bursts from the source, it was identified as a neutron star low-mass X-ray binary. The results of data analyses of the MAXI GSC and the Swift XRT follow-up observations suggest that the spectral hardness remained unchanged during the first two weeks of the outburst. All the XRT spectra in the 0.5-10 keV band can be well explained by thermal Comptonization of multi-color disk blackbody emission. The photon index of the Comptonized component is $\approx$ 2, which is typical of low-mass X-ray binaries in the low/hard state. Since X-ray bursts have a maximum peak luminosity, it is possible to estimate the (maximum) distance from its observed peak flux. The peak flux of the second X-ray burst, which was observed by the GSC, is about 5 photons cm$^{-2}$ s$^{-1}$. By assuming a blackbody spectrum of 2.5 keV, the maximum distance to the source is estimated as 7 kpc. The position of this source is contained by the large error regions of two bright X-ray sources detected with Orbiting Solar Observatory-7 (OSO-7) in 1970s. Besides this, no past activities at the XRT position are reported in the literature. If MAXI J1421-613 is the same source as (one of) them, the outburst observed with MAXI may have occurred after the quiescence of 30-40 years.

Multi-wavelength observations of the black-widow binary system 2FGL J2339.6-0532 are reported. The Fermi gamma-ray source 2FGL J2339.6-0532 was recently categorized as a black widow in which a recycled millisecond pulsar (MSP) is evaporating up the companion star with its powerful pulsar wind. Our optical observations show clear sinusoidal light curves due to the asymmetric temperature distribution of the companion star. Assuming a simple geometry, we constrained the range of the inclination angle of the binary system to 52$^{\circ}$ < i < 59$^{\circ}$, which enables us to discuss the interaction between the pulsar wind and the companion in detail. The X-ray spectrum consists of two components: a soft, steady component that seems to originate from the surface of the MSP, and a hard variable component from the wind-termination shock near the companion star. The measured X-ray luminosity is comparable to the bolometric luminosity of the companion, meaning that the heating efficiency is less than 0.5. In the companion orbit, 10$^{11}$ cm from the pulsar, the pulsar wind is already in particle dominant-stage, with a magnetization parameter of $\sigma$ < 0.1. In addition, we precisely investigated the time variations of the X-ray periodograms and detected a weakening of orbital modulation. The observed phenomenon may be related to an unstable pulsar-wind activity or a weak mass accretion, both of which can result in the temporal extinction of radio-pulse.

We report a new observation of the Jupiter-family comet 209P/LINEAR during its 2014 return. The comet is recognized as a dust source of a new meteor shower, the May Camelopardalids. 209P/LINEAR was apparently inactive at a heliocentric distance rh = 1.6 au and showed weak activity at rh < 1.4 au. We found an active region of <0.001% of the entire nuclear surface during the comet's dormant phase. An edge-on image suggests that particles up to 1 cm in size (with an uncertainty of factor 3-5) were ejected following a differential power-law size distribution with index q=-3.25+-0.10. We derived a mass loss rate of 2-10 kg/s during the active phase and a total mass of ~5x10^7 kg during the 2014 return. The ejection terminal velocity of millimeter- to centimeter-sized particles was 1-4 m/s, which is comparable to the escape velocity from the nucleus (1.4 m/s). These results imply that such large meteoric particles marginally escaped from the highly dormant comet nucleus via the gas drag force only within a few months of the perihelion passage.

We have made near-infrared (JHKs) imaging polarimetry of a bright-rimmed cloud (SFO 74). The polarization vector maps clearly show that the magnetic field in the layer just behind the bright rim is running along the rim, quite different from its ambient magnetic field. The direction of the magnetic field just behind the tip rim is almost perpendicular to that of the incident UV radiation, and the magnetic field configuration appears to be symmetric as a whole with respect to the cloud symmetry axis. We estimated the column and number densities in the two regions (just inside and far inside the tip rim), and then derived the magnetic field strength, applying the Chandrasekhar-Fermi method. The estimated magnetic field strength just inside the tip rim, ~90 uG, is stronger than that far inside, ~30 uG. This suggests that the magnetic field strength just inside the tip rim is enhanced by the UV radiation induced shock. The shock increases the density within the top layer around the tip, and thus increases the strength of the magnetic field. The magnetic pressure seems to be comparable to the turbulent one just inside the tip rim, implying a significant contribution of the magnetic field to the total internal pressure. The mass-to-flux ratio was estimated to be close to the critical value just inside the tip rim. We speculate that the flat-topped bright rim of SFO 74 could be formed by the magnetic field effect.

Our preliminary results from laboratory experiments studying white dwarf (WD) photospheres show a systematic difference between experimental plasma conditions inferred from measured H$\beta$ absorption line profiles versus those from H$\gamma$. One hypothesis for this discrepancy is an inaccuracy in the relative theoretical line profiles of these two transitions. This is intriguing because atmospheric parameters inferred from H Balmer lines in observed WD spectra show systematic trends such that inferred surface gravities decrease with increasing principal quantum number, $n$. If conditions inferred from lower-$n$ Balmer lines are indeed more accurate, this suggests that spectroscopically determined DA WD masses may be greater than previously thought and in better agreement with the mean mass determined from gravitational redshifts.

We have conducted a systematic near-infrared circular polarization (CP) survey in star-forming regions, covering high-mass, intermediate-mass, and low-mass young stellar objects. All the observations were made using the SIRPOL imaging polarimeter on the Infrared Survey Facility 1.4 m telescope at the South African Astronomical Observatory. We present the polarization properties of ten sub-regions in six star-forming regions. The polarization patterns, extents, and maximum degrees of linear and circular polarizations are used to determine the prevalence and origin of CP in the star-forming regions. Our results show that the CP pattern is quadrupolar in general, the CP regions are extensive, up to 0.65 pc, the CP degrees are high, up to 20 %, and the CP degrees decrease systematically from high- to low-mass young stellar objects. The results are consistent with dichroic extinction mechanisms generating the high degrees of CP in star forming regions.

We report our measurement of the trigonometric distance and proper motion IRAS 20056+3350, obtained from the annual parallax of H2O masers. Our distance of D = 4.69 +0.65-0.51 kpc, which is more than two times larger than the near kinematic distance adopted in the literature, places IRAS 20056+3350 at the leading tip of the Local arm, and proximal to the Solar circle. Using our distance we re-evaluate past observations to reveal IRAS 20056+3350 as a site of massive star formation at a young stage of evolution. This result is consistent with the spectral energy distribution of the source evaluated with published photometric data from UKIDSS, WISE, AKARI, IRAS and sub-millimetre continuum. Both analytical approaches reveal the luminosity of the region to be 2.4 x 10^4 Lo, and suggest that IRAS 20056+3350 is forming an embedded star of 16 Mo. We estimated the proper motion of IRAS 20056+3350 to be ($\mu_{\alpha}\cos\delta$, $\mu_{\delta}$) = ($-2.62\pm0.33$, $-5.65\pm0.52$) mas yr$^{-1}$ from the group motion of H$_{2}$O masers, and use our results to estimate the angular velocity of Galactic rotation at the Galactocentric distance of the Sun, $\Omega_{0} = 28.71\pm0.63$ km s$^{-1}$ kpc$^{-1}$, which is consistent with the values obtained for other tangent point and Solar circle objects.

We present results of optical polarimetric and multi-band photometric observations for gamma-ray loud narrow-line Seyfert 1 galaxy 1H 0323+342. This object has been monitored by 1.5 m Kanata telescope since 2012 September but following a gamma-ray flux enhancement detected by Fermi-LAT on MJD 56483 (2013 July 10) dense follow-up was performed by ten 0.5-2.0 m telescopes in Japan over one week. The 2-year R_C-band light curve showed clear brightening corresponding to the gamma-ray flux increase and then decayed gradually. The high state as a whole lasted for ~20 days, during which we clearly detected optical polarization from this object. The polarization degree (PD) of the source increased from 0-1% in quiescence to ~3% at maximum and then declined to the quiescent level, with the duration of the enhancement of less than 10 days. The moderate PD around the peak allowed us to precisely measure the daily polarization angle (PA). As a result, we found that the daily PAs were almost constant and aligned to the jet axis, suggesting that the magnetic field direction at the emission region is transverse to the jet. This implies either a presence of helical/toroidal magnetic field or transverse magnetic field compressed by shock(s). We also found small-amplitude intra-night variability during the 2-hour continuous exposure on a single night. We discuss these findings based on the turbulent multi-zone model recently advocated by Marscher (2014). Optical to ultraviolet (UV) spectrum showed a rising shape in the higher frequency and the UV magnitude measured by Swift/UVOT was steady even during the flaring state, suggesting that thermal emission from accretion disk is dominant in that band.

Experimental tests are in progress to evaluate the accuracy of the modeled iron opacity at solar interior conditions [J.E. Bailey et al., Phys. Plasmas 16, 058101 (2009)]. The iron sample is placed on top of the Sandia National Laboratories z-pinch dynamic hohlraum (ZPDH) radiation source. The samples are heated to 150 - 200 eV electron temperatures and 7e21 - 4e22 e/cc electron densities by the ZPDH radiation and backlit at its stagnation [T. Nagayama et al., Phys. Plasmas 21, 056502 (2014)]. The backlighter attenuated by the heated sample plasma is measured by four spectrometers along +/- 9 degree with respect to the z-pinch axis to infer the sample iron opacity. Here we describe measurements of the source-to-sample distance that exploit the parallax of spectrometers that view the half-moon-shaped sample from +/-9 degree. The measured sample temperature decreases with increased source-to-sample distance. This distance must be taken into account for understanding the sample heating.

We investigated the magnitude-phase relation of (162173) 1999 JU3, a target asteroid for the JAXA Hayabusa 2 sample return mission. We initially employed the international Astronomical Union's H-G formalism but found that it fits less well using a single set of parameters. To improve the inadequate fit, we employed two photometric functions, the Shevchenko and Hapke functions. With the Shevchenko function, we found that the magnitude-phase relation exhibits linear behavior in a wide phase angle range (alpha = 5-75 deg) and shows weak nonlinear opposition brightening at alpha< 5 deg, providing a more reliable absolute magnitude of Hv = 19.25 +- 0.03. The phase slope (0.039 +- 0.001 mag/deg) and opposition effect amplitude (parameterized by the ratio of intensity at alpha=0.3 deg to that at alpha=5 deg, I(0.3)/I(5)=1.31+-0.05) are consistent with those of typical C-type asteroids. We also attempted to determine the parameters for the Hapke model, which are applicable for constructing the surface reflectance map with the Hayabusa 2 onboard cameras. Although we could not constrain the full set of Hapke parameters, we obtained possible values, w=0.041, g=-0.38, B0=1.43, and h=0.050, assuming a surface roughness parameter theta=20 deg. By combining our photometric study with a thermal model of the asteroid (Mueller et al. in preparation), we obtained a geometric albedo of pv = 0.047 +- 0.003, phase integral q = 0.32 +- 0.03, and Bond albedo AB = 0.014 +- 0.002, which are commensurate with the values for common C-type asteroids.

WASP-80b is a warm Jupiter transiting a bright late-K/early-M dwarf, providing a good opportunity to extend the atmospheric study of hot Jupiters toward the lower temperature regime. We report multi-band, multi-epoch transit observations of WASP-80b by using three ground-based telescopes covering from optical (g', Rc, and Ic bands) to near-infrared (NIR; J, H, and Ks bands) wavelengths. We observe 5 primary transits, each of which in 3 or 4 different bands simultaneously, obtaining 17 independent transit light curves. Combining them with results from previous works, we find that the observed transmission spectrum is largely consistent with both a solar abundance and thick cloud atmospheric models at 1.7$\sigma$ discrepancy level. On the other hand, we find a marginal spectral rise in optical region compared to the NIR region at 2.9$\sigma$ level, which possibly indicates the existence of haze in the atmosphere. We simulate theoretical transmission spectra for a solar abundance but hazy atmosphere, finding that a model with equilibrium temperature of 600 K can explain the observed data well, having a discrepancy level of 1.0$\sigma$. We also search for transit timing variations, but find no timing excess larger than 50 s from a linear ephemeris. In addition, we conduct 43 day long photometric monitoring of the host star in the optical bands, finding no significant variation in the stellar brightness. Combined with the fact that no spot-crossing event is observed in the five transits, our results confirm previous findings that the host star appears quiet for spot activities, despite the indications of strong chromospheric activities.

We report the results from an X-ray and near-infrared observation of the Galactic black hole binary 4U 1630--47 in the very high state, performed with {\it Suzaku} and IRSF around the peak of the 2012 September-October outburst. The X-ray spectrum is approximated by a steep power law, with photon index of 3.2, identifying the source as being in the very high state. A more detailed fit shows that the X-ray continuum is well described by a multi-color disc, together with thermal and non-thermal Comptonization. The inner disc appears slightly truncated by comparison with a previous high/soft state of this source, even taking into account energetic coupling between the disc and corona, although there are uncertainties due to the dust scattering correction. The near-infrared fluxes are higher than the extrapolated disc model, showing that there is a contribution from irradiation in the outer disk and/or the companion star at these wavelengths. Our X-ray spectra do not show the Doppler shifted iron emission lines indicating a baryonic jet which were seen four days previously in an XMM-Newton observation, despite the source being in a similar state. There are also no significant absorption lines from highly ionized irons as are seen in the previous high/soft state data. We show that the increased source luminosity is not enough on its own to make the wind so highly ionized as to be undetectable. This shows that the disc wind has changed in terms of its launch radius and/or density compared to the high/soft state.

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.

Context. The Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey is a variability survey of the Milky Way bulge and an adjacent section of the disk carried out from 2010 on ESO Visible and Infrared Survey Telescope for Astronomy (VISTA). VVV will eventually deliver a deep near-IR atlas with photometry and positions in five passbands (ZYJHK_S) and a catalogue of 1-10 million variable point sources - mostly unknown - which require classifications. Aims. The main goal of the VVV Templates Project, that we introduce in this work, is to develop and test the machine-learning algorithms for the automated classification of the VVV light-curves. As VVV is the first massive, multi-epoch survey of stellar variability in the near-infrared, the template light-curves that are required for training the classification algorithms are not available. In the first paper of the series we describe the construction of this comprehensive database of infrared stellar variability. Methods. First we performed a systematic search in the literature and public data archives, second, we coordinated a worldwide observational campaign, and third we exploited the VVV variability database itself on (optically) well-known stars to gather high-quality infrared light-curves of several hundreds of variable stars. Results. We have now collected a significant (and still increasing) number of infrared template light-curves. This database will be used as a training-set for the machine-learning algorithms that will automatically classify the light-curves produced by VVV. The results of such an automated classification will be covered in forthcoming papers of the series.

We report on the results from Suzaku observations of the Galactic black hole X-ray binary H 1743-322 in the low/hard state during its outburst in 2012 October. We appropriately take into account the effects of dust-scattering to accurately analyze the X-ray spectra. The time-averaged spectra in the 1-200 keV band are dominated by a hard power-law component of a photon index of \approx 1.6 with a high-energy cutoff at \approx 60 keV, which is well described with Comptonization of the disk emission by hot corona. We estimate the inner disk radius from the multi-color disk component, and find that it is 1.3-2.3 times larger than the radius in the high/soft state. This suggests that the standard disk was not extended to the innermost stable circular orbit. A reflection component from the disk is detected with R = \Omega/2\pi \approx 0.6 (\Omega is the solid angle). We also successfully estimate the stable disk component in a way independent of the time-averaged spectral modeling, by analyzing short-term spectral variability on the \sim 1-sec timescale. A weak low-frequency quasi-periodic oscillation (LF QPO) at 0.1-0.2 Hz is detected, whose frequency is found to correlate with the X-ray luminosity and photon index. This result may be explained by the evolution of the disk truncation radius.

Experimental tests are in progress to evaluate the accuracy of the modeled iron opacity at solar interior conditions, in particular to better constrain the solar abundance problem [S. Basu and H.M. Antia, Physics Reports 457, 217 (2008)]. Here we describe measurements addressing three of the key requirements for reliable opacity experiments: control of sample conditions, independent sample condition diagnostics, and verification of sample condition uniformity. The opacity samples consist of iron/magnesium layers tamped by plastic. By changing the plastic thicknesses, we have controlled the iron plasma conditions to reach i) Te=167+/-3 eV and ne=(7.1+/-1.5)e21 e/cc, ii) Te=170+/-2 eV and ne=(2.0+/-0.2)e22 e/cc, and iii) Te=196+/-6 eV and ne=(3.8+/-0.8)e22 e/cc, which were measured by magnesium tracer K-shell spectroscopy. The opacity sample non-uniformity was directly measured by a separate experiment where Al is mixed into the side of the sample facing the radiation source and Mg into the other side. The iron condition was confirmed to be uniform within their measurement uncertainties by Al and Mg K-shell spectroscopy. The conditions are suitable for testing opacity calculations needed for modeling the solar interior, other stars, and high energy density plasmas.

We report the results of VLBI observations of H$_{2}$O masers in the IRAS 20143+3634 star forming region using VERA (VLBI Exploration of Radio Astronomy). By tracking masers for a period of over two years we measured a trigonometric parallax of $\pi = 0.367 \pm 0.037$ mas, corresponding to a source distance of $D = 2.72 ^{+0.31}_{-0.25}$ kpc and placing it in the Local spiral arm. Our trigonometric distance is just 60% of the previous estimation based on radial velocity, significantly impacting the astrophysics of the source. We measured proper motions of $-2.99 \pm 0.16$ mas yr$^{-1}$ and $-4.37 \pm 0.43$ mas yr$^{-1}$ in R.A. and Decl. respectively, which were used to estimate the peculiar motion of the source as $(U_{s},V_{s},W_{s}) = (-0.9 \pm 2.9, -8.5 \pm 1.6, +8.0 \pm 4.3)$ km s$^{-1}$ for $R_0=8$ kpc and $\Theta_0=221$ km s$^{-1}$, and $(U_{s},V_{s},W_{s}) = (-1.0 \pm 2.9, -9.3 \pm 1.5, +8.0 \pm 4.3)$ km s$^{-1}$ for $R_0=8.5$ kpc and $\Theta_0=235$ km s$^{-1}$. IRAS 20143+3634 was found to be located near the tangent point in the Cygnus direction. Using our observations we derived the angular velocity of Galactic rotation of the local standard of rest (LSR), $\Omega_{0} = 27.3 \pm 1.6$ km s$^{-1}$ kpc$^{-1}$, which is consistent with previous values derived using VLBI astrometry of SFRs at the tangent points and Solar circle. It is higher than the value recommended by the IAU of $\Omega_{0} = 25.9$ km s$^{-1}$ kpc$^{-1}$ which was calculated using the Galactocentric distance of the Sun and circular velocity of the LSR.

Aims. Eta Car's ultra-violet, optical, and X-ray light curves and its spectrum suggest a physical change in its stellar wind over the last decade. It was proposed that the mass-loss rate decreased by a factor of about 2 in the last 15 years. We complement these recent results by investigating the past evolution and the current state of eta Car in the near-infrared (IR). Methods. We present JHKL photometry of eta Car obtained at SAAO Sutherland from 2004-2013 with the Mk II photometer at the 0.75-m telescope and JHKs photometry with SIRIUS at the 1.4-m IRSF telescope from 2012-2013. The near-IR light curves since 1972 are analyzed. Results. The long-term brightening trends in eta Car's JHKL light curves were discontinuous around the 1998 periastron passage. After 1998, the star shows excess emission above the extrapolated trend from earlier dates, foremost in J and H, and the blueward, cyclical progression in its near-IR colors is accelerated. The near-IR color evolution is strongly correlated with the periastron passages. After correcting for the secular trend we find that the color evolution matches an apparent increase in blackbody temperature of an optically thick near-IR emitting plasma component from about 3500 to 6000 K over the last 20 years. Conclusions. We suggest that the changing near-IR emission may be caused by variability in optically thick bremsstrahlung emission. Periastron passages play a key role in the observed excess near-IR emission after 1998 and the long-term color evolution. We thus propose as a hypothesis that angular momentum transfer (via tidal acceleration) during periastron passages leads to sudden changes in eta Car's atmosphere resulting in a long-term decrease in the mass-loss rate.

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

We report the results from X-ray studies of the newly discovered black hole candidate MAXI J1305-704 based on Suzaku and Swift observations in the low/hard and high/soft states, respectively. The long Suzaku observation shows two types of clear absorption dips, both of which recur on a dip interval of 9.74 +- 0.04 hours, which we identify with the orbital period. There is also partially ionized absorption in the non-dip (persistent) emission in both the high/soft state and, very unusually, the low/hard state. However, this absorption (in both states) has substantially lower ionization than that seen in other high inclination systems, where the material forms a homogeneous disk wind. Here instead the absorption is most probably associated with clumpy, compact structures associated with the dipping material, which we see uniquely in this source likely because we view it at a very large inclination angle. A large inclination angle is also favored, together with a low black hole mass, to explain the high disk temperature seen in the fairly low luminosity high/soft state, as Doppler boosting enhances the disk temperature at high inclination. The disk radius inferred from these data is significantly smaller than that of the soft component seen in the low/hard state, supporting models where the disk is truncated at low luminosities. We find, however, that the lack of variability power on time scales of ~50 sec in the Suzaku low/hard state data is difficult to explain, even with a low mass black hole.

We present the first measurement of the absolute proper motions of IRAS 00259+5625 (CB3, LBN594) associated with the HI loop called the "NGC281 superbubble" that extends from the Galactic plane over ~300 pc toward decreasing galactic latitude. The proper motion components measured with VERA are (mu_alpha cos(delta), mu_delta) = (-2.48 +/- 0.32, -2.85 +/- 0.65) mas yr^{-1}, converted into (mu_l cos(b), mu_b) = (-2.72 +/- 0.32, -2.62 +/- 0.65) mas yr^{-1} in the Galactic coordinates. The measured proper motion perpendicular to the Galactic plane (mu_b) shows vertical motion away from the Galactic plane with a significance of about ~4-sigma. As for the source distance, the distance measured with VERA is marginal, 2.4^{+1.0}_{-0.6} kpc. Using the distance, an absolute vertical motion (v_{b}) of -17.9 +/- 12.2 km s^{-1} is determined with ~1.5-sigma significance. The tendency of the large vertical motion is consistent with previous VLBI results for NGC 281 associated with the same superbubble. Thus, our VLBI results indicate the superbubble expansion motion whose origin is believed to be sequential supernova explosions.

We present 5 new transit light curves of GJ 1214b taken in BJHKs-bands. Two transits were observed in B-band using the Suprime-Cam and the FOCAS instruments onboard the Subaru 8.2m telescope, and one transit was done in JHKs-bands simultaneously with the SIRIUS camera on the IRSF 1.4m telescope. MCMC analyses show that the planet-to-star radius ratios are, Rp/Rs = 0.11651 \pm 0.00065 (B-band, Subaru/Suprime-Cam), Rp/Rs = 0.11601 \pm 0.00117 (B-band, Subaru/FOCAS), Rp/Rs = 0.11654 \pm 0.00080 (J-band, IRSF/SIRIUS), Rp/Rs = 0.11550 ^{+0.00142}_{-0.00153} (H-band, IRSF/SIRIUS), and Rp/Rs = 0.11547 \pm 0.00127 (Ks-band, IRSF/SIRIUS). The Subaru Suprime-Cam transit photometry shows a possible spot-crossing feature. Comparisons of the new transit depths and those from previous studies with the theoretical models by Howe & Burrows (2012) suggest that the high molecular weight atmosphere (e.g., 1% H$_2$O + 99% N$_2$) models are most likely, however, the low molecular weight (hydrogen dominated) atmospheres with extensive clouds are still not excluded. We also report a long-term monitoring of the stellar brightness variability of GJ 1214 observed with the MITSuME 50cm telescope in g'-, Rc-, and Ic-bands simultaneously. The monitoring was conducted for 32 nights spanning 78 nights in 2012, and we find a periodic brightness variation with a period of Ps = 44.3 \pm 1.2 days and semi-amplitudes of 2.1% \pm 0.4% in g'-band, 0.56% \pm 0.08% in Rc-band, and 0.32% \pm 0.04% in Ic-band.

We present multi-color light curves and optical spectra of V1280 Scorpii obtained from 2007 to 2012. It is shown that V1280 Sco is the extremely slow nova and the mass of white dwarf appears to be $\sim$ 0.6 M$\odot$ or lower. Blue-shifted multiple absorption lines of Na {\sc i} D, Ca {\sc ii} HK, and He {\sc i*} are detected on high-resolution spectra. We also discuss that an approach using metastable He absorption lines is useful to investigate structures of nova shells.

CTA 102, classified as a flat spectrum radio quasar at z=1.037, produced exceptionally bright optical flare in 2012 September. Following Fermi-LAT detection of enhanced gamma-ray activity, we densely monitored this source in the optical and near-infrared bands for the subsequent ten nights using twelve telescopes in Japan and South-Africa. On MJD 56197 (2012 September 27, 4-5 days after the peak of bright gamma-ray flare), polarized flux showed a transient increase, while total flux and polarization angle remained almost constant during the "orphan polarized-flux flare". We also detected an intra-night and prominent flare on MJD 56202. The total and polarized fluxes showed quite similar temporal variations, but PA again remained constant during the flare. Interestingly, the polarization angles during the two flares were significantly different from the jet direction. Emergence of a new emission component with high polarization degree (PD) up to 40% would be responsible for the observed two flares, and such a high PD indicates a presence of highly ordered magnetic field at the emission site. We discuss that the well-ordered magnetic field and even the observed directions of polarization angle which is grossly perpendicular to the jet are reasonably accounted for by transverse shock(s) propagating down the jet.

We present the result of near-infrared (near-IR) [Fe II] line mapping of the supernova remnant IC443 with the IRSF/SIRIUS, using the two narrow-band filters tuned for the [Fe II] 1.257 micron and [Fe II] 1.644 micron lines. Covering a large area of 30' x 35', our observations reveal that [Fe II] filamentary structures exist all over the remnant, not only in an ionic shock shell, but also in a molecular shock shell and a central region inside the shells. With the two [Fe II] lines, we performed corrections for dust extinction to derive the intrinsic line intensities. We also obtained the intensities of thermal emission from the warm dust associated with IC443, using the far- and mid-IR images taken with AKARI and Spitzer, respectively. As a result, we find that the [Fe II] line emission relative to the dust emission notably enhances in the inner central region. We discuss causes of the enhanced [Fe II] line emission, estimating the Fe+ and dust masses.

We present results from deep imaging linear and circular polarimetry of the massive star-forming region NGC 6334- V. These observations show high degrees of circular polarization (CP), as much as 22 % in the Ks band, in the infrared nebula associated with the outflow. The CP has an asymmetric positive/negative pattern and is very extended (~80" or 0.65 pc). Both the high CP and its extended size are larger than those seen in the Orion CP region. Three-dimensional Monte Carlo light-scattering models are used to show that the high CP may be produced by scattering from the infrared nebula followed by dichroic extinction by an optically thick foreground cloud containing aligned dust grains. Our results show not only the magnetic field orientation of around young stellar objects but also the structure of circumstellar matter such as outflow regions and their parent molecular cloud along the line of sight. The detection of the large and extended CP in this source and the Orion nebula may imply the CP origin of the biological homochirality on Earth.

We present an analysis of the anomalous microlensing event, MOA-2010-BLG-073, announced by the Microlensing Observations in Astrophysics survey on 2010-03-18. This event was remarkable because the source was previously known to be photometrically variable. Analyzing the pre-event source lightcurve, we demonstrate that it is an irregular variable over time scales >200d. Its dereddened color, $(V-I)_{S,0}$, is 1.221$\pm$0.051mag and from our lens model we derive a source radius of 14.7$\pm$1.3 $R_{\odot}$, suggesting that it is a red giant star. We initially explored a number of purely microlensing models for the event but found a residual gradient in the data taken prior to and after the event. This is likely to be due to the variability of the source rather than part of the lensing event, so we incorporated a slope parameter in our model in order to derive the true parameters of the lensing system. We find that the lensing system has a mass ratio of q=0.0654$\pm$0.0006. The Einstein crossing time of the event, $T_{\rm{E}}=44.3$\pm$0.1d, was sufficiently long that the lightcurve exhibited parallax effects. In addition, the source trajectory relative to the large caustic structure allowed the orbital motion of the lens system to be detected. Combining the parallax with the Einstein radius, we were able to derive the distance to the lens, $D_L$=2.8$\pm$0.4kpc, and the masses of the lensing objects. The primary of the lens is an M-dwarf with $M_{L,p}$=0.16$\pm0.03M_{\odot}$ while the companion has $M_{L,s}$=11.0$\pm2.0M_{\rm{J}}$ putting it in the boundary zone between planets and brown dwarfs.

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 report the result of our near-infrared survey of short-period variable stars (P<60d) in a field-of-view of 20'x30' towards the Galactic Centre. Forty-five variables are discovered and we classify the variables based on their light curve shapes and other evidence. In addition to 3 classical Cepheids reported previously, we find 16 type II Cepheids, 24 eclipsing binaries, one pulsating star with P=0.265d (RR Lyr or delta Sct) and one Cepheid-like variable whose nature is uncertain. Eclipsing binaries are separated into the foreground objects and those significantly obscured by interstellar extinction. One of the reddened binaries contains an O-type supergiant and its light curve indicates an eccentric orbit. We discuss the nature and distribution of type II Cepheids as well as the distance to the Galactic Centre based on these Cepheids and other distance indicators. The estimates of R0(GC) we obtained based on photometric data agree with previous results obtained with kinematics of objects around the GC. Furthermore, our result gives a support to the reddening law obtained by Nishiyama and collaborators, A(Ks)/E(H-Ks)=1.44, because a different reddening law would result in a rather different distance estimate.

In this Letter, we report detections of SiO v=3 J=1--0 maser emission in very long baseline interferometric (VLBI) observations towards 4 out of 12 long-period variable stars: WX Psc, R Leo, W Hya, and T Cep. The detections towards WX Psc and T Cep are new ones. We also present successful astrometric observations of SiO v=2 and v=3 J=1--0 maser emissions associated with two stars: WX Psc and W Hya and their position-reference continuum sources: J010746.0+131205 and J135146.8-291218 with the VLBI Exploration of Radio Astrometry (VERA). The relative coordinates of the position-reference continuum source and SiO v=3 maser spots were measured with respect to those of an SiO v=2 maser spot adopted as fringe-phase reference. Thus the faint continuum sources were inversely phase-referenced to the bright maser sources. It implies possible registration of multiple SiO maser line maps onto a common coordinate system with 10 microarcsecond-level accuracy.

We report high precision transit photometry of GJ1214b in JHKs bands taken simultaneously with the SIRIUS camera on the IRSF 1.4m telescope at Sutherland, South Africa. Our MCMC analyses show that the observed planet-to-star radius ratios in JHKs bands are R_{\rm p}/R_{\rm s,J} = 0.11833 \pm 0.00077, R_{\rm p}/R_{\rm s,H} = 0.11522 \pm 0.00079, R_{\rm p}/R_{\rm s,Ks} = 0.11459 \pm 0.00099, respectively. The radius ratios are well consistent with the previous studies by Bean et al. (2011) within 1\sigma, while our result in Ks band is shallower than and inconsistent at 4\sigma\ level with the previous measurements in the same band by Croll et al. (2011). We have no good explanation for this discrepancy at this point. Our overall results support a flat transmission spectrum in the observed bands, which can be explained by a water-dominated atmosphere or an atmosphere with extensive high-altitude clouds or haze. To solve the discrepancy of the radius ratios and to discriminate a definitive atmosphere model for GJ1214b in the future, further transit observations around Ks band would be especially important.

We report the extremely high magnification (A > 1000) binary microlensing event OGLE-2007-BLG-514. We obtained good coverage around the double peak structure in the light curve via follow-up observations from different observatories. The binary lens model that includes the effects of parallax (known orbital motion of the Earth) and orbital motion of the lens yields a binary lens mass ratio of q = 0.321 +/- 0.007 and a projected separation of s = 0.072 +/- 0.001$ in units of the Einstein radius. The parallax parameters allow us to determine the lens distance D_L = 3.11 +/- 0.39 kpc and total mass M_L=1.40 +/- 0.18 M_sun; this leads to the primary and secondary components having masses of M_1 = 1.06 +/- 0.13 M_sun and M_2 = 0.34 +/- 0.04 M_sun, respectively. The parallax model indicates that the binary lens system is likely constructed by the main sequence stars. On the other hand, we used a Bayesian analysis to estimate probability distributions by the model that includes the effects of xallarap (possible orbital motion of the source around a companion) and parallax (q = 0.270 +/- 0.005, s = 0.083 +/- 0.001). The primary component of the binary lens is relatively massive with M_1 = 0.9_{-0.3}^{+4.6} M_sun and it is at a distance of D_L = 2.6_{-0.9}^{+3.8} kpc. Given the secure mass ratio measurement, the companion mass is therefore M_2 = 0.2_{-0.1}^{+1.2} M_sun. The xallarap model implies that the primary lens is likely a stellar remnant, such as a white dwarf, a neutron star or a black hole.

Microlensing detections of cool planets are important for the construction of an unbiased sample to estimate the frequency of planets beyond the snow line, which is where giant planets are thought to form according to the core accretion theory of planet formation. In this paper, we report the discovery of a giant planet detected from the analysis of the light curve of a high-magnification microlensing event MOA-2010-BLG-477. The measured planet-star mass ratio is $q=(2.181\pm0.004)\times 10^{-3}$ and the projected separation is $s=1.1228\pm0.0006$ in units of the Einstein radius. The angular Einstein radius is unusually large $\theta_{\rm E}=1.38\pm 0.11$ mas. Combining this measurement with constraints on the "microlens parallax" and the lens flux, we can only limit the host mass to the range $0.13<M/M_\odot<1.0$. In this particular case, the strong degeneracy between microlensing parallax and planet orbital motion prevents us from measuring more accurate host and planet masses. However, we find that adding Bayesian priors from two effects (Galactic model and Keplerian orbit) each independently favors the upper end of this mass range, yielding star and planet masses of $M_*=0.67^{+0.33}_{-0.13}\ M_\odot$ and $m_p=1.5^{+0.8}_{-0.3}\ M_{\rm JUP}$ at a distance of $D=2.3\pm0.6$ kpc, and with a semi-major axis of $a=2^{+3}_{-1}$ AU. Finally, we show that the lens mass can be determined from future high-resolution near-IR adaptive optics observations independently from two effects, photometric and astrometric.

We observed the microlensing event MOA-2007-BLG-192 at high angular resolution in JHKs with the NACO adaptive optics system on the VLT while the object was still amplified by a factor 1.23 and then at baseline 18 months later. We analyzed and calibrated the NACO photometry in the standard 2MASS system in order to accurately constrain the source and the lens star fluxes. We detect light from the host star of MOA-2007-BLG-192Lb, which significantly reduces the uncertainties in its char- acteristics as compared to earlier analyses. We find that MOA-2007-BLG-192L is most likely a very low mass late type M-dwarf (0.084 [+0.015] [-0.012] M\odot) at a distance of 660 [+100] [-70] pc orbited by a 3.2 [+5.2] [-1.8] M\oplus super-Earth at 0.66 [+0.51] [-0.22] AU. We then discuss the properties of this cold planetary system.

We present optical ($B$, $V$, $R_{\rm c}$, $I_{\rm c}$ and $y$) and near infrared ($J$, $H$ and $K_{\rm s}$) photometric and spectroscopic observations of a classical nova V1280 Scorpii for five years from 2007 to 2011. Our photometric observations show a declining event in optical bands shortly after the maximum light which continues $\sim$ 250 days. The event is most probably caused by a dust formation. The event is accompanied by a short ($\sim$ 30 days) re-brightening episode ($\sim$ 2.5 mag in $V$), which suggests a re-ignition of the surface nuclear burning. After 2008, the $y$ band observations show a very long plateau at around $y$ = 10.5 for more than 1000 days until April 2011 ($\sim$ 1500 days after the maximum light). The nova had taken a very long time ($\sim$ 50 months) before entering the nebular phase (clear detection of both [\ion{O}{iii}] 4959 and 5007) and is still continuing to generate the wind caused by H-burning. The finding suggests that V1280 Sco is going through the historically slowest evolution. The interval from the maximum light (2007 February 16) to the beginning of the nebular phase is longer than any previously known slow novae: V723 Cas (18 months), RR Pic (10 months), or HR Del (8 months). It suggests that the mass of a white dwarf in the V1280 Sco system might be 0.6 $M_\mathrm{\sun}$ or smaller. The distance, based on our measurements of the expansion velocity combined with the directly measured size of the dust shell, is estimated to be 1.1 $\pm$ 0.5 kpc.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

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.

The nuclear bulge is a region with a radius of about 200 parsecs around the centre of the Milky Way. It contains stars with ages ranging from a few million years to over a billion years, yet its star-formation history and the triggering process for star formation remain to be resolved. Recently, episodic star formation, powered by changes in the gas content, has been suggested. Classical Cepheid variable stars have pulsation periods that decrease with increasing age, so it is possible to probe the star-formation history on the basis of the distribution of their periods. Here we report the presence of three classical Cepheids in the nuclear bulge with pulsation periods of approximately 20 days, within 40 parsecs (projected distance) of the central black hole. No Cepheids with longer or shorter periods were found. We infer that there was a period about 25 million years ago, and possibly lasting until recently, in which star formation increased relative to the period of 30-70 million years ago.

Near-infrared (IR) imaging polarimetry in the J, H, and Ks bands has been carried out for the protostellar cluster region around NGC 2264 IRS 2 in the Monoceros OB1 molecular cloud. Various infrared reflection nebulae clusters (IRNCs) associated with NGC 2264 IRS 2 and IRAS 12 S1 core were detected as well as local infrared reflection nebulae (IRNe). The illuminating sources of the IRNe were identified with known or new near- and mid-IR sources. In addition, 314 point-like sources were detected in all three bands and their aperture polarimetry was studied. Using a color-color diagram, reddened field stars and diskless pre-main sequence stars were selected to trace the magnetic field (MF) structure of the molecular cloud. The mean polarization position angle of the point-like sources is 81 \pm 29 degree in the cluster core, and 58 \pm 24 degree in the perimeter of the cluster core, which is interpreted as the projected direction on the sky of the MF in the observed region of the cloud. The Chandrasekhar-Fermi method gives a rough estimate of the MF strength to be about 100 {\mu}G. A comparison with recent numerical simulations of the cluster formation implies that the cloud dynamics is controlled by the relatively strong MF. The local MF direction is well associated with that of CO outflow for IRAS 12 S1 and consistent with that inferred from submillimeter polarimetry. In contrast, the local MF direction runs roughly perpendicular to the Galactic MF direction.

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.

We report the discovery of two field brown dwarfs, ULAS J0128-0041 and ULAS J0321+0051, and the rediscovery of ULAS J0226+0051 (IfA 0230-Z1), in the Sloan Digital Sky Survey (SDSS) southern equatorial stripe. They are found in the course of our follow-up observation program of 1 um excess sources in the United Kingdom Infrared Telescope Infrared Deep Sky Survey. The Gemini Multi-Object Spectrographs spectra at red optical wavelengths (6500-10500 A) are presented, which reveal that they are early-T dwarfs. The classification is also supported by their optical to near-infrared colors. It is noted that ULAS J0321+0051 is one of the faintest currently known T dwarfs. The estimated distances to the three objects are 50-110 pc, thus they are among the most distant field T dwarfs known. Dense temporal coverage of the target fields achieved by the SDSS-II Supernova Survey allows us to perform a simple time-series analysis, which leads to the finding of significant proper motions of 150-290 mas/yr or the transverse velocities of 40-100 km/s for ULAS J0128-0041 and ULAS J0226+0051. We also find that there are no detectable, long-term (a-few-year) brightness variations above a few times 0.1 mag for the two brown dwarfs.

This paper presents the analysis of a deep near-infrared J,H,Ks-imaging survey (37.5 sq deg) aimed at tracing the galaxy distribution of the Great Attractor (GA) in the Zone of Avoidance along the so-called Norma Wall. The resulting galaxy catalog is complete to extinction-corrected magnitudes Ks^o = 14.8 mag for extinctions less than A_K = 1.0 mag and star densities below log N(Ks<14.0) < 4.72. Of the 4360 cataloged galaxies, 99.2% lie in the hereby constrained 89.5% of the survey area. Although the analyzed galaxy distribution reveals no new major galaxy clusters at the GA distance (albeit some more distant ones), the overall number counts and luminosity density indicate a clear and surprisingly smooth overdensity at the GA distance that extends over the whole surveyed region. A mass estimate of the Norma Wall overdensity derived from (a) galaxy number counts and (b) photometric redshift distribution gives a lower value compared to the original prediction by Lynden-Bell et al. 1988 (~14%), but is consistent with more recent independent assessments.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

XSS J12270-4859 (J12270) is an enigmatic source of unknown nature. Previous studies revealed that the source has unusual X-ray temporal characteristics, including repetitive short-term flares followed by spectral hardening, non-periodic dips, and dichotomy in activity; i.e. intervals filled with flares and those without. Together with a power-law X-ray spectrum, it is suggested to be a low-mass X-ray binary (LMXB). In order to better understand the object, we present the results of our near-infrared (NIR) photometry and linear polarimetry observations as well as X-ray spectroscopy observations, which overlap with each other partially in time, taken respectively with the InfraRed Survey Facility (IRSF) and the Rossi X-ray Timing Explorer (RXTE). We detected several simultaneous NIR and X-ray flares for the first time. No significant NIR polarization was obtained. We assembled data taken with IRSF, RXTE, Suzaku, Swift, and other missions in the literature and compared the flare profile and the spectral energy distribution (SED) with some representative high-energy sources. Based on some similarities of the repetitive NIR and X-ray flaring characteristics and the broad SED, we argue that J12270 is reminiscent of microquasars with a synchrotron jet, which is at a very low luminosity state of ~1e-4 Eddington luminosity for a stellar mass black hole or neutron star at a reference distance of 1 kpc.

X-ray and near-infrared ($J$-$H$-$K_{\rm s}$) observations of the Galactic black hole binary GX 339--4 in the low/hard state were performed with Suzaku and IRSF in 2009 March. The spectrum in the 0.5--300 keV band is dominated by thermal Comptonization of multicolor disk photons, with a small contribution from a direct disk component, indicating that the inner disk is almost fully covered by hot corona with an electron temperature of $\approx$175 keV. The Comptonizing corona has at least two optical depths, $\tau \approx 1,0.4$. Analysis of the iron-K line profile yields an inner disk radius of $(13.3^{+6.4}_{-6.0}) R_{\rm g}$ ($R_{\rm g} $ represents the gravitational radius $GM/c^2$), with the best-fit inclination angle of $\approx50^\circ$. This radius is consistent with that estimated from the continuum fit by assuming the conservation of photon numbers in Comptonization. Our results suggest that the standard disk of GX 339--4 is likely truncated before reaching the innermost stable circular orbit (for a non rotating black hole) in the low/hard state at $\sim$1% of the Eddington luminosity. The one-day averaged near-infrared light curves are found to be correlated with hard X-ray flux with $F_{\rm Ks} \propto F_{\rm X}^{0.45}$. The flatter near infrared $\nu F_{\nu}$ spectrum than the radio one suggests that the optically thin synchrotron radiation from the compact jets dominates the near-infrared flux. Based on a simple analysis, we estimate the magnetic field and size of the jet base to be $5\times10^4$ G and $6\times 10^8$ cm, respectively. The synchrotron self Compton component is estimated to be approximately 0.4% of the total X-ray flux.

We propose a method to determine the galactic constants R_0 (distance to the Galactic Center) and V_0 (rotation velocity of the Sun) from measurements of distances, radial velocities and proper motions of objects near the solar circle. This is a modification of the solar-circle method to a more practical observational method. We apply the method to determine R_0 using data from the literature with known distances and radial velocities, and obtain R_0 = 7.54 +/- 0.77 kpc.

The Serpens South embedded cluster, which is located at the constricted part in a long filamentary infrared dark cloud, is believed to be in very early stage of cluster formation. We present results of near-infrared (JHKs) polarization observations toward the filamentary cloud. Our polarization measurements of near-infrared point sources indicate a well-ordered global magnetic field that is perpendicular to the main filament, implying that the magnetic field is likely to have controlled the formation of the main filament. On the other hand, the sub-filaments, which converge on the central part of the cluster, tend to run along the magnetic field. The global magnetic field appears to be curved in the southern part of the main filament. Such morphology is consistent with the idea that the global magnetic field is distorted by gravitational contraction along the main filament toward the northern part that contains larger mass. Applying the Chandrasekhar-Fermi method, the magnetic field strength is roughly estimated to be a few x 100 microgauss, suggesting that the filamentary cloud is close to magnetically critical as a whole.

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 first example of binary microlensing for which the parameter measurements can be verified (or contradicted) by future Doppler observations. This test is made possible by a confluence of two relatively unusual circumstances. First, the binary lens is bright enough (I=15.6) to permit Doppler measurements. Second, we measure not only the usual 7 binary-lens parameters, but also the 'microlens parallax' (which yields the binary mass) and two components of the instantaneous orbital velocity. Thus we measure, effectively, 6 'Kepler+1' parameters (two instantaneous positions, two instantaneous velocities, the binary total mass, and the mass ratio). Since Doppler observations of the brighter binary component determine 5 Kepler parameters (period, velocity amplitude, eccentricity, phase, and position of periapsis), while the same spectroscopy yields the mass of the primary, the combined Doppler + microlensing observations would be overconstrained by 6 + (5 + 1) - (7 + 1) = 4 degrees of freedom. This makes possible an extremely strong test of the microlensing solution. We also introduce a uniform microlensing notation for single and binary lenses, we define conventions, summarize all known microlensing degeneracies and extend a set of parameters to describe full Keplerian motion of the binary lenses.

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 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 report the gravitational microlensing discovery of a sub-Saturn mass planet, MOA-2009-BLG-319Lb, orbiting a K or M-dwarf star in the inner Galactic disk or Galactic bulge. The high cadence observations of the MOA-II survey discovered this microlensing event and enabled its identification as a high magnification event approximately 24 hours prior to peak magnification. As a result, the planetary signal at the peak of this light curve was observed by 20 different telescopes, which is the largest number of telescopes to contribute to a planetary discovery to date. The microlensing model for this event indicates a planet-star mass ratio of q = (3.95 +/- 0.02) x 10^{-4} and a separation of d = 0.97537 +/- 0.00007 in units of the Einstein radius. A Bayesian analysis based on the measured Einstein radius crossing time, t_E, and angular Einstein radius, \theta_E, along with a standard Galactic model indicates a host star mass of M_L = 0.38^{+0.34}_{-0.18} M_{Sun} and a planet mass of M_p = 50^{+44}_{-24} M_{Earth}, which is half the mass of Saturn. This analysis also yields a planet-star three-dimensional separation of a = 2.4^{+1.2}_{-0.6} AU and a distance to the planetary system of D_L = 6.1^{+1.1}_{-1.2} kpc. This separation is ~ 2 times the distance of the snow line, a separation similar to most of the other planets discovered by microlensing.

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 present a large-scale view of the magnetic field in the central 2deg * 2deg region of our Galaxy. The polarization of point sources has been measured in the J, H, and Ks bands using the near-infrared polarimetric camera SIRPOL on the 1.4 m telescope IRSF. Comparing the Stokes parameters between high extinction stars and relatively low extinction ones, we obtain polarization originating from magnetically aligned dust grains in the central few-hundred pc of our Galaxy. We find that near the Galactic plane, the magnetic field is almost parallel to the Galactic plane (i.e., toroidal configuration) but at high Galactic latitudes (| b | > 0.4deg), the field is nearly perpendicular to the plane (i.e., poloidal configuration). This is the first detection of a smooth transition of the large-scale magnetic field configuration in this region.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

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

From multi-epoch adaptive optics imaging and integral field unit spectroscopy we report the discovery of an expanding and narrowly confined bipolar shell surrounding the helium nova V445 Puppis (Nova Puppis 2000). An equatorial dust disc obscures the nova remnant, and the outflow is characterised by a large polar outflow velocity of 6720 +/- 650 km/s and knots moving at even larger velocities of 8450 +/- 570 km/s. We derive an expansion parallax distance of 8.2 +/- 0.5 kpc and deduce a pre-outburst luminosity of the underlying binary of log L/L_Sun = 4.34 +/- 0.36. The derived luminosity suggests that V445 Puppis probably contains a massive white dwarf accreting at high rate from a helium star companion making it part of a population of binary stars that potentially lead to supernova Ia explosions due to accumulation of helium-rich material on the surface of a massive white dwarf.

We report the results of a near-infrared survey for long-period variables in a field of view of 20 arcmin by 30 arcmin towards the Galactic Centre (GC). We have detected 1364 variables, of which 348 are identified with those reported in Glass et al. (2001). We present a catalogue and photometric measurements for the detected variables and discuss their nature. We also establish a method for the simultaneous estimation of distances and extinctions using the period-luminosity relations for the JHKs bands. Our method is applicable to Miras with periods in the range 100--350 d and mean magnitudes available in two or more filter bands. While J-band means are often unavailable for our objects because of the large extinction, we estimated distances and extinctions for 143 Miras whose H- and Ks-band mean magnitudes are obtained. We find that most are located at the same distance to within our accuracy. Assuming that the barycentre of these Miras corresponds to the GC, we estimate its distance modulus to be 14.58+-0.02(stat.)+-0.11(syst.) mag, corresponding to 8.24+-0.08(stat.)+-0.42(syst.) kpc. We have assumed the distance modulus to the LMC to be 18.45 mag, and the uncertainty in this quantity is included in the systematic error above. We also discuss the large and highly variable extinction. Its value ranges from 1.5 mag to larger than 4 mag in A(Ks) except towards the thicker dark nebulae and it varies in a complicated way with the line of sight. We have identified mid-infrared counterparts in the Spitzer/IRAC catalogue of Ramirez et al. (2008) for most of our variables and find that they follow rather narrow period-luminosity relations in the 3.6 to 8.0 micrometre wavelength range.

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.

The detection of GeV photons from gamma-ray bursts (GRBs) has important consequences for the interpretation and modelling of these most-energetic cosmological explosions. The full exploitation of the high-energy measurements relies, however, on the accurate knowledge of the distance to the events. Here we report on the discovery of the afterglow and subsequent redshift determination of GRB 080916C, the first GRB detected by the Fermi Gamma-Ray Space Telescope with high significance detection of photons at >0.1 GeV. Observations were done with 7-channel imager GROND at the 2.2m MPI/ESO telescope, the SIRIUS instrument at the Nagoya-SAAO 1.4m telescope in South Africa, and the GMOS instrument at Gemini-S. The afterglow photometric redshift of z=4.35+-0.15, based on simultaneous 7-filter observations with the Gamma-Ray Optical and Near-infrared Detector (GROND), places GRB 080916C among the top 5% most distant GRBs, and makes it the most energetic GRB known to date. The detection of GeV photons from such a distant event is rather surprising. The observed gamma-ray variability in the prompt emission together with the redshift suggests a lower limit for the Lorentz factor of the ultra-relativistic ejecta of Gamma > 1090. This value rivals any previous measurements of Gamma in GRBs and strengthens the extreme nature of GRB 080916C.

In the J, H, and Ks bands survey of the the Galactic Center region over an area of 2deg x 5deg, we have found many dark clouds, among which a distinguished chain of dark clouds can be identified with a quiescent CO cloud. The distances of the clouds is estimated to be 3.2-4.2 kpc, corresponding to the Norma arm by our new method to determine distance to dark clouds using the cumulative number of stars against J-Ks colors. Adopting these estimated distances, the size is about 70 pc in length and the total mass of the cloud is 6x10^4 M_solar. Three compact HII regions harbor in the cloud, indicating that star forming activities are going on at the cores of the quiescent CO cloud on the spiral arm.

We present the results of wide-field JHKs polarimetry toward the HII region S106 using the IRSF (Infrared Survey Facility) telescope. Our polarimetry data revealed an extended (up to ~ 5') polarized nebula over S106. We confirmed the position of the illuminating source of most of the nebula as consistent with S106 IRS4 through an analysis of polarization vectors. The bright portion of the polarized intensity is consistent with the red wing component of the molecular gas. Diffuse polarized intensity emission is distributed along the north--south molecular gas lanes. We found the interaction region between the radiation from S106 IRS4 and the dense gas. In addition, we also discovered two small polarization nebulae, SIRN1 and SIRN2, associated with a young stellar objects (YSO). Aperture polarimetry of point-like sources in this region was carried out for the first time. The regional magnetic field structures were derived using point-like source aperture polarimetry, and the magnetic field structure position angle around the cluster region in S106 was found to be ~ 120$\arcdeg$. The magnetic fields in the cluster region, however, have three type position angles: ~ 20$\arcdeg$, ~ 80$\arcdeg$, and ~ 120$\arcdeg$. The present magnetic field structures are consistent with results obtained by submillimeter continuum observations. We found that the magnetic field direction in the dense gas region is not consistent with that of the low density gas region.

We present a polarimetric map of a 20'x20' area toward the Galactic center. The polarization of point sources has been measured in the J, H, and Ks bands using the near-infrared polarimetric camera SIRPOL on the 1.4 m telescope IRSF. One percent or better accuracy of polarization degree is achieved for sources with J<14.5, H<13.5, and Ks<12.0. Comparing the Stokes parameters between high extinction stars and relatively low extinction ones, we have obtained a polarization originating from magnetically aligned dust grains at the central region of our Galaxy of at most 1-2 kpc. The distribution of the position angles shows a peak at about 20 deg, nearly parallel to the Galactic plane, suggesting a toroidal magnetic configuration. The derived direction of the magnetic field is in good agreement with that obtained from far-infrared/submillimeter observations, which detect polarized thermal emission from dust in the molecular clouds at the Galactic center. Our results show that by subtracting foreground components, near-infrared polarimetry allows investigation of the magnetic field structure "at" the Galactic center.

HIZOA J0836-43 is an extreme gas-rich ($M_{\rm{HI}}$=7.5$\times10^{10} M_{\sun}$) disk galaxy which lies hidden behind the strongly obscuring Vela region of the Milky Way. Utilizing observations from the {\it Spitzer Space Telescope}, we have found it to be a luminous infrared starburst galaxy with a star formation rate of $\sim 21 M_{\sun} \rm{yr^{-1}}$, arising from exceptionally strong molecular PAH emission ($L_{7.7\micron} = 1.50 \times 10^{9} L_{\odot}$) and far-infrared emission from cold dust. The galaxy exhibits a weak mid-infrared continuum compared to other starforming galaxies and U/LIRGs. This relative lack of emission from small grains suggests atypical interstellar medium conditions compared to other starbursts. We do not detect significant $[$Ne {\sc v}$]$ or $[$O {\sc iv}$]$, which implies an absent or very weak AGN. The galaxy possesses a prominent bulge of evolved stars and a stellar mass of 4.4($\pm$1.4)$\times10^{10} M_{\sun}$. With its plentiful gas supply and current star formation rate, a doubling of stellar mass would occur on a timescale of $\sim$2 Gyr. Compared to local galaxies, HIZOA J0836-43 appears to be a "scaled-up" spiral undergoing inside-out formation, possibly resembling stellar disk building processes at intermediate redshifts.

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 conducted aperture polarimetry of ~500 stars of the Orion Nebula Cluster (ONC) in M42 based on our wide-field (~8'\times 8') $JHKs$ band polarimetry. Most of the near-infrared (NIR) polarizations are dichroic, with position angles of polarization agreeing, both globally and locally, with previous far-infrared (FIR) and submillimeter observations, having taken into account the 90$^\circ $ difference in angles between dichroic absorption and emission. This is consistent with the idea that both NIR dichroic polarizations and FIR/submillimeter thermal polarizations trace the magnetic fields in the OMC-1 region. The magnetic fields inferred from these observations show a pinch at scales less than 0.5 pc with a centroid near IRc2. The hourglass-shaped magnetic field pattern is explained by the models in which the magnetic field lines are dragged along with the contracting gas and then wound up by rotation in a disk. The highly polarized region to the northwest of IRc2 and the low-polarized region near the bright bar are also common among NIR and FIR/submillimeter data, although a few regions of discrepancy exist. We have also discerned ~50 possible highly polarized sources whose polarizations are more likely to be intrinsic rather than dichroic. Their polarization efficiencies ($P(H)/A(H)$) are too large to be explained by the interstellar polarization. These include 10 young brown dwarfs that suggest a higher polarization efficiency, which may present geometrical evidence for (unresolved) circumstellar structures around young brown dwarfs.

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 optical to near-infrared (IR) spectrum of the galaxy TSPS J1329-0957, a red and bright member of the class of extremely red objects (EROs) at z = 1.26. This galaxy was found in the course of the Tokyo-Stromlo Photometry Survey (TSPS) which we are conducting in the southern sky. The spectroscopic observations were carried out with the Gemini Multi-Object Spectrograph (GMOS) and the Gemini Near Infra-Red Spectrograph (GNIRS) mounted on the Gemini-South telescope. The wide wavelength coverage of 0.6 - 2.3 um provides useful clues as to the nature of EROs while most published spectra are limited to a narrower spectral range which is dictated by the need for efficient redshift determination in a large survey. We compare our spectrum with several optical composite spectra obtained in recent large surveys, and with stellar population synthesis models. The effectiveness of using near-IR broad-band data, instead of the spectral data, in deriving the galaxy properties are also investigated. We find that TSPS J1329-0957 formed when the universe was 2 - 3 Gyr old, and subsequently evolved passively to become one of the most massive galaxies found in the z = 1 - 2 universe. Its early type and estimated stellar mass of M* = 10^{11.5} Msun clearly point to this galaxy being a direct ancestor of the brightest elliptical and spheroidal galaxies in the local universe.

We have carried out near-IR/optical observations to examine star formation toward a bright-rimmed cometary globule (BRC37) facing the exciting star(s) of an HII region (IC1396) containing an IRAS source, which is considered to be an intermediate-mass protostar. With slit-less spectroscopy we detected ten H_alpha emission stars around the globule, six of which are near the tip of the globule and are aligned along the direction to the exciting stars. There is evidence that this alignment was originally towards an O9.5 star, but has evolved to align towards a younger O6 star when that formed. Near-IR and optical photometry suggests that four of these six stars are low-mass young stellar objects (YSOs) with masses of ~0.4 M_sun. Their estimated ages of ~1 Myr indicate that they were formed at the tip in advance of the formation of the IRAS source. Therefore, it is likely that sequential star formation has been taking place along the direction from the exciting stars towards the IRAS source, due to the UV impact of the exciting star(s). Interestingly, one faint, H_alpha emission star, which is the closest to the exciting star(s), seems to be a young brown dwarf that was formed by the UV impact in advance of the formation of other YSOs at the tip.

We report the detection of eighteen infrared reflection nebulae (IRNe) in the $J$, $H$, & $Ks$ linear polarimetric observations of the NGC 6334 massive star-formation complex, of which 16 IRNe are new discoveries. Our images cover $\sim$180 square arcminutes, one of the widest near-infrared polarization data in star-formation regions so far. These IRNe are most likely associated with embedded young OB stars at different evolutionary phases, showing a variety of sizes, morphologies, and polarization properties, which can be divided into four categories. We argue the different nebula characteristics to be a possible evolutionary sequence of circumstellar structures around young massive stars.

We present proper motions of water masers toward the Onsala 1 star forming region, observed with the Japanese VLBI network at three epochs spanning 290 days. We found that there are two water maser clusters (WMC1 and WMC2) separated from each other by 1".6 (2900 AU at a distance of 1.8 kpc). The proper motion measurement reveals that WMC1 is associated with a bipolar outflow elongated in the east-west direction with an expansion velocity of 69+-11 km/s. WMC1 and WMC2 are associated with two 345 GHz continuum dust emission sources, and located 2" (3600 AU) east from the core of an ultracompact HII region traced by 8.4 GHz radio continuum emission. This indicates that star formation activity of Onsala 1 could move from the west side of ultracompact HII region to the east side of two young stellar objects associated with the water masers. We also find that WMC1 and UC HII region could be gravitationally bound. Their relative velocity along the line of sight is 3 km/s, and total mass is 37 Mo. Onsala 1 seems to harbor a binary star at different evolutionary stage.

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.

We present a map of the major part of the central molecular zone (CMZ) of simultaneous observations in the NH3 (J,K) = (1,1) and (2,2) lines using the Kagoshima 6-m telescope. The mapped area is -1.000 < l < 1.625 deg, -0.375 < b < +0.250 deg. The kinetic temperatures derived from the (2,2) to (1,1) intensity ratios are 20--80 K or exceed 80 K. The gases corresponding to temperature of 20--80 K and > 80 K contain 75% and 25% of the total NH3 flux, respectively. These temperatures indicate that the dense molecular gas in the CMZ is dominated by gas that is warmer than the majority of the dust present there. A comparison with the CO survey by Sawada et al. (2001) shows that the NH3 emitting region is surrounded by a high pressure region on the l-v plane. Although NH3 emission traces dense gas, it is not extended into a high pressure region. Therefore, the high pressure region is less dense and has to be hotter. This indicates that the molecular cloud complex in the Galactic center region has a ``core'' of dense and warm clouds which are traced by the NH3 emission, and an ``envelope'' of less dense and hotter gas clouds. Besides heating by ambipolar diffusion, the hot plasma gas emitting the X-ray emission may heat the hot ``envelope''.

We carried out deep and wide (about 8' x 8') JHKs imaging polarimetry in the southern region of the Eagle Nebula (M16). The polarization intensity map reveals that two YSOs with near-IR reflection nebulae are located at the tips of two famous molecular pillars (Pillars 1 and 2) facing toward the exciting stars of M16. The centrosymmetric polarization pattern are consistent with those around class I objects having circumstellar envelopes, confirming that star formation is now taking place at the two tips of the pillars under the influence of UV radiation from the exciting stars. Polarization measurements of point sources show that magnetic fields are aligned along some of the pillars but in a direction that is quite different to the global structure in M16.

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 present the results of a high-resolution X-ray imaging study of the stellar population in the Galactic massive star-forming region RCW49 and its central OB association Westerlund 2. We obtained a 40 ks X-ray image of a 17'x17' field using the Chandra X-ray Observatory and deep NIR images using the Infrared Survey Facility in a concentric 8'3x8'3 region. We detected 468 X-ray sources and identified optical, NIR, and Spitzer Space Telescope MIR counterparts for 379 of them. The unprecedented spatial resolution and sensitivity of the X-ray image, enhanced by optical and infrared imaging data, yielded the following results: (1) The central OB association Westerlund 2 is resolved for the first time in the X-ray band. X-ray emission is detected from all spectroscopically-identified early-type stars in this region. (2) Most (86%) X-ray sources with optical or infrared identifications are cluster members in comparison with a control field in the Galactic Plane. (3) A loose constraint (2--5 kpc) for the distance to RCW49 is derived from the mean X-ray luminosity of T Tauri stars. (4) The cluster X-ray population consists of low-mass pre--main-sequence and early-type stars as obtained from X-ray and NIR photometry. About 30 new OB star candidates are identified. (5) We estimate a cluster radius of 6'--7' based on the X-ray surface number density profiles. (6) A large fraction (90%) of cluster members are identified individually using complimentary X-ray and MIR excess emission. (7) The brightest five X-ray sources, two Wolf-Rayet stars and three O stars, have hard thermal spectra.

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 conducted wide-field JHKs imaging polarimetry toward NGC 2024. We found a prominent and extended polarized nebula over NGC 2024, and constrained the location of illuminating source of the nebula through the analysis of polarization vectors. A massive star, IRS 2b, is located in the center of the symmetric vector pattern. Five small polarized nebulae associated with YSOs are discovered on our polarization images. These nebulae are responsible for the structures of circumstellar matter that produce strongly polarized light through dust scattering. For the point-like sources, we performed software aperture polarimetry in order to measure integrated polarizations, and found five young brown dwarfs with highly polarized integrated emission. These sources serve as direct evidence for the existence of disk/envelope system around brown dwarfs. We investigated the magnetic field structure of NGC 2024 through the measurements of dichroic polarization. The average position angle of projected magnetic fields across the region is found to be 110 degrees. We found a good consistency in magnetic field structures obtained using near-infrared dichroic polarization and sub-mm/far-infrared dust emission polarization, indicating that the dichroic polarizations at near-infrared wavelengths trace magnetic field structures inside dense molecular clouds.

We have conducted deep JHKs imaging polarimetry of a ~8' x 8' area of the NGC 2071 star forming region. Our polarization data have revealed various infrared reflection nebulae (IRNe) associated with the central IR young star cluster NGC2071IR and identified their illuminating sources. There are at least 4 IRNe in NGC2071IR and several additional IRNe are identified around nearby young stars in the same field-of-view. Each illuminating source coincides with a known near-IR source except for IRS3, which is only a part of IRN2 and is illuminated by the radio source 1c. Aperture polarimetry of each cluster source is used to detect unresolved circumstellar disk/outflow systems. Aperture polarimetry of the other point-like sources within the field is made in this region for the first time. The magnetic field structures (from ~1 pc down to \~0.1 pc) are derived using both aperture polarimetry of the point-like sources and imaging polarimetry of the shocked H2 emission that is seen as the dominant knotty nebulae in the Ks band image; they are both of dichroic origin and the derived field directions are consistent with each other. The magnetic field direction projected on the sky is also consistent with that inferred from the 850 micron thermal continuum emission polarimetry of the central 0.2 pc region, but running roughly perpendicular (~75 degrees) to the direction of the large scale outflow. We argue that the field strength is too weak to align the outflow in the large scale field direction via magnetic braking.

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

We present the polarization images in the $J$, $H$, & $Ks$ bands of the Orion Molecular Cloud 1 South region. The polarization images clearly show at least six infrared reflection nebulae (IRNe) which are barely seen or invisible in the intensity images. Our polarization vector images also identify the illuminating sources of the nebulae: IRN 1 & 2, IRN 3, 4, & 5, and IRN 6 are illuminated by three IR sources, Source 144-351, Source 145-356, and Source 136-355, respectively. Moreover, our polarization images suggest the candidate driving sources of the optical Herbig-Haro objects for the first time; HH529, a pair of HH202 and HH528 or HH 203/204, HH 530 and HH269 are originated from Source 144-351, Source 145-356, and Source 136-355, respectively.

Wide-field (~8' x 8') and deep near-infrared (JHKs bands) polarization images of the Orion nebulae (IRNe) around young stellar objects (YSOs), both massive and low-mass. We found the IRNe around both IRc2 and BN to be very extensive, suggesting that there might be two extended (>0.7 pc) bipolar/monopolar IRNe in these sources. We discovered at least 13 smaller-scale (~0.01-0.1 pc) IRNe around less-massive YSOs including the famous source theta^2 Ori C. We also suggest the presence of many unresolved (<690 AU) systems around low-mass YSOs and young brown dwarfs showing possible intrinsic polarizations. Wide-field infrared polarimetry is thus demonstrated to be a powerful technique in revealing IRNe and hence potential disk/outflow systems among high-mass to substellar YSOs.

On the basis of the near infrared observations of bulge red clump stars near the Galactic center, we have determined the galactocentric distance to be R_0 = 7.52 +- 0.10 (stat) +- 0.35 (sys) kpc. We observed the red clump stars at |l| < 1.0 deg and 0.7 deg < |b| < 1.0 deg with the IRSF 1.4 m telescope and the SIRIUS camera in the H and Ks bands. After extinction and population corrections, we obtained (m - M)_0 = 14.38 +- 0.03 (stat) +- 0.10 (sys). The statistical error is dominated by the uncertainty of the intrinsic local red clump stars' luminosity. The systematic error is estimated to be +- 0.10 including uncertainties in extinction and population correction, zero-point of photometry, and the fitting of the luminosity function of the red clump stars. Our result, R_0 = 7.52 kpc, is in excellent agreement with the distance determined geometrically with the star orbiting the massive black hole in the Galactic center. The recent result based on the spatial distribution of globular clusters is also consistent with our result. In addition, our study exhibits that the distance determination to the Galactic center with the red clump stars, even if the error of the population correction is taken into account, can achieve an uncertainty of about 5%, which is almost the same level as that in recent geometrical determinations.

An extensive photometric monitoring of KH 15D, an enigmatic variable in the young star cluster NGC 2264, has been conducted. Simultaneous and accurate near-infrared (JHKs-bands) photometry is presented between 2003 December and 2005 March covering most of the variable phase. The infrared variability is characterized by large-amplitude and long-lasting eclipse, as observed at optical. The period of variability is 48.3 +/- 0.2 days, the maximum photometric amplitude of variability is ~4.2 mag, and the eclipse duration is \~0.5 in phase units. These are consistent with the most recent period, amplitude, and duration at optical. The blueing of the J-H color (~0.16 mag) during the eclipse, which has been suggested before, is unambiguously confirmed; a similar blueing at H-Ks is less clear but is probably present at a similar level. The overall shape of the JHKs light curves is very similar to the optical one, including a fair time-symmetry and a less stable flux during the eclipse with a slight hump near the zero phase. Most of these variability features of KH 15D observed at near-infrared wavelengths can be explained with the recent model employing an eclipse by the inclined, precessing disk and an outer scattering region around a pre-main-sequence binary.

On the basis of near-infrared imaging observations, we derived visual extinction (Av) distribution toward ten Bok globules through measurements of both the color excess (E_{H-K}) and the stellar density at J, H, and Ks (star count). Radial column density profiles for each globule were analyzed with the Bonnor-Ebert sphere model. Using the data of our ten globules and four globules in the literature, we investigated the stability of globules on the basis of xi_max, which characterizes the Bonnor-Ebert sphere as well as the stability of the equilibrium state against the gravitational collapse. We found that more than half of starless globules are located near the critical state (xi_max = 6.5 +/- 2). Thus, we suggest that a nearly critical Bonnor-Ebert sphere characterizes the typical density structure of starless globules. Remaining starless globules show clearly unstable states (xi_max > 10). Since unstable equilibrium states are not long maintained, we expect that these globules are on the way to gravitational collapse or that they are stabilized by non-thermal support. It was also found that all the star-forming globules show unstable solutions of xi_max >10, which is consistent with the fact that they have started gravitational collapse. We investigated the evolution of a collapsing gas sphere whose initial condition is a nearly critical Bonnor-Ebert sphere. We found that the column density profiles of the collapsing sphere mimic those of the static Bonnor-Ebert spheres in unstable equilibrium. The collapsing gas sphere resembles marginally unstable Bonnor-Ebert spheres for a long time. We found that the frequency distribution of xi_max for the observed starless globules is consistent with that from model calculations of the collapsing sphere.

We present the result of a near-infrared (J H Ks) survey along the Galactic plane, -10.5deg < l < +10.5deg and b=+1.0deg, with the IRSF 1.4m telescope and the SIRIUS camera. Ks vs. H-Ks color-magnitude diagrams reveal a well-defined population of red clump (RC) stars whose apparent magnitude peak changes continuously along the Galactic plane, from Ks=13.4 at l=-10deg to Ks=12.2 at l=+10deg after dereddening. This variation can be explained by the bar-like structure found in previous studies, but we find an additional inner structure at |l| < 4deg, where the longitude - apparent magnitude relation is distinct from the outer bar, and the apparent magnitude peak changes by only 0.1 mag over the central 8deg. The exact nature of this inner structure is as yet uncertain.

We present a study of the newly discovered McNeil's nebula in Orion using the JHKs-band simultaneous observations with the near-infrared (NIR) camera SIRIUS on the IRSF 1.4m telescope. The cometary infrared nebula is clearly seen extending toward north and south from the NIR source (V1647 Orionis) that illuminates McNeil's nebula. The compact nebula has an apparent diameter of about 70 arcsec. The nebula is blue (bright in J) and has a cavity structure with two rims extending toward north-east and north-west. The north-east rim is brighter and sharp, while the north-west rim is diffuse. The north-east rim can be traced out to ~ 40 arcsec from the location of the NIR source. In contrast, no cavity structure is seen toward the south, although diffuse nebula is extended out to ~ 20 arcsec. New NIR photometric data show a significant variation in the magnitudes (> 0.15 mag) of the source of McNeil's nebula within a period of one week, that is possibly under the phase of eruptive variables like FUors or EXors.

We present sub-arcsecond (FWHM ~ 0".7), NIR JHKs-band images and a high sensitivity radio continuum image at 1280 MHz, using SIRIUS on UH 88-inch telescope and GMRT. The NIR survey covers an area of ~ 24 arcmin^2 with 10-sigma limiting mags of ~ 19.5, 18.4, and 17.3 in J, H, and Ks-band, respectively. Our NIR images are deeper than any JHK surveys to date for the larger area of NGC 7538 star forming region. We construct JHK CC and J-H/J and H-K/K CM diagrams to identify YSOs and to estimate their masses. Based on these CC and CM diagrams, we identified a rich population of YSOs (Class I and Class II), associated with the NGC 7538 region. A large number of red sources (H-K > 2) have also been detected around NGC 7538. We argue that these red stars are most probably PMS stars with intrinsic color excesses. Most of YSOs in NGC 7538 are arranged from the N-W toward S-E regions, forming a sequence in age: the diffuse H II region (N-W, oldest: where most of the Class II and Class I sources are detected); the compact IR core (center); and the regions with the extensive IR reflection nebula and a cluster of red young stars (S-E and S). We find that the slope of the KLF of NGC 7538 is lower than the typical values reported for the young embedded clusters, although equally low values have also been reported in the W3 Main star forming region. From the slope of the KLF and the analysis by Megeath et al. (1996), we infer that the embedded stellar population is comprised of YSOs with an age of ~ 1 Myr. Based on the comparison between models of PMS stars with the observed CM diagram we find that the stellar population in NGC 7538 is primarily composed of low mass PMS stars similar to those observed in the W3 Main star forming region.

We present deep near-infrared (J, H, Ks) images toward an embedded cluster which lies in a C18O clump in the cloud C of the Vela Molecular Ridge. This cluster has at least ~ 350 members and a radius of ~ 0.5 pc. The stellar surface number density is approximately 3000 stars pc^-2 in the central 0.1 pc x 0.1 pc region of the cluster. This is much higher than most of the young clusters within 1 kpc of the Sun. From the comparison of the luminosity function and near-infrared excess fraction with those of other embedded clusters, we estimate that the age of this cluster is approximately 2-3 Myr. This cluster exhibits an excess of brighter stars in its central region, from which we conclude that the more massive stars are located near the cluster center.

The data of 8,852 and 2,927 variable stars detected by OGLE survey in the Large and Small Magellanic Clouds are presented. They are cross-identified with the SIRIUS JHK survey data, and their infrared properties are discussed. Variable red giants are well separated on the period-J - K plane, suggesting that it could be a good tool to distinguish their pulsation mode and type.

We present a deep JHKs-band imaging survey of the W3 Main star forming region, using the near-infrared camera, SIRIUS, mounted on the University of Hawaii 2.2m telescope. The near-infrared survey covers an area of ~ 24 sq. arcmin with 10 sigma limiting magnitudes of ~ 19.0, 18.1, and 17.3 in J, H, and Ks-band, respectively. We construct JHK color-color and J/J-H and K/H-K color-magnitude diagrams to identify young stellar objects and estimate their masses. Based on these color-color and color-magnitude diagrams, a rich population of YSOs is identified which is associated with the W3 Main region. A large number of previously unreported red sources (H-K > 2) have also been detected around W3 Main. We argue that these red stars are most probably pre-main sequence stars with intrinsic color excesses. We find that the slope of the Ks-band luminosity function of W3 Main is lower than the typical values reported for the young embedded clusters. The derived slope of the KLF is the same as that found by Megeath et al. (1996), from which analysis by Megeath et al. indicates that the W3 Main region has an age in the range of 0.3--1 Myr. Based on the comparison between models of pre-main sequence stars with the observed color-magnitude diagram we find that the stellar population in W3 Main is primarily composed of low mass pre-main sequence stars. We also report the detection of isolated young stars with large infrared excesses which are most probably in their earliest evolutionary phases.

We report on the discovery of an infrared counterpart to the X-ray transient XTE J1720-318 on 2003 January 18, nine days after an X-ray outburst, and the infrared light curve during the first 130 days after the outburst. The infrared light curve shows a decline of about 1.2 mag from the peak magnitude of Ks about 15.3 over the observation period, and a secondary maximum, about 40 days after the outburst. Another small increase in the flux was also recorded about 20 days after the outburst. These increases were also detected in the X-ray light curve. The J H Ks colors are consistent with an X-ray irradiated accretion disk suffering an extinction of Av about 8, which is also inferred from its X-ray spectrum and the extinction map constructed from far-infrared dust emission of this line of sight. These J, H, and Ks observations demonstrate that useful data can be obtained even for such an object, which suffers heavy optical extinction, possibly located beyond the Galactic center.

We have performed a cross-identification between OGLE-II data and single-epoch SIRIUS JHK survey data in the LMC and SMC. After eliminating obvious spurious variables, we determined the pulsation periods for 9,681 and 2,927 variables in the LMC and SMC, respectively. Based on these homogeneous data, we studied the pulsation properties and metallicity effects on period-K magnitude (PK) relations by comparing the variable stars in the LMC and SMC. The sample analyzed here is much larger, and we found the following new features: (1) variable red giants in the SMC form parallel sequences on the PK plane, just like those found by Wood (2000) in the LMC; (2) both of the sequences A and B of Wood (2000) have discontinuities, and they occur at the K-band luminosity of the TRGB; (3) the sequence B of Wood (2000) separates into three independent sequences B+- and C'; (4) comparison between the theoretical pulsation models (Wood et al. 1996) and observational data suggests that the variable red giants on sequences C and newly discovered C' are pulsating in the fundamental and first overtone mode, respectively; (5) the theory can not explain the pulsation mode of sequences A+- and B+-, and they are unlikely to be the sequences for the first and second overtone pulsators, as was previously suggested; (6) the zero points of PK relations of Cepheids in the metal deficient SMC are fainter than those of LMC ones by ~0.1 mag but those of SMC Miras are brighter than those of LMC ones by ~0.13 mag, which are probably due to metallicity effects.

We have carried out a brown dwarf search over an area of 14' x 23' near the central portion of the Pleiades open cluster in five optical and near infrared bands (i', Z, J, H, Ks) with 10 sigma detection limits of i' = 22.0, J = 20.0 and Ks = 18.5 mag. The surveyed area has large extinction in excess of Av = 3 in the Pleiades region. We detected four new brown dwarf candidates from the colour-colour (J-K, i'-J) and the colour-magnitude (J, i'-K) diagrams. We estimated their masses as 0.046 Msun down to 0.028 Msun. The least massive one is estimated to have a mass smaller than Roque 25 or int-pl-IZ-69, and maybe the lowest mass object found so far in the Pleiades cluster.

We study the evolution of elliptical galaxies at z >~ 1 on the basis of a sample of 247 Extremely Red Objects (EROs) with R-Ks >= 3.35 (AB) and Ks <= 22.1 (AB) constructed from BVRi'z'JHKs multicolor data of a 114 arcmin2 area in the Subaru/XMM Deep Survey Field. By fitting template spectra of old ellipticals (OEs) and young, dusty starbursts (DSs) to the multicolor data, we classify EROs into these two classes and estimate their redshifts. We find that 58% of the EROs in our sample belong to the OE class and that 24% of the OEs are fit by a spectrum having a disk component with the B-band bulge-to-total luminosity ratio of <= 0.9. We compare the observed surface density of OEs with predictions from passive evolution models. We then derive rest-frame B-band luminosity functions (LFs) of OEs in our sample at z=1-1.5 and 1.5-2.5. We find that the LF at z=1-1.5 roughly agrees with the LF of local ellipticals if a dimming of 1.3 mag from z=1.25 to the present epoch is assumed. On the other hand, the amplitude of the LF at z=1.5-2.5 is found to be lower than that of the local LF by a factor of ~ 3 over the whole range of magnitude observed. Taking account of a strong decrease in the number density of morphologically classified early-type galaxies at z >~ 1.5 found by several authors, we conclude that the majority of ellipticals seen at present have not established either a red color or a smooth 1/4-law profile before z ~ 1.5$. (abridged)

The first results of our ongoing near-infrared (NIR) survey of the variable red giants in the Large Magellanic Cloud, using the Infrared Survey Facility (IRSF) and the SIRIUS infrared camera, are presented. Many very red stars were detected and found that most of them are variables. In the observed colour-magnitude diagram ($J-K, K$) and the stellar $K$ magnitude distribution, the tip of the first giant branch (TRGB), where helium burning in the core starts, is clearly seen. Apart from the genuine AGB variables, we found many variable stars at luminosities around the TRGB. From this result, we infer that a substantial fraction of them are RGB variables.