We report an automated morphological classification of galaxies into S-wise spirals, Z-wise spirals, and non-spirals using big image data taken from Subaru/Hyper Suprime-Cam (HSC) Survey and a convolutional neural network(CNN)-based deep learning technique. The HSC i-band images are about 25 times deeper than those from the Sloan Digital Sky Survey (SDSS) and have a two times higher spatial resolution, allowing us to identify substructures such as spiral arms and bars in galaxies at z>0.1. We train CNN classifiers by using HSC images of 1447 S-spirals, 1382 Z-spirals, and 51,650 non-spirals. As the number of images in each class is unbalanced, we augment the data of spiral galaxies by horizontal flipping, rotation, and rescaling of images to make the numbers of three classes similar. The trained CNN models correctly classify 97.5% of the validation data, which is not used for training. We apply the CNNs to HSC images of a half million galaxies with an i-band magnitude of i<20 over an area of 320 deg^2. 37,917 S-spirals and 38,718 Z-spirals are identified, indicating no significant difference between the numbers of two classes. Among a total of 76,635 spiral galaxies, 48,576 are located at z>0.2, where we are hardly able to identify spiral arms in the SDSS images. Our attempt demonstrates that a combination of the HSC big data and CNNs has a large potential to classify various types of morphology such as bars, mergers and strongly-lensed objects.

We analyzed the young (2.8-Myr-old) binary system FS Tau A using near-infrared (H-band) high-contrast polarimetry data from Subaru/HiCIAO and sub-millimeter CO (J=2-1) line emission data from ALMA. Both the near-infrared and sub-millimeter observations reveal several clear structures extending to $\sim$240 AU from the stars. Based on these observations at different wavelengths, we report the following discoveries. One arm-like structure detected in the near-infrared band initially extends from the south of the binary with a subsequent turn to the northeast, corresponding to two bar-like structures detected in ALMA observations with an LSRK velocity of 1.19-5.64 km/s. Another feature detected in the near-infrared band extends initially from the north of the binary, relating to an arm-like structure detected in ALMA observations with an LSRK velocity of 8.17-16.43 km/s. From their shapes and velocities, we suggest that these structures can mostly be explained by two streamers that connect the outer circumbinary disk and the central binary components. These discoveries will be helpful for understanding the evolution of streamers and circumstellar disks in young binary systems.

The SR24 multi-star system hosts both circumprimary and circumsecondary disks, which are strongly misaligned from each other. The circumsecondary disk is circumbinary in nature. Interestingly, both disks are interacting, and they possibly rotate in opposite directions. To investigate the nature of this unique twin disk system, we present 0.''1 resolution near-infrared polarized intensity images of the circumstellar structures around SR24, obtained with HiCIAO mounted on the Subaru 8.2 m telescope. Both the circumprimary disk and the circumsecondary disk are resolved and have elongated features. While the position angle of the major axis and radius of the NIR polarization disk around SR24S are 55$^{\circ}$ and 137 au, respectively, those around SR24N are 110$^{\circ}$ and 34 au, respectively. With regard to overall morphology, the circumprimary disk around SR24S shows strong asymmetry, whereas the circumsecondary disk around SR24N shows relatively strong symmetry. Our NIR observations confirm the previous claim that the circumprimary and circumsecondary disks are misaligned from each other. Both the circumprimary and circumsecondary disks show similar structures in $^{12}$CO observations in terms of its size and elongation direction. This consistency is because both NIR and $^{12}$CO are tracing surface layers of the flared disks. As the radius of the polarization disk around SR24N is roughly consistent with the size of the outer Roche lobe, it is natural to interpret the polarization disk around SR24N as a circumbinary disk surrounding the SR24Nb-Nc system.

Whether the spiral structure of galaxies is trailing or leading has been a subject of debate. We present a new spin parity catalog of 146 spiral galaxies that lists the following three pieces of information: whether the spiral structure observed on the sky is S-wise or Z-wise; which side of the minor axis of the galaxy is darker and redder, based on examination of Pan-STARRS and/or ESO/DSS2 red image archives, and which side of the major axis of the galaxy is approaching to us based on the published literature. This paper confirms that all of the spiral galaxies in the catalog show a consistent relationship among these three parameters, without any confirmed counterexamples, that supports the generally accepted interpretation that all the spiral galaxies are trailing and that the darker/redder side of the galactic disk is closer to us. Although the results of this paper may not be surprising, they provide a rationale for analyzing the S/Z winding distribution of spiral galaxies, using the large and uniform image data bases available now and in the near future, to study the spin vorticity distribution of galaxies in order to constrain the formation scenarios of galaxies and the large-scale structure of the Universe.

We present H-band polarized scattered light imagery and JHK high-contrast spectroscopy of the protoplanetary disk around HD 163296 observed with the HiCIAO and SCExAO/CHARIS instruments at Subaru Observatory. The polarimetric imagery resolve a broken ring structure surrounding HD 163296 that peaks at a distance along the major axis of 0.65 (66 AU) and extends out to 0.98 (100 AU) along the major axis. Our 2011 H-band data exhibit clear axisymmetry, with the NW- and SE- side of the disk exhibiting similar intensities. Our data are clearly different than 2016 epoch H-band observations from VLT/SPHERE that found a strong 2.7x asymmetry between the NW- and SE-side of the disk. Collectively, these results indicate the presence of time variable, non-azimuthally symmetric illumination of the outer disk. Based on our 3D-MCRT modeling of contemporaneous IR spectroscopic and H-band polarized intensity imagery of the system, we suggest that while the system could plausibly host an inclined inner disk component, such a component is unlikely to be responsible for producing the observed time-dependent azimuthal variations in the outer scattered light disk of the system. While our SCExAO/CHARIS data are sensitive enough to recover the planet candidate identified from NIRC2 in the thermal IR, we fail to detect an object with a corresponding JHK brightness estimated from the atmospheric models of Baraffe et al. 2003. This suggests that the candidate is either fainter in JHK bands than model predictions, possibly due to extinction from the disk or atmospheric dust/clouds, or that it is an artifact of the dataset/data processing. Our SCExAO/CHARIS data lower the IR mass limits for planets inferred at larger stellocentric separations; however, these ALMA-predicted protoplanet candidates are currently still consistent with direct imaging constraints.

We report an orbital characterization of GJ1108Aab that is a low-mass binary system in pre-main-sequence phase. Via the combination of astrometry using adaptive optics and radial velocity measurements, an eccentric orbital solution of $e$=0.63 is obtained, which might be induced by the Kozai-Lidov mechanism with a widely separated GJ1108B system. Combined with several observed properties, we confirm the system is indeed young. Columba is the most probable moving group, to which the GJ1108A system belongs, although its membership to the group has not been established. If the age of Columba is assumed for GJ1108A, the dynamical masses of both GJ1108Aa and GJ1108Ab ($M_{\rm dynamical,GJ1108Aa}=0.72\pm0.04 M_{\odot}$ and $M_{\rm dynamical,GJ1108Ab}=0.30\pm0.03 M_{\odot}$) are more massive than what an evolutionary model predicts based on the age and luminosities. We consider the discrepancy in mass comparison can attribute to an age uncertainty; the system is likely older than stars in Columba, and effects that are not implemented in classical models such as accretion history and magnetic activity are not preferred to explain the mass discrepancy. We also discuss the performance of the evolutionary model by compiling similar low-mass objects in evolutionary state based on the literature. Consequently, it is suggested that the current model on average reproduces the mass of resolved low-mass binaries without any significant offsets.

We present $H$- and $K_{\rm s}$-bands observations of the LkH$\alpha$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investigated the azimuthal profiles at the ring and the outer-disk regions as well as radial profiles in the directions of the spiral-like structures and semi-major axis. Azimuthal analysis shows a large variety in wavelength and implies that the disk has non-axisymmetric dust distributions. The radial profiles in the major-axis direction (PA=$271^\circ$) suggest that the outer region ($r\geq0\farcs25$) may be influenced by shadows of the inner region of the disk. The spiral-like directions (PA=10$^\circ$ and 230$^\circ$) show different radial profiles, which suggests that the surfaces of the spiral-like structures are highly flared and/or have different dust properties. Finally, a color-map of the disk shows a lack of an outer eastern region in the $H$-band disk, which may hint the presence of an inner object that casts a directional shadow onto the disk.

We conducted high-contrast polarimetry observations of T Tau in the H-band, using the HiCIAO instrument mounted on the Subaru Telescope, revealing structures as near as 0.$\arcsec$1 from the stars T Tau N and T Tau S. The whole T Tau system is found to be surrounded by nebula-like envelopes, and several outflow-related structures are detected in these envelopes. We analyzed the detailed polarization patterns of the circumstellar structures near each component of this triple young star system and determined constraints on the circumstellar disks and outflow structures. We suggest that the nearly face-on circumstellar disk of T Tau N is no larger than 0.$\arcsec$8, or 117 AU, in the northwest, based on the existence of a hole in this direction, and no larger than 0.$\arcsec$27, or 40 AU, in the south. A new structure "N5" extends to about 0.$\arcsec$42, or 59 AU, on the southwest of the star, believed to be part of the disk. We suggest that T Tau S is surrounded by a highly inclined circumbinary disk with a radius of about 0.$\arcsec$3, or 44 AU, with a position angle of about 30$^\circ$, that is misaligned with the orbit of the T Tau S binary. After analyzing the positions and polarization vector patterns of the outflow-related structures, we suggest that T Tau S should trigger the well-known E-W outflow, and is also likely to be responsible for a southwest precessing outflow "coil" and a possible south outflow.

Large exoplanet surveys have successfully detected thousands of exoplanets to-date. Utilizing these detections and non-detections to constrain our understanding of the formation and evolution of planetary systems also requires a detailed understanding of the basic properties of their host stars. We have determined the basic stellar properties of F, K, and G stars in the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS) survey from echelle spectra taken at the Apache Point Observatory's 3.5m telescope. Using ROBOSPECT to extract line equivalent widths and TGVIT to calculate the fundamental parameters, we have computed Teff, log(g), vt, [Fe/H], chromospheric activity, and the age for our sample. Our methodology was calibrated against previously published results for a portion of our sample. The distribution of [Fe/H] in our sample is consistent with that typical of the Solar neighborhood. Additionally, we find the ages of most of our sample are $< 500 Myrs$, but note that we cannot determine robust ages from significantly older stars via chromospheric activity age indicators. The future meta-analysis of the frequency of wide stellar and sub-stellar companions imaged via the SEEDS survey will utilize our results to constrain the occurrence of detected co-moving companions with the properties of their host stars.

We detect 20 $z=7.0$ Ly$\alpha$ emitter (LAE) candidates to $L({\rm Ly}\alpha) \geq 2 \times 10^{42}$ erg s$^{-1}$ or $0.3$ $L^*_{z=7}$ and in $6.1\times 10^5$ Mpc$^3$ volume in the Subaru Deep Field and the Subaru/XMM-Newton Deep Survey field by 82 and 37 hours of Subaru Suprime-Cam narrowband NB973 and reddest optical $y$-band imaging. We compare their Ly$\alpha$ and UV luminosity functions (LFs) and densities and Ly$\alpha$ equivalent widths (EWs) to those of $z=5.7$, 6.6 and 7.3 LAEs from previous Suprime-Cam surveys. The Ly$\alpha$ LF (density) rapidly declines by a factor of $\times$1.5 (1.9) in $L({\rm Ly}\alpha)$ at $z=5.7-6.6$ (160 Myr), $\times$1.5 (1.6) at $z=6.6-7.0$ (60 Myr) at the faint end and $\times$2.0 (3.8) at $z=7.0-7.3$ (40 Myr). Also, in addition to the systematic decrease in EW at $z=5.7-6.6$ previously found, 2/3 of the $z=7.0$ LAEs detected in the UV continuum exhibit lower EWs than the $z=6.6$ ones. Moreover, while the UV LF and density do not evolve at $z=5.7-6.6$, they modestly decline at $z=6.6-7.0$, implying galaxy evolution contributing to the decline of the Ly$\alpha$ LF. Comparison of the $z=7.0$ Ly$\alpha$ LF to the one predicted by an LAE evolution model further reveals that galaxy evolution alone cannot explain all the decline of Ly$\alpha$ LF. If we attribute the discrepancy to Ly$\alpha$ attenuation by neutral hydrogen, the intergalactic medium transmission of Ly$\alpha$ photons at $z=7.0$ would be $T_{{\rm Ly}\alpha}^{\rm IGM} \leq 0.6-0.7$. It is lower (higher) than $T_{{\rm Ly}\alpha}^{\rm IGM}$ at $z=6.6$ (7.3) derived by previous studies, suggesting rapid increase in neutral fraction at $z > 6$.

We present high-contrast observations of 68 young stellar objects (YSOs) explored as part of the SEEDS survey on the Subaru telescope. Our targets are very young ($<$10 Myr) stars, which often harbor protoplanetary disks where planets may be forming. We achieve a typical contrast of $\sim$$10^{-4}$--$10^{-5.5}$ at an angular distance of 1\arcsec\ from the central star, corresponding to typical mass sensitivities (assuming hot-start evolutionary models) of $\sim$10 ${\rm M_J}$ at 70 AU and $\sim$6 ${\rm M_J}$ at 140 AU. We detected a new stellar companion to HIP 79462 and confirmed the substellar objects GQ Lup b and ROXs 42B b. An additional six companion candidates await follow-up observations to check for common proper motion. Our SEEDS YSO observations probe the population of planets and brown dwarfs at the very youngest ages; these may be compared to the results of surveys targeting somewhat older stars. Our sample and the associated observational results will help enable detailed statistical analyses of giant planet formation.

By performing non-masked polarization imaging with Subaru/HiCIAO, polarized scattered light from the inner region of the disk around the GG Tau A system was successfully detected in the $H$ band with a spatial resolution of approximately 0.07$\arcsec$, revealing the complicated inner disk structures around this young binary. This paper reports the observation of an arc-like structure to the north of GG Tau Ab and part of a circumstellar structure that is noticeable around GG Tau Aa extending to a distance of approximately 28 AU from the primary star. The speckle noise around GG Tau Ab constrains its disk radius to <13 AU. Based on the size of the circumbinary ring and the circumstellar disk around GG Tau Aa, the semi-major axis of the binary's orbit is likely to be 62 AU. A comparison of the present observations with previous ALMA and near-infrared (NIR) H$_2$ emission observations suggests that the north arc could be part of a large streamer flowing from the circumbinary ring to sustain the circumstellar disks. According to the previous studies, the circumstellar disk around GG Tau Aa has enough mass and can sustain itself for a duration sufficient for planet formation; thus, our study indicates that planets can form within close (separation $\lesssim$ 100 AU) young binary systems.

We present high-contrast H-band polarized intensity (PI) images of the transitional disk around the young solar-like star GM Aur. The near-infrared direct imaging of the disk was derived by polarimetric differential imaging using the Subaru 8.2-m Telescope and HiCIAO. An angular resolution and an inner working angle of 0."07 and r~0."05, respectively, were obtained. We clearly resolved a large inner cavity, with a measured radius of 18+/-2 au, which is smaller than that of a submillimeter interferometric image (28 au). This discrepancy in the cavity radii at near-infrared and submillimeter wavelengths may be caused by a 3-4M_Jup planet about 20 au away from the star, near the edge of the cavity. The presence of a near-infrared inner is a strong constraint on hypotheses for inner cavity formation in a transitional disk. A dust filtration mechanism has been proposed to explain the large cavity in the submillimeter image, but our results suggest that this mechanism must be combined with an additional process. We found that the PI slope of the outer disk is significantly different from the intensity slope obtained from HST/NICMOS, and this difference may indicate the grain growth process in the disk.

We present high signal-to-noise ratio, precise $YJH$ photometry and $Y$ band (\gpiwave~$\mu$m) spectroscopy of HD 1160 B, a young substellar companion discovered from the Gemini NICI Planet Finding Campaign, using the Subaru Coronagraphic Extreme Adaptive Optics instrument and the Gemini Planet Imager. HD 1160 B has typical mid-M dwarf-like infrared colors and a spectral type of M5.5$^{+1.0}_{-0.5}$, where the blue edge of our $Y$ band spectrum rules out earlier spectral types. Atmospheric modeling suggests HD 1160 B having an effective temperature of 3000--3100 $K$, a surface gravity of log $g$ = 4--4.5, a radius of~\bestfitradius~$R_{\rm J}$, and a luminosity of log $L$/$L_{\odot} = -2.76 \pm 0.05$. Neither the primary's Hertzspring-Russell diagram position nor atmospheric modeling of HD 1160 B show evidence for a sub-solar metallicity. The interpretation of the HD 1160 B depends on which stellar system components are used to estimate an age. Considering HD 1160 A, B and C jointly, we derive an age of 80--125 Myr, implying that HD 1160 B straddles the hydrogen-burning limit (70--90 $M_{\rm J}$). If we consider HD 1160 A alone, younger ages (20--125 Myr) and a brown dwarf-like mass (35--90 $M_{\rm J}$) are possible. Interferometric measurements of the primary, a precise GAIA parallax, and moderate resolution spectroscopy can better constrain the system's age and how HD 1160 B fits within the context of (sub)stellar evolution.

