Coalescence of neutron stars gives rise to kilonova, thermal emission powered by radioactive decays of freshly synthesized r-process nuclei. Although observational properties are largely affected by bound-bound opacities of r-process elements, available atomic data have been limited. In this paper, we study element-to-element variation of the opacities in the ejecta of neutron star mergers by performing systematic atomic structure calculations of r-process elements for the first time. We show that the distributions of energy levels tend to be higher as electron occupation increases for each electron shell due to the larger energy spacing caused by larger effects of spin-orbit and electron-electron interactions. As a result, elements with a fewer number of electrons in the outermost shells tend to give larger contributions to the bound-bound opacities. This implies that Fe is not representative for the opacities of light r-process elements. The average opacities for the mixture of r-process elements are found to be kappa ~ 20-30 cm^2 g^{-1} for the electron fraction of Ye < 0.20, kappa ~ 3-5 cm^2 g^{-1} for Ye = 0.25-0.35, and kappa ~ 1 cm^2 g^{-1} for Ye = 0.40 at T = 5,000-10,000 K, and they steeply decrease at lower temperature. We show that, even with the same abundance or Ye, the opacity in the ejecta changes with time by one order of magnitude from 1 to 10 days after the merger. Our radiative transfer simulations with the new opacity data confirm that ejecta with a high electron fraction (Ye >~ 0.25, with no lanthanide) are needed to explain the early, blue emission in GW170817/AT2017gfo while lanthanide-rich ejecta (with a mass fraction of lanthanides ~ 5 x 10^{-3}) reproduce the long-lasting near-infrared emission.

Lanthanide elements play important roles as an opacity source in the ejected material from neutron star mergers. Accurate and complete atomic data are necessary to evaluate the opacities and to analyze the observed data. In this paper, we perform extended, ab-initio atomic calculations from Pr II (Z=59) to Gd II (Z=64). By using multiconfiguration Dirac-Hartree-Fock and relativistic configuration-interaction methods, implemented in the general-purpose relativistic atomic structure package GRASP2K, we calculate the energy levels and transition data of electric dipole transitions. These computations are based on strategies (with small variations) of Nd II published by Gaigalas et al. (2019). Accuracy of data is evaluated by comparing computed energy levels with the NIST database or other works. For the energy levels, we obtain the average relative accuracy of 8%, 12%, 6%, 8%, and 7% for Pr II, Pm II, Sm II, Eu II, and Gd II ions, respectively as compared with the NIST data. Accuracy of energy transfer to the wavelength as 3%, 14% and 11% for Pr II, Eu II and Gd II. Our computed E1 type transition probabilities are in good agreement with experimental values presented by other authors especially for strong transitions.

The main aim of this paper is to present accurate energy levels of the ground [Xe]$4f^{12}$ and first excited [Xe]$4f^{11}5d$ configurations of Er$^{2+}$. The energy level structure of the Er$^{2+}$ ion was computed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction (RCI) methods, as implemented in the GRASP2018 program package. The Breit interaction, self-energy and vacuum polarization corrections were included in the RCI computations. The zero-first-order approach was used in the computations. Energy levels with the identification in $LS$ coupling for all (399) states belonging to the [Xe]$4f^{12}$ and [Xe]$4f^{11}5d$ configurations are presented. Electric dipole (E1) transition data between the levels of these two configurations are computed. The accuracy of the these data are evaluated by studying the behaviour of the transition rates as functions of the gauge parameter as well as by evaluating the cancellation factors. The core electron correlations were studied using different strategies. Root-mean-square deviations obtained in this study for states of the ground and excited configurations from the available experimental or semi-empirical data are 649 cm$^{-1}$, and 747 cm$^{-1}$, respectively.

Coalescence of binary neutron star give rise to electromagnetic emission, kilonova, powered by radioactive decays of r-process nuclei. Observations of kilonova associated with GW170817 provided unique opportunity to study the heavy element synthesis in the Universe. However, atomic data of r-process elements to decipher the light curves and spectral features of kilonova are not fully constructed yet. In this paper, we perform extended atomic calculations of neodymium (Nd, Z=60) to study the impact of accuracies in atomic calculations to the astrophysical opacities. By employing multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods, we calculate energy levels and transition data of electric dipole transitions for Nd II, Nd III, and Nd IV ions. Compared with previous calculations, our new results provide better agreement with the experimental data. The accuracy of energy levels was achieved in the present work 10 %, 3 % and 11 % for Nd II, Nd III and Nd IV, respectively, comparing with the NIST database. We confirm that the overall properties of the opacity are not significantly affected by the accuracies of the atomic calculations. The impact to the Planck mean opacity is up to a factor of 1.5, which affects the timescale of kilonova at most 20 %. However, we find that the wavelength dependent features in the opacity are affected by the accuracies of the calculations. We emphasize that accurate atomic calculations, in particular for low-lying energy levels, are important to provide predictions of kilonova light curves and spectra.

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.

