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.

Luminous infrared galaxies (LIRGs) are enshrouded by a large amount of dust, produced by their active star formation, and it is difficult to measure their activity in the optical wavelength. We have carried out Pa$\alpha$ narrow-band imaging observations of 38 nearby star-forming galaxies including 33 LIRGs listed in $IRAS$ RBGS catalog with the Atacama Near InfraRed camera (ANIR) on the University of Tokyo Atacama Observatory (TAO) 1.0 m telescope (miniTAO). Star formation rates (SFRs) estimated from the Pa$\alpha$ fluxes, corrected for dust extinction using the Balmer Decrement Method (typically $A_V$ $\sim$ 4.3 mag), show a good correlation with those from the bolometric infrared luminosity of $IRAS$ data within a scatter of 0.27 dex. This suggests that the correction of dust extinction for Pa$\alpha$ flux is sufficient in our sample. We measure the physical sizes and the surface density of infrared luminosities ($\Sigma_{L(\mathrm{IR})}$) and $SFR$ ($\Sigma_{SFR}$) of star-forming region for individual galaxies, and find that most of the galaxies follow a sequence of local ultra luminous or luminous infrared galaxies (U/LIRGs) on the $L(\mathrm{IR})$-$\Sigma_{L(\mathrm{IR})}$ and $SFR$-$\Sigma_{SFR}$ plane. We confirm that a transition of the sequence from normal galaxies to U/LIRGs is seen at $L(\mathrm{IR})=8\times10^{10}$ $L_{\odot}$. Also, we find that there is a large scatter in physical size, different from those of normal galaxies or ULIRGs. Considering the fact that most of U/LIRGs are merging or interacting galaxies, this scatter may be caused by strong external factors or differences of their merging stage.

The properties of K-band selected galaxies (K_AB<24) in the z = 3.09 SSA22 protocluster field are studied. 430 galaxies at 2.6 < z_phot < 3.6 are selected as potential protocluster members in a 112 arcmin^2 area based on their photometric redshifts. We find that \approx 20% of the massive galaxies with stellar masses >10^11 M_sun at z_phot \sim 3.1 have colors consistent with those of quiescent galaxies with ages > 0.5 Gyr. This fraction increases to \approx 50% after correcting for unrelated foreground/background objects. We also find that 30% of the massive galaxies are heavily reddened dusty star-forming galaxies. Few such quiescent galaxies at similar redshifts are seen in typical survey fields. An excess surface density of 24\mu m sources at z_phot \sim 3.1 is also observed, implying the presence of dusty star-formation activity in the protocluster. Cross-correlation with the X-ray data indicates that the fraction of K-band selected protocluster galaxies hosting active galactic nuclei (AGN) is also high compared with the field. The sky distribution of the quiescent galaxies, the 24\mu m sources, and the X-ray AGNs show clustering around a density peak of z=3.1 Ly\alpha emitters (LAEs). A significant fraction of the massive galaxies have already become quiescent, while the dusty star-formation is still active in the SSA22 protocluster. These findings indicate that we are witnessing the formation epoch of massive early-type galaxies at the center of predecessors to present-day rich galaxy clusters.

We present the results of wide-field deep JHK imaging of the SSA22 field using MOIRCS instrument equipped with Subaru telescope. The observed field is 112 arcmin^2 in area, which covers the z=3.1 protocluster characterized by the overdensities of Ly Alpha emitters (LAEs) and Ly Alpha Blobs (LABs). The 5 sigma limiting magnitude is K_{AB} = 24.3. We extract the potential protocluster members from the K-selected sample by using the multi-band photometric-redshift selection as well as the simple color cut for distant red galaxies (DRGs; J-K_{AB}>1.4). The surface number density of DRGs in our observed fields shows clear excess compared with those in the blank fields, and the location of the densest area whose projected overdensity is twice the average coincides with the large-scale density peak of LAEs. We also found that K-band counterparts with z_{phot} = 3.1 are detected for 75% (15/20) of the LABs within their Ly Alpha halo, and the 40 % (8/20) of LABs have multiple components, which gives a direct evidence of the hierarchical multiple merging in galaxy formation. The stellar mass ofLABs correlates with their luminosity, isophotal area, and the Ly Alpha velocity widths, implying that the physical scale and the dynamical motion of Ly Alpha emission are closely related to their previous star-formation activities. Highly dust-obscured galaxies such as hyper extremely red objects (HEROs; J-K_{AB}>2.1) and plausible K-band counterparts of submillimeter sources are also populated in the high density region.