We present H-band (1.6 {\mu}m) scattered light observations of the transitional disk RX J1615.3-3255, located in the ~1 Myr old Lupus association. From a polarized intensity image, taken with the HiCIAO instrument of the Subaru Telescope, we deduce the position angle and the inclination angle of the disk. The disk is found to extend out to 68 $\pm$ 12 AU in scattered light and no clear structure is observed. Our inner working angle of 24 AU does not allow us to detect a central decrease in intensity similar to that seen at 30 AU in the 880 {\mu}m continuum observations. We compare the observations with multiple disk models based on the Spectral Energy Distribution (SED) and submm interferometry and find that an inner rim of the outer disk at 30 AU containing small silicate grains produces a polarized intensity signal which is an order of magnitude larger than observed. We show that a model in which the small dust grains extend smoothly into the cavity found for large grains is closer to the actual H-band observations. A comparison of models with different dust size distributions suggests that the dust in the disk might have undergone significant processing compared to the interstellar medium.

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and radial velocity (RV) measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be $m_1=1.141_{-0.091}^{+0.037}$ and $m_2=0.279^{+0.023}_{-0.020}$ M$_\odot$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. The Hipparcos time-series photometry shows a periodicity of $P=5.743$ d, which is also seen in SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is $380^{+220}_{-100}$ Myr old, consistent with an isochrone analysis ($220^{+2120}_{-90}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs and exoplanets by combination of RV and direct imaging data.

We present high-contrast H-band polarized intensity images of the transitional disk around the young solar-like star LkCa 15. By utilizing Subaru/HiCIAO for polarimetric differential imaging, both the angular resolution and the inner working angle reach 0.07" and r=0.1", respectively. We obtained a clearly resolved gap (width <~ 27 AU) at ~ 48 AU from the central star. This gap is consistent with images reported in previous studies. We also confirmed the existence of a bright inner disk with a misaligned position angle of 13+/-4 degree with respect to that of the outer disk, i.e., the inner disk is possibly warped. The large gap and the warped inner disk both point to the existence of a multiple planetary system with a mass of <~1Mjup.

We find a new substellar companion to the Pleiades member star, Pleiades HII 3441, using the Subaru telescope with adaptive optics. The discovery is made as part of the high-contrast imaging survey to search for planetary-mass and substellar companions in the Pleiades and young moving groups. The companion has a projected separation of 0".49 +/- 0".02 (66 +/- 2 AU) and a mass of 68 +/- 5 M_J based on three observations in the J-, H-, and K_S-band. The spectral type is estimated to be M7 (~2700 K), and thus no methane absorption is detected in the H band. Our Pleiades observations result in the detection of two substellar companions including one previously reported among 20 observed Pleiades stars, and indicate that the fraction of substellar companions in the Pleiades is about 10.0 +26.1/-8.8 %. This is consistent with multiplicity studies of both the Pleiades stars and other open clusters.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

We present new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We used a Monte Carlo, radiative transfer code to simultaneously model the system's SED and H-band polarized intensity imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aur's SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aur's spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and polarized intensity imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8 year baseline. If such structures are caused by disk-planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be >47 AU.

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

We present high-contrast angular differential imaging (ADI) observations of the debris disk around HD 32297 in H-band, as well as the first polarimetric images for this system in polarized differential imaging (PDI) mode with Subaru/HICIAO. In ADI, we detect the nearly edge-on disk at >5sigma levels from ~0.45 arcsec to ~1.7 arcsec (50-192 AU) from the star and recover the spine deviation from the midplane already found in previous works. We also find for the first time imaging and surface brightness (SB) indications for the presence of a gapped structure on both sides of the disk at distances of ~0.75 arcsec (NE side) and ~0.65 arcsec (SW side). Global forward-modelling work delivers a best-fit model disk and well-fitting parameter intervals that essentially match previous results, with high-forward scattering grains and a ring located at 110 AU. However, this single ring model cannot account for the gapped structure seen in our SB profiles. We create simple double ring models and achieve a satisfactory fit with two rings located at 60 and 95 AU, respectively, low-forward scattering grains and very sharp inner slopes. In polarized light we retrieve the disk extending from ~0.25-1.6 arcsec, although the central region is quite noisy and high S/N are only found in the range ~0.75-1.2 arcsec. The disk is polarized in the azimuthal direction, as expected, and the departure from the midplane is also clearly observed. Evidence for a gapped scenario is not found in the PDI data. We obtain a linear polarization degree of the grains that increases from ~10% at 0.55 arcsec to ~25% at 1.6 arcsec. The maximum is found at scattering angles of ~90degrees, either from the main components of the disk or from dust grains blown out to larger radii.

We present the first near-infrared scattered-light detection of the transitional disk around V1247 Ori, which was obtained using high-resolution polarimetric differential imaging observations with Subaru/HiCIAO. Our imaging in the H band reveals the disk morphology at separations of ~0.14"-0.86" (54-330 au) from the central star. The polarized intensity (PI) image shows a remarkable arc-like structure toward the southeast of the star, whereas the fainter northwest region does not exhibit any notable features. The shape of the arm is consistent with an arc of 0.28" $\pm$ 0.09" in radius (108 au from the star), although the possibility of a spiral arm with a small pitch angle cannot be excluded. V1247 Ori features an exceptionally large azimuthal contrast in scattered, polarized light; the radial peak of the southeastern arc is about three times brighter than the northwestern disk measured at the same distance from the star. Combined with the previous indication of an inhomogeneous density distribution in the gap at $\lesssim$46 au, the notable asymmetry in the outer disk suggests the presence of unseen companions and/or planet-forming processes ongoing in the arc.

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $\gamma$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($\iota$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $\iota$ Dra b, HD 5608 b, and HD 14067 b.

We present the results of an imaging observation campaign conducted with the Subaru Telescope adaptive optics system (IRCS+AO188) on 28 gravitationally lensed quasars (23 doubles, 1 quad, and 1 possible triple, and 3 candidates) from the SDSS Quasar Lens Search. We develop a novel modelling technique that fits analytical and hybrid point spread functions (PSFs), while simultaneously measuring the relative astrometry, photometry, as well as the lens galaxy morphology. We account for systematics by simulating the observed systems using separately observed PSF stars. The measured relative astrometry is comparable with that typically achieved with the Hubble Space Telescope, even after marginalizing over the PSF uncertainty. We model for the first time the quasar host galaxies in 5 systems, without a-priory knowledge of the PSF, and show that their luminosities follow the known correlation with the mass of the supermassive black hole. For each system, we obtain mass models far more accurate than those previously published from low-resolution data, and we show that in our sample of lensing galaxies the observed light profile is more elliptical than the mass, for ellipticity > 0.25. We also identify eight doubles for which the sources of external and internal shear are more reliably separated, and should therefore be prioritized in monitoring campaigns aimed at measuring time-delays in order to infer the Hubble constant.

A high angular resolution near-infrared polarized-intensity image of the GG Tau A binary system was obtained with the Subaru Telescope. The image shows the circumbinary disk scattering the light from the central binary. The azimuthal profile of the polarized intensity of the circumbinary disk is roughly reproduced by a simple disk model with the Henyey-Greenstein function and the Rayleigh function, indicating small dust grains at the surface of the disk. Combined with a previous observation of the circumbinary disk, our image indicates that the gap structure in the circumbinary disk orbits anti-clockwise, while material in the disk orbit clockwise. We propose a shadow of material located between the central binary and the circumbinary disk. The separations and position angles of the stellar components of the binary in the past 20 years are consistent with the binary orbit with a = 33.4 AU and e = 0.34.

We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 transitional disk. Our two epochs of imagery detect the scattered light disk from our effective inner working angle of 0.10" (13 AU) out to 0.50" (65 AU). This inner working angle is interior to the location of the system's gap inferred by previous studies using SED modeling (15 AU). We detected a candidate point source companion 1.08" northwest of the system; however, our second epoch of imagery strongly suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer models of the system, and our best fit models utilize a modestly inclined (50 deg), 0.01 Msun disk that has a partially depleted inner gap from the dust sublimation radius out to ~8 AU. Subtracting this best fit, axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the disk, which could be attributable to variety of mechanisms.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We report the late-time evolution of Type IIb Supernova (SN IIb) 2013df. SN 2013df showed a dramatic change in its spectral features at ~1 year after the explosion. Early on it showed typical characteristics shared by SNe IIb/Ib/Ic dominated by metal emission lines, while later on it was dominated by broad and flat-topped Halpha and He I emissions. The late-time spectra are strikingly similar to SN IIb 1993J, which is the only previous example clearly showing the same transition. This late-time evolution is fully explained by a change in the energy input from the $^{56}$Co decay to the interaction between the SN ejecta and dense circumstellar matter (CSM). The mass loss rate is derived to be ~(5.4 +- 3.2) x 10^{-5} Msun/yr (for the wind velocity of ~20 km/s), similar to SN 1993J but larger than SN IIb 2011dh by an order of magnitude. The striking similarity between SNe 2013df and 1993J in the (candidate) progenitors and the CSM environments, and the contrast in these natures to SN 2011dh, infer that there is a link between the natures of the progenitor and the mass loss: SNe IIb with a more extended progenitor have experienced a much stronger mass loss in the final centuries toward the explosion. It might indicate that SNe IIb from a more extended progenitor are the explosions during a strong binary interaction phase, while those from a less extended progenitor have a delay between the strong binary interaction and the explosion.

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at \lambda = 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

We present new high-resolution ($\sim$0\farcs09) $H$-band imaging observations of the circumstellar disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the presence of scattered light as close as 0\farcs15 ($\sim$20 AU) to the star. Within our image, we identify bright emission associated with a disk with a minimum radius of $\sim$90 AU, an inclination of $\sim$35\degr from the plane of the sky, and an approximate P.A. of 15\degr for the major axis. We find a brightness asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk structure. We also identify a pair of asymmetric tail structures extending east and west from the disk. The western tail extends at least 2\farcs5 (350 AU) from the star, and is probably associated with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern tail extends at least 1\arcsec (140 AU) at the present signal-to-noise. These tails are likely due to an encounter with an unseen brown dwarf, but our results do not exclude the explanation that these tails are outflow cavities or jets.

We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely.

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$\mu$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

We compared the number of faint stars detected in deep survey fields with the current stellar distribution model of the Galaxy and found that the detected number in the H band is significantly smaller than the predicted number. This indicates that M-dwarfs, the major component, are fewer in the halo and the thick disk. We used archived data of several surveys in both the north and south field of GOODS (Great Observatories Origins Deep Survey), MODS in GOODS-N, and ERS and CANDELS in GOODS-S. The number density of M-dwarfs in the halo has to be 20+/-13% relative to that in the solar vicinity, in order for the detected number of stars fainter than 20.5 mag in the H band to match with the predicted value from the model. In the thick disk, the number density of M-dwarfs must be reduced (52+/-13%) or the scale height must be decreased (~600 pc). Alternatively, overall fractions of the halo and thick disks can be significantly reduced to achieve the same effect, because our sample mainly consists of faint M-dwarfs. Our results imply that the M-dwarf population in regions distant from the Galactic plane is significantly smaller than previously thought. We then discussed the implications this has on the suitability of the model predictions for the prediction of non-companion faint stars in direct imaging extrasolar planet surveys by using the best-fit number densities.

We have obtained deep J & Ks-band images centered on a bright radio quiet QSO UM402 (z_{em}=2.856) using IRCS camera and AO systems on Subaru Telescope, as well as retrieved WFC3/F140W archive images. A faint galaxy (m_{k}=23.32 +/- 0.05 in the Vega system) that lies ~ 2.4" north of the QSO sightline has been clearly resolved in all three deep high resolution datasets, and appears as an irregular galaxy with two close components in the Ks-band images (separation ~ 0.3"). Given the small impact parameter (b=19.6 kpc, at z_{lls}=2.531), as well as the red color of (J-Ks)_{vega} ~ 1.6, it might be a candidate galaxy giving rise to the Lyman Limit system absorption at z_{abs}=2.531 seen in the QSO spectrum. After carefully subtracting the PSF from the QSO images, the host galaxy of this bright radio quiet QSO at z ~ 3 was marginally revealled. We placed a low limit of the host component of m_{k} ~ 23.3 according to our analyses.

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections ($\kappa$ And b, two $\sim$60 M$_{\rm J}$ brown dwarf companions in the Pleiades, PZ Tel B, and CD$-$35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of $\sim$30--100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at $\sim$5 M$_{\rm Jup}$, with a single power law distribution. We find that $p(M, a) \propto M^{-0.65\pm0.60} a^{-0.85\pm0.39}$ (1$\sigma$ errors) provides an adequate fit to our data, with 1.0--3.1\% (68\% confidence) of stars hosting 5--70 $M_{\rm Jup}$ companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs.

We present the ultra-deep Subaru narrowband imaging survey for Lya emitters (LAEs) at $z=7.3$ in SXDS and COSMOS fields with a total integration time of 106 hours. Exploiting our new sharp bandwidth filter, NB101, installed on Suprime-Cam, we have reached $L(Lya)=2.4\times10^{42} \ erg \ s^{-1}$ ($5\sigma$) for $z=7.3$ LAEs, about 4 times deeper than previous Subaru $z \gtrsim 7$ studies, which allows us to reliably investigate evolution of Lya luminosity function (LF), for the first time, down to the luminosity limit same as those of Subaru $z=3.1-6.6$ LAE samples. Surprisingly, we only find three and four LAEs in SXDS and COSMOS fields, respectively, while one expects a total of $\sim 65$ LAEs by our survey in the case of no Lya LF evolution from $z=6.6$ to $7.3$.We identify a decrease of Lya LF from $z=6.6$ to $7.3$ at the $>90\%$ confidence level from our $z=7.3$ Lya LF.Moreover, the evolution of Lya LF is clearly accelerated at $z>6.6$ beyond the measurement uncertainties including cosmic variance. Because no such accelerated evolution of UV-continuum LF or cosmic star-formation rate (SFR) is found at $z\sim 7$, but suggested only at $z>8$ (Oesch et al. 2013, Bouwens et al. 2014), this accelerated Lya LF evolution is explained by physical mechanisms different from a pure SFR decrease but related to Lya production and escape in the process of cosmic reionization. Because a simple accelerating increase of IGM neutral hydrogen absorbing Lya would not reconcile with Thomson scattering optical depth measurements from WMAP and Planck, our findings may support new physical pictures suggested by recent theoretical studies, such as the existence of HI clumpy clouds within cosmic ionized bubbles selectively absorbing Lya and the large ionizing photon escape fraction of galaxies making weak Lya emission.

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

We present the result of Subaru Telescope multi-band adaptive optics observations of the complex gravitationally lensed quasar SDSS J1405+0959, which is produced by two lensing galaxies. These observations reveal dramatically enhanced morphological detail, leading to the discovery of an additional object 0. 26'' from the secondary lensing galaxy, as well as three collinear clumps located in between the two lensing galaxies. The new object is likely to be the third quasar image, although the possibility that it is a galaxy cannot be entirely excluded. If confirmed via future observations, it would be the first three image lensed quasar produced by two galaxy lenses. In either case, we show based on gravitational lensing models and photometric redshift that the collinear clumps represent merging images of a portion of the quasar host galaxy, with a magnification factor of 15 - 20, depending on the model.

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

We present results from the first three years of observations of moving group targets in the SEEDS high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of ~10^5 at 1" and ~10^6 beyond 2" around 63 proposed members of nearby kinematic moving groups. We review each of the kinematic associations to which our targets belong, concluding that five, \beta Pictoris (~20 Myr), AB Doradus (~100 Myr), Columba (~30 Myr), Tucana-Horogium (~30 Myr), TW Hydrae (~10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these moving groups. For all of our targets, we combine proposed moving group membership with other age indicators where available, including Ca II HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, \kappa And, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper.

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array ($\sim1\arcsec$--3$\arcsec$ resolution), and high-resolution imaging of polarized intensity at the $K_s$-band by using the HiCIAO instrument on the Subaru Telescope ($0\farcs25$ resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H$_2$ mass of $2.4\times10^{-3}$ $M_\sun$ in the cold ($T<$30 K) outer part at $65<r<170$ AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ($>3\times10^{-9}$ $M_\sun$) of hot ($T\sim$180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

Most exoplanets detected by direct imaging so far have been characterized by relatively hot (> ~1000 K) and cloudy atmospheres. A surprising feature in some of their atmospheres has been a distinct lack of methane, possibly implying non-equilibrium chemistry. Recently, we reported the discovery of a planetary companion to the Sun-like star GJ 504 using Subaru/HiCIAO within the SEEDS survey. The planet is substantially colder (<600 K) than previously imaged planets, and has indications of fewer clouds, which implies that it represents a new class of planetary atmospheres with expected similarities to late T-type brown dwarfs in the same temperature range. If so, one might also expect the presence of significant methane absorption, which is characteristic of such objects. Here, we report the detection of deep methane absorption in the atmosphere of GJ 504 b, using the Spectral Differential Imaging mode of HiCIAO to distinguish the absorption feature around 1.6 um. We also report updated JHK photometry based on new Ks-band data and a re-analysis of the existing data. The results support the notion that GJ 504 b has atmospheric properties distinct from other imaged exoplanets, and will become a useful reference object for future planets in the same temperature range.

We report a discovery of a companion candidate around one of {\it Kepler} Objects of Interest (KOIs), KOI-94, and results of our quantitative investigation of the possibility that planetary candidates around KOI-94 are false positives. KOI-94 has a planetary system in which four planetary detections have been reported by {\it Kepler}, suggesting that this system is intriguing to study the dynamical evolutions of planets. However, while two of those detections (KOI-94.01 and 03) have been made robust by previous observations, the others (KOI-94.02 and 04) are marginal detections, for which future confirmations with various techniques are required. We have conducted high-contrast direct imaging observations with Subaru/HiCIAO in $H$ band and detected a faint object located at a separation of $\sim0.6''$ from KOI-94. The object has a contrast of $\sim 1\times 10^{-3}$ in $H$ band, and corresponds to an M type star on the assumption that the object is at the same distance of KOI-94. Based on our analysis, KOI-94.02 is likely to be a real planet because of its transit depth, while KOI-94.04 can be a false positive due to the companion candidate. The success in detecting the companion candidate suggests that high-contrast direct imaging observations are important keys to examine false positives of KOIs. On the other hand, our transit light curve reanalyses lead to a better period estimate of KOI-94.04 than that on the KOI catalogue and show that the planetary candidate has the same limb darkening parameter value as the other planetary candidates in the KOI-94 system, suggesting that KOI-94.04 is also a real planet in the system.

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star {\kappa} Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. We present new angular differential imaging (ADI) images of the Kappa Andromedae system at 2.146 (Ks), 3.776 (L'), 4.052 (NB 4.05) and 4.78 {\mu}m (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We derive a more accurate J = 15.86 +- 0.21, H = 14.95 +- 0.13, Ks = 14.32 +- 0.09 mag for {\kappa} And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at Ks and L' band. We derive NB 4.05 = 13.0 +- 0.2 and M' = 13.3 +- 0.3 mag and estimate Log10(L/Lsun) = -3.76 +- 0.06. We build the 1-5 microns spectral energy distribution of the companion and compare it to seven PHOENIX-based atmospheric models in order to derive Teff = 1900+100-200 K. Models do not set constrains on the surface gravity. ``Hot-start" evolutionary models predict masses of 14+25-2 MJup based on the luminosity and temperature estimates, and considering a conservative age range for the system (30+120-10 Myr). ``warm-start" evolutionary tracks constrain the mass to M >= 11 MJup. Therefore, the mass of {\kappa} Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of {\kappa} And b.

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures of 800--1800 K and very red colors (J - H > 0.5 mag), implying thick cloud covers. Furthermore, substantial model uncertainties exist at these young ages due to the unknown initial conditions at formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as a result, its estimated mass remains in the planetary regime independent of uncertainties related to choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling model, and given the system age of 160 [+350, -60] Myr, GJ 504 b has an estimated mass of 4 [+4.5, -1.0] Jupiter masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds the typical outer boundary of ~30 AU predicted for the core accretion mechanism. GJ 504 b is also significantly cooler (510 [+30, -20] K) and has a bluer color (J-H = -0.23 mag) than previously imaged exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their atmospheric properties.

Supernovae (SNe) have been proposed to be the main production sites of dust grains in the Universe. Our knowledge on their importance to dust production is, however, limited by observationally poor constraints on the nature and amount of dust particles produced by individual SNe. In this paper, we present a spectrum covering optical through near-Infrared (NIR) light of the luminous Type IIn supernova (SN IIn) 2010jl around one and half years after the explosion. This unique data set reveals multiple signatures of newly formed dust particles. The NIR portion of the spectrum provides a rare example where thermal emission from newly formed hot dust grains is clearly detected. We determine the main population of the dust species to be carbon grains at a temperature of ~1,350 - 1,450K at this epoch. The mass of the dust grains is derived to be ~(7.5 - 8.5) x 10^{-4} Msun. Hydrogen emission lines show wavelength-dependent absorption, which provides a good estimate on the typical size of the newly formed dust grains (~0.1 micron, and most likely <~0.01 micron). We attribute the dust grains to have been formed in a dense cooling shell as a result of a strong SN-circumstellar media (CSM) interaction. The dust grains occupy ~10% of the emitting volume, suggesting an inhomogeneous, clumpy structure. The average CSM density is required to be >~3 x 10^{7} cm^{-3}, corresponding to a mass loss rate of >~0.02 Msun yr^{-1} (for a mass loss wind velocity of ~100 km s^{-1}). This strongly supports a scenario that SN 2010jl and probably other luminous SNe IIn are powered by strong interactions within very dense CSM, perhaps created by Luminous Blue Variable (LBV)-like eruptions within the last century before the explosion.

We present near-infrared coronagraphic imaging polarimetry of RY Tau. The scattered light in the circumstellar environment was imaged at H-band at a high resolution (~0".05) for the first time, using Subaru-HiCIAO. The observed polarized intensity (PI) distribution shows a butterfly-like distribution of bright emission with an angular scale similar to the disk observed at millimeter wavelengths. This distribution is offset toward the blueshifted jet, indicating the presence of a geometrically thick disk or a remnant envelope, and therefore the earliest stage of the Class II evolutionary phase. We perform comparisons between the observed PI distribution and disk models with: (1) full radiative transfer code, using the spectral energy distribution (SED) to constrain the disk parameters; and (2) monochromatic simulations of scattered light which explore a wide range of parameters space to constrain the disk and dust parameters. We show that these models cannot consistently explain the observed PI distribution, SED, and the viewing angle inferred by millimeter interferometry. We suggest that the scattered light in the near-infrared is associated with an optically thin and geometrically thick layer above the disk surface, with the surface responsible for the infrared SED. Half of the scattered light and thermal radiation in this layer illuminates the disk surface, and this process may significantly affect the thermal structure of the disk.

We carried out an imaging survey for extrasolar planets around stars in the Pleiades (125 Myr, 135 pc) in the $H$ and $K_{S}$ bands using HiCIAO combined with the adaptive optics, AO188, on the Subaru telescope. We found 13 companion candidates fainter than 14.5 mag in the $H$ band around 9 stars. Five of these 13 were confirmed to be background stars by measurement of their proper motion. One was not found in the second epoch observation, and thus was not a background or companion object. One had multi-epoch image, but the precision of its proper motion was not sufficient to conclude whether it was background object. Four other candidates are waiting for second epoch observations to determine their proper motion. Finally, the remaining 2 were confirmed to be 60 $M_{J}$ brown dwarf companions orbiting around HD 23514 (G0) and HII 1348 (K5) respectively, as had been reported in previous studies. In our observations, the average detection limit for a point source was 20.3 mag in the $H$ band beyond 1''.5 from the central star. On the basis of this detection limit, we calculated the detection efficiency to be 90% for a planet with 6 to 12 Jovian masses and a semi-major axis of 50--1000 AU. For this we extrapolated the distribution of planet mass and semi-major axis derived from RV observations and adopted the planet evolution model of Baraffe et al. (2003). As there was no detection of a planet, we estimated the frequency of such planets to be less than 17.9% ($2\sigma$) around one star of the Pleiades cluster.

Debris disks around young main-sequence stars often have gaps and cavities which for a long time have been interpreted as possibly being caused by planets. In recent years, several giant planet discoveries have been made in systems hosting disks of precisely this nature, further implying that interactions with planets could be a common cause of such disk structures. As part of the SEEDS high-contrast imaging survey, we are surveying a population of debris disk-hosting stars with gaps and cavities implied by their spectral energy distributions, in order to attempt to spatially resolve the disk as well as to detect any planets that may be responsible for the disk structure. Here we report on intermediate results from this survey. Five debris disks have been spatially resolved, and a number of faint point sources have been discovered, most of which have been tested for common proper motion, which in each case has excluded physical companionship with the target stars. From the detection limits of the 50 targets that have been observed, we find that beta Pic b-like planets (~10 Mjup planets around G--A-type stars) near the gap edges are less frequent than 15--30%, implying that if giant planets are the dominant cause of these wide (27 AU on average) gaps, they are generally less massive than beta Pic b.

We report the discovery of a pair of quasars at $z=1.487$, with a separation of $8\farcs585\pm0\farcs002$. Subaru Telescope infrared imaging reveals the presence of an elliptical and a disk-like galaxy located almost symmetrically between the quasars, creating a cross-like configuration. Based on absorption lines in the quasar spectra and the colors of the galaxies, we estimate that both galaxies are located at redshift $z=0.899$. This, as well as the similarity of the quasar spectra, suggests that the system is a single quasar multiply imaged by a galaxy group or cluster acting as a gravitational lens, although the possibility of a binary quasar cannot be fully excluded. We show that the gravitational lensing hypothesis implies these galaxies are not isolated, but must be embedded in a dark matter halo of virial mass $\sim 4 \times 10^{14}\ h_{70}^{-1}\ {M}_\odot$ assuming an NFW model with a concentration parameter of $c_{vir}=6$, or a singular isothermal sphere profile with a velocity dispersion of $\sim 670$ km s$^{-1}$. We place constraints on the location of the dark matter halo, as well as the velocity dispersions of the galaxies. In addition, we discuss the influence of differential reddening, microlensing and intrinsic variability on the quasar spectra and broadband photometry.

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.

We describe Algorithms for Calibration, Optimized Registration, and Nulling the Star in Angular Differential Imaging (ACORNS-ADI), a new, parallelized software package to reduce high-contrast imaging data, and its application to data from the SEEDS survey. We implement several new algorithms, including a method to register saturated images, a trimmed mean for combining an image sequence that reduces noise by up to ~20%, and a robust and computationally fast method to compute the sensitivity of a high-contrast observation everywhere on the field-of-view without introducing artificial sources. We also include a description of image processing steps to remove electronic artifacts specific to Hawaii2-RG detectors like the one used for SEEDS, and a detailed analysis of the Locally Optimized Combination of Images (LOCI) algorithm commonly used to reduce high-contrast imaging data. ACORNS-ADI is written in python. It is efficient and open-source, and includes several optional features which may improve performance on data from other instruments. ACORNS-ADI requires minimal modification to reduce data from instruments other than HiCIAO. It is freely available for download at www.github.com/t-brandt/acorns-adi under a BSD license.

We present Subaru/HiCIAO H-band high-contrast images of the debris disk around HIP 79977, whose pres- ence was recently inferred from an infrared excess. Our images resolve the disk for the first time, allowing characterization of its shape, size, and dust grain properties. We use angular differential imaging (ADI) to reveal the disk geometry in unpolarized light out to a radius of ~2", as well as polarized differential imaging (PDI) to measure the degree of scattering polarization out to ~1.5". In order to strike a favorable balance between suppression of the stellar halo and conservation of disk flux, we explore the application of principal component analysis (PCA) to both ADI and reference star subtraction. This allows accurate forward modeling of the effects of data reduction on simulated disk images, and thus direct comparison with the imaged disk. The resulting best-fit values and well-fitting intervals for the model parameters are a surface brightness power-law slope of S_out = -3.2 [-3.6,-2.9], an inclination of i = 84{\deg} [81{\deg},86{\deg}], a high Henyey-Greenstein forward-scattering parameter of g = 0.45 [0.35, 0.60], and a non-significant disk-star offset of u = 3.0 [-1.5, 7.5] AU = 24 [-13, 61] mas along the line of nodes. Furthermore, the tangential linear polarization along the disk rises from ~10% at 0.5" to ~45% at 1.5". These measurements paint a consistent picture of a disk of dust grains produced by collisional cascades and blown out to larger radii by stellar radiation pressure.

We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micron HST/NICMOS data. While sub-millimeter studies suggested there is a dust-depleted cavity with r=0.35, we find scattered light as close as 0.1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks. We find two small-scaled spiral structures which asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h ~ 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5+3,-4 Mj, in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks, we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0.5. We reach 5-sigma contrasts limiting companions to planetary masses, 3-4 MJ at 1.0 and 2 MJ at 1.55 using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.

We present the direct imaging discovery of an extrasolar planet, or possible low-mass brown dwarf, at a projected separation of 55 +/- 2 AU (1.058 +/- 0.007 arcsec) from the B9-type star Kappa And. The planet was detected with Subaru/HiCIAO during the SEEDS survey, and confirmed as a bound companion via common proper motion measurements. Observed near-infrared magnitudes of J = 16.3 +/- 0.3, H = 15.2 +/- 0.2, Ks = 14.6 +/- 0.4, and L' = 13.12 +/- 0.09 indicate a temperature of ~1700 K. The galactic kinematics of the host star are consistent with membership in the Columba association, implying a corresponding age of 30 +20 -10 Myr. The system age, combined with the companion photometry, points to a model-dependent companion mass ~12.8 MJup. The host star's estimated mass of 2.4-2.5 Msun places it among the most massive stars ever known to harbor an extrasolar planet or low-mass brown dwarf. While the mass of the companion is close to the deuterium burning limit, its mass ratio, orbital separation, and likely planet-like formation scenario imply that it may be best defined as a `Super-Jupiter' with properties similar to other recently discovered companions to massive stars.

We report high-resolution (0.07 arcsec) near-infrared polarized intensity images of the circumstellar disk around the star 2MASS J16042165-2130284 obtained with HiCIAO mounted on the Subaru 8.2 m telescope. We present our $H$-band data, which clearly exhibits a resolved, face-on disk with a large inner hole for the first time at infrared wavelengths. We detect the centrosymmetric polarization pattern in the circumstellar material as has been observed in other disks. Elliptical fitting gives the semimajor axis, semiminor axis, and position angle (P.A.) of the disk as 63 AU, 62 AU, and -14 $^{\circ}$, respectively. The disk is asymmetric, with one dip located at P.A.s of $\sim85^{\circ}$. Our observed disk size agrees well with a previous study of dust and CO emission at submillimeter wavelength with Submillimeter Array. Hence, the near-infrared light is interpreted as scattered light reflected from the inner edge of the disk. Our observations also detect an elongated arc (50 AU) extending over the disk inner hole. It emanates at the inner edge of the western side of the disk, extending inward first, then curving to the northeast. We discuss the possibility that the inner hole, the dip, and the arc that we have observed may be related to the existence of unseen bodies within the disk.

We present the near-infrared images and spectra of four silhouette disks in the Orion Nebula Cluster (ONC; M42) and M43 using the Subaru Adaptive Optics system. While d053-717 and d141-1952 show no water ice feature at 3.1 micron, a moderately deep (tau~0.7) water ice absorption is detected toward d132-1832 and d216-0939. Taking into account the water ice so far detected in the silhouette disks, the critical inclination angle to produce a water ice absorption feature is confirmed to be 65-75deg. As for d216-0939, the crystallized water ice profile is exactly the same as in the previous observations taken 3.63 years ago. If the water ice material is located at 30AU, then the observations suggest it is uniform at a scale of about 3.5AU.

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is ~70 AU. Our data show that the geometric center of the disk shifts by ~6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at ~30 to ~50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius, the sample shows the wedge pattern expected from a radial shell, which is detected clearly here for the first time. This confirms predictions from an earlier model of formation of the GSS, which proposed that the W Shelf is a shell from the third orbital wrap of the same tidal debris stream that produces the GSS, with the main body of the progenitor lying in the second wrap. We calculate the distortions in the shelf wedge pattern expected from its outward expansion and angular momentum, and show that these effects are echoed in the data. In addition, a hot, relatively smooth spheroid population is clearly present. We construct a bulge-disk-halo N-body model that agrees with surface brightness and kinematic constraints, and combine it with a simulation of the GSS. From the contrasting kinematic signatures of the hot spheroid and shelf components, we decompose the observed stellar metallicity distribution into contributions from each component using a non-parametric mixture model. The shelf component's metallicity distribution matches previous observations of the GSS superbly, further strengthening the evidence they are connected and bolstering the case for a massive progenitor of this stream.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

We derived the luminosity function (LF) of dwarf galaxies in the Coma Cluster down to M_R=-10 at three fields located at the center, intermediate, and outskirt. The LF (-19<M_R<-10) shows no significant differences among the three fields. It shows a clear dip at M_R\sim-13, and is composed of two distinct components of different slopes; the bright component with -19<M_R<-13 has a flatter slope than the faint component with -13<M_R<-10 which has a steep slope. The bright component (-19<M_R<-13) consists of mostly red extended galaxies including few blue galaxies whose colors are typical of late-type galaxies. On the other hand, the faint component (-13<M_R<-10) consists of largely PSF-like compact galaxies. We found that both these compact galaxies and some extended galaxies are present in the center while only compact galaxies are seen in the outskirt. In the faint component, the fraction of blue galaxies is larger in the outskirt than in the center. We suggest that the dwarf galaxies in the Coma Cluster, which make up the two components in the LF, are heterogeneous with some different origins.

We study the multi-dimensional geometry of supernova (SN) explosions by means of spectropolarimetric observations of stripped-envelope SNe, i.e., SNe without a H-rich layer. We perform spectropolarimetric observations of 2 stripped-envelope SNe, the Type Ib SN 2009jf and the Type Ic SN 2009mi. Both objects show non-zero polarization at the wavelength of the strong lines. They also show a loop in the Stokes Q-U diagram, which indicates a non-axisymmetric, three-dimensional ion distribution in the ejecta. We show that five out of six stripped-envelope SNe which have been observed spectropolarimetrically so far show such a loop. This implies that a three-dimensional geometry is common in stripped-envelope SNe. We find that stronger lines tend to show higher polarization. This effect is not related to the geometry, and must be corrected to compare the polarization of different lines or different objects. Even after the correction, however, there remains a dispersion of polarization degree among different objects. Such a dispersion might be caused by three-dimensional clumpy ion distributions viewed from different directions.

One of the key predictions of modeling from the IR excess of Herbig Ae stars is that for protoplanetary disks, where significant grain growth and settling has occurred, the dust disk has flattened to the point that it can be partially or largely shadowed by the innermost material at or near the dust sublimation radius. When the self-shadowing has already started, the outer disk is expected to be detected in scattered light only in the exceptional cases that the scale height of the dust disk at the sublimation radius is smaller than usual. High-contrast imaging combined with the IR spectral energy distribution allow us to measure the degree of flattening of the disk, as well as to determine the properties of the outer disk. We present polarimetric differential imaging in $H$ band obtained with Subaru/HiCIAO of one such system, MWC 480. The HiCIAO data were obtained at a historic minimum of the NIR excess. The disk is detected in scattered light from 0\farcs2-1\farcs0 (27.4-137AU). Together with the marginal detection of the disk from 1998 February 24 by HST/NICMOS, our data constrain the opening half angle for the disk to lie between 1.3$\leq\theta\leq 2.2^\circ$. When compared with similar measures in CO for the gas disk from the literature, the dust disk subtends only $\sim$30% of the gas disk scale height (H/R$\sim$0.03). Such a dust disk is a factor of 5-7 flatter than transitional disks, which have structural signatures that giant planets have formed.

We searched for z=7.3 Lya emitters (LAEs) behind two lensing clusters, Abell 2390 and CL 0024, with the Subaru Telescope Suprime-Cam and a narrowband NB1006 (FWHM ~ 21 nm centered at 1005 nm). We investigated if there exist objects consistent with the color of z=7.3 LAEs behind the clusters but could not detect any LAEs to the unlensed line limit F(Lya) ~ 6.9 x 10^{-18} erg/s/cm^2. Using several z=7 Lya luminosity functions (LFs) from the literature, we estimated and compared the expected detection numbers of z ~ 7 LAEs in lensing and blank field surveys in the case of using an 8m class ground based telescope. Given the steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient in detecting z ~ 7 LAEs than imaging a blank field. However, the gain is expected to be modest, a factor of two at most and likely much less depending on the adopted LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain in depth and not necessarily in detection efficiency. For much larger detectors, the lensing effect becomes negligible and the efficiency of LAE detection is proportional to the instrumental FOV. We also inspected NB1006 images of three z ~ 7 z-dropouts previously detected in Abell 2390 and found that none of them are detected in NB1006. Two of them are consistent with predictions from the previous studies that they would be at lower redshifts. The other one has a photometric redshift of z ~ 7.3, and if it is at z=7.3, its unlensed Lya line flux would be very faint: F(Lya) < 4.4 x 10^{-18} erg/s/cm^2 (1 sigma upper limit) or rest frame equivalent width of W(Lya) < 26A. Its Lya emission might be attenuated by neutral hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying strong Lya emission is lower at z ~ 7 than at z <~ 6.

We have performed deep imaging surveys for LyA emitters (LAEs) at redshift ~7.3 in two blank fields, the Subaru Deep Field (SDF) and the Subaru/XMM-Newton Deep survey Field (SXDF), using the Subaru/Suprime-Cam equipped with new red-sensitive CCDs and a new narrow-band filter, NB1006 (lambda_c=10052 Ang, FWHM=214 Ang). We identified four objects as LAE candidates that exhibit luminosity excess in NB1006. By carrying out deep follow-up spectroscopy for three of them using Subaru/FOCAS and Keck/DEIMOS, a definitively asymmetric emission line is detected for one of them, SXDF-NB1006-2. Assuming this line is LyA, this object is a LAE at z=7.215 which has luminosity of 1.2^{+1.5}_{-0.6} x 10^43 [erg s-1] and a weighted skewness S_w=4.90+-0.86. Another object, SDF-NB1006-2, shows variable photometry and is thus probably a quasar (QSO) or an active galactic nucleus (AGN). It shows an asymmetric emission line at 10076 Ang, which may be due to either LyA at z=7.288 or [OII] at z=1.703. The third object, SDF-NB1006-1, is likely a galaxy with temporal luminosity enhancement associated with a supernova explosion, as the brightness of this object varies between the observed epochs. Its spectrum does not show any emission lines. The inferred decrease in the number density of LAEs toward higher redshift is n_LyA(z=7.3)/n_LyA(z=5.7) = 0.05^+0.11_-0.05 from z=5.7 to 7.3 down to L(LyA)=1.0 x 10^43 [erg s-1]. The present result is consistent with the interpretation in previous studies that the neutral hydrogen fraction is rapidly increasing from z=5.7 to 7.3.

Recent observations of z >~ 7 Ly{\alpha} emitters (LAEs) have derived a variety of Ly{\alpha} luminosity functions (LFs) with contradictory results, evolution or non-evolution from z <~ 6, the epoch after reionization. This could be because most of z >~ 7 LFs comprise photometric candidates and might include some contaminations. We conducted the Subaru Telescope Faint Object Camera And Spectrograph narrowband NB980 ({\lambda}c ~ 9800A, FWHM ~ 100A) imaging and spectroscopy survey of z=7-7.1 LAEs to compare its "contamination-free" result with z >~ 7 photometric Ly{\alpha} LFs previously derived. We imaged the Subaru Deep Field and the sky around a cluster MS 1520.1+3002 and found one LAE candidate, but spectroscopy did not reveal Ly{\alpha} though deep enough to detect it. We calculated the expected number of LAEs in our survey, using five z=7 and three z=7.7 Ly{\alpha} LFs from recent surveys. Seven of them are consistent with null detection (0.1^{+1.8}_{-0.1}-1.1^{+2.2}_{-1.0} LAEs) within errors including Poisson statistics and cosmic variance, but average values (0.7-1.1 LAEs) predicted from one z=7 and two z=7.7 LFs among the seven indicate nearly a single detection. The remaining one z=7 LF predicts 3.0^{+3.2}_{-2.0} LAEs. As to z=7, the discrepancy likely comes from different LAE selection criteria. For z=7.7, there are two possibilities; (1) If z=7.7 LAEs are somehow brighter in Ly{\alpha} luminosity than lower redshift LAEs, z=7.7 LF is observed to be similar to or higher than lower redshift LFs even if attenuated by neutral hydrogen. (2) All/most of the z=7.7 candidates are not LAEs. This supports the decline of LF from z ~ 6 to 7.7 and reionization at z ~ 6-7.7.

We present a panoramic narrow-band imaging survey of [OII] emitters in and around the ClG J0218.3-0510 cluster at z=1.62 with Suprime-Cam on Subaru telescope. 352 [OII] emitters were identified on the basis of narrow-band excesses and photometric redshifts. We discovered a huge filamentary structure with some clumps traced by [OII] emitters and found that the ClG J0218.3-0510 cluster is embedded in an even larger super-structure than the one reported previously. 31 [OII] emitters were spectroscopically confirmed with the detection of H-alpha and/or [OIII] emission lines by FMOS observations. In the high density regions such as cluster core and clumps, star-forming [OII] emitters show a high overdensity by a factor of more than 10 compared to the field region. Although the star formation activity is very high even in the cluster core, some massive quiescent galaxies also exits at the same time. Furthermore, the properties of the individual [OII] emitters, such as star formation rates, stellar masses and specific star formation rates, do not show a significant dependence on the local density, either. Such lack of environmental dependence is consistent with our earlier result by Hayashi et al. (2011) on a z=1.5 cluster and its surrounding region. The fact that the star-forming activity of galaxies in the cluster core is as high as that in the field at z~1.6 may suggest that the star-forming galaxies are probably just in a transition phase from a starburst mode to a quiescent mode, and are thus showing comparable level of star formation rates to those in lower density environments. We may be witnessing the start of the reversal of the local SFR--density relation due to the "biased" galaxy formation and evolution in high density regions at high this redshift, beyond which massive galaxies would be forming vigorously in a more biased way in proto-cluster cores.

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

We are developing the Laser Guide Star Adaptive Optics (LGS/AO188) system for the Subaru Telescope at Mauna Kea, Hawaii. This system utilizes a combination of an all-solid-state mode-locked sum-frequency generation (SFG) laser (1.7-GHz bandwidth, 0.7-ns pulse width) as a light source and a single-mode optical fiber for beam transference. However, optical fibers induce nonlinear effects, especially self-phase modulation (SPM). We studied SPM in our photonic crystal fiber (PCF). SPM broadens the spectrum of a laser beam and decrease the efficiency of bright laser guide star generation. We measured the spectrum width using a spectrum analyzer. We found a spectrum width of 8.4 GHz at full width at half maximum (FWHM). The original FWHM of our laser spectrum was 1.4 GHz. This was equivalent to a 70 % loss in laser energy. We also measured the brightness of the sodium cell and evaluated its performance as a function of laser wavelength. The cell's brightness showed a peculiar tendency; specifically, it did not extinguish even though the wavelength varied by more than 5 pm. To reduce the impact of SPM, we developed an optical system that divides one laser pulse into four lower-power pulses. The laser peak power after passing through the new optical system was decreased to one-fourth the original, reducing the impact of SPM on the sodium cell. An actual laser guide star created with the new system was 0.41 mag brighter than the laser guide star created with the original system. We achieved brighter laser guide star generation by dividing a laser pulse to reduce its peak intensity. This is an effective method for laser relay using optical fiber.

Using deep narrow-band and broad-band imaging, we identify 401 z~0.40 and 249 z~0.49 H-alpha line-emitting galaxies in the Subaru Deep Field. Compared to other H-alpha surveys at similar redshifts, our samples are unique since they probe lower H-alpha luminosities, are augmented with multi-wavelength (rest-frame 1000AA--1.5 microns) coverage, and a large fraction (20%) of our samples has already been spectroscopically confirmed. Our spectra allow us to measure the Balmer decrement for nearly 60 galaxies with H-beta detected above 5-sigma. The Balmer decrements indicate an average extinction of A(H-alpha)=0.7^{+1.4}_{-0.7} mag. We find that the Balmer decrement systematically increases with higher H-alpha luminosities and with larger stellar masses, in agreement with previous studies with sparser samples. We find that the SFRs estimated from modeling the spectral energy distribution (SED) is reliable---we derived an "intrinsic" H-alpha luminosity which is then reddened assuming the color excess from SED modeling. The SED-predicted H-alpha luminosity agrees with H-alpha narrow-band measurements over 3 dex (rms of 0.25 dex). We then use the SED SFRs to test different statistically-based dust corrections for H-alpha and find that adopting one magnitude of extinction is inappropriate: galaxies with lower luminosities are less reddened. We find that the luminosity-dependent dust correction of Hopkins et al. yields consistent results over 3 dex (rms of 0.3 dex). Our comparisons are only possible by assuming that stellar reddening is roughly half of nebular reddening. The strong correspondence argue that with SED modeling, we can derive reliable intrinsic SFRs even in the absence of H-alpha measurements at z~0.5.

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (~9 Jupiter masses at 0.5", ~3 Jupiter masses at 1").

An overview on the current status of the census of the early universe population is given. Observational surveys of high redshift objects provide direct opportunities to study the early epoch of the Universe. The target population included are Lyman Alpha Emitters (LAE), Lyman Break Galaxies (LBG), gravitationally lensed galaxies, quasars and gamma-ray bursts (GRB). The basic properties of these objects and the methods used to study them are reviewed. The present paper highlights the fact that the Subaru Telescope group made significant contributions in this field of science to elucidate the epoch of the cosmic dawn and to improve the understanding of how and when infant galaxies evolve into mature ones.

We report a search for z=7 Ly-alpha emitters (LAEs) using a custom-made Narrow-Band filter, centered at 9755 Angstroms, with the instrument Suprime-Cam installed at the Subaru telescope. We observed two different fields and obtained two sample of 7 Ly-alpha emitters of which 4 are robust in each field. We are covering the luminosity range of 9.10^{42} - 2.10^{43} erg/s in comoving volumes of ~ 4.10^{5} and 4.3.10^{5} Mpc^{3}. From this result, we derived possible z~7 Ly-alpha luminosity functions for the full samples and for a subsample of 4 objects in each field. We do not observe, in each case, any strong evolution between the z=6.5 and z~7 Ly-alpha luminosity functions. Spectroscopic confirmation for these candidate samples is required to establish a definitive measure of the luminosity function at z~7.

The quasar SDSS J133401.39+331534.3 at z = 2.426 is found to be a two-image gravitationally lensed quasar with the image separation of 0.833. The object is first identified as a lensed quasar candidate in the Sloan Digital Sky Survey Quasar Lens Search, and then confirmed as a lensed system from follow-up observations at the Subaru and University of Hawaii 2.2-meter telescopes. We estimate the redshift of the lensing galaxy to be 0.557 based on absorption lines in the quasar spectra as well as the color of the galaxy. In particular, we observe the system with the Subaru Telescope AO188 adaptive optics with laser guide star, in order to derive accurate astrometry, which well demonstrates the usefulness of the laser guide star adaptive optics imaging for studying strong lens systems. Our mass modeling with improved astrometry implies that a nearby bright galaxy $\sim 4"$ apart from the lensing galaxy is likely to affect the lens potential.

We carried out extended spectroscopic confirmations of Ly-alpha emitters (LAEs) at z=6.5 and 5.7 in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 at z=6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z=6.5 (5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Ly-alpha luminosity, both photometrically and spectroscopically, to accurately determine the Ly-alpha and rest-UV luminosity functions (LFs). The substantially improved evaluation of the Ly-alpha LF at z=6.5 shows an apparent deficit from z=5.7 at least at the bright end, and a possible decline even at the faint end, though small uncertainties remain. The rest-UV LFs at z=6.5 and 5.7 are in good agreement, at least at the bright end, in clear contrast to the differences seen in the Ly-alpha LF. These results imply an increase in the neutral fraction of the intergalactic medium from z=5.7 to 6.5. The rest-frame equivalent width (EW_0) distribution at z=6.5 seems to be systematically smaller than z=5.7, and it shows an extended tail toward larger EW_0. The bright end of the rest-UV LF can be reproduced from the observed Ly-alpha LF and a reasonable EW_0-UV luminosity relation. Integrating this rest-UV LF provides the first measurement of the contribution of LAEs to the photon budget required for reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon budget among Lyman break galaxies significantly increases towards faint magnitudes. Low-luminosity LAEs could dominate the ionizing photon budget, though this inference depends strongly on the uncertain faint-end slope of the Ly-alpha LF.

Faint Object Camera and Spectrograph (FOCAS) is a versatile common-use optical instrument for the 8.2m Subaru Telescope, offering imaging and spectroscopic observations. FOCAS employs grisms with resolving powers ranging from 280 to 8200 as dispersive optical elements. A grism is a direct-vision grating composed of a transmission grating and prism(s). FOCAS has five grisms with replica surface-relief gratings including an echelle-type grism, and eight grisms with volume-phase holographic (VPH) gratings. The size of these grisms is 110 mm X 106 mm in aperture with a maximum thickness of 110 mm. We employ not only the dichromated gelatin, but also the hologram resin as a recording material for VPH gratings. We discuss the performance of these FOCAS grisms measured in the laboratory, and verify it by test observations, and show examples of astronomical spectroscopic observations.

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.

GJ 758 B is a cold (~600K) companion to a Sun-like star at 29 AU projected separation, which was recently detected with high-contrast imaging. Here we present photometry of the companion in seven photometric bands from Subaru/HiCIAO, Gemini/NIRI and Keck/NIRC2, providing a rich sampling of the spectral energy distribution in the 1-5 micron wavelength range. A clear detection at 1.58 micron combined with an upper limit at 1.69 micron shows methane absorption in the atmosphere of the companion. The mass of the companion remains uncertain, but an updated age estimate indicates that the most likely mass range is ~30-40 Mjup. In addition, we present an updated astrometric analysis that imposes tighter constraints on GJ 758 B's orbit and identifies the proposed second candidate companion, "GJ 758 C", as a background star.

We present Lya luminosity function (LF), clustering measurements, and Lya line profiles based on the largest sample, to date, of 207 Lya emitters (LAEs) at z=6.6 on the 1-deg^2 sky of Subaru/XMM-Newton Deep Survey (SXDS) field. Our z=6.6 Lya LF including cosmic variance estimates yields the best-fit Schechter parameters of phi*=8.5 +3.0/-2.2 x10^(-4) Mpc^(-3) and L*(Lya)=4.4 +/-0.6 x10^42 erg s^(-1) with a fixed alpha=-1.5, and indicates a decrease from z=5.7 at the >~90% confidence level. However, this decrease is not large, only =~30% in Lya luminosity, which is too small to be identified in the previous studies. A clustering signal of z=6.6 LAEs is detected for the first time. We obtain the correlation length of r_0=2-5 h^(-1) Mpc and bias of b=3-6, and find no significant boost of clustering amplitude by reionization at z=6.6. The average hosting dark halo mass inferred from clustering is 10^10-10^11 Mo, and duty cycle of LAE population is roughly ~1% albeit with large uncertainties. The average of our high-quality Keck/DEIMOS spectra shows an FWHM velocity width of 251 +/-16 km s^(-1). We find no large evolution of Lya line profile from z=5.7 to 6.6, and no anti-correlation between Lya luminosity and line width at z=6.6. The combination of various reionization models and our observational results about the LF, clustering, and line profile indicates that there would exist a small decrease of IGM's Lya transmission owing to reionization, but that the hydrogen IGM is not highly neutral at z=6.6. Our neutral-hydrogen fraction constraint implies that the major reionization process took place at z>~7.

We obtained mid-infrared 3.6 and 4.5 micron imaging of a z=6.96 Lyman alpha emitter (LAE) IOK-1 discovered in the Subaru Deep Field, using Spitzer Space Telescope Infrared Array Camera observations. After removal of a nearby bright source, we find that IOK-1 is not significantly detected in any of these infrared bands to m_3.6 ~ 24.00 and m_4.5 ~ 23.54 at 3 sigma. Fitting population synthesis models to the spectral energy distribution consisting of the upper limit fluxes of the optical to infrared non-detection images and fluxes in detection images, we constrain the stellar mass M* of IOK-1. This LAE could have either a mass as low as M* <~ 2-9 x 10^8 Msun for the young age (<~ 10 Myr) and the low dust reddening (A_V ~ 0) or a mass as large as M* <~ 1-4 x 10^{10} Msun for either the old age (> 100 Myr) or the high dust reddening (A_V ~ 1.5). This would be within the range of masses of z ~ 3-6.6 LAEs studied to date, ~ 10^6-10^{10} Msun. Hence, IOK-1 is not a particularly unique galaxy with extremely high mass or low mass but is similar to one of the LAEs seen at the later epochs.

We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts.

We investigate the sky distribution of z ~ 6 Lyman break galaxies selected as i'-dropouts having i' - z' > 1.45 down to z' < 26.5 in the Subaru Deep Field (SDF). We discover 37 i'-dropouts clustered in a projected comoving 21.6 x 21.6 Mpc^2 region at z = 6, showing a local density excess. Carrying out follow-up spectroscopy, we identify four of them as Lyman-alpha emitters at z = 5.92, 6.01, 6.03 and 6.03 (spread over a distance of 46.6 Mpc). The number density of the cluster itself in SDF is ~ 2.2 x 10^{-7} Mpc^{-3}, smaller than those of protoclusters (i.e., forming galaxy clusters) at z ~ 2-5.7. Also, the structure shows ~4-21 times larger galaxy number density than those of z ~ 6 galaxies in a general field. It has a mass of M ~ 1.5^{+1.8}_{-0.5} x 10^{15}M_sun, comparable to those of z ~ 0-5 protoclusters. Since the contamination of our sample by interlopers is estimated to be quite low, 5.9%, most of the other unconfirmed i'-dropouts in the overdense region can be also z ~ 6 galaxies. Hence, it could be a candidate forming cluster at z ~ 6, representing a progenitor of galaxy clusters seen in the recent-day Universe.

We present optical and near infrared observations of GRB 080325 classified as a "Dark GRB". Near-infrared observations with Subaru/MOIRCS provided a clear detection of afterglow in Ks band, although no optical counterpart was reported. The flux ratio of rest-wavelength optical to X-ray bands of the afterglow indicates that the dust extinction along the line of sight to the afterglow is Av = 2.7 - 10 mag. This large extinction is probably the major reason for optical faintness of GRB 080325. The J - Ks color of the host galaxy, (J - Ks = 1.3 in AB magnitude), is significantly redder than those for typical GRB hosts previously identified. In addition to J and Ks bands, optical images in B, Rc, i', and z' bands with Subaru/Suprime-Cam were obtained at about one year after the burst, and a photometric redshift of the host is estimated to be z_{photo} = 1.9. The host luminosity is comparable to L^{*} at z \sim 2 in contrast to the sub-L^{*} property of typical GRB hosts at lower redshifts. The best-fit stellar population synthesis model for the host shows that a large dust extinction (Av = 0.8 mag) attributes to the red nature of the host and that the host galaxy is massive (M_{*} = 7.0 \times 10^{10} Msun) which is one of the most massive GRB hosts previously identified. By assuming that the mass-metallicity relation for star-forming galaxies at z \sim 2 is applicable for the GRB host, this large stellar mass suggests the high metallicity environment around GRB 080325, consistent with inferred large extinction.

We present the first spectropolarimetric observations of a candidate of super-Chandrasekhar mass Type Ia supernova (SN): SN 2009dc. The observations were performed at 5.6 and 89.5 days after the B-band maximum. The data taken at the later epoch are used to determine the interstellar polarization. Continuum polarization is found to be small (<0.3 %), indicating that the explosion is nearly spherically symmetric. This fact suggests that a very aspherical explosion is not a likely scenario for SN 2009dc. Polarization at the Si II and Ca II lines clearly shows a loop in the Q-U plane, indicating a non-axisymmetric, clumpy distribution of intermediate-mass elements. The degree of line polarization at the Si and Ca line is moderate (0.5% +- 0.1% and 0.7% +- 0.1%, respectively), but it is higher than expected from the trend of other Type Ia SNe. This may suggest that there are thick enough, clumpy Si-rich layers above the thick 56Ni-rich layers (>~ 1.2 Msun). The observed spectropolarimetric properties, combined with the photometric and spectroscopic properties, suggest that the progenitor of SN 2009dc has a super-Chandrasekhar mass, and that the explosion geometry is globally spherically symmetric, with clumpy distribution of intermediate-mass elements.

We present a photometric survey of the stellar halo of the Andromeda galaxy, using Suprime-Cam on the Subaru Telescope. A detailed analysis of VI color-magnitude diagrams of the resolved stellar population is used to measure properties such as line-of-sight distance, surface brightness, metallicity, and age, and these are used to isolate and characterize different components of the M31 halo: (1) several substructures, and (2) the smooth halo. First, we study M31's halo substructure along the north-west/south-east minor axis out to R ~ 100 kpc and the south-west major axis region at R ~ 60 kpc. We confirm two substructures in the south-east halo reported by Ibata et al. (2007) and discover two overdense substructures in the north-west halo. We investigate the properties of these four substructures as well as other structures including the western shelf and find that differences in stellar populations among these systems, thereby suggesting each has a different origin. Our statistical analysis implies that the M31 halo as a whole may contain at least 16 substructures, each with a different origin. Second, we investigate the properties of an underlying, smooth and extended halo component out to R > 100 kpc. We find that the surface density of this smooth halo can be fitted to a Hernquist model of scale radius ~ 17 kpc or a power-law profile with ~ R^{-2.17 +/- 0.15}. In contrast to the relative smoothness of the halo density profile, its metallicity distribution appears to be spatially non-uniform with non-monotonic variations with radius, suggesting that the halo population has not had sufficient time to dynamically homogenize the accreted populations. Further implications for the formation of the M31 halo are discussed.

National Astronomical Observatory is an inter-university institute serving as the national center for ground based astronomy offering observational facilities covering the optical, infrared, radio wavelength domain. NAOJ also has solar physics and geo-lunar science groups collaborating with JAXA for space missions and a theoretical group with computer simulation facilities. The outline of NAOJ, its various unique facilities, and some highlights of recent science achievements are reviewed.

We have obtained deep images in the BVRIJHKs bands of the field centered on QSO 1508+5714 (z_{em} =4.28) with the Suprime camera, FOCAS and MOIRCS cameras on Subaru telescope. We report here the detection of a B-dropout galaxy, which is $3\secpoint 5$ north-west of the QSO sightline. A photometric redshift analysis is presented to complement the color selection. Given the photometric properties of this object ($M = -22.2$, making $ L\approx 3 L^{\ast}$, if placed at its photometric redshift $z\sim 3.5$), as well as the S$\acute{e}$rsic index ($ n \sim 1$) derived from a 2-D imaging decomposition of the HST WFPC2 image taken in the $I_{F814}$ filter, the identified system is consistent with a massive disk galaxy at z>3. If confirmed, it would be one of the most distant massive disk galaxies known so far.

We present an optical spectropolarimetric observation of the unique Type Ib supernova (SN) 2005bf at 8 days after the second maximum. The data, combined with the polarization spectrum taken at 6 days before the second maximum (Maund et al. 2007a), enable us to closely examine the intrinsic properties of the SN. The polarization percentage is smaller at the later epoch over a wide wavelength range, while the position angle is similar at the two epochs. We find that an assumption of complete depolarization of strong lines at the emission peak is not necessarily correct. The intrinsic polarization of the SN is larger, and thus, the ejecta of SN 2005bf would be more asymmetric than previously expected. The axis ratio of the photosphere projected on the sky deviates from unity by at least 20 %. If the position angle of interstellar polarization is aligned with the spiral arm of the host galaxy, the deviation is larger than 25 %. The line polarization at the He I, Ca II and Fe II lines is also smaller at the later epoch. The changes in the position angle across these lines, which were observed at the earlier epoch, are still marginally present at the later epoch. The similar polarimetric behavior suggests that the distributions of these ions are correlated. Properties of polarization, as well as the light curve and the spectra both at photospheric and nebular phases, can be explained by an aspherical, possibly unipolar explosion of a WN star in which the blob of 56Ni penetrates C+O core and stops within the He layer.

We present new measurements of the diameter of o Cet (Mira) as a function of wavelength in the 2.2 micron atmospheric window using the adaptive optics system and the infrared camera and spectrograph mounted on the Subaru Telescope. We found that the angular size of the star at the wavelengths of CO and H2O absorption lines were up to twice as large as the continuum photosphere. This size difference is attributable to the optically thick CO and H2O molecular layers surrounding the photosphere. This measurement is the first direct differential spectroscopic imaging of stellar extension that resolves individual molecular lines with high spectral-resolution observations. This observation technique is extremely sensitive to differences in spatial profiles at different wavelengths; we show that a difference in diameter much smaller than the point spread function can be measured.

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

We examine variability of absorption line strength of intervening systems along the line of sight to GRB 060206 at $z = 4.05$, by the low-resolution optical spectra obtained by the Subaru telescope from six to ten hours after the burst. Strong variabilities of Fe\emissiontype{II} and Mg\emissiontype{II} lines at $z=1.48$ during $t=$5--8 hours have been reported for this GRB \citep{Hao07}, and this has been used to support the idea of clumpy Mg\emissiontype{II} cloudlets that was originally proposed to explain the anomalously high incidence of Mg\emissiontype{II} absorbers in GRB spectra compared with quasars . However, our spectra with higher signal-to-noise ratio do not show any evidence for variability in $t=$6-10 hours. There is a clear discrepancy between our data and Hao et al. data in the overlapping time interval. Furthermore, the line strengths in our data are in good agreement with those observed at $t \sim$ 2 hours by \citet{Tho08}. We also detected Fe\emissiontype{II} and Mg\emissiontype{II} absorption lines for a system at $z = 2.26$, and these lines do not show evidence for variability either. Therefore we conclude that there is no strong evidence for variability of intervening absorption lines toward GRB 060206, significantly weakening the support to the Mg\emissiontype{II} cloudlet hypothe sis by the GRB 060206 data.

A brief overview on the current status of the census of the early universe population is given. Observational surveys of high redshift galaxies provide direct opportunities to witness the cosmic dawn and to have better understanding of how and when infant galaxies evolve into mature ones. It is a much more astronomical approach in contrast to the physical approach of to study the spatial fluctuation of cosmic microwave radiation. Recent findings in these two areas greatly advanced our understanding of the early Universe. I will describe the basic properties of several target objects we are looking for and the concrete methods astronomers are using to discover those objects in early Universe. My talk starts with Lyman alpha emitters and Lyman break galaxies, then introduces a clever approach to use gravitational lensing effect of clusters of galaxies to detect distant faint galaxies behind the clusters. Finally I will touch on the status and prospects of surveys for quasars and gamma-ray bursts.

Using the HST ACS, we have obtained deep optical images reaching well below the oldest main sequence turnoff in fields on the southeast minor-axis of the Andromeda Galaxy, 35 kpc from the nucleus. These data probe the star formation history in the extended halo of Andromeda -- that region beyond 30 kpc that appears both chemically and morphologically distinct from the metal-rich, highly-disturbed inner spheroid. The present data, together with our previous data for fields at 11 and 21 kpc, do not show a simple trend toward older ages and lower metallicities, as one might expect for populations further removed from the obvious disturbances of the inner spheroid. Specifically, the mean ages and [Fe/H] values at 11 kpc, 21 kpc, and 35 kpc are 9.7 Gyr and -0.65, 11.0 Gyr and -0.87, and 10.5 Gyr and -0.98, respectively. In the best-fit model of the 35 kpc population, one third of the stars are younger than 10 Gyr, while only ~10% of the stars are truly ancient and metal-poor. The extended halo thus exhibits clear evidence of its hierarchical assembly, and the contribution from any classical halo formed via early monolithic collapse must be small.

We present luminosity functions (LFs) and various properties of Lya emitters (LAEs) at z=3.1, 3.7, and 5.7, in a 1 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field. We obtain a photometric sample of 858 LAE candidates based on deep Subaru/Suprime-Cam imaging data, and a spectroscopic sample of 84 confirmed LAEs from Subaru/FOCAS and VLT/VIMOS spectroscopy in a survey volume of ~10^6 Mpc^3 with a limiting Lya luminosity of ~3x10^42 erg/s. We derive the LFs of Lya and UV-continuum (~1500 \AA) for each redshift, taking into account the statistical error and the field-to-field variation. We find that the apparent Lya LF shows no significant evolution between z=3.1 and 5.7 within factors of 1.8 and 2.7 in L* and phi*, respectively. On the other hand, the UV LF of LAEs increases from z=3.1 to 5.7, indicating that galaxies with Lya emission are more common at earlier epochs. We identify six LAEs with AGN activities from our spectra combined with VLA, Spitzer, and XMM-Newton data. Among the photometrically selected LAEs at z=3.1 and 3.7, only ~1 % show AGN activities, while the brightest LAEs with logL(Lya) >~ 43.4-43.6 erg/s appear to always host AGNs. Our LAEs are bluer in UV-continuum color than dropout galaxies, suggesting lower extinction and/or younger stellar populations. Our stacking analyses provide upper limits to the radio luminosity and the f(HeII)/f(Lya) line fraction, and constrain the hidden star formation (+low-luminosity AGN) and the primordial population in LAEs.

We have examined in detail the morphologies of seven z~1.5 passively evolving luminous red galaxies using high resolution HST NICMOS and ACS imaging data. Almost all of these galaxies appear to be relaxed systems, with smooth morphologies at both rest-frame UV and visible wavelengths. Previous results from spectral synthesis modeling favor a single burst of star formation more than 1 Gyr before the observed epoch. The prevalence of old stellar populations, however, does not correlate exclusively with early-type morphologies as it does in the local universe; the light profiles for some of these galaxies appear to be dominated by massive exponential disks. This evidence for massive old disks, along with the apparent uniformity of stellar age across the disk, suggests formation by a mechanism better described as a form of monolithic collapse than as a hierarchical merger. These galaxies could not have undergone a single major merging event since the bulk of their stars were formed, more than 1 Gyr earlier. There is at least one case, however, that appears to be undergoing a "dry merger", which may be an example of the process that converts these unusual galaxies into the familiar spheroids that dominate galaxies comprising old stellar populations at the present epoch.

We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7.

Core-collapse supernovae (CC-SNe) are the explosions that announce the death of massive stars. Some CC-SNe are linked to long-duration gamma-ray bursts (GRBs) and are highly aspherical. One important question is to what extent asphericity is common to all CC-SNe. Here we present late-time spectra for a number of CC-SNe from stripped-envelope stars, and use them to explore any asphericity generated in the inner part of the exploding star, near the site of collapse. A range of oxygen emission-line profiles is observed, including a high incidence of double-peaked profiles, a distinct signature of an aspherical explosion. Our results suggest that all CC-SNe from stripped-envelope stars are aspherical explosions and that SNe accompanied by GRBs exhibit the highest degree of asphericity.

We present multi-waveband optical imaging data obtained from observations of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, centered at R.A.=02:18:00, decl.=-05:00:00, has been the focus of a wide range of multi-wavelength observing programs spanning from X-ray to radio wavelengths. A large part of the optical imaging observations are carried out with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Subaru Telescope Observatory Projects. This paper describes our optical observations, data reduction and analysis procedures employed, and the characteristics of the data products. A total area of 1.22 sqdeg is covered in five contiguous sub-fields, each of which corresponds to a single Suprime-Cam field of view (34'x27'), in five broad-band filters B, V, Rc, i', z' to the depths of B=28.4, V=27.8, Rc=27.7, i'=27.7 and z'=26.6 (AB, 3-sigma, 2-arcsec aperture). The data are reduced and compiled into five multi-waveband photometric catalogs, separately for each Suprime-Cam pointing. The i'-band catalogs contain about 900,000 objects, making the SXDS catalogs one of the largest multi-waveband catalogs in corresponding depth and area coverage. The SXDS catalogs can be used for an extensive range of astronomical applications such as the number density of the Galactic halo stars to the large scale structures at the distant universe. The number counts of galaxies are derived and compared with those of existing deep extragalactic surveys. The optical data, the source catalogs, and configuration files used to create the catalogs are publicly available via the SXDS web page (http://www.naoj.org/Science/SubaruProject/SXDS/index.html)

The results of a world-wide coordinated observational campaign on the broad-lined Type Ic SN 2003jd are presented. In total, 74 photometric data points and 26 spectra were collected using 11 different telescopes. SN 2003jd is one of the most luminous SN Ic ever observed. A comparison with other Type Ic supernovae (SNe Ic) confirms that SN 2003jd represents an intermediate case between broad-line events (2002ap, 2006aj), and highly energetic SNe (1997ef, 1998bw, 2003dh, 2003lw), with an ejected mass of M_{ej} = 3.0 +/- 1 Mo and a kinetic energy of E_{k}(tot) = 7_{-2}^{+3} 10^{51} erg. SN 2003jd is similar to SN 1998bw in terms of overall luminosity, but it is closer to SNe 2006aj and 2002ap in terms of light-curve shape and spectral evolution. The comparison with other SNe Ic, suggests that the V-band light curves of SNe Ic can be partially homogenized by introducing a time stretch factor. Finally, due to the similarity of SN 2003jd to the SN 2006aj/XRF 060218 event, we discuss the possible connection of SN 2003jd with a GRB.

We present the results of NICMOS imaging of two massive galaxies photometrically selected to have old stellar populations at z ~ 2.5. Both galaxies are dominated by apparent disks of old stars, although one of them also has a small bulge comprising about 1/3 of the light at rest-frame 4800 A. The presence of massive disks of old stars at high redshift means that at least some massive galaxies in the early universe have formed directly from the dissipative collapse of a large mass of gas. The stars formed in disks like these may have made significant contributions to the stellar populations of massive spheroids at the present epoch.

The nebular spectra of the broad-lined, SN 1998bw-like Type Ic SN 2002ap are studied by means of synthetic spectra. Two different modelling techniques are employed. In one technique, the SN ejecta are treated as a single zone, while in the other a density and abundance distribution in velocity is used from an explosion model. In both cases, heating caused by gamma-ray and positron deposition is computed (in the latter case using a Monte Carlo technique to describe the propagation of gamma-rays and positrons), as is cooling via forbidden-line emission. The results are compared, and although general agreement is found, the stratified models are shown to reproduce the observed line profiles much more accurately than the single-zone model. The explosion produced ~ 0.1 Msun of 56Ni. The distribution in velocity of the various elements is in agreement with that obtained from the early-time models, which indicated an ejected mass of ~ 2.5 Msun with a kinetic energy of 4 x 10^{51} erg. Nebular spectroscopy confirms that most of the ejected mass (~ 1.2 Msun) was oxygen. The presence of an oxygen-rich inner core, combined with that of 56Ni at high velocities as deduced from early-time models, suggests that the explosion was asymmetric, especially in the inner part.

We list and analyze H I absorption lines at redshifts 2 < z < 4 with column density (12 < log(N_HI) < 19) in 40 high-resolutional (FWHM = 8.0 km/s) quasar spectra obtained with the Keck+HIRES. We de-blend and fit all H I lines within 1,000 km/s of 86 strong H I lines whose column densities are log(N_HI/[cm^-2]) > 15. Unlike most prior studies, we use not only Lya but also all visible higher Lyman series lines to improve the fitting accuracy. This reveals components near to higher column density systems that can not be seen in Lya. We list the Voigt profile fits to the 1339 H I components that we found. We examined physical properties of H I lines after separating them into several sub-samples according to their velocity separation from the quasars, their redshift, column density and the S/N ratio of the spectrum. We found two interesting trends for lines with 12 < log(N_HI) < 15 which are within 200-1000 km/s of systems with log(N_HI) > 15. First, their column density distribution becomes steeper, meaning relatively fewer high column density lines, at z < 2.9. Second, their column density distribution also becomes steeper and their line width becomes broader by about 2-3 km/s when they are within 5,000 km/s of their quasar.

We report the serendipitous findings of 13 faint meteors and 44 artificial space objects by Subaru SuprimeCam imaging observations during 11-16 August 2004. The meteors, at about 100km altitude, and artificial satellites/debris in orbit, at 500km altitude or higher, were clearly discriminated by their apparent defocused image sizes. CCD photometry of the 13 meteors, including 1 Perseid, 1 Aquarid, and 11 sporadic meteors, was performed. We defined a peak video-rate magnitude by comparing the integrated photon counts from the brightest portion of the track traversed within 33ms to those from a 0-mag star during the same time duration. This definition gives magnitudes in the range 4.0< V_{vr} <6.4 and 4.1< I_{vr}<5.9 for these 13 meteors. The corresponding magnitude for virtual naked-eye observers could be somewhat fainter especially for the V-band observation, in which the [OI] 5577 line lasting about 1 sec as an afterglow could contribute to the integrated flux of the present 5-10 min CCD exposures. Although the spatial resolution is insufficient to resolve the source size of anything smaller than about 1 m, we developed a new estimate of the collisionally excited column diameter of these meteors. A diameter as small as a few mm was derived from their collisionally excited photon rates, meteor speed, and the volume density of the oxygen atoms at the 100km altitude. The actual column diameter of the radiating zone, however, could be as large as few 100m because the excited atoms travel that distance before they emit forbidden lines in 0.7 sec of its average lifetime. Among the 44 artificial space objects, we confirmed that 17 were cataloged satellites/space debris.

We have used Suprime-Cam on the Subaru Telescope to conduct a V- and I-band imaging survey of fields sampling the spheroid of the Andromeda galaxy along its south-east minor axis. Our photometric data are deep enough to resolve stars down to the red clump. Based on a large and reliable sample of red giant stars available from this deep wide-field imager, we have derived metallicity distributions vs. radius and a surface brightness profile over projected distances of R=23-66 kpc from the galaxy's center. The metallicity distributions across this region shows a clear high mean metallicity and a broad distribution ([Fe/H] ~ -0.6 +/- 0.5), and indicates no metallicity gradient within our observed range. The surface brightness profile at R>40 kpc is found to be flatter than previously thought. It is conceivable that this part of the halo samples as yet unidentified, metal-rich substructure.

We study the galaxy population in the central region and a region about 0.6 Mpc away from the center of the Hydra I cluster in B- and Rc-bands down to M~-10 using the Subaru Suprime-Cam photometry. We find that the luminosity function of the entire population has a slightly steeper slope (alpha~-1.6) in the range of -20<M<-10 than those reported for other clusters in slightly brighter ranges. The slope appears to be steeper in poorer clusters. The number of faint galaxies (M>-14) increases in the cluster center as well as bright galaxies. The Hydra I cluster is dominated by red dwarfs and the luminosity function shows a significant upturn at M~-16 as is seen in several other nearby clusters, but not in the field. This upturn and the variation in the faint-end slope of the luminosity function may be caused by the cluster environment or the evolution history of individual clusters.

We present 3 um spectroscopy of the carbon-rich proto-planetary nebulae IRAS 04296+3429 and IRAS 05341+0852 conducted with the adaptive optics system at the Subaru Telescope. We utilize the nearly diffraction-limited spectroscopy to probe the spatial extent of the hydrocarbon dust emitting zone. We find a hydrocarbon emission core extending up to 100--160 mas from the center of IRAS 04296+3429, corresponding to a physical diameter of 400--640 AU, assuming a distance of 4 kpc. On the other hand, we find that IRAS 05341+0852 is not spatially resolved with this instrumentation. The physical extent of these proto-planetary nebulae, along with the reanalyzed data of IRAS 22272+5435 published previously, suggests a correlation between the physical extent of the hydrocarbon dust emission and the spectral evolution of the aliphatic to aromatic features in these post-AGB stars. These measurements represent the first direct test of the proposed chemical synthesis route of carbonaceous dust in the circumstellar environment of evolved stars.

The gamma-ray burst (GRB) 050904 at z = 6.3 provides the first opportunity of probing the intergalactic medium (IGM) by GRBs at the epoch of the reionization. Here we present a spectral modeling analysis of the optical afterglow spectrum taken by the Subaru Telescope, aiming to constrain the reionization history. The spectrum shows a clear damping wing at wavelengths redward of the Lyman break, and the wing shape can be fit either by a damped Ly alpha system with a column density of log (N_HI/cm^{-2}) ~ 21.6 at a redshift close to the detected metal absorption lines (z_metal = 6.295), or by almost neutral IGM extending to a slightly higher redshift of z_{IGM,u} ~ 6.36. In the latter case, the difference from z_metal may be explained by acceleration of metal absorbing shells by the activities of the GRB or its progenitor. However, we exclude this possibility by using the light transmission feature around the Ly beta resonance, leading to a firm upper limit of z_{IGM,u} < 6.314. We then show an evidence that the IGM was largely ionized already at z=6.3, with the best-fit neutral fraction of IGM, x_HI = N_HI/N_H = 0.00, and upper limits of x_HI < 0.17 and 0.60 at 68 and 95% C.L., respectively. This is the first direct and quantitative upper limit on x_HI at z > 6. Various systematic uncertainties are examined, but none of them appears large enough to change this conclusion. To get further information on the reionization, it is important to increase the sample size of z >~ 6 GRBs, in order to find GRBs with low column densities (log N_HI <~ 20) within their host galaxies, and for statistical studies of Ly alpha line emission from host galaxies.

Abridged: Subaru Deep Field line-emitting galaxies in four narrow-band filters at low and intermediate redshifts are presented. Broad-band colors, follow-up optical spectroscopy, and multiple NB filters are used to distinguish Ha, [O II], and [O III] emitters between redshifts of 0.07 and 1.47 to construct their luminosity functions (LFs). These LFs are derived down to faint magnitudes, which allows for a more accurate determination of the faint end slope. With a large (N~200-900) sample for each redshift interval, a Schechter profile is fitted to each LF. Prior to dust extinction corrections, the [O III] and [O II] LFs reported in this paper agree reasonably well with those of Hippelein et al. The z=0.08 Ha LF, which reaches two orders of magnitude fainter than Gallego et al., is steeper by 25%. This indicates that there are more low luminosity star-forming galaxies for z<0.1. The faint end slope \alpha and \phi_{\star} show a strong evolution with redshift while L_{\star} show little evolution. The evolution in \alpha indicates that low-luminosity galaxies have a stronger evolution compared to brighter ones. Integrated star formation rate densities are derived via Ha, [O III], and [O II] for 0.07<z<1.47. A steep increase in the star-formation rate density, as a function of redshift, is seen for 0.4<z<0.9. For z>1, the star-formation rate densities are more or less constant. The latter is consistent with previous UV and [O II] measurements. Below z<0.4, the SFR densities are consistent with several Ha, [O II], and UV measurements, but other measurements are a factor of two higher. For example, the z=0.066-0.092 Ha LF agrees with those of Jones & Bland-Hawthorn, but at z=0.24 and 0.40, their number density is higher by a factor of two. This discrepancy can be explained by cosmic variance.

When galaxy formation started in the history of the Universe remains unclear. Studies of the cosmic microwave background indicate that the Universe, after initial cooling (following the Big Bang), was reheated and reionized by hot stars in newborn galaxies at a redshift in the range 6 < z < 14 (ref. 1). Though several candidate galaxies at redshift z > 7 have been identified photometrically (refs 2,3), galaxies with spectroscopically confirmed redshifts have been confined to z < 6.6 (refs. 4-8). Here we report a spectroscopic redshift of z = 6.96 (corresponding to just 750 Myr after the Big Bang) for a galaxy whose spectrum clearly shows Lyman-alpha emission at 9,682 A, indicating active star formation at a rate of about 10 M_sun/yr, where M_sun us the mass of the Sun. This demonstrates that galaxy formation was under way when the Universe was only about 6 per cent of its present age. The number density of galaxies at z = 7 seems to be only 18-36 per cent of the density at z = 6.6.

We report an extensive search for Lyman-alpha emitters (LAEs) at z=6.5 in the Subaru Deep Field. Subsequent spectroscopy with Subaru and Keck identified eight more LAEs, giving a total of 17 spectroscopically confirmed LAEs at z=6.5. Based on this spectroscopic sample of 17, complemented by a photometric sample of 58 LAEs, we have derived a more accurate Lyman-alpha luminosity function of LAEs at z=6.5, which reveals an apparent deficit at the bright end of ~0.75 mag fainter L*, compared with that observed at z=5.7. The difference in the LAE luminosity functions between z=5.7 and 6.5 is significant at the 3-sigma level, which is reduced to 2-sigma when cosmic variance is taken into account. This result may imply that the reionization of the universe has not been completed at z=6.5. We found that the spatial distribution of LAEs at z=6.5 was homogeneous over the field. We discuss the implications of these results for the reionization of the universe.

We present new results of [Fe II] 1.644-micron spectroscopy toward the jets from HL Tau and RW Aur carried out with the Subaru Telescope combined with the adaptive optics system. We observed the regions within 2" - 3" from the stars with the sub-arcsecond resolutions of 0."5 and 0."2 for HL Tau and RW Aur, respectively. In addition to the strong, high velocity emission extended along each jet, we detected a blueshifted low velocity emission feature seen as a wing or shoulder of the high velocity emission at each stellar position. Detailed analysis shows that the position-velocity diagrams (PVDs) of HL Tau and RW Aur show a characteristic similar to those of the cold disk wind and X-wind models in that the [Fe II] line width is broad in the vicinity of the stellar position and is narrower at the extended jet. A closer comparison suggests, however, that the disk wind model tends to have too large line width at the jet while the X-wind model has excess emission on the redshifted side at the stellar position. The narrow velocity width with symmetric line profiles of the observed high velocity emission supports an X-wind type model where the launching region is localized in a small radial range, while the low velocity emission located away from the star favors the presence of a disk wind. The [Fe II] emission from the HL Tau jet shows a gap of 0."8 between the redshifted jet and the star, indicating the presence of an optically thick disk of ~ 160 AU in radius. The [Fe II] emission from the RW Aur jet shows a marked drop from the redshifted peak at Y ~ -0."2 toward the star, suggesting that its disk radius is smaller than 40 AU.

The Telescope to Observe Planetary Systems (TOPS) is a proposed space mission to image in the visible (0.4-0.9 micron) planetary systems of nearby stars simultaneously in 16 spectral bands (resolution R~20). For the ~10 most favorable stars, it will have the sensitivity to discover 2 R_E rocky planets within habitable zones and characterize their surfaces or atmospheres through spectrophotometry. Many more massive planets and debris discs will be imaged and characterized for the first time. With a 1.2m visible telescope, the proposed mission achieves its power by exploiting the most efficient and robust coronagraphic and wavefront control techniques. The Phase-Induced Amplitude Apodization (PIAA) coronagraph used by TOPS allows planet detection at 2 lambda/d with nearly 100% throughput and preserves the telescope angular resolution. An efficient focal plane wavefront sensing scheme accurately measures wavefront aberrations which are fed back to the telescope active primary mirror. Fine wavefront control is also performed independently in each of 4 spectral channels, resulting in a system that is robust to wavefront chromaticity.

A Near-infrared (1.18-1.35 micron) high-resolution spectrum of the gravitationally-lensed QSO APM 08279+5255 was obtained with the IRCS mounted on the Subaru Telescope using the AO system. We detected strong NaI D 5891,5897 doublet absorption in high-redshift DLAs at z=1.062 and 1.181, confirming the presence of NaI, which was first reported for the rest-frame UV NaI 3303.3,3303.9 doublet by Petitjean et al. This is the first detection of NaI D absorption in a high-redshift (z>1) DLA. In addition, we detected a new NaI component in the z=1.062 DLA and four new components in the z=1.181 DLA. Using an empirical relationship between NaI and HI column density, we found that all "components" have large HI column density, so that each component is classified as DLA absorption. We also detected strong NaI D absorption associated with a MgII system at z=1.173. Because no other metal absorption lines were detected in this system at the velocity of the NaI absorption in previously reported optical spectra (observed 3.6 years ago), we interpret this NaI absorption cloud probably appeared in the line of sight toward the QSO after the optical observation. This newly found cloud is likely to be a DLA based upon its large estimated HI column density. We found that the N(NaI)/N(CaII) ratios in these DLAs are systematically smaller than those observed in the Galaxy; they are more consistent with the ratios seen in the Large Magellanic Cloud. This is consistent with dust depletion generally being smaller in lower metallicity environments. However, all five clouds of the z=1.181 system have a high N(NaI)/N(CaII) ratio, which is characteristic of cold dense gas. We tentatively suggest that the host galaxy of this system may be the most significant contributor to the gravitational-lens toward APM 08279+5255.

We present the properties of Ly Alpha emitters (LAEs) at z = 5.7 in the Subaru Deep Field. A photometric sample of 89 LAE candidates is constructed from narrow-band (NB816) data down to NB816 = 26.0 (AB) in a continuous 725 arcmin^2 area. Spectra of 39 objects satisfying the photometric selection criteria for LAEs were obtained with Subaru and Keck II Telescopes, among which 28 were confirmed LAEs, one was a nearby galaxy, and eight were unclassified. We also obtained spectra of another 24 NB816-excess objects in the field, identifying six additional LAEs. We find that the Ly Alpha luminosity function derived from the photometric sample is reproduced well by a Schechter function with L* = (7.9+3.0-2.2) x 10^42 erg/s and phi* = (6.3+3.0-2.0) x 10^-4 Mpc^-3 for alpha = -1.5 (fixed) over the whole luminosity range of L ~= 3x10^42 - 3x10^43 erg/s. We then measure rest-frame Ly Alpha equivalent widths (EWs) for the confirmed LAEs, to find that the median among the 28 objects satisfying the photometric selection criteria is W_0^i = 233 A. We infer that 30% - 40% of LAEs at z=5.7 exceed W_0^i = 240 A. These large-EW objects probably cannot be accounted for by ordinary star-forming populations with a Salpeter IMF. We also find that LAEs with fainter far-UV luminosities have larger EWs. Finally, we derive the far-UV luminosity function of LAEs down to M_UV ~= -19.6 using the photometric sample, and compare it with that of Lyman-break galaxies (LBGs). We find that as high as about 80% of LBGs at z ~ 6 have W_0^i >= 100 A, in sharp contrast to lower-z counterparts.

We explored the clustering properties of Lyman Break Galaxies (LBGs) at z=4 and 5 with an angular two-point correlation function on the basis of the very deep and wide Subaru Deep Field data. We found an apparent dependence of the correlation function slope on UV luminosity for LBGs at both z=4 and 5. More luminous LBGs have a steeper correlation function. To compare these observational results, we constructed numerical mock LBG catalogs based on a semianalytic model of hierarchical clustering combined with high-resolution N-body simulation, carefully mimicking the observational selection effects. The luminosity functions for LBGs predicted by this mock catalog were found to be almost consistent with the observation. Moreover, the overall correlation functions of LBGs were reproduced reasonably well. The observed dependence of the clustering on UV luminosity was not reproduced by the model, unless subsamples of distinct halo mass were considered. That is, LBGs belonging to more massive dark haloes had steeper and larger-amplitude correlation functions. With this model, we found that LBG multiplicity in massive dark halos amplifies the clustering strength at small scales, which steepens the slope of the correlation function. The hierarchical clustering model could therefore be reconciled with the observed luminosity-dependence of the angular correlation function, if there is a tight correlation between UV luminosity and halo mass. Our finding that the slope of the correlation function depends on luminosity could be an indication that massive dark halos hosted multiple bright LBGs (abridged).

Gamma-ray bursts (GRBs) and their afterglows have been proposed as an excellent probe to study the evolution of cosmic star formation, the reionization of the intergalactic medium, and the metal enrichment history of the universe, since the prompt gamma-ray emission of GRBs should be detectable out to distances z>10. Hitherto, the highest measured redshift for a GRB has been z=4.50. Here we report the optical spectrum of the afterglow of GRB 050904 obtained 3.4 days after the burst. The spectrum shows a clear continuum at the long wavelength end of the spectrum with a sharp cutoff at around 9000 A due to Ly alpha absorption at a redshift of 6.3 with a damping wing. Little flux is present in the waveband shortward of the Ly alpha break. A system of absorption lines of heavy elements at redshift z=6.295 +- 0.002 were also detected, yielding a precise measurement of the largest known redshift of a GRB. Analysis of the Si II fine structure lines suggest a dense metal-enriched environment around the GRB progenitor, providing unique information on the properties of the gas in a galaxy when the universe was younger than one billion years.

We report angular correlation function (ACF) of Lyman Break Galaxies (LBGs) with unprecedented statistical quality on the basis of 16,920 LBGs at z=4 detected in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field. The ACF significantly departs from a power law, and shows an excess on small scale. Particularly, the ACF of LBGs with i'<27.5 have a clear break between the small and large-scale regimes at the angular separation of ~7'' whose projected length corresponds to the virial radius of dark halos with a mass of 10^11-12 Mo, indicating multiple LBGs residing in a single dark halo. Both on small (2''<theta<3'') and large (40''<theta<400'') scales, clustering amplitudes monotonically increase with luminosity for the magnitude range of i'=24.5-27.5, and the small-scale clustering shows a stronger luminosity dependence than the large-scale clustering. The small-scale bias reaches b~10-50, and the outskirts of small-scale excess extend to a larger angular separation for brighter LBGs. The ACF and number density of LBGs can be explained by the cold dark matter model.

Observations and modeling for the light curve (LC) and spectra of supernova (SN) 2005bf are reported. This SN showed unique features: the LC had two maxima, and declined rapidly after the second maximum, while the spectra showed strengthening He lines whose velocity increased with time. The double-peaked LC can be reproduced by a double-peaked $^{56}$Ni distribution, with most $^{56}$Ni at low velocity and a small amount at high velocity. The rapid post-maximum decline requires a large fraction of the $\gamma$-rays to escape from the $^{56}$Ni-dominated region, possibly because of low-density ``holes''. The presence of Balmer lines in the spectrum suggests that the He layer of the progenitor was substantially intact. Increasing $\gamma$-ray deposition in the He layer due to enhanced $\gamma$-ray escape from the $^{56}$Ni-dominated region may explain both the delayed strengthening and the increasing velocity of the He lines. The SN has massive ejecta ($\sim6-7\Msun$), normal kinetic energy ($\sim 1.0-1.5\times 10^{51}$ ergs), high peak bolometric luminosity ($\sim 5\times 10^{42}$ erg s$^{-1}$) for an epoch as late as $\sim$ 40 days, and a large $^{56}$Ni mass ($\sim0.32\Msun$). These properties, and the presence of a small amount of H suggest that the progenitor was initially massive (M$\sim 25-30 \Msun$) and had lost most of its H envelope, and was possibly a WN star. The double-peaked $^{56}$Ni distribution suggests that the explosion may have formed jets that did not reach the He layer. The properties of SN 2005bf resemble those of the explosion of Cassiopeia A.

We carried out a Subaru and UKIRT near infrared imaging survey for H-alpha emitting galaxies around two pair quasar systems (Q0301-005/Q0302-003 and Q2343+125/Q2344+125), and a triple quasar system (KP76/KP77/KP78). Narrow band near infrared filters covering the H-alpha emission expected for galaxies at the confirmed C IV absorption redshift toward the quasar systems were used for this survey. These quasar pairs or triplet are separated at most by 17 arcmins (~5 h^-1 Mpc in proper distance) from each other on the sky, and have common C IV absorption lines at almost identical redshifts at z=2.24-2.43, which suggests there could be a Mpc-scale absorbing systems such as a cluster, or a group, of galaxies that cover all the line-of-sights to the pair/triple quasars. Using narrow-band deep images, we detected five candidates for H-alpha emitting galaxies around two of the six fields, Q2343+125 and Q2344+125, whose apparent star formation rates are, extremely high, 20-466 M_solar/year. However, all or most of them are not likely to be galaxies at the absorption redshift but galaxies at lower redshift, because of their extreme brightness. In the fields of the other quasars, we detected no star-forming galaxies, nor did we find any number excess of galaxy counts around them. This no-detection results could be because the luminosities and star formation rates of galaxies are lower than the detection limits of our observations (K' > 21 and SFR < 1.8-240 h^-2 M_solar/year). They could be located outside of the observed fields by chance. Otherwise, most C IV absorption lines could be ascribed not to cluster of galaxies, but to isolated star forming pockets far from bright galaxies and could be analogous objects to weak Mg II absorbers.

We obtained the number counts and the rest-frame B-band luminosity function of the color-selected old passively-evolving galaxies (OPEGs) at z=1 with very high statistical accuracy using a large and homogeneous sample of about 4000 such objects with z' <25 detected in the area of 1.03 deg^2 in the Subaru/XMM-Newton Deep Survey (SXDS) field. Our selection criteria are defined on the i'-z' and R-z' color-magnitude plane so that OPEGs at z=0.9-1.1 with formation redshift z_f=2-10 are properly sampled. The limiting magnitude corresponds to the luminosity of galaxies with M_*+3 at z=0. We made a pilot redshift observations for 99 OPEG candidates with 19 < z' < 22 and found that at least 78% (73/93) of the entire sample, or 95% (73/77) of these whose redshifts were obtained are indeed lie between z=0.87 and 1.12 and the most of their spectra show the continuum break and strong Ca H and K lines, indicating that these objects are indeed dominated by the old stellar populations. We then compare our results with the luminosity functions of the color- or the morphologically-selected early type galaxies at z=0 taking the evolutionary factor into account and found that the number density of old passive galaxies with sim M_* magnitude at z~1 averaged over the SXDS area is 40-60% of the equivalently red galaxies and 60-85% of the morphologically-selected E/S0 galaxies at z=0 depending on their luminosity evolution. It is revealed that more than half, but not all, of the present-day early-type galaxies had already been formed into quiescent passive galaxies at z=1.

We report new spectroscopy of two star-forming galaxies with strong Ly_alpha emission at z=6.03 and z=6.04 in the Subaru Deep Field. These two objects are originally selected as i'-dropouts (i'-z' > 1.5) showing an interesting photometric property, the ``NB921 depression''. The NB921-band (centered at 9196A) magnitude is significantly depressed with respect to the z'-band magnitude. The optical spectra of these two objects exhibit asymmetric emission-lines at lambda_obs ~ 8540A and ~ 8560A, suggesting that these objects are Ly_alpha emitters at z~6. The rest-frame equivalent widths of the Ly_alpha emission of the two objects are 94A and 236A; the latter one is the Ly_alpha emitter with the largest Ly_alpha equivalent width at z > 6 ever spectroscopically confirmed. The spectroscopically measured Ly_alpha fluxes of these two objects are consistent with the interpretation that the NB921 depression is caused by the contribution of the strong Ly_alpha emission to the z'-band flux. Most of the NB921-depressed i'-dropout objects are thought to be strong Ly_alpha emitters at 6.0 < z < 6.5; Galactic L and T dwarfs and NB921-dropout galaxies at z > 6.6 do not dominate the NB921-depressed i'-dropout sample. Thus the NB921-depression method is very useful for finding high-z Ly_alpha emitters with a large Ly_alpha equivalent width over a large redshift range, 6.0 < z < 6.5. Although the broadband-selected sample at z ~ 3 contains only a small fraction of objects with a Ly_alpha equivalent width larger than 100A, the i'-dropout sample of the Subaru Deep Field contains a much larger fraction of such strong Ly_alpha emitters. This may imply a strong evolution of the Ly_alpha equivalent width from z > 6 to z ~ 3.

We present a deep $K^{\prime}$-band (2.12$\mu$m) imaging of 1\arcmin\ $\times$ 1\arcmin Subaru Super Deep Field (SSDF) taken with the Subaru adaptive optics (AO) system. Total integration time of 26.8 hours results in the limiting magnitude of $K^{\prime} \sim 24.7$ (5$\sigma$, 0\farcs2 aperture) for point sources and $K^{\prime} \sim 23.5$ (5$\sigma$, 0\farcs6 aperture) for galaxies, which is the deepest limit ever achieved in the $K^{\prime}$ band. The average stellar FWHM of the co-added image is 0\farcs18. Based on the photometric measurements of detected galaxies, we obtained the differential galaxy number counts, for the first time, down to $K^{\prime} \sim 25$, which is more than 0.5 mag deeper than the previous data. We found that the number count slope $d\log N/dm$ is about 0.15 at $22 < K^{\prime} < 25$, which is flatter than the previous data. Therefore, detected galaxies in the SSDF have only negligible contribution to the near-infrared extragalactic background light (EBL), and the discrepancy claimed so far between the diffuse EBL measurements and the estimated EBL from galaxy count integration has become more serious . The size distribution of detected galaxies was obtained down to the area size of less than 0.1 arcsec$^2$, which is less than a half of the previous data in the $K^{\prime}$ band. We compared the observed size-magnitude relation with a simple pure luminosity evolution model allowing for intrinsic size evolution, and found that a model with no size evolution gives the best fit to the data. It implies that the surface brightness of galaxies at high redshift is not much different from that expected from the size-luminosity relation of present-day galaxies.

Type Ic supernovae, the explosions following the core collapse of massive stars that have previously lost their hydrogen and helium envelopes, are particularly interesting because of the link with long-duration gamma-ray bursts. Although indications exist that these explosions are aspherical, direct evidence has still been missing. Late-time observations of SN 2003jd, a luminous Type Ic supernova, provide such evidence. Recent Subaru and Keck spectra reveal double-peaked profiles in the nebular lines of neutral oxygen and magnesium. These profiles are different from those of known Type Ic supernovae, with or without a gamma-ray burst, and they can be understood if SN 2003jd was an aspherical, axisymmetric explosion viewed from near the equatorial plane. If SN 2003jd was associated with a gamma-ray burst, we missed the burst as it was pointing away from us.

A new optical design concept of telescopes to provide an aberration-free, wide field, unvignetted flat focal plane is described. The system employs three aspheric mirrors to remove aberrations, and provides a semi-circular field of view with minimum vignetting. The third mirror reimages the intermediate image made by the first two-mirror system with a magnification factor on the order of unity. The present system contrasts with the Korsch system where the magnification factor of the third mirror is usually much larger than unity. Two separate optical trains can be deployed to cover the entire circular field, if necessary.

We present simulations of quasar colors, magnitudes, and numbers at redshifts 5<z<10 based on our discovery of ten new high-redshift quasars and the cloning of lower redshift Sloan Digital Sky Survey (SDSS) quasars. The ten quasars have redshifts ranging from z=4.7 to z=5.3 and i-magnitudes of 20.21 to 20.94. The natural diversity of spectral features in the cloned sample allows more realistic simulation of the quasar locus width than previously possible with synthetic template spectra. Colors are generated for the z>6 epoch taking advantage of the new UKIDSS near-infrared filter set, and we examine the redshift intervals of maximum productivity, discussing color selection and survey depth issues. On the basis of the SDSS sample, we find that the surface density of z>4.7 quasars increases by a factor of 3X by extending 0.7 i-magnitudes deeper than the SDSS spectroscopic survey limit of i=20.2 -- correspondingly we predict a total of ~400 faint quasars in the SDSS main area that have redshift z>4.7 and magnitudes $<20.9.

We report the discovery of primeval large-scale structures (LSSs) including two proto-clusters in a forming phase at z=5.7. We carried out extensive deep narrow-band imaging in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field, and obtained a cosmic map of 515 Ly alpha emitters (LAEs) in a volume with a transverse dimension of 180 Mpc x 180 Mpc and a depth of ~40 Mpc in comoving units. This cosmic map shows filamentary LSSs, including clusters and surrounding 10-40 Mpc scale voids, similar to the present-day LSSs. Our spectroscopic follow-up observations identify overdense regions in which two dense clumps of LAEs with a sphere of 1-Mpc diameter in physical units are included. These clumps show about 130 times higher star-formation rate density, mainly due to a large overdensity, ~80, of LAEs. These clumps would be clusters in a formation phase involving a burst of galaxy formation.

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy survey over a blank field as large as 34'x27'. The program consists of very deep multi-band optical imaging, near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific goals of the project are to construct large samples of Lyman-break galaxies at z~4-5 and Lyman alpha emitters at z~5.7 and 6.6, and to make detailed studies these very high-redshift galaxy populations. In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked images and object catalogs in seven bandpasses.The optical imaging was made through five broad-band filters, B, V, R, i', z', and two narrow-band filters, NB816 (lambda_c=8150A) and NB921 (lambda_c=9196A) with almost 10 hours long integrations for each band. The limiting magnitudes measured at 3-sigma on a 2" aperture are B=28.45, V=27.74, R=27.80, i'=27.43, z'=26.62, NB816=26.63, and NB921=26.54 in the AB system. The object catalog constructed for each of the seven bands contains more than 10^5 objects. The galaxy number counts corrected for detection incompleteness and star count contribution are found to be consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/.

Using H I absorption alone, we attempt to separate H I absorption lines in quasar spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density Lines), and we discuss the difference in their physical properties. We deblend and fit all H I lines with Voigt profiles, and make an unbiased sample of H I lines covering a wide column density range (12 < log N_HI < 19 cm^-2). To reduce the influence of line blending, we simultaneously fit several Lyman series lines. As a result of a two-point correlation analysis, we found that higher column density H I lines are clustering at $\Delta$ v < 200 km/s, while lower ones at $\Delta$ v < 100 km/s. We define HDLs as H I lines with 15 < log N_HI < 19 cm^-2 and all H I lines within +/- 200 km/s of a line with log N_HI > 15 cm^-2, and LDLs as others with 12 < log N_HI < 15 cm^-2. We found that the HDLs have smaller minimum b-values for a given column density than the LDLs. This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler dense gas near to galaxies.

We report on the discovery of a star-forming galaxy at z=6.33 in the Subaru Deep Field. This object is selected as a candidate of an i'-dropout, high-redshift galaxy around z=6 because of its red i'-z' color in our deep optical imaging survey in the Subaru Deep Field. Our follow up optical spectroscopy reveals that this object is a strong Ly-alpha emitter with only very faint ultraviolet continuum. The rest-frame equivalent width of the detected Ly-alpha emission is as much as 130 A. Thus the light detected in our z' image is largely attributed to the Ly-alpha emission, i.e., ~40% of the z'-band flux is the strong Ly-alpha emission, giving a very red i'-z' color. This is consistent with the photometric property of this object because the narrow-band data obtained with the NB921 filter shows a significant depression, z'-NB921 = -0.54 mag. By using the photometric data, we show that some other objects among the 48 i'-dropout high-redshift galaxy candidates found in the Subaru Deep Field also show a significant NB921 depression. We briefly discuss the nature of these NB921-depressed objects.

An overview of the current status of the 8.2 m Subaru Telescope constructed and operated at Mauna Kea, Hawaii, by the National Astronomical Observatory of Japan is presented. The basic design concept and the verified performance of the telescope system are described. Also given are the status of the instrument package offered to the astronomical community, the status of operation, and some of the future plans. The status of the telescope reported in a number of SPIE papers as of the summer of 2002 are incorporated with some updates included as of 2004 February. However, readers are encouraged to check the most updated status of the telescope through the home page, http://subarutelescope.org/index.html, and/or the direct contact with the observatory staff.

Optical spectroscopy of the Type Ia supernova SN 2002ic obtained on 2003 June 27.6 UT, i.e., about 222 rest-frame days after explosion, is presented. Strong H emission indicates an interaction between the expanding SN ejecta and an H-rich circumstellar medium (CSM). The spectrum of SN 2002ic resembles those of SNe 1997cy and 1999E. The three SNe also have similar luminosities, suggesting that they are the same phenomenon and that the CSM is also similar. We propose a new classification, Type IIa SNe, for these events. The observed line profiles and line ratios are measured and discussed within the ejecta-CSM interaction scenario. The emission in H Balmer, [O III], and He I lines, and in permitted Fe II blends, resembles the spectra of the Type IIn SN 1987F and of Seyfert 1 galaxies. A high-density, clumpy CSM is inferred. Strong, very broad [Ca II]/Ca II and [O I]/O I missions imply that not all the outer SN ejecta were decelerated in the interaction, suggesting that the CSM is aspherical.

An H$_2$ emission filament is found in close proximity to the unique object KH 15D using the adaptive optics system of the Subaru Telescope. The morphology of the filament, the presence of spectroscopic outflow signatures observed by Hamilton et al., and the detection of extended H$_2$ emission from KH 15D by Deming, Charbonneau, & Harrington suggest that this filament arises from shocked H$_2$ in an outflow. The filament extends about 15" to the north of KH 15D.

Adaptive optics imaging of the bright visual T Tauri binary AS 353 with the Subaru Telescope shows that it is a hierarchical triple system. The secondary component, located 5.6" south of AS 353A, is resolved into a subarcsecond binary, AS 353Ba and Bb, separated by 0.24". Resolved spectroscopy of the two close components shows that both have nearly identical spectral types of about M1.5. Whereas AS 353A and Ba show clear evidence for an infrared excess, AS 353Bb does not. We discuss the possible role of multiplicity in launching the large Herbig-Haro flow associated with AS 353A.

We report here the results of deep optical spectroscopy of the very extended emission-line region (VEELR) found serendipitously around the Seyfert 2 galaxy NGC 4388 in the Virgo cluster using the Subaru Telescope. The H-alpha recession velocities of most of the filaments of the region observed are highly blue-shifted with respect to the systemic velocity of the galaxy. The velocity field is complicated, and from the kinematic and morphological points of view, there seem to be several streams of filaments: low velocity (v = -100 km/s) filaments, high velocity (v = -300 km/s) filaments, and a very high velocity (v > -500 km/s) cloud complex. The emission-line ratios of the VEELR filaments are well explained by power-law photoionization models with solar abundances, suggesting that the Seyfert nucleus of NGC 4388 is the dominant ionization source of the VEELR and that the VEELR gas has moderate metallicity. In addition to photoionization, shock heating probably contributes to the ionization of the gas. In particular, the filaments outside the ionization cone of the Seyfert nucleus are mainly excited by shocks. We conclude that the VEELR was formerly the disk gas of NGC 4388, which has been stripped by ram pressure due to the interaction between the hot intra-cluster medium (ICM) and the galaxy. The velocity field and the morphology of the VEELR closely resemble snapshots from some numerical simulations of this process. In the case of NGC 4388, the ram pressure-stripped gas, which is normally seen as extended HI filaments, happens to be exposed and ionized by the radiation from the AGN, and so can be seen as optical emission-line gas.

We performed deep K'-band imaging observations of 2 massive clusters, MS 0451.6-0305 at z = 0.55 and MS 0440.5+0204 at z = 0.19, for searching counterparts of the faint sub-mm sources behind these clusters, which would provide one of the deepest extremely red object(ERO) samples. Comparing our near-infrared images with optical images taken by the Hubble Space Telescope and by the Subaru Telescope, we identified 13 EROs in these fields. The sky distributions of EROs are consistent with the previous results, that there is a sign of strong clustering among detected EROs. Also, the surface density with corrected lensing amplification factors in both clusters are in good agreement with that derived from previous surveys. We found 7 EROs and 3 additional very red objects in a small area (\sim 0.6 arcmin^2) of the MS 0451.6-0305 field around an extended SCUBA source. Many of their optical and near-infrared colors are consistent with dusty star-forming galaxies at high redshifts(z \sim 1.0-4.0), and they may be constituting a cluster of dusty starburst galaxies and/or lensed star-forming galaxies at high redshift. Their red J-K' colors and faint optical magnitudes suggest they are relatively old massive stellar systems with ages(>300 Mega years) suffering from dust obscuration. We also found a surface-density enhancement of EROs around the SCUBA source in the MS 0440.5+0204 field.

Spectroscopic and spectropolarimetric observations of SN 2003dh/GRB 030329 obtained in 2003 May using the Subaru 8.2 m telescope are presented. The properties of the SN are investigated through a comparison with spectra of the Type Ic hypernovae SNe 1997ef and 1998bw. (Hypernovae being a tentatively defined class of SNe with very broad absorption features: these features suggest a large velocity of the ejected material and possibly a large explosion kinetic energy.) Comparison with spectra of other hypernovae shows that the spectrum of SN 2003dh obtained on 2003 May 8 and 9, i.e., 34-35 rest-frame days after the GRB (for z=0.1685), are similar to those of SN 1997ef obtained ~34-42 days after the fiducial time of explosion of that SN. The match with SN 1998bw spectra is not as good (at rest 7300-8000 A, but again spectra obtained ~33-43 days after GRB 980425 are preferred. This indicates that the SN may have intermediate properties between SNe 1997ef and 1998bw. Based on the analogy with the other hypernovae, the time of explosion of SN 2003dh is then constrained to be between -8 and +2 days of the GRB. The Si and O P-Cygni lines of SN 2003dh seem comparable to those of SN 1997ef, which suggests that the ejected mass in SN 2003dh may match that in SN 1997ef. Polarization was marginally detected at optical wavelengths. This is consistent with measurements of the late afterglow, implying that it mostly originated in the interstellar medium of the host galaxy.

Several high-z (z > 5.7) quasars have been found in the course of Sloan Digital Sky Survey. The presence of such very high-z quasars is expected to give constraints on early structure formation. On one hand, it is suggested that these most luminous objects at high redshift are biased toward the highly magnified objects by gravitational lensing. To clarify the effect of gravitational lensing on the high-z quasars, we began the imaging survey of intervening lensing galaxies. Indeed our previous optical image showed that SDSSp J104433.04+012502.2 at z=5.74 is gravitationally magnified by a factor 2. In this paper, we report our new optical imaging of other two high-z quasars, SDSSp J103027.10+052455.0 at z=6.28 and SDSSp J130608.26+035626.3 at z=5.99. Since we find neither intervening galaxy nor counter image with i^{\prime} < 25.4-25.8 around each quasar, we conclude that they are not strongly magnified regardless that a lens galaxy is dusty.

We present deep optical spectroscopy of an extended Ly$\alpha$ emission-line blob located in an over-dense region at redshift $z \approx 3.1$; `blob 1' of Steidel et al. (2000). The origin of such Ly$\alpha$ blobs has been debated for some time; two of the most plausible models are (1) that it comes from a dust-enshrouded, extreme starburst galaxy with a large-scale galactic outflow (superwind/hyperwind) or (2) that it is the cooling radiation of proto-galaxies in dark matter halos. Examination of the kinematic properties of the Ly$\alpha$ emission-line gas should allow us to determine its nature. With this motivation, we performed optical spectroscopy of `blob 1' using the Subaru Telescope, and found that its kinematic properties can be well explained in terms of superwind activity.

We present further observations of the extremely red object ERO J094258+4659.2, identified by \citet{iye00} as ERO R1 in their deep images of the cluster A851. We estimate its redshift independently by eight-band photometric redshift determination and cross-correlation of a new H-band spectrum with the optical spectra of local E/S0 galaxies, and conclude that it lies at $z \sim 1.5$. Although its colors are consistent both with an elliptical galaxy and an S0 galaxy at that redshift, its elongated shape and exponential luminosity profile suggest the presence of an evolved stellar disk component. We rule out the possibility that these properties are strongly influenced by gravitational lensing by the foreground cluster, and therefore conclude that this object is more likely to be an S0-like galaxy, rather than a lensed elliptical. The H-band spectrum does not show strong H$\alpha$ emission and the star formation rate therefore appears to be very modest. The presence of such a galaxy with an apparently relaxed disk of stars at this high redshift provides a new and strong constraint on theoretical models which aim to explain the formation and evolution of galaxies.

We present medium-resolution 3 um spectroscopy of the carbon-rich proto-planetary nebula IRAS 22272+5435. Spectroscopy with the Subaru Telescope adaptive optics system revealed a spatial variation of hydrocarbon molecules and dust surrounding the star. The ro-vibrational bands of acetylene (C2H2) and hydrogen cyanide (HCN) at 3.0 um are evident in the central star spectra. The molecules are concentrated in the compact region near the center. The 3.3 and 3.4 um emission of aromatic and aliphatic hydrocarbons is detected at 600--1300 AU from the central star. The separation of spatial distribution between gas and dust suggests that the small hydrocarbon molecules are indeed the source of solid material, and that the gas leftover from the grain formation is being observed near the central star. The intensity of aliphatic hydrocarbon emission relative to the aromatic hydrocarbon emission decreases with distance from the central star. The spectral variation is well matched to that of a laboratory analog thermally annealed with different temperatures. We suggest that either the thermal process after the formation of a grain or the variation in the temperature in the dust-forming region over time determines the chemical composition of the hydrocarbon dust around the proto-planetary nebula.

We present a high resolution spectrum of the quasar, HS1603+3820 (z_em=2.542), observed with the High Dispersion Spectrograph (HDS) on Subaru Telescope. This quasar, first discovered in the Hamburg/CfA Quasar Survey, has 11 C IV lines at 1.96 < z_abs < 2.55. Our spectrum covers 8 of the 11 C IV lines at z_abs > 2.29 and resolves some of them into multiple narrow components with b < 25 km/s because of the high spectral resolution R=45000, while other lines show broad profiles (b > 65 km/s). We use three properties of C IV lines, specifically, time variability, covering factor, and absorption line profile, to classify them into quasar intrinsic absorption lines (QIALs) and spatially intervening absorption lines (SIALs). The C IV lines at 2.42 < z_abs < 2.45 are classified as QIALs in spite of their large velocity shifts from the quasar. Perhaps they are produced by gas clouds ejected from the quasar with the velocity of v_ej = 8000 km/s -- 10000 km/s. On the other hand, three C IV lines at 2.48 < z_abs < 2.55 are classified as SIALs, which suggests there exist intervening absorbers near the quasar. We, however, cannot rule out QIALs for the two lines at z_abs ~ 2.54 and 2.55, because their velocity shifts, 430 km/s blueward and 950 km/s redward of the quasar, are very small. The C IV line at z_abs ~ 2.48 consists of many narrow components, and has also corresponding low-ionization metal lines (Al II, Si II, and Fe II). The velocity distribution of these low-ionization ions is concentrated at the center of the system compared to that of high-ionization C IV ion. Therefore we ascribe this system of absorption lines to an intervening galaxy.

We present spectropolarimetry of the Type Ic supernova SN 2002ap and give a preliminary analysis: the data were taken at two epochs, close to and one month later than the visual maximum (2002 February 8). In addition we present June 9 spectropolarimetry without analysis. The data show the development of linear polarization. Distinct polarization profiles were seen only in the O I \lambda 7773 multiplet/Ca II IR triplet absorption trough at maximum light and in the Ca II IR triplet absorption trough a month later, with the latter showing a peak polarization as high as ~2 %. The intrinsic polarization shows three clear position angles: 80 degs for the February continuum, 120 degs for the February line feature, and 150 degs for the March data. We conclude that there are multiple asymmetric components in the ejecta. We suggest that the supernova has a bulk asymmetry with an axial ratio projected on the sky that is different from 1 by of order 10 %. Furthermore, we suggest very speculatively that a high velocity ejecta component moving faster than ~0.115c (e.g., a jet) contributes to polarization in the February epoch.

We have constructed a very deep K'-selected multicolor BVRIz'JK' sample of 439 field galaxies. Photometric redshifts for sample galaxies were estimated. The overall redshift distribution N(z) for the K'<21.0 sample is consistent with previous observations, and for the first time we derive the N(z) down to K'=24.0. After taking account of the dust extinction and selection effects of the sample, the observed N(z) distribution is well described with the PLE model, while the hierarchical galaxy formation model shows an apparent deficiency of galaxies especially at z>2. The photometric redshift and the best-fit SED model evaluations allow the derivation of the rest-frame K', B, and UV(2000A) luminosity functions (LFs) and their evolutions. The rest-K' LF shows almost no evolution up to z=3, while the rest B LF shows mild luminosity evolution, and the rest UV LF shows strong luminosity evolution. This trend corresponds to the evolution of the rest UV-K' color, which gets bluer with increasing redshift. We also found that more massive galaxies are redder in this rest-frame color in every epoch. The rest-UV LF of our K'-selected galaxies shows a much shallower faint end slope at z=3 than that of previous estimations for rest-UV selected Lyman break galaxies. As a consequence, the contribution to the global star formation rate of our K'-selected galaxies is about 42% of that derived from the integration of LF of Lyman break galaxies at z=3. This result suggests that a large fraction of the star formation rate density at z>1.5 comes from the contribution from the faint (M_2000A>-20) blue galaxy population at high redshift universe that have not yet obviously been identified. (abridged)

We present new optical images ($B$, $V$, and H$\alpha$) of the archetypical starburst/superwind galaxy M82 obtained with the 8.2 m Subaru Telescope to reveal new detailed structures of the superwind-driven nebula and the high-latitude dark lanes. The emission-line nebula is decomposed into (1) a ridge-dominated component comprising numerous filament/loop sub-structures whose overall morphology appears as a pair of narrow cylinders, and (2) a diffuse component extended over much wider opening angle from the nucleus. We suggest that these two components have different origins. The ridge-dominated component appears as a pair of cylinders rather than a pair of cones. Since this morphological property is similar to that of hot plasma probed by soft X-ray, this component seems to surround the hot plasma. On the other hand, the diffuse component may arise from dust grains which scatter stellar light from the galaxy. Since inner region of this component is seen over the prominent ^^ ^^ X"-shaped dark lanes streaming out from the nuclear region and they can be reproduced as a conical distribution of dust grains, there seems to be a dusty cold outflow as well as the hot one probed by soft X-ray and shock-excited optical emission lines. If this is the case, the presence of such high-latitude dust grains implies that neutral gaseous matter is also blown out during the course of the superwind activity.

During the course of our deep optical imaging survey for Ly alpha emitters at z approximately 5.7 in the field around the z=5.74 quasar SDSSp J104433.04-012502.2, we have found a candidate strong emission-line source. Follow-up optical spectroscopy shows that the emission line profile of this object is asymmetric, showing excess red-wing emission. These properties are consistent with an identification of Ly alpha emission at a redshift of z=5.687 +/- 0.002. The observed broad line width, Delta V_{FWHM} ~= 340 km s^{-1} and excess red-wing emission also suggest that this object hosts a galactic superwind.

We studied the relation between the ratio of rotational velocity to velocity dispersion and the metallicity (<v_{rot}>/\sigma_{v}-metallicity relation) of globular cluster systems (GCS) of disk galaxies by comparing the relation predicted from simple chemo-dynamical models for the formation and evolution of disk galaxies with the observed kinematical and chemical properties of their GCSs. We conclude that proto disk galaxies underwent a slow initial collapse that was followed by a rapid contraction and derive that the ratio of the initial collapse time scale to the active star formation time scale is \sim 6 for our Galaxy and \sim 15 for M31. The fundamental formation process of disk galaxies was simulated based on simple chemo-dynamical models assuming the conservation of their angular momentum. We suggest that there is a typical universal pattern in the <v_{rot}>/\sigma_{v}-metallicity relation of the GCS of disk galaxies. This picture is supported by the observed properties of GCSs in the Galaxy and in M31. This relation would deviate from the universal pattern, however, if large-scale merging events took major role in chemo-dynamical evolution of galaxies and will reflect the epoch of such merging events. We discuss the properties of the GCS of M81 and suggest the presence of past major merging event.

Photometric and spectroscopic data of the energetic Type Ic supernova (SN) 2002ap are presented, and the properties of the SN are investigated through models of its spectral evolution and its light curve. The SN is spectroscopically similar to the "hypernova" SN 1997ef. However, its kinetic energy [$\sim (4-10) \times 10^{51}$ erg] and the mass ejected (2.5-5 $M_{\odot}$) are smaller, resulting in a faster-evolving light curve. The SN synthesized $\sim 0.07 M_{\odot}$ of $^{56}$Ni, and its peak luminosity was similar to that of normal SNe. Brightness alone should not be used to define a hypernova, whose defining character, namely very broad spectral features, is the result of a high kinetic energy. The likely main-sequence mass of the progenitor star was 20-25 $M_{\odot}$, which is also lower than that of both hypernovae SNe 1997ef and 1998bw. SN 2002ap appears to lie at the low-energy and low-mass end of the hypernova sequence as it is known so far. Observations of the nebular spectrum, which is expected to dominate by summer 2002, are necessary to confirm these values.

We have obtained a high resolution (R~90,000) spectrum of the carbon- and s-process-element-rich, very metal-poor ([Fe/H]=-2.7) subgiant LP625-44, as well as that of HD140283 (a metal-poor subgiant with normal abundance ratio) for comparison, with the High Dispersion Spectrograph (HDS) for the Subaru Telescope for detailed abundance study. The oxygen abundance derived from the O I triplet around 7770A is uncertain, but the excess of oxygen in LP625-44 seems remarkable (perhaps by nearly a factor 10), in comparison with that of HD140283 derived from the same lines. The Na enhancement in LP625-44 is by about a factor 50, suggesting hydrogen burning in the 22Ne-rich layer in an asymptotic giant branch star which produces the abundance pattern of this object. In our new spectrum of LP625-44, the Pb I lambda 3683A line has been detected, as well as the Pb I lambda 4057A line which has already been studied, confirming the Pb abundance (log epsilon(Pb)~1.9) derived by the previous work. The abundance ratio of s-process elements at the second peak (e.g., La, Ce and Nd) to that at the third peak (Pb) in LP 625-44 is significantly higher (by a factor 5) than that in other three s-process element-rich objects recently studied by van Eck et al.. Recent theoretical works have modeled the s-process nucleosynthesis in the radiative layer of asymptotic giant branch stars in the inter-pulse phase, and the above results means that these processes produced a large scatter in the abundance ratios. Another possiblity is that different processes (e.g., s-process nucleosynthesis during thermal pulses) have contributed to heavy elements in the early Galaxy.

We present new optical spectroscopy of the lens elliptical galax in the ``Einstein Cross'' lens system HST 14176+5226, using the Faint Object Camera and Spectrograph (FOCAS) of the Subaru t Our spectroscopic observations are aimed at measuring the stella dispersion of the lens galaxy, located at high redshift of z_L= as an important component to lens models. We have measured this 230 +- 14 km s^{-1} (1 sigma) inside 0.35 effective radi based on the comparison between the observed galaxy spectrum and templates of three G-K giants by means of the Fourier cross-corr To extract the significance of this information on the geometry universe which also affects the lensing of the background image, to fit three different lens models to the available data of the Provided that the lens galaxy has the structural and dynamical p (i.e., its radial density profile, core radius, and velocity ani similar to those of local elliptical galaxies, we calculate the function for the simultaneous reproduction of both the observed and newly measured velocity dispersion of the lens. Although the interval depends rather sensitively on the adopted lens models o parameters, our experiments suggest the larger likelihood for a cosmological constant, Omega_Lambda: formal 1 sigma lower Omega_Lambda in the flat universe ranges 0.73 to 0.97, where lower limit is basically unavailable. This method for determinin model is thus dependent on lens models but is insensitive to oth ambiguities, such as the dust absorption or the evolutionary eff galaxies. Exploring spectroscopic observations of more lens gala redshift may minimize the model uncertainties and thus place a m constraint on Omega_Lambda.

We present near-infrared spectroscopy of low-mass companions in a nearby triple system HD 130948 (Gliese 564, HR 5534). Adaptive optics on the Subaru Telescope allowed spectroscopy of the individual components of the 0".13 binary system. Based on a direct comparison with a series of template spectra, we determined the spectral types of HD 130948B and C to be L4 +- 1. If we take the young age of the primary star into account (0.3-0.8 Gyr), HD 130948B and C most likely are a binary brown dwarf system.

We found a very large, 35 kpc, emission-line region around the Seyfert type 2 galaxy NGC 4388, using deep narrow-band imaging with the prime focus camera (Suprime-Cam) of the Subaru telescope. This region consists of many faint gas clouds or filaments, and extends northeastwards from the galaxy. The typical H-alpha luminosity L(H-alpha) of the clouds is 10^37 erg s^-1, and the total L(H-alpha) of the region within 10 kpc from the nucleus is 2 x 10^38 erg s^-1, which corresponds to an ionized gas mass of 10^5 M_solar. The map of the emission-line intensity ratio I([OIII])/I(H-alpha) indicates that the inner (r < 12 kpc) region of the VEELR may be excited by nuclear ionizing radiation. The excitation mechanism of the outer (r > 12 kpc) region is unclear, but it is likely that the nuclear radiation is also a dominant source of its ionization. We discuss the origin of the ionized gas. Two plausible origins of the gas in the VEELR are (i) the tidal debris resulting from a past interaction with a gas-rich dwarf galaxy, i.e., a minor merger, or (ii) the interstellar medium of NGC 4388, stripped by the ram pressure of the hot intracluster medium of the Virgo cluster.

We investigate the radial dependence of the spectroscopic properties, in particular the Mg2, <Fe> and H beta spectroscopic indices, in a sample of galaxies spanning a wide range of absolute luminosity in the Coma cluster. After allowing for the magnitude dependence of these indices, we find a significant gradient in Mg2, in the sense that galaxies in the core of the cluster have stronger Mg2. We find only weak gradients in <Fe> and H beta. Using the model grids presented in an earlier paper in this series, we attribute the Mg2 gradient to changes in metal abundance. One possible mechanism to create this abundance gradient is pressure confinement by the intracluster medium of material from Supernova driven winds early in the history of the galaxies.

We discuss the properties of galaxies around the radio galaxy 3C 324 at z=1.2 based on BVRIK' multi-band imaging data. We have applied a photometric-redshift technique to objects in the 3C 324 field, and identified 35 objects as plausible cluster members. We have found that red and luminous members are concentrated in a small region enclosed by a circle of 40" radius (0.33 Mpc at z=1.2 for \Omega_0=0.3, \lambda_0=0.7, H_0=70km/s/Mpc cosmology) from the 3C 324 galaxy. The 3C 324 cluster is probably much more compact in size compared with the local clusters. We constructed a K'-band luminosity function of the cluster members and fit a Schechter function, and found the characteristic magnitude to be K*_{AB}=20.2+-0.6. This value is consistent with the extrapolation of the pure passive evolution seen for z<1 clusters. We have identified eight bright galaxies which form a red color-magnitude sequence. The slope of the sequence is consistent with the passive evolution model down to K'_{AB}<22; we also found that there is no clear age variation in these bright red galaxies. However, seven out of these eight galaxies exhibit a significant excess in the rest UV light with respect to the passive evolution model. This may suggest that the massive early-type galaxies in this high-redshift cluster are still forming stars to some extent. We have confirmed a truncation of the color-magnitude sequence at K'_{AB}\sim22; faint passively-evolving galaxies may not yet be present in this cluster at z\sim1.2. The overall color distribution of the cluster members, selected by the photometric redshift technique, is found to be very broad. We derived the fraction of blue galaxies in this cluster following a

(abr.) We present the analysis of the spectroscopic catalog of galaxies in the Coma cluster from Mobasher et al. This catalog comprises ~300 spectra of cluster members with magnitudes M_B=-20.5 to -14 in two areas of ~1 X 1.5 Mpc towards the center and the SW region of the cluster. In this paper the analysis is restricted to the 257 galaxies with no emission lines in their spectra. The strength of the age-sensitive indices is found to correlate with galaxy magnitude over the whole magnitude range. Similarly, the metallicity-sensitive indices anticorrelate with magnitude. By comparing the observed indices with model grids based on the Padova isochrones, we derive luminosity-weighted ages and metallicities. We present the distributions of ages and metallicities for galaxies in various magnitude bins. The mean metallicity decreases with galaxy magnitude and, at a given luminosity, appears to be generally lower for galaxies in the SW region of Coma as compared to the center of the cluster. A broad range of ages, from younger than 3 Gyr to older than 9 Gyr, is found in galaxies of any magnitude. However, systematic trends of age with luminosity are present among galaxies in the central field. In the central Mpc of Coma, a large fraction of galaxies at any luminosity (50-60% of the giants, >30 % of the dwarfs) show no evidence in their central regions of star formation occurred at redshift z<2, while the proportion of galaxies with significant star formation occurring at intermediate (0.35<z<2) and low (z<0.35) redshifts is found to depend on galaxy luminosity. Faint galaxies with young luminosity-weighted ages appear to have a bimodal metallicity distribution that, if confirmed, would point to a composite formation scenario. Finally, we present an interpretation of the index-magnitude relations observed.

The ages of stellar populations in 52 elliptical and S0 galaxies in the Coma cluster are investigated using a new spectroscopic sample. More than 40% of the S0s are found to have experienced recent star formation in their central regions during the last ~5 Gyrs, while such activity is absent in the ellipticals. Galaxies in this sample have absolute magnitudes in the range -20.5 < M_B < -17.5, and the fraction of S0 galaxies with recent star formation is higher at fainter luminosities. The observed luminosity range of S0 galaxies with signs of recent star formation activity is consistent with them being the descendants of typical star-forming spirals at intermediate redshift whose star formation has been halted as a consequence of the dense environment.

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

This is the second paper in a series studying the photometric and spectroscopic properties of galaxies of different luminosities in the Coma cluster. The sample selection, spectroscopic observations and completeness functions are presented here. To study the spectral properties of galaxies as a function of their local environment, two fields were selected for spectroscopic observations to cover both the core (Coma1) and outskirts (ie. south-west of the core and centered on NGC4839)- (Coma3) of the cluster. Medium resolution spectroscopy (6-9 \AA) was carried out for a total of 490 galaxies in both fields (302 in Coma1 and 188 in Coma3), using the WYFFOS multi-fiber spectrograph on the William Herschel Telescope. The galaxies cover a range of $12 < R < 20$, corresponding to -23 < M_R < -15 (H0=65 km/sec/Mpc). The redshifts are measured with an accuracy of 100 km/sec. The spectral line strengths and equivalent widths are also measured for the same galaxies and analysed in Poggianti et al (2001- paper III). A total of 189 (Coma1) and 90 (Coma3) galaxies are identified as members of the Coma cluster. An analysis of the colors show that only two members of the Coma cluster in our sample have B-R > 2. The completeness functions for the spectroscopic sample is presented.

We present low-resolution, near-IR JHK spectra of the weak z=2.39 radio galaxy 53W002, obtained with the OH-airglow Suppressor spectrograph (OHS) and Cooled Infrared Spectrograph and Camera for OHS (CISCO) on the Subaru Telescope. They cover rest-frame wavelengths of 3400-7200 A, and the emission lines of [O II]3727, Hb, [O III]4959, 5007, Ha, [N II]6548, 6583 and [S II]6716, 6731 were detected. Using the Ha/Hb line ratio, we find an extinction of E(B-V)=0.14. The emission-line ratios are reproduced by a cloud of electron density n_e=1x10^{3-4}(/cm3) with solar metallicity, ionized by an alpha=-0.7 power-law continuum with ionizing parameter U=1x10^-3. In addition to these emission lines, we make the first spectroscopic confirmation of the Balmer discontinuity in a high-z radio galaxy. Together with rest-frame UV photometry from the literature, we show that at least 1/3 of the present stellar mass was formed in the current starburst. The stellar mass was estimated to be (1-1.4)x10^11 M_sol by one-component model fitting, which is smaller than that of typical z~1 B2/6C radio galaxies. We suggest that 53W002 is currently assembling a large part of its stellar mass through merger events with the surrounding sub-galactic clumps, some of which can be identified with the Lya emitters detected in narrow-band imaging. After a few such events over the next few Gyr, 53W002 will evolve into a massive elliptical galaxy.

Deep near-infrared images of a blank 2'x2' section of sky near the Galactic north pole taken by Subaru Telescope are presented. The total integration times of the J and K' bands are 12.1 hours and 9.7 hours, resulting in 5-sigma limiting magnitudes of 25.1 and 23.5 mag, respectively. The numbers of sources within these limiting magnitudes found with an automated detection procedure are 385 in the J band and 350 in K'. Based on photometric measurements of these sources, we present number count vs. magnitude relations, color vs. magnitude diagrams, size vs. color relationships, etc. The slope of the galaxy number count plotted against the AB magnitude scale is about 0.23 in the 22 to 26 AB magnitude range of both bands. The spatial number density of galaxies as well as the slopes in the faint-end region given by the Subaru Deep Field (SDF) survey is consistent with those given by HST-NICMOS surveys as expressed on the AB magnitude diagram. Several sources having very large J-K' color are found including a few K' objects without detection at J. In addition, a number of faint Galactic stars are also detected, most of which are assigned to M-subdwarfs, together with a few brown dwarf candidates.

We construct infrared star count models of the Galaxy applicable at faint magnitudes. In addition to dwarfs and subdwarfs, we also take into account L dwarfs, T dwarfs, and white dwarfs of both disk and halo. As a test case of the models, we analyze the infrared imaging data at J and K' obtained during the Subaru Deep Field Survey to study stellar objects. Infrared star count models we have obtained predict that the Next Generation Space Telescope will primarily see T dwarfs, M subdwarfs, and old halo white dwarfs at faint magnitudes.

Results of a narrow-band imaging for the redshifted [OII]3727 emission around a quasar at z=4.7 obtained with the Subaru telescope and CISCO (a Cassegrain near infrared camera) are presented. A significant emission line is detected in the narrow-band H_2 (v=1-0 S(1)) filter at a location 2.4" northwest from the quasar, where the presence of a companion has been reported in Lyman alpha emission and the rest-frame UV continuum. We identify this line as [OII]3727 emission and confirm that the source really is a companion at z=4.7. The [OII]3727 flux from the companion is estimated to be 2.5 x 10^{-17} erg s^{-1} cm^{-2}. If the companion is a star forming object, the inferred star formation rate is as high as 45-230 M_{solar} yr^{-1} even without assuming the extinction correction. This value is higher than those derived from the Lyman alpha emission or from the UV continuum. Thus, provided that the difference is caused by dust extinction, the extinction corrected star formation rate is calculated to be 45 to 2300 M_{solar} yr^{-1} depending on the assuming extinction curves.

We have obtained high-resolution K'-band images of the powerful z=1.206 radio galaxy 3C 324 with the Subaru telescope under seeing conditions of 0.3--0.4 arcsec. We clearly resolved the galaxy and directly compared it to the optical images obtained with the Hubble Space Telescope. The host galaxy of 3C 324 is revealed to be a moderately luminous elliptical galaxy with a smooth light profile. The effective radius of the galaxy, as determined by profile fitting, is 1.3+-0.1 arcsec (1.2 kpc), which is significantly smaller than the value of 2.2 arcsec, published in Best et al. (1998, MNRAS, 292, 758). The peak of the K'-band light coincides with the position of the radio core, which implies that the powerful AGN lies at the nucleus of the host galaxy. The peak also coincides with the gap in the optical knotty structures which may be a dust lane hiding the UV-optical emission of the AGN from our line of sight; it is very likely that we are seeing the obscuring structure almost edge-on. We clearly detected the `aligned component' in the K'-band image by subtracting a model elliptical galaxy from the observed image. The red R_F702W-K color of the outer region of the galaxy avoiding the aligned component indicates that the near infrared light of the host galaxy is dominated by an old stellar population.

We investigate the $K$-band luminosity distribution of galaxies in the region of clusters at $z\sim1.2$ near to the radio galaxy 3C 324. The imaging data were obtained during the commissioning period of the Subaru telescope. There is a significant excess of the surface number density of the galaxies with $K =$ 17--20 mag in the region within $\sim$ 40'' from 3C 324. At this bright end, the measured luminosity distribution shows a drop, which can be represented by the exponential cut off of the Schechter-function formula; the best-fitted value of the characteristic magnitude, $K^{*}$, is $\sim 18.4\pm0.8$. This measurement follows the evolutionary trend of the $K^*$ of the rich clusters observed at an intermediate redshift, which is consistent with passive evolution models with a formation redshift $z_f \gtsim 2$. At $K \gtsim 20$ mag, however, the excess of the galaxy surface density in the region of the clusters decreases abruptly, which may imply that the luminosity function of the cluster galaxies has a negative slope at the faint end. This may imply strong luminosity segregation between the inner and outer parts of the clusters, or some deficit of faint galaxies in the cluster central region of the cluster.

We have investigated the optical and near-infrared colors of K'-selected galaxies in clusters at z ~ 1.2 near to the radio galaxy 3C 324 using images obtained with the Subaru telescope and archival HST data. The distribution of colors of the galaxies in the cluster region is found to be fairly broad, and it may imply significant scatter in their star-formation histories, although the effect of contamination of field galaxies is uncertain. The red sequence of galaxies whose R-K colors are consistent with passive evolution models for old galaxies is found to be truncated at K' ~ 20 mag, and there are few fainter galaxies with similar red colors in the cluster region. We find that the bulge-dominated galaxies selected by quantitative morphological classification form a broad sequence in the color-magnitude diagram, whose slope is much steeper than that expected from metallicity variations within a passively evolving coeval galaxy population. We argue that the observed color-magnitude sequence can be explained by metallicity and age variations, and the fainter galaxies with K' > 20 mag may be 1-2 Gyr younger than the brighter galaxies. Some spatial segregation of the color and K'-band luminosity is seen in the sky distribution; the redder and the brighter objects tend to be located near 3C 324.

We have performed a long-slit K band spectroscopic observation of the luminous infrared galaxy NGC 6240. The peak position of the H_2 v=1-0 S(1) emission in the slit is located ~0.3" - 0.4" north of the southern nucleus. It is almost the midpoint between the southern nucleus and the peak position of the ^12CO J=1-0 emission. Based on the line-ratio analyses, we suggest the excitation mechanism of H_2 is pure thermal at most positions. In the southern region we find the following three velocity components in the H_2 emission: the blueshifted shell component (~-250 km s^-1 with respect to V_sys) which is recognized as a distinct C-shape distortion in the velocity field around the southern nucleus, the high-velocity blueshifted ``wing'' component (~-1000 km s^-1 with respect to V_sys), and the component indicating possible line splitting of ~500 km s^-1. The latter two components are extended to the south from the southern nucleus. We show that the kinematic properties of these three components can be reproduced by expanding motion of a shell-like structure around the southern nucleus. The offset peak position of the H_2 emission can be understood if we assume that the shell expanding to the north interacts with the extragalactic molecular gas. At the interface between the shell and the molecular gas concentration the cloud-crushing mechanism proposed by Cowie et al. (1981) may work efficiently, and the intense H_2 emission is thus expected there. All these findings lead us to propose a model that the most H_2 emission is attributed to the shock excitation driven by the superwind activity of the southern nucleus.

We present the 2.12~$\mu$m narrow-band image of the Hubble Deep Field North taken with the near-infrared camera (CISCO) on the Subaru telescope. Among five targets whose H$\alpha$ or [O~{\sc iii}] emission lines are redshifted into our narrow-band range expected from their spectroscopic redshift, four of them have strong emission lines, especially for the two [O~{\sc iii}] emission-line objects. The remaining one target shows no H$\alpha$ emission in spite of its bright rest-UV luminosity, indicating that this object is already under the post-starburst phase. The volume-averaged $SFR$ derived from the detected two H$\alpha$ emission is roughly consistent with that evaluated from the rest-UV continuum.

We report on infrared imaging observations of the z=2.429 radio galaxy B3 0731+438 with the Subaru telescope. The images were taken with the K'-band filter and the 2.25 um narrow-band filter to examine the structure and properties of the Ha+[N II] 6548,6583 emission-line components. The Ha+[N II] emission-line image shows biconical lobes with an extent of 40 kpc, which are aligned with the radio axis. The rest-frame equivalent widths of the emission lines at these cones are as large as 1100 AA, and can be well explained by a gas-cloud model photoionized by power-law continuum radiation. The isotropic ionizing photon luminosity necessary to ionize the hydrogen gas in these cones amounts to 1e57(photons/s), which is larger than that in the majority of radio-loud QSOs. From these results, we propose that the Ha alignment effect in this object is produced by biconical gas clouds, which are swept up by the passage of radio jets, and are ionized by strong UV radiation from a hidden AGN. The continuum image consists of two components, a stellar-like point source and an extended diffuse galaxy. These are supposed to be a type-2 AGN and its host galaxy. The SED is fitted by a combination of spectra of a reddened dust-scattered AGN and an instantaneous starburst population of 500 Myr old. The stellar mass of the galaxy is estimated to be 3e11 M_solar, which is as large as that of typical 3C radio galaxies at z=1.