We present laboratory spectra of the $3p$--$3d$ transitions in Fe$^{14+}$ and Fe$^{15+}$ excited with a mono-energetic electron beam. In the energy dependent spectra obtained by sweeping the electron energy, resonant excitation is confirmed as an intensity enhancement at specific electron energies. The experimental results are compared with theoretical cross sections calculated based on fully relativistic wave functions and the distorted-wave approximation. Comparisons between the experimental and theoretical results show good agreement for the resonance strength. A significant discrepancy is, however, found for the non-resonant cross section in Fe$^{14+}$. %, which can be considered as a fundamental cause of the line intensity ratio problem that has often been found in both observatory and laboratory measurements. This discrepancy is considered to be the fundamental cause of the previously reported inconsistency of the model with the observed intensity ratio between the $^3\!P_2$ -- $^3\!D_3$ and $^1\!P_1$ -- $^1\!D_2$ transitions.

Ejected material from neutron star mergers give rise to electromagnetic emission powered by radioactive decays of r-process nuclei, which is so called kilonova or macronova. While properties of the emission are largely affected by opacities in the ejected material, available atomic data for r-process elements are still limited. We perform atomic structure calculations for r-process elements: Se (Z=34), Ru (Z=44), Te (Z=52), Ba (Z=56), Nd (Z=60), and Er (Z=68). We confirm that the opacities from bound-bound transitions of open f-shell, Lanthanide elements (Nd and Er) are higher than those of the other elements over a wide wavelength range. The opacities of open s-shell (Ba), p-shell (Se and Te), and d-shell (Ru) elements are lower than those of open f-shell elements and their transitions are concentrated in the ultraviolet wavelengths. We show that the optical brightness can be different by >2 mag depending on the element abundances in the ejecta such that post-merger, Lanthanide-free ejecta produce brighter and bluer optical emission. Such blue emission from post-merger ejecta can be observed from the polar directions if the mass of the preceding dynamical ejecta in these regions is small. For the ejecta mass of 0.01 Msun, observed magnitudes of the blue emission will reach 21.0 mag (100 Mpc) and 22.5 mag (200 Mpc) in g and r bands within a few days after the merger, which are detectable with 1m or 2m-class telescopes.

We present a far-infrared all-sky atlas from a sensitive all-sky survey using the Japanese $AKARI$ satellite. The survey covers $> 99$% of the sky in four photometric bands centred at 65 $\mu$m, 90 $\mu$m, 140 $\mu$m, and 160 $\mu$m with spatial resolutions ranging from 1 to 1.5 arcmin. These data provide crucial information for the investigation and characterisation of the properties of dusty material in the Interstellar Medium (ISM), since significant portion of its energy is emitted between $\sim$50 and 200 $\mu$m. The large-scale distribution of interstellar clouds, their thermal dust temperatures and column densities, can be investigated with the improved spatial resolution compared to earlier all-sky survey observations. In addition to the point source distribution, the large-scale distribution of ISM cirrus emission, and its filamentary structure, are well traced. We have made the first public release of the full-sky data to provide a legacy data set for use by the astronomical community.

We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of Tokyo Atacama Observatory 1.0m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5640 m above sea level) in northern Chile. The camera provides a field of view of 5'.1 $\times$ 5'.1 with a spatial resolution of 0".298 /pixel in the wavelength range of 0.95 to 2.4 $\mu$m. Taking advantage of the dry site, the camera is capable of hydrogen Paschen-$\alpha$ (Pa$\alpha$, $\lambda=$1.8751 $\mu$m in air) narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pa$\alpha$ images of Galactic objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band filters ($N1875$, $N191$) including the atmospheric absorption show larger dispersion (~10%) than those at broad-band filters (a few %), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor (PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, and derive that the median and the dispersion of the distribution of the PWV are 0.40+/-0.30 mm for $N1875$ and 0.37+/-0.21 mm for $N191$, which are remarkably smaller (49+/-38% for $N1875$ and 59+/-26% for $N191$) than radiometry measurements at the base of Cerro Chajnantor (5100 m alt.). The decrease in PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor is approximated at an exponential profile with scale heights within 0.3-1.9 km (previously observed values at night). We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability for a "ground-based" Pa$\alpha$ observation.