We have developed high dispersion VPH (volume phase holographic) grisms with zinc selenide (ZnSe) prisms for the cryogenic optical system of MOIRCS (Multi-Object near InfraRed Camera and Spectrograph) for Y-, J-, H- and K- band observations. We fabricated the VPH gratings using a hologram resin. After several heat cycles at between room temperature and 120 K, the VPH gratings were assembled to grisms by gluing with two ZnSe prisms. Several heat cycles were also carried out for the grisms before being installed into MOIRCS. We measured the efficiencies of the VPH grisms in a laboratory, and found them to be 70% - 82%. The performances obtained by observations of MOIRCS with the 8.2 m Subaru Telescope have been found to be very consistent with the results in the laboratory test. This is the first astronomical application of cryogenic VPH grisms.

We present deep J-, H-, and Ks-band imaging data of the MOIRCS Deep Survey (MODS), which was carried out with Multi-Object Infrared Camera and Spectrograph (MOIRCS) mounted on the Subaru telescope in the GOODS-North region. The data reach 5sigma total limiting magnitudes for point sources of J=23.9, H=22.8, and Ks=22.8 (Vega magnitude) over 103 arcmin^2 (wide field). In 28 arcmin^2 of the survey area, which is ultra deep field of the MODS (deep field), the data reach the 5sigma depths of J=24.8, H=23.4, and Ks=23.8. The spatial resolutions of the combined images are FWHM ~ 0.6 arcsec and ~ 0.5 arcsec for the wide and deep fields in all bands, respectively. Combining the MODS data with the multi-wavelength public data taken with the HST, Spitzer, and other ground-based telescopes in the GOODS field, we construct a multi-wavelength photometric catalog of Ks-selected sources. Using the catalog, we present Ks-band number counts and near-infrared color distribution of the detected objects, and demonstrate some selection techniques with the NIR colors for high redshift galaxies. These data and catalog are publicly available via internet.

We present results from Submillimeter Array (SMA) 860-micron sub-arcsec astrometry and multiwavelength observations of the brightest millimeter (S_1.1mm = 8.4 mJy) source, SSA22-AzTEC1, found near the core of the SSA22 protocluster that is traced by Ly\alpha emitting galaxies at z = 3.09. We identify a 860-micron counterpart with a flux density of S_860um = 12.2 +/- 2.3 mJy and absolute positional accuracy that is better than 0.3". At the SMA position, we find radio to mid-infrared counterparts, whilst no object is found in Subaru optical and near-infrared deep images at wavelengths \le 1 micron (J > 25.4 in AB, 2\sigma). The photometric redshift estimate, using flux densities at \ge 24 microns, indicates z_phot = 3.19^{+0.26}_{-0.35}, consistent with the protocluster redshift. We then model the near-to-mid-infrared spectral energy distribution (SED) of SSA22-AzTEC1, and find that the SED modeling requires a large extinction (A_V \approx 3.4 mag) of starlight from a stellar component with M_star ~ 10^{10.9} M_sun, assuming z = 3.1. Additionally, we find a significant X-ray counterpart with a very hard spectrum (Gamma_eff = -0.34 ^{+0.57}_{-0.61}), strongly suggesting that SSA22-AzTEC1 harbors a luminous AGN (L_X ~ 3*10^{44} ergs s^{-1}) behind a large hydrogen column (N_H ~ 10^{24} cm^{-2}). The AGN, however, is responsible for only ~10% of the bolometric luminosity of the host galaxy, and therefore the star-formation activity likely dominates the submillimeter emission. It is possible that SSA22-AzTEC1 is the first example of a protoquasar growing at the bottom of the gravitational potential underlying the SSA22 protocluster.

