Search

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., Comment: Accepted in ApJ

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., Comment: 50 pages, 6 Tables, 13 Figures; Accepted for publication in The Astrophysical Journal

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., Comment: 15 pages, 9 figures, 7 tables, to appear in ApJ

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., Comment: 10 pages, 4 figures, AJ accepted for publication

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: Accepted for publication in ApJ (May 7, 2015). 14 pages, 6 figures

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., Comment: 8 pages, 2 figures, ApJL accepted for publication

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., Comment: ApJ in press, replaced with accepted version. Methodology clarified, conclusions mostly unchanged

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., Comment: 9 pages, 4 figures, accepted for publication in ApJ Letters

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., Comment: 20 pages, 16 figures, accepted for publication in Astronomy and Astrophysics on August 6, 2013

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., Comment: 18 pages, 11 figures. Accepted by ApJ on 30 July 2013. The accepted version was submitted on 8 July 2013, and the original version was submitted on 3 March 2013

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., Comment: 20 pages, 12 figures, Accepted for publication in ApJ. Minor updates from the version 1

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., Comment: 26 pages, 4 figures, 6 tables. Replaced with published ApJ version

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., Comment: 6 pages, 4 figures, 1 table. Accepted for publication in ApJ Letters

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., Comment: 45 pages, 9 figures; ApJ in press

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., Comment: Accepted for publication in Astrophysical Journal Letters. 25 pages, 4 figures, 1 table

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., Comment: 21 pages, 3 figures, published 2012 November 7 by ApJL, typo corrected

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., Comment: 28 pages (single column), 7 figures, 1 table, ApJ accepted

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., Comment: 15 pages, 9 figures, accepted to ApJ. Replaced with accepted version; mostly minor changes. Software updated

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., Comment: 20 pages, 8 figures, and 1 table. accepted to PASJ

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., Comment: 21 pages, 6 figures, 1 table, ApJ accepted 2012-05-06

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., Comment: 41 pages (single column), 1 table, 10 figures, ApJ accepted

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., Comment: 8 pages, 5 figures, ApJL in press, typo corrected

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")., Comment: 6 pages, 4 figures, 1 table; accepted for publication in ApJ Letters

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., Comment: 10 pages, 4 figures, accepted for publication in ApJ Letters. New version: Corrected a few numbers in the astrometry section (which were already correct in the print version, but were based on an outdated simulation in the astro-ph version)

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., Comment: 22 pages, 10 figures, accepted for publication in The Astrophysical Journal

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., Comment: 40 pages, 44 figures, accepted for publication in the ApJ. High resolution version available at http://www.astr.tohoku.ac.jp/~mikito/AndromedaHalo.pdf

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), Comment: 55 pages, 22 figures, 9 tables, accepted by ApJS, a higher-resolution version is available at http://step.mtk.nao.ac.jp/sxds/

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., Comment: 19 pages, 19 figures, Accepted for publication in MNRAS

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., Comment: 24 pages, 6 figures, 2 Tables. Accepted by the Astrophysical Journal

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., Comment: 32 pages, 14 figures, accepted for publication in the Astronomical Journal. A complete version with all tables and figures is available at http://www.astro.psu.edu/users/misawa/pub/Paper/40hires.ps.gz

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., Comment: 14 pages, 13 figures, 5 tables, submitted to PASJ

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-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., Comment: 21 pages, 6 figures, submitted to AJ

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., Comment: 6 pages, 4 figures, Accepted for publication in the Astrophysical Journal

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., Comment: To appear in APJ vol.648. Only minor corrections have been made. Black&White version is available at http://zone.mtk.nao.ac.jp/~kashik/sdf/z6p5lae/paper/sdf_z6p5lae_bw.pdf

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., Comment: 23 pages, 20 figs, accepted for PASJ, a high resolution version of Figs.7,8 is available at http://hikari.astron.s.u-tokyo.ac.jp/~shima/z57LAE/

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., Comment: 12 pages, 1 figure (figure only in PostScript format, not in PDF)

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., Comment: 10 pages, 4 figures. Figure of UBVr'i'JHKs light curves is revised. Accepted for publication in the Astrophysical Journal (Letters) on 5 October 2005

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., Comment: 21 pages, 6 figures, to appear in The Astrophysical Journal

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

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., Comment: To be published in Science on 27 May 2005. 14 pages, 4 figures

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., Comment: 8 pages, 4 figures, to be published in ApJL. The paper with high resolution figures is available at http://www-int.stsci.edu/~ouchi/work/astroph/sxds_LSS/

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/., Comment: 18 pages, 10figures, accepted for publication in PASJ, and see also http://step.mtk.nao.ac.jp/sdf/ [Japan site], http://soaps.naoj.org/sdf/ [US site]

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., Comment: 22 pages, 7 figures, accepted for publication in the Astronomical Journal

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., Comment: 13 pages including 3 figures, to appear in The Astrophysical Journal Letters

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., Comment: 13 pages including 4 figures, aastex preprint style. Minor changes. Accepted by ApJL

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., Comment: 20 pages, LaTeX, Accepted by Astron. J

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., Comment: 7 pages, 4 figures, uses emulateapj5.sty, accepted for puclication in ApJL

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., Comment: 12 pages, including 4 figures, accepted for publication in ApJL

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., Comment: 14 pages, 7 figures, Accepted for publication in the Astrophyical Journal