The evolution of galaxies and the history of star formation in the Universe are among the most important topics in today's astrophysics. Especially, the role of small, irregular galaxies in the star-formation history of the Universe is not yet clear. Using the data from the AKARI IRC survey of the Large Magellanic Cloud at 3.2, 7, 11, 15, and 24 {\mu}m wavelengths, i.e., at the mid- and near-infrared, we have constructed a multiwavelength catalog containing data from a cross-correlation with a number of other databases at different wavelengths. We present the separation of different classes of stars in the LMC in color-color, and color-magnitude, diagrams, and analyze their contribution to the total LMC flux, related to point sources at different infrared wavelengths.

Near-infrared polarimetric imaging observations toward the Galactic center have been carried out to examine the efficiency and wavelength dependence of interstellar polarization. A total area of about 5.7 deg$^2$ is covered in the $J$, $H$, and $K_S$ bands. We examined the polarization efficiency, defined as the ratio of degree of polarization to color excess. The interstellar medium between the Galactic center and us shows the polarization efficiency lower than that in the Galactic disk by a factor of three. Moreover we investigated the spatial variation of the polarization efficiency by comparing it with those of color excess, degree of polarization, and position angle. The spatial variations of color excess and degree of polarization depend on the Galactic latitude, while the polarization efficiency varies independently of the Galactic structure. Position angles are nearly parallel to the Galactic plane, indicating the longitudinal magnetic field configuration between the Galactic center and us. The polarization efficiency anticorrelates with dispersions of position angles. The low polarization efficiency and its spatial variation can be explained by the differences of the magnetic field directions along the line-of-sight. From the lower polarization efficiency, we suggest a higher strength of a random component relative to a uniform component of the magnetic field between the Galactic center and us. We also derived the ratios of degree of polarization $p_H/p_J$ = 0.581 $\pm$ 0.004 and $p_{K_S}/p_H$ = 0.620 $\pm$ 0.002. The power law indices of the wavelength dependence of polarization are $\beta_{JH}$ = 2.08 $\pm$ 0.02 and $\beta_{HK_S}$ = 1.76 $\pm$ 0.01. Therefore the wavelength dependence of interstellar polarization exhibits flattening toward longer wavelengths in the range of 1.25$-$2.14 $\micron$. The flattening would be caused by aligned large-size dust grains.

We present a near- to mid-infrared point source catalog of 5 photometric bands at 3.2, 7, 11, 15 and 24 um for a 10 deg2 area of the Large Magellanic Cloud (LMC) obtained with the Infrared Camera (IRC) onboard the AKARI satellite. To cover the survey area the observations were carried out at 3 separate seasons from 2006 May to June, 2006 October to December, and 2007 March to July. The 10-sigma limiting magnitudes of the present survey are 17.9, 13.8, 12.4, 9.9, and 8.6 mag at 3.2, 7, 11, 15 and 24 um, respectively. The photometric accuracy is estimated to be about 0.1 mag at 3.2 um and 0.06--0.07 mag in the other bands. The position accuracy is 0.3" at 3.2, 7 and 11um and 1.0" at 15 and 24 um. The sensitivities at 3.2, 7, and 24 um are roughly comparable to those of the Spitzer SAGE LMC point source catalog, while the AKARI catalog provides the data at 11 and 15 um, covering the mid-infrared spectral range contiguously. Two types of catalog are provided: a Catalog and an Archive. The Archive contains all the detected sources, while the Catalog only includes the sources that have a counterpart in the Spitzer SAGE point source catalog. The Archive contains about 650,000, 140,000, 97,000, 43,000, and 52,000 sources at 3.2, 7, 11, 15, and 24 um, respectively. Based on the catalog, we discuss the luminosity functions at each band, the color-color diagram, and the color-magnitude diagram using the 3.2, 7, and 11 um band data. Stars without circumstellar envelopes, dusty C-rich and O-rich stars, young stellar objects, and background galaxies are located at distinct regions in the diagrams, suggesting that the present catalog is useful for the classification of objects towards the LMC.