We investigate rest-frame near-infrared (NIR) morphologies of a sample of 139 galaxies with M_{s} >= 1 x 10^{10} M_{sun} at z=0.8-1.2 in the GOODS-North field using our deep NIR imaging data (MOIRCS Deep Survey, MODS). We focus on Luminous Infrared Galaxies (LIRGs), which dominate high star formation rate (SFR) density at z~1, in the sample identified by cross-correlating with the Spitzer/MIPS 24um source catalog. We perform two-dimensional light profile fitting of the z~1 galaxies in the Ks-band (rest-frame J-band) with a single component Sersic model. We find that at z~1, ~90% of LIRGs have low Sersic indices (n<2.5, similar to disk-like galaxies) in the Ks-band, and those disk-like LIRGs consist of ~60% of the whole disk-like sample above M_{s} >= 3 x 10^{10} M_{sun}. The z~1 disk-like LIRGs are comparable or ~20% small at a maximum in size compared to local disk-like galaxies in the same stellar mass range. If we examine rest-frame UV-optical morphologies using the HST/ACS images, the rest-frame B-band sizes of the z~1 disk-like galaxies are comparable to those of the local disk-like galaxies as reported by previous studies on size evolution of disk-like galaxies in the rest-frame optical band. Measuring color gradients (galaxy sizes as a function of wavelength) of the z~1 and local disk-like galaxies, we find that the z~1 disk-like galaxies have 3-5 times steeper color gradient than the local ones. Our results indicate that (i) more than a half of relatively massive disk-like galaxies at z~1 are in violent star formation epochs observed as LIRGs, and also (ii) most of those LIRGs are constructing their fundamental disk structure vigorously. The high SFR density in the universe at z~1 may be dominated by such star formation in disk region in massive galaxies.

We study the evolution of star formation activity of galaxies at 0.5<z<3.5 as a function of stellar mass, using very deep NIR data taken with Multi-Object Infrared Camera and Spectrograph (MOIRCS) on the Subaru telescope in the GOODS-North region. The NIR imaging data reach K ~ 23-24 Vega magnitude and they allow us to construct a nearly stellar mass-limited sample down to ~ 10^{9.5-10} Msun even at z~3. We estimated star formation rates (SFRs) of the sample with two indicators, namely, the Spitzer/MIPS 24um flux and the rest-frame 2800A luminosity. The SFR distribution at a fixed Mstar shifts to higher values with increasing redshift at 0.5<z<3.5. More massive galaxies show stronger evolution of SFR at z>~1. We found galaxies at 2.5<z<3.5 show a bimodality in their SSFR distribution, which can be divided into two populations by a constant SSFR of ~2 Gyr^{-1}. Galaxies in the low-SSFR group have SSFRs of ~ 0.5-1.0 Gyr^{-1}, while the high-SSFR population shows ~10 Gyr^{-1}. The cosmic SFRD is dominated by galaxies with Mstar = 10^{10-11} Msun at 0.5<z<3.5, while the contribution of massive galaxies with Mstar = 10^{11-11.5} Msun shows a strong evolution at z>1 and becomes significant at z~3, especially in the case with the SFR based on MIPS 24um. In galaxies with Mstar = 10^{10-11.5} Msun, those with a relatively narrow range of SSFR (<~1 dex) dominates the cosmic SFRD at 0.5<z<3.5. The SSFR of galaxies which dominate the SFRD systematically increases with redshift. At 2.5<z<3.5, the high-SSFR population, which is relatively small in number, dominates the SFRD. Major star formation in the universe at higher redshift seems to be associated with a more rapid growth of stellar mass of galaxies.

We present near-infrared (NIR) spectroscopy for 37 BzK-color-selected star-forming galaxies with MOIRCS on Subaru. The sample is drawn from the Ks-band selected catalog of the MOIRCS Deep Survey (MODS) in the GOODS-N region. About half of our samples are selected from the MIPS 24um-source catalog. Ha emission lines are detected from 23 galaxies, of which the median redshift is 2.12. We derived the star formation rates (SFRs) from extinction-corrected Ha luminosities. The extinction correction is estimated from the SED fitting of multi-band photometric data covering UV to NIR wavelengths. The Balmer decrement of the stacked emission lines shows that the amount of extinction for the ionized gas is larger than that for the stellar continuum. From a comparison of the extinction corrected Ha luminosity and other SFR indicators we found that the relation between the dust properties of stellar continuum and ionized gas is different depending on the intrinsic SFR. The comparison between the Ha SFR and stellar mass estimated from SED fitting shows no correlation between them. Some galaxies with stellar mass smaller than ~10^10 Msun show SFRs higher than ~100 Msun/yr. The specific SFRs (SSFRs) of these galaxies are remarkably high. The average stellar-population age of these high-SSFR galaxies is younger than 100 Myr, which is consistent with the implied high SSFR. The large SFR implies the possibility that the high SSFR galaxies significantly contribute to the cosmic SFR density of the universe at z~2. The total SFR density estimated from the Ha emission line galaxies is 0.089-0.136 Msun/yr/Mpc^3, which is consistent with the total SFR densities in the literature. The metallicity of the high-SSFR galaxies is larger than that expected from the mass-metallicity relation of UV-selected galaxies at z~2 by Erb et al. (2006a). [Abridged]

We examine whether the spectral energy distribution of optical continuum emission of active galactic nuclei (AGNs) changes during flux variation, based on accurate and frequent monitoring observations of 11 nearby Seyfert galaxies and QSOs carried out in the B, V, and I bands for seven years by the MAGNUM telescope. The multi-epoch flux data in any two different bands obtained on the same night show a very tight linear flux to flux relationship for all target AGNs. The flux of the host galaxy within the photometric aperture is carefully estimated by surface brightness fitting to available high-resolution HST images and MAGNUM images. The flux of narrow emission lines in the photometric bands is also estimated from available spectroscopic data. We find that the non-variable component of the host galaxy plus narrow emission lines for all target AGNs is located on the fainter extension of the linear regression line of multi-epoch flux data in the flux to flux diagram. This result strongly indicates that the spectral shape of AGN continuum emission in the optical region does not systematically change during flux variation. The trend of spectral hardening that optical continuum emission becomes bluer as it becomes brighter, which has been reported by many studies, is therefore interpreted as the domination of the variable component of the nearly constant spectral shape of an AGN as it brightens over the non-variable component of the host galaxy plus narrow lines, which is usually redder than AGN continuum emission.

We use very deep near-infrared (NIR) imaging data obtained in MOIRCS Deep Survey (MODS) to investigate the evolution of the galaxy stellar mass function back to z~3. The MODS data reach J=24.2, H=23.1, K=23.1 (5sigma, Vega magnitude) over 103 arcmin^2 (wide) and J=25.1, H=23.7, K=24.1 over 28 arcmin^2 (deep) in the GOODS-North region. The wide and very deep NIR data allow us to measure the number density of galaxies down to low stellar mass (10^9-10^10 Msun) even at high redshift with high statistical accuracy. The normalization of the mass function decreases with redshift and the integrated stellar mass density becomes ~ 8-18% of the local value at z~2 and ~ 4-9% at z~3, which are consistent with results of previous studies in general fields. Furthermore, we found that the low-mass slope becomes steeper with redshift from alpha ~- 1.3 at z~1 to alpha ~- 1.6 at z~3, and that the evolution of the number density of low-mass (10^9-10^10 Msun) galaxies is weaker than that of M* (~10^11 Msun) galaxies. This indicates that the contribution of low-mass galaxies to the total stellar mass density has been significant at high redshift. The steepening of the low-mass slope with redshift is opposite trend expected from the stellar mass dependence of the specific star formation rate reported in previous studies. The present result suggests that the hierarchical merging process overwhelmed the effect of the stellar mass growth by star formation and was very important for the stellar mass assembly of these galaxies at 1<~z<~3.

We investigate the X-ray properties of the K-band-selected galaxies at redshift 2 < z < 4 by using our deep near-infrared images obtained in the MOIRCS Deep Survey project and the published Chandra X-ray source catalog. 61 X-ray sources with the 2-10 keV luminosity L_X = 10^{42}-10^{44} erg/s are identified with the K-selected galaxies and we found that they are exclusively (90%) associated with the massive objects with stellar mass larger than 10^{10.5} Msun. Our results are consistent with the idea that the M_BH/M_str ratio of the galaxies at z=2-4 is similar to the present-day value. On the other hand, the AGN detection rate among the very massive galaxies with the stellar mass larger than 10^{11} Msun is high, 33% (26/78). They are active objects in the sense that the black-hole mass accretion rate is ~ 1-50% of the Eddington limit if they indeed have similar M_BH/M_str ratio with those observed in the local universe. The active duration in the AGN duty cycle of the high-redshift massive galaxies seems large.

We present spectroscopy of 15 star-forming BzK galaxies (sBzKs) with K(AB)<23 in the Subaru Deep Field, for which Halpha and some other emission lines are detected in 0.9 to 2.3 micron spectra with a resolution of R=500. Using Halpha luminosities, we obtain star formation rates (SFRs), and then specific SFRs (SSFRs) dividing SFRs by stellar masses, which are derived from SED fitting to BVRi'z'K photometry. It is found that sBzKs with higher stellar masses have larger SFRs. A negative correlation is seen between stellar mass and SSFR, which is consistent with the previous results for z~2 galaxies. This implies that a larger growth of stellar mass occurs in less massive galaxies. In addition, gas-phase oxygen abundances, 12+log(O/H), are derived from the ratio of NII(lambda 6584) to Halpha using the N2 index method. We have found a correlation between stellar mass and oxygen abundance in the sense that more massive sBzKs tend to be more metal rich, which is qualitatively consistent with the relation for UV-selected z~2 galaxies. However, the metallicity of the sBzKs is ~0.2 dex higher than that of UV-selected galaxies with similar stellar masses, which is significant considering the small uncertainties. The sBzKs in our sample have redder R-K colors than the UV-selected galaxies. This galaxy color-dependence in the oxygen abundance may be caused by older or dustier galaxies having higher metallicities at z~2.

We present the results of deep near-infrared imaging observations of the z=3.1 proto-cluster region in the SSA22a field taken by MOIRCS mounted on the Subaru Telescope. We observed a 21.7 arcmin^2 field to the depths of J=24.5, H=24.3, and K=23.9 (5 sigma). We examine the distribution of the K-selected galaxies at z~3 by using the simple color cut for distant red galaxies (DRGs) as well as the photometric-redshift selection technique. The marginal density excess of DRGs and the photo-z selected objects are found around the two most luminous Ly alpha blobs (LABs). We investigate the correlation between the K-selected objects and the LABs, and find that several galaxies with stellar mass M_* = 10^9-10^11 M_solar exist in vicinity of LABs, especially around the two most luminous ones. We also find that 7 of the 8 LABs in the field have plausible K_s-band counterparts and the sum of the stellar mass possibly associated with LABs correlates with the luminosity and surface brightness of them, which implies that the origin of Ly alpha emission may be closely correlated with their stellar mass or their previous star formation phenomena.

We use very deep near-infrared imaging data taken with Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru Telescope to investigate the number counts of Distant Red Galaxies (DRGs). We have observed a 4x7 arcmin^2 field in the Great Observatories Origins Deep Survey North (GOODS-N), and our data reach J=24.6 and K=23.2 (5sigma, Vega magnitude). The surface density of DRGs selected by J-K>2.3 is 2.35+-0.31 arcmin^-2 at K<22 and 3.54+-0.38 arcmin^-2 at K<23, respectively. These values are consistent with those in the GOODS-South and FIRES. Our deep and wide data suggest that the number counts of DRGs turn over at K~22, and the surface density of the faint DRGs with K>22 is smaller than that expected from the number counts at the brighter magnitude. The result indicates that while there are many bright galaxies at 2<z<4 with the relatively old stellar population and/or heavy dust extinction, the number of the faint galaxies with the similar red color is relatively small. Different behaviors of the number counts of the DRGs and bluer galaxies with 2<z_phot<4 at K>22 suggest that the mass-dependent color distribution, where most of low-mass galaxies are blue while more massive galaxies tend to have redder colors, had already been established at that epoch.