Far-infrared observations provide crucial data for the investigation and characterisation of the properties of dusty material in the Interstellar Medium (ISM), since most of its energy is emitted between ~100 and 200 um. We present the first all-sky image from a sensitive all-sky survey using the Japanese AKARI satellite, in the wavelength range 50 -- 180 um. Covering >99% of the sky in four photometric bands with four filters centred at 65 um, 90 um, 140 um, and 160 um wavelengths, this achieved spatial resolutions from 1 to 2 arcmin and a detection limit of <10 MJy sr-1, with absolute and relative photometric accuracies of <20%. All-sky images of the Galactic dust continuum emission enable astronomers to map the large-scale distribution of the diffuse ISM cirrus, to study its thermal dust temperature, emissivity and column density, and to measure the interaction of the Galactic radiation field and embedded objects with the surrounding ISM. In addition to the point source population of stars, protostars, star-forming regions, and galaxies, the high Galactic latitude sky is shown to be covered with a diffuse filamentary-web of dusty emission that traces the potential sites of high latitude star formation. We show that the temperature of dust particles in thermal equilibrium with the ambient interstellar radiation field can be estimated by using 90 um, 140 um, and 160 um data. The FIR AKARI full-sky maps provide a rich new data set within which astronomers can investigate the distribution of interstellar matter throughout our Galaxy, and beyond.

Mid-infrared images frequently suffer artifacts and extended point spread functions (PSFs). We investigate the characteristics of the artifacts and the PSFs in images obtained with the Infrared Camera (IRC) onboard AKARI at four mid-infrared bands of the S7 (7{\mu}m), S11 (11{\mu}m), L15 (15{\mu}m), and L24 (24 {\mu}m). Removal of the artifacts significantly improves the reliability of the ref- erence data for flat-fielding at the L15 and L24 bands. A set of models of the IRC PSFs is also constructed from on-orbit data. These PSFs have extended components that come from diffraction and scattering within the detector arrays. We estimate the aperture correction factors for point sources and the surface brightness correction factors for diffuse sources. We conclude that the surface brightness correction factors range from 0.95 to 0.8, taking account of the extended component of the PSFs. To correct for the extended PSF effects for the study of faint structures, we also develop an image reconstruction method, which consists of the deconvolution with the PSF and the convolution with an appropriate Gaussian. The appropriate removal of the artifacts, improved flat-fielding, and image reconstruction with the extended PSFs enable us to investigate de- tailed structures of extended sources in IRC mid-infrared images.

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.

In this paper, we derive the period-luminosity (P-L) relation for Large Magellanic Cloud (LMC) Cepheids based on mid-infrared AKARI observations. AKARI's IRC sources were matched to the OGLE-III LMC Cepheid catalog. Together with the available I band light curves from the OGLE-III catalog, potential false matches were removed from the sample. This procedure excluded most of the sources in the S7 and S11 bands: hence only the P-L relation in the N3 band was derived in this paper. Random-phase corrections were included in deriving the P-L relation for the single epoch AKARI data, even though the derived P-L relation is consistent with the P-L relation without random-phase correction, though there is a \sim 7 per-cent improvement in the dispersion of the P-L relation. The final adopted N3 band P-L relation is N3 = -3.246 log(P) + 15.844, with a dispersion of 0.149.

The aim of this study is to understand the chemical conditions of ices around embedded young stellar objects (YSOs) in the metal-poor Large Magellanic Cloud (LMC). We performed near-infrared (2.5-5 micron) spectroscopic observations toward 12 massive embedded YSOs and their candidates in the LMC using the Infrared Camera (IRC) onboard AKARI. We estimated the column densities of the H2O, CO2, and CO ices based on their 3.05, 4.27, and 4.67 micron absorption features, and we investigated the correlation between ice abundances and physical properties of YSOs.The ice absorption features of H2O, CO2, 13CO2, CO, CH3OH, and possibly XCN are detected in the spectra. In addition, hydrogen recombination lines and PAH emission bands are detected toward the majority of the targets. The derived typical CO2/H2O ice ratio of our samples (~0.36 +- 0.09) is greater than that of Galactic massive YSOs (~0.17 +- 0.03), while the CO/H2O ice ratio is comparable. It is shown that the CO2 ice abundance does not correlate with the observed characteristics of YSOs; the strength of hydrogen recombination line and the total luminosity. Likewise, clear no correlation is seen between the CO ice abundance and YSO characteristics, but it is suggested that the CO ice abundance of luminous samples is significantly lower than in other samples.The systematic difference in the CO2 ice abundance around the LMC's massive YSOs, which was suggested by previous studies, is confirmed with the new near-infrared data. We suggest that the strong ultraviolet radiation field and/or the high dust temperature in the LMC are responsible for the observed high abundance of the CO2 ice. It is suggested that the internal stellar radiation does not play an important role in the evolution of the CO2 ice around a massive YSO, while more volatile molecules like CO are susceptible to the effect of the stellar radiation.

We carried out a near- to mid-infrared imaging and spectroscopic observations of the patchy areas in the Small Magellanic Cloud using the Infrared Camera on board AKARI. Two 100 arcmin2 areas were imaged in 3.2, 4.1, 7, 11, 15, and 24 um and also spectroscopically observed in the wavelength range continuously from 2.5 to 13.4 um. The spectral resolving power (lambda/Delta lambda) is about 20, 50, and 50 at 3.5, 6.6 and 10.6 um, respectively. Other than the two 100 arcmin2 areas, some patchy areas were imaged and/or spectroscopically observed as well. In this paper, we overview the observations and present a list of near- to mid-infrared photometric results, which lists ~ 12,000 near-infrared and ~ 1,800 mid-infrared bright point sources detected in the observed areas. The 10 sigma limits are 16.50, 16.12, 13.28, 11.26, 9.62, and 8.76 in Vega magnitudes at 3.2, 4.1, 7, 11, 15, and 24 um bands, respectively.

We have determined interstellar extinction law toward the Galactic center (GC) at the wavelength from 1.2 to 8.0 micron, using point sources detected in the IRSF/SIRIUS near-infrared survey and those in the 2MASS and Spitzer/IRAC/GLIMPSE II catalogs. The central region |l| < 3deg and |b| < 1deg has been surveyed in the J, H and Ks bands with the IRSF telescope and the SIRIUS camera whose filters are similar to the Mauna Kea Observatories (MKO) near-infrared photometric system. Combined with the GLIMPSE II point source catalog, we made Ks versus (Ks - lambda) color-magnitude diagrams where lambda = 3.6, 4.5, 5.8, and 8.0 micron. The Ks magnitudes of bulge red clump stars and the (Ks - lambda) colors of red giant branches are used as a tracer of the reddening vector in the color-magnitude diagrams. From these magnitudes and colors, we have obtained the ratios of total to selective extinction A(Ks)/E(Ks-lambda) for the four IRAC bands. Combined with A(lambda)/A(Ks) for the J and H bands derived by Nishiyama et al., we obtain A(J):A(H):A(Ks):A([3.6]):A([4.5]):A([5.8]):A([8.0])=3.02:1.73:1:0.50:0.39:0.36:0.43 for the line of sight toward the GC. This confirms the flattening of the extinction curve at lambda > 3 micron from a simple extrapolation of the power-law extinction at shorter wavelengths, in accordance with recent studies. The extinction law in the 2MASS JHKs bands has also been calculated, and a good agreement with that in the MKO system is found. In nearby molecular clouds and diffuse interstellar medium, the lack of reliable measurements of the total to selective extinction ratios hampers unambiguous determination of the extinction law; however, observational results toward these lines of sight cannot be reconciled with a single extinction law.

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 the first results of AKARI Infrared Camera near-infrared spec- troscopic survey of the Large Magellanic Cloud (LMC). We detected absorption features of the H2O ice 3.05 um and the CO2 ice 4.27 um stretching mode toward seven massive young stellar objects (YSOs). These samples are for the first time spectroscopically confirmed to be YSOs. We used a curve-of-growth method to evaluate the column densities of the ices and derived the CO2/H2O ratio to be 0.45 pm 0.17. This is clearly higher than that seen in Galactic massive YSOs (0.17 pm 0.03). We suggest that the strong ultraviolet radiation field and/or the high dust temperature in the LMC may be responsible for the observed high CO2 ice abundance.

We observed an area of 10 deg^2 of the Large Magellanic Cloud using the Infrared Camera on board AKARI. The observations were carried out using five imaging filters (3, 7, 11, 15, and 24 micron) and a dispersion prism (2 -- 5 micron, $\lambda / \Delta\lambda$ $\sim$ 20) equipped in the IRC. This paper describes the outline of our survey project and presents some initial results using the imaging data that detected over 5.9x10^5 near-infrared and 6.4x10^4 mid-infrared point sources. The 10 $\sigma$ detection limits of our survey are about 16.5, 14.0, 12.3, 10.8, and 9.2 in Vega-magnitude at 3, 7, 11, 15, and 24 micron, respectively. The 11 and 15 micron data, which are unique to AKARI IRC, allow us to construct color-magnitude diagrams that are useful to identify stars with circumstellar dust. We found a new sequence in the color-magnitude diagram, which is attributed to red giants with luminosity fainter than that of the tip of the first red giant branch. We suggest that this sequence is likely to be related to the broad emission feature of aluminium oxide at 11.5 micron. The 11 and 15 micron data also indicate that the ([11] - [15]) micron color of both oxygen-rich and carbon-rich red giants once becomes blue and then turns red again in the course of their evolution, probably due to the change in the flux ratio of the silicate or silicon carbide emission feature at 10 or 11.3 micron to the 15 micron flux.

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 have found Herbig Ae/Be star candidates in the western region of the Magellanic Bridge. Using the near infrared camera SIRIUS and the 1.4 m telescope IRSF, we surveyed about 3.0 deg x 1.3 deg (24 deg < RA < 36 deg, -75 deg < Dec. < -73.7 deg) in the J, H, and Ks bands. On the basis of colors and magnitudes, about 200 Herbig Ae/Be star candidates are selected. Considering the contaminations by miscellaneous sources such as foreground stars and early-type dwarfs in the Magellanic Bridge, we estimate that about 80 (about 40%) of the candidates are likely to be Herbig Ae/Be stars. We also found one concentration of the candidates at the young star cluster NGC 796, strongly suggesting the existence of pre-main-sequence (PMS) stars in the Magellanic Bridge. This is the first detection of PMS star candidates in the Magellanic Bridge, and if they are genuine PMS stars, this could be direct evidence of recent star formation. However, the estimate of the number of Herbig Ae/Be stars depends on the fraction of classical Be stars, and thus a more precise determination of the Be star fraction or observations to differentiate between the Herbig Ae/Be stars and classical Be stars are required.

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.

We report the result of our near-infrared observations (JHKs) for type II Cepheids (including possible RV Tau stars) in galactic globular clusters. We detected variations of 46 variables in 26 clusters (10 new discoveries in seven clusters) and present their light curves. Their periods range from 1.2 d to over 80 d. They show a well-defined period-luminosity relation at each wavelength. Two type II Cepheids in NGC6441 also obey the relation if we assume the horizontal branch stars in NGC6441 are as bright as those in metal-poor globular clusters in spite of the high metallicity of the cluster. This result supports the high luminosity which has been suggested for the RR Lyr variables in this cluster. The period-luminosity relation can be reproduced using the pulsation equation (P sqrt(rho)=Q) assuming that all the stars have the same mass. Cluster RR Lyr variables were found to lie on an extrapolation of the period-luminosity relation. These results provide important constraints on the parameters of the variable stars. Using Two Micron All-Sky Survey (2MASS) data, we show that the type II Cepheids in the Large Magellanic Cloud (LMC) fit our period-luminosity relation within the expected scatter at the shorter periods. However, at long periods ($P>40$ d, i.e. in the RV Tau star range) the LMC field variables are brighter by about one magnitude than those of similar periods in galactic globular clusters. The long-period cluster stars also differ from both these LMC stars and galactic field RV Tau stars in a colour-colour diagram. The reasons for these differences are discussed.

We have determined the ratios of total to selective extinction in the near-infrared bands (J, H, Ks) toward the Galactic center from the observations of the region |l| < 2.0deg and 0.5deg < |b| < 1.0deg with the IRSF telescope and the SIRIUS camera. Using the positions of red clump stars in color-magnitude diagrams as a tracer of the extinction and reddening, we determine the average of the ratios of total to selective extinction to be A(Ks)/E(H-Ks) = 1.44+-0.01, A(Ks)/E(J-Ks) = 0.494+-0.006, and A(H)/E(J-H) = 1.42+-0.02, which are significantly smaller than those obtained in previous studies. From these ratios, we estimate that A(J) : A(H) : A(Ks) = 1 : 0.573+-0.009 : 0.331+-0.004 and E(J-H)/E(H-Ks) = 1.72+-0.04, and we find that the power law A(lambda) \propto lambda^{-1.99+-0.02} is a good approximation over these wavelengths. Moreover, we find a small variation in A(Ks)/E(H-Ks) across our survey. This suggests that the infrared extinction law changes from one line of sight to another, and the so-called ``universality'' does not necessarily hold in the infrared wavelengths.

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

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.