We present a description of the second data release for the photometric redshift (photo-$z$) of the Subaru Strategic Program for the Hyper-Suprime Cam survey. Our photo-$z$ products for the entire area in the Data Release 2 are publicly available, and both our point estimate catalog products and full PDFs can be retrieved from the data release site, \url{https://hsc-release.mtk.nao.ac.jp/}.

We perform a high-cadence transient survey with Subaru Hyper Suprime-Cam (HSC), which we call the Subaru HSC survey Optimized for Optical Transients (SHOOT). We conduct HSC imaging observations with time intervals of about one hour on two successive nights, and spectroscopic and photometric follow-up observations. A rapidly declining blue transient SHOOT14di at $z=0.4229$ is found in observations on two successive nights with an image subtraction technique. The rate of brightness change is $+1.28^{+0.40}_{-0.27}~{\rm mag~day^{-1}}$ ($+1.83^{+0.57}_{-0.39}~{\rm mag~day^{-1}}$) in the observer (rest) frame and the rest-frame color between $3400$ and $4400~\unicode[.8,0]{x212B}$ is $M_{\rm 3400\unicode[.8,0]{x212B}}-M_{\rm 4400\unicode[.8,0]{x212B}}=-0.4$. The nature of the object is investigated by comparing its peak luminosity, decline rate, and color with those of transients and variables previously observed, and those of theoretical models. None of the transients or variables share the same properties as SHOOT14di. Comparisons with theoretical models demonstrate that, while the emission from the cooling envelope of a Type IIb supernova shows a slower decline rate than SHOOT14di, and the explosion of a red supergiant star with a dense circumstellar wind shows a redder color than SHOOT14di, the shock breakout at the stellar surface of the explosion of a $25M_{\odot}$ red supergiant star with a small explosion energy of $\leq0.4\times10^{51}$ erg reproduces the multicolor light curve of SHOOT14di. This discovery shows that a high-cadence, multicolor optical transient survey at intervals of about one hour, and continuous and immediate follow-up observations, is important for studies of normal core-collapse supernovae at high redshifts.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

We develop an infinite mixture model of Ornstein-Uhlenbeck(OU) processes for describing the optical variability of QSOs based on treating the variability as a stochastic process. This enables us to get the parameters of the power spectral densities(PSDs) on their brightness variations by providing more flexible description of PSDs than the models based on single OU process(damped random walk). We apply this model to 67,507 variable objects extracted from SDSS Stripe82 photometric data and succeed in showing very high precision in identifying QSOs (~99% levels in completeness and purity) among variable objects based only on their variability, by investigating on 9,855 spectroscopically confirmed objects(7,714 QSOs and 2,141 stars) in the data of SDSS Data Release 12(DR12), with sufficient and accurate multiple measurements of their brightness. By comparing our results with the values based on other models that are used in previous research, it is revealed that our model can be used as the most effective method for selecting QSOs from variable object catalog, especially regarding completeness and purity. The main reason of improved identification rates are the ability of our model to separate clearly QSOs and stars, especially on the small fraction of QSOs with variabilities which can be described better than simple damped random walk model.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We present the deepest optical images of the COSMOS field based on a joint dataset taken with Hyper Suprime-Cam (HSC) by the HSC Subaru Strategic Program (SSP) team and the University of Hawaii (UH). The COSMOS field is one of the key extragalactic fields with a wealth of deep, multi-wavelength data. However, the current optical data are not sufficiently deep to match with, e.g., the UltraVista data in the near-infrared. The SSP team and UH have joined forces to produce very deep optical images of the COSMOS field by combining data from both teams. The coadd images reach depths of g=27.8, r=27.7, i=27.6, z=26.8, and y=26.2 mag at 5 sigma for point sources based on flux uncertainties quoted by the pipeline and they cover essentially the entire COSMOS 2 square degree field. The seeing is between 0.6 and 0.9 arcsec on the coadds. We perform several quality checks and confirm that the data are of science quality; ~2% photometry and 30 mas astrometry. This accuracy is identical to the Public Data Release 1 from HSC-SSP. We make the joint dataset including fully calibrated catalogs of detected objects available to the community at https://hsc-release.mtk.nao.ac.jp/.

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies.

We present clustering analysis results from 10,381 Lyman break galaxies (LBGs) at z~ 4-7, identified in the Hubble legacy deep imaging and new complimentary large-area Subaru/Hyper Suprime-Cam data. We measure the angular correlation functions (ACFs) of these LBGs at z~4, 5, 6, and 7, and fit these measurements using halo occupation distribution (HOD) models that provide an estimate of halo masses, M_h~(1-20)x10^11 Msun. Our M_h estimates agree with those obtained by previous clustering studies in a UV-magnitude vs. M_h plane, and allow us to calculate stellar-to-halo mass ratios (SHMRs) of LBGs. By comparison with the z~0 SHMR, we identify evolution of the SHMR from z~0 to z~4, and z~4 to z~7 at the >98% confidence levels. The SHMR decreases by a factor of ~2 from z~0 to 4, and increases by a factor of ~4 from z~4 to 7. We compare our SHMRs with results of a hydrodynamic simulation and a semi-analytic model, and find that these theoretical studies do not predict the SHMR increase from z~4 to 7. We obtain the baryon conversion efficiency (BCE) of LBGs at z~4, and find that the BCE increases with increasing dark matter halo mass. Finally, we compare our clustering+HOD estimates with results from abundance matching techniques, and conclude that the M_h estimates of the clustering+HOD analyses agree with those of the simple abundance matching within a factor of 3, and that the agreement improves when using more sophisticated abundance matching techniques that include subhalos, incompleteness, and/or evolution in the star formation and stellar mass functions.

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.

We report the multi-wavelength identification of the X-ray sources found in the Subaru-XMM-Newton Deep Survey (SXDS) using deep imaging data covering the wavelength range between the far-UV to the mid-IR. We select a primary counterpart of each X-ray source by applying the likelihood ratio method to R-band, 3.6micron, near-UV, and 24micron source catalogs as well as matching catalogs of AGN candidates selected in 1.4GHz radio and i'-band variability surveys. Once candidates of Galactic stars, ultra-luminous X-ray sources in a nearby galaxy, and clusters of galaxies are removed there are 896 AGN candidates in the sample. We conduct spectroscopic observations of the primary counterparts with multi-object spectrographs in the optical and NIR; 65\% of the X-ray AGN candidates are spectroscopically-identified. For the remaining X-ray AGN candidates, we evaluate their photometric redshift with photometric data in 15 bands. Utilising the multi-wavelength photometric data of the large sample of X-ray selected AGNs, we evaluate the stellar masses, M*, of the host galaxies of the narrow-line AGNs. The distribution of the stellar mass is remarkably constant from z=0.1 to 4.0. The relation between M* and 2--10 keV luminosity can be explained with strong cosmological evolution of the relationship between the black hole mass and M*. We also evaluate the scatter of the UV-MIR spectral energy distribution (SED) of the X-ray AGNs as a function of X-ray luminosity and absorption to the nucleus. The scatter is compared with galaxies which have redshift and stellar mass distribution matched with the X-ray AGN. The UV-NIR SEDs of obscured X-ray AGNs are similar to those of the galaxies in the matched sample. In the NIR-MIR range, the median SEDs of X-ray AGNs are redder, but the scatter of the SEDs of the X-ray AGN broadly overlaps that of the galaxies in the matched sample.

We present the results from a 1.1 mm imaging survey of the SSA22 field, known for having an overdensity of z=3.1 Lyman-alpha emitting galaxies (LAEs), taken with the AzTEC camera on the Atacama Submillimeter Telescope Experiment (ASTE). We imaged a 950 arcmin$^2$ field down to a 1 sigma sensitivity of 0.7-1.3 mJy/beam to find 125 submillimeter galaxies (SMGs) with a signal to noise ratio >= 3.5. Counterpart identification using radio and near/mid-infrared data was performed and one or more counterpart candidates were found for 59 SMGs. Photometric redshifts based on optical to near-infrared images were evaluated for 45 SMGs of these SMGs with Spitzer/IRAC data, and the median value is found to be z=2.4. By combining these estimation with estimates from the literature we determined that 10 SMGs might lie within the large-scale structure at z=3.1. The two-point angular cross-correlation function between LAEs and SMGs indicates that the positions of the SMGs are correlated with the z=3.1 protocluster. These results suggest that the SMGs were formed and evolved selectively in the high dense environment of the high redshift universe. This picture is consistent with the predictions of the standard model of hierarchical structure formation.

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 report the detection of an extremely bright ($\sim$37 mJy at 1100 $\mu$m and $\sim$91 mJy at 880 $\mu$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $\mu$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $\mu$m and SMA 880 $\mu$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $\sigma$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs.

We investigate mass-dependent galaxy evolution based on a large sample of (more than 50,000) K-band selected galaxies in a multi-wavelength catalog of the Subaru/XMM-Newton Deep Survey (SXDS) and the UKIRT Infrared Deep Sky Survey (UKIDSS)/Ultra Deep Survey (UDS). We employ the optical to near-infrared photometry to determine photometric redshifts of these galaxies. Then, we estimate the stellar mass of our sample galaxies using a standard fitting procedure. From the sample galaxies, we obtain the stellar mass function of galaxies and the cosmic stellar mass density up to z<4. Our results are consistent with previous studies and we find a considerable number of low-mass galaxies (M<10^{10.5}) at the redshift range 3<z<4. The stellar-mass dependent correlation functions of our sample galaxies show clear evolution and they connect to that in the local universe consistently. Also, the massive galaxies show strong clustering throughout our studied redshift range. The correlation length of massive galaxies rapidly decreases from z=4 to 2. We also find some high mass density regions of massive galaxies at 1.4<z<2.5 in our sample, which may be candidate progenitors of the present-day clusters of galaxies. At this redshift range, massive star-forming galaxies are the dominant population making up the structures and the passively evolving galaxies show stronger clustering and they may have formed earlier than those star-forming galaxies.

We present the result of projected cross correlation analysis of AGNs and galaxies at redshifts from 0.3 to 3.0. The Japanese Virtual Observatory (JVO) was used to obtain the Subaru Suprime-Cam images and UKIDSS catalog data around AGNs. We investigated 1,809 AGNs, which is about ten times larger a sample than that used in previous studies on AGN-galaxy clustering at redshifts larger than 0.6. 90% of the AGN samples are optically-selected AGN from the SDSS and 2dF catalogs. The galaxy samples at low redshift includes many redder objects from UKIDSS survey, while at higher redshift they are mainly blue galaxies from Suprime-Cam. We found significant excess of galaxies around the AGNs at redshifts from 0.3 to 1.8. For the low redshift samples ($z<0.9$), we obtained correlation length of $r_{0} = $5--6 $h^{-1}$Mpc ($\gamma = 1.8$), which indicates that the AGNs at this redshift range reside in a similar environment around typical local galaxies. We also found that AGNs at higher redshift ranges reside in a denser environment than lower redshift AGNs; For $z=1.3 \sim 1.8$ AGNs, the cross correlation length was measured as 11$^{+6}_{-3}$ $h^{-1}$Mpc ($\gamma=1.8$). Considering that our galaxies sample is based on optical observations with Suprime-Cam at the redshift range, it is expected that blue star-forming galaxies comprise the majority of objects that are observed to be clustered around the AGNs. It is successfully demonstrated that the use of the archive through the Virtual Observatory system can provide a powerful tool for investigating the small scale environment of the intermediate redshift AGNs.

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

To investigate the present situation of the merging in the southern outer region of Abell 85, we carried out long (~100 ks) observations with Suzaku, and produced an X-ray hardness ratio map. We found a high hardness ratio peak in the east side of a subcluster located in the south of the cluster; an X-ray spectrum of the region including this peak indicates a high temperature of ~8.5 keV. This hot spot has not been reported so far. We consider that this hot spot is a postshock region produced by the infall of the subcluster from the southwest. By using the Rankine--Hugoniot jump conditions for shocks, the Mach number and the infall velocity of the subcluster are obtained as 1.5 +/- 0.2 and 1950^{+290}_{-280} km s^{-1}, respectively, in the case of merging with the subcluster from the southwest direction. By using the redshift difference between the A 85 and the subcluster obtained from optical observations, the angle between the line of sight and the direction of the motion of the subcluster is estimated to be 75^{+7}_{-8} degrees. We estimate the kinetic energy of the subcluster and the energy used for intracluster medium (ICM) heating to be ~10^{63} and \lesssim 8 \times 10^{60} erg, respectively. This shows that the deceleration of the subcluster by ICM heating has been negligibly small.

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.

Lyman-alpha emitters are thought to be young, low-mass galaxies with ages of ~10^8 yr. An overdensity of them in one region of the sky (the SSA 22 field) traces out a filamentary structure in the early Universe at a redshift of z = 3.1 (equivalent to 15 per cent of the age of the Universe) and is believed to mark a forming protocluster. Galaxies that are bright at (sub)millimetre wavelengths are undergoing violent episodes of star formation, and there is evidence that they are preferentially associated with high-redshift radio galaxies, so the question of whether they are also associated with the most significant large-scale structure growing at high redshift (as outlined by Lyman-alpha emitters) naturally arises. Here we report an imaging survey of 1,100-um emission in the SSA 22 region. We find an enhancement of submillimetre galaxies near the core of the protocluster, and a large-scale correlation between the submillimetre galaxies and the low-mass Lyman-alpha emitters, suggesting synchronous formation of the two very different types of star-forming galaxy within the same structure at high redshift. These results are in general agreement with our understanding of the formation of cosmic structure.

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.

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 estimates of the photometric redshifts, stellar masses and star formation histories of sources in the SCUBA HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ~320 square arcmin. Using photometry spanning the B band to 8um, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift but the exact ratio of early to late mass is quite sensitive to the way extinction is treated in the modelling. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of SED types, only two out of the 51 sources we have obtained photometric redshifts for are best fit by a quasar-like spectral energy distribution, with approximately 80 per cent of the sources being best fit with late-type spectra (Sc, Im and starburst). By including photometry at 850um, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15-65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright (L_* +2 mag < L_optical < L_* -1 mag) galaxies in the field over the interval 0 < z < 3 will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around z=2.

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)

We have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

We describe deep radio imaging at 1.4-GHz of the 1.3 square degree Subaru/XMM-Newton Deep Field (SXDF), made with the Very Large Array in B and C configurations. We present a radio map of the entire field, and a catalogue of 505 sources covering 0.8 square degrees to a peak flux density limit of 100 microJy. Robust optical identifications are provided for 90% of the sources, and suggested IDs are presented for all but 14 (of which 7 are optically blank, and 7 are close to bright contaminating objects). We show that the optical properties of the radio sources do not change with flux density, suggesting that AGNs continue to contribute significantly at faint flux densities. We test this assertion by cross-correlating our radio catalogue with the X-ray source catalogue and conclude that radio-quiet AGNs become a significant population at flux densities below 300 microJy, and may dominate the population responsible for the flattening of the radio source counts if a significant fraction of them are Compton-thick.

We report the results of a systematic near-infrared spectroscopic survey using the Subaru, VLT and Keck Telescopes of a sample of high redshift Ultra-luminous Infrared Galaxies (ULIRGs) mainly composed of submillimeter-selected galaxies. Our observations span the restframe optical range containing nebular emission lines such as H_beta, [OIII], and [OII], which are essential for making robust diagnostics of the physical properties of these ULIRGs. Using the H_alpha/H_beta emission line ratios, we derive internal extinction estimates for these galaxies similar to those of local ULIRGs:A_v~2.9+-0.5. Correcting the H_alpha estimates of the star formation rate for dust extinction using the Balmer decrement, results in rates which are consistent with those estimated from the far-infrared luminosity. The majority (>60%) of our sample show spectral features characteristic of AGN, with ~65% exhibiting broad Balmer emission lines. A proportion of these sources show relatively low [OIII]/H_beta line ratios, which are similar to those of Narrow Line Seyfert 1 galaxies suggesting small mass black holes which are rapidly growing. In the subsample of our survey with both [OIII] and hard X-ray coverage, at least ~60% show an excess of [OIII] emission, by a factor of 5-10 times, relative to the hard X-ray luminosity compared to the correlation between these two properties seen in Seyferts and QSOs locally. From our spectral diagnostics, we propose that the strong [OIII] emission in these galaxies arises from shocks in dense gaseous regions in this vigorously star-forming population.

We present low-resolution ($R\sim 90$) and medium-resolution ($R\sim 2500$) spectropolarimetry of Nova V475 Sct with the HBS instrument, mounted on the 0.91-m telescope at the Okayama Astrophysical Observatory, and with FOCAS, mounted on the 8.2-m Subaru telescope. We estimated the interstellar polarization toward the nova from the steady continuum polarization components and H$\alpha$ line emission components. After subtracting the interstellar polarization component from the observations, we found that the H$\alpha$ emission seen on 2003 October 7 was clearly polarized. In the polarized flux spectrum, the H$\alpha$ emission had a distinct red wing extending to $\sim +4900$ km s$^{-1}$ and a shoulder around $+3500$ km s$^{-1}$, showing a constant position angle of linear polarization $\theta_{\rm *}\simeq 155\arcdeg\pm 15\arcdeg$. This suggests that the nova had an asymmetric outflow with a velocity of $v_{\rm wind}\simeq 3500$ km s$^{-1}$ or more, which is six times higher than the expansion velocity of the ionized shell at the same epoch. Such a high-velocity component has not previously been reported for a nova in the `moderately fast' speed class. Our observations suggest the occurrence of violent mass-loss activity in the nova binary system even during the common-envelope phase. The position angle of the polarization in the H$\alpha$ wing is in good agreement with that of the continuum polarization found on 2003 September 26 ($p_{\rm *}\simeq 0.4$--0.6 %), which disappeared within the following 2 d. The uniformity of the PA between the continuum polarization and the wing polarization on October 7 suggests that the axis of the circumstellar asymmetry remained nearly constant during the period of our observations.

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 present the results of a deep X-ray survey conducted with XMM-Newton, centred on the UK ROSAT 13H deep field area. This region covers 0.18 deg^2 and is the first of two areas covered with XMM-Newton as part of an extensive multi-wavelength survey designed to study the nature and evolution of the faint X-ray source population. We have produced detailed Monte-Carlo simulations to obtain a quantitative characterisation of the source detection procedure and to assess the reliability of the resultant sourcelist. We use the simulations to establish a likelihood threshold above which we expect less than 7 (3%) of our sources to be spurious. We present the final catalogue of 225 sources. Within the central 9 arcmin, 68 per cent of source positions are accurate to 2 arcsec, making optical follow-up relatively straightforward. We construct the N(>S) relation in four energy bands: 0.2-0.5 keV, 0.5-2 keV, 2-5 keV and 5-10 keV. In all but our highest energy band we find that the source counts can be represented by a double powerlaw with a bright end slope consistent with the Euclidean case and a break around 10^-14 cgs. Below this flux the counts exhibit a flattening. Our source counts reach densities of 700, 1300, 900 and 300 deg^-2 at fluxes of 4.1x10^-16, 4.5x10^-16, 1.1x10^-15 and 5.3x10^-15 cgs in the 0.2-0.5, 0.5-2, 2-5 and 5-10 keV energy bands respectively. We have compared our source counts with those in the two Chandra deep fields and Lockman hole and find our source counts to be amongst the highest of these fields in all energy bands. We resolve >51% (>50%) of the X-ray background emission in the 1-2 keV (2-5 keV) energy bands.

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.

We have carried out near-infrared JHK spectroscopy of a gravitationally lensed submillimeter galaxy SMM J14011+0252 at z=2.565, using OHS and CISCO on the Subaru telescope. This object consists of two optical components, J1 and J2, which are lensed by the cluster Abell 1835. J1 suffers additional strong lensing by a foreground galaxy at z=0.25 in the cluster. The rest-optical H-alpha, H-beta, and [O II]3727 lines are detected in both J1 and J2, and [N II]6548,6583 lines are also detected in J1. A diagnosis of emission-line ratios shows that the excitation source of J1 is stellar origin, consistent with previous X-ray observations. The continua of J1 and J2 show breaks at rest 4000A, indicating relatively young age. Combined with optical photometry, we have carried out model spectrum fitting of J2 and find that it is a very young (~50 Myr) galaxy of rather small mass (~10e8 M_sol) which suffers some amount of dust extinction. A new gravitational lensing model is constructed to assess both magnification factor and contamination from the lensing galaxy of the component J1, using HST-F702W image. We have found that J1 suffers strong lensing with magnification of ~30, and its stellar mass is estimated to be < 10e9 M_sol. These results suggest that SMM J14011+0252 is a major merger system at high redshift that undergoes intense star formation, but not a formation site of a giant elliptical. Still having plenty of gas, it will transform most of the gas into stars and will evolve into a galaxy of < 10e10 M_sol. Therefore, this system is possibly an ancestor of a less massive galaxy such as a mid-sized elliptical or a spiral at the present.

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.

We report on our new analysis of the spatial and kinematical distribution of warm and cold molecular gas in NGC 6240, which was undertaken to explore the origin of its unusually luminous H$_2$ emission. By comparing three-dimensional emission-line data (in space and velocity) of CO (J=2-1) in the radio and H$_2$ in the near infrared, we are able to study the H$_2$ emitting efficiency, defined in terms of the intensity ratio of H$_2$ to CO [$I$(H$_2$)/$I$(CO)], as a function of velocity. The integrated H$_2$ emitting efficiency is calculated by integrating the velocity profile of H$_2$ emitting efficiency in blue, red, and total (blue + red) velocity regions of the profile. We find that (1) both the total H$_2$ emitting efficiency and the blue-to-red ratio of the efficiency are larger in regions surrounding the CO and H$_2$ intensity peaks, and (2) the H$_2$ emitting efficiency and the kinematical conditions in the warm molecular gas are closely related to each other. A collision between the molecular gas concentration and the external superwind outflow from the southern nucleus seems plausible to explain these characteristics, since it can reproduce the enhanced emitting efficiency of blueshifted H$_2$ around the molecular gas concentration, if we assume that the superwind blows from the southern nucleus toward us, hitting the entire gas concentration from behind. In this model, internal cloud-cloud collisions within the molecular gas concentration are enhanced by the interaction with the superwind outflow, and efficient and intense shock-excited H$_2$ emission is expected as a result of the cloud-crushing mechanism.

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.

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 optical and near-IR spectroscopy of a z=2.38 hyperluminous IR galaxy, covering the restframe wavelength range from 1000-5000A. It appears to comprise two components separated by less than 1" on the sky (<8kpc); one component (B) is blue, the other (P) is red in restframe UV-optical colours. The combined system has an optical luminosity of ~8L* and its restframe optical spectrum is characteristic of a Seyfert AGN. However, its restframe UV spectrum exhibits striking features associated with young stars, including P-Cygni lines from stellar winds and blue-shifted interstellar absorption lines indicative of a galactic outflow. Redshifts are derived from stellar photospheric lines in the UV and from narrow emission lines in the restframe optical, and these are compared to that measured for the molecular gas recently detected with the IRAM interferometer. The offsets indicate that the far-IR emission is most likely associated with the near-IR source P, which hosts the Seyfert nucleus, while the UV-bright component B is blueshifted by 400km/s. This suggests that the two components are probably merging and the resulting gravitational interactions have triggered the hyperluminous activity. Modelling of the UV spectral features implies that the starburst within the UV component of this system has been going on for at least ~10Myrs. Assuming that the bolometrically-dominant obscured component has a similar lifetime, we estimate that it has so far formed a total stellar mass of ~10^11 Mo. If this star formation continues at its present level for substantially longer, or if this activity is repeated, then the present-day descendant of N2 850.4 will be a very luminous galaxy.

We present the results of a Chandra ACIS-I survey of a high latitude region at 13h +38 which was earlier observed with ROSAT (McHardy et al. 1998) and which has recently been observed by XMM-Newton for 200ksec. XMM-Newton will provide good quality X-ray spectra for over 200 sources with fluxes around the knee of the logN/logS, which are responsible for the bulk of the X-ray background (XRB). The main aim of the Chandra observations is to provide arcsecond, or better, positions, and hence reliable identifications, for the XMM-Newton sources. The ACIS-I observations were arranged in a mosaic of four 30ksec pointings, covering almost all of the 15' radius XMM-Newton/ROSAT field. We detect 214 Chandra sources above 5 sigma significance, to a limiting flux of \~1.3e-15 erg/cm2/s (0.5-7 keV). Optical counterparts are derived from a Subaru SuprimeCam image reaching to R~27. The very large majority of the Chandra sources have an optical counterpart, with the distribution peaking at 23<R<24, although 14 have no counterpart to R=27. The fraction of X-ray sources with no identification brighter than R=27 is similar to that found in deeper Chandra surveys. The majority of the identifications are with galaxies. As found in other Chandra surveys, there is a very wide range of optical magnitude for given X-ray flux, implying a range of emission mechanisms, and many sources have high L_X/L_opt ratios, implying absorption at moderate redshift. Comparison with the earlier ROSAT survey shows that the accuracy of the ROSAT positions agrees very well with the predictions from simulations in McHardy et al. (1998) and that the large majority of the identifications were correct (abridged).

The 13hr XMM-Newton/Chandra deep survey is the first of two extremely deep XMM-Newton fields observed by the XMM-OM consortium. A 120 ks Chandra mosaic, covering 0.2 square degrees, provides sensitive, confusion-free point source detection with sub-arcsecond positions, while the 200 ks XMM-Newton observation provides high quality X-ray spectroscopy over the same sky area. We have optical spectroscopic identifications for 70 X-ray sources. Of these, 42 are broad emission-line AGN with a wide range of redshifts. The optical counterparts of a further 23 sources are narrow emission line galaxies and absorption line galaxies. These 23 sources all lie at z < 1 and typically have lower X-ray luminosities than the broad-line AGN. About half of them show significant X-ray absorption and are almost certainly intrinsically absorbed AGN. However some of them have unabsorbed, AGN-like, power-law components in their X-ray spectra, but do not show broad emission lines in their optical spectra. These sources may be weak, unobscured AGN in bright galaxies and their existence at low redshifts could be a consequence of the strong cosmological evolution of AGN characteristic luminosities.

Our VLA observations of the XMM-Newton/Chandra 13hr deep survey field (see Page et al., this proceedings) result in one of the two deepest 1.4GHz radio maps ever made. Within the 15' radius field covered by the deep X-ray data (0.19 sq deg), a total of 556 radio sources are detected, down to a 4 sigma flux limit of 28uJy. Of the 214 Chandra sources, 55 have radio counterparts. The sub-arcsecond accuracy of the VLA and Chandra positions enable us to determine with high confidence the sources common to both surveys. Here we present the relationship between the X-ray and radio source populations at the faintest radio flux limits yet probed by such a study. We discuss how the X-ray/radio relationship differs as a function of optical morphology, ie between unresolved `stellar' objects and well resolved galaxies. We then discuss the origin of the X-ray and radio emission, ie AGN, starburst or a mixture of both, in these two classes of object.

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

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

The results of deep-imaging observations of a candidate of an absorbed QSO at z=0.653, AX J131831+3341, are presented. AX J131831+3341 was found during the course of optical follow-up observations of the ASCA Large Sky Survey, and has an X-ray luminosity of 10^{45} erg s^{-1} (2-10keV), which corresponds to those of QSOs. Its optical spectrum shows no significant broad H_beta emission line, suggesting that the object is an absorbed QSO. Deep R and V band images reveal the presence of a point-like nucleus and an asymmetric extended component. The nuclear component has a blue color, and the optical magnitude is much fainter than that expected from the observed X-ray flux for typical type-1 AGNs. These photometric properties and the presence of broad MgII 2800A emission can be explained simultaneously if the observed nuclear light is dominated by scattered nuclear light, though there is a possibility that the nuclear component is a slightly absorbed nucleus if its intrinsic X-ray to optical flux ratio is the largest among X-ray selected AGNs. The size of the extended component, which is thought to be the host galaxy of the QSO, is larger than those of normal disk galaxies at z=0-0.75, and the absolute magnitude is similar to those of the brightest host galaxies of QSOs at redshifts smaller than 0.5. The V-R and R-I colors of the component are consistent with a 1 Gyr-old stellar population model without absorption.

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 investigated the central gas systems of E/S0 galaxies by making use of the WFPC2 images of the Hubble Space Telescope archive. We searched the gas systems that were traced by the dust with a new method of making color excess images in F555W - F814W (V-I). Out of 25 sample galaxies, we detected gas system in 14 galaxies. The dust was newly detected in two galaxies that were thought to contain no dust based on single band, pre-refurbishment data. The full extents of the gas systems are 0.1 to 3.5 kpc, and the masses of the gas, log M(gas) [M(solar)], are 4.2 to 7.2. The AGN activity is well correlated with existence of the gas systems. None of galaxies without the gas systems show the AGN activity. On the other hand, some galaxies with the gas systems show the AGN activity; optical AGN activities are shown in 5 out of 11 galaxies of which AGNs are optically studied, and radio activities are shown in 6 out of 14 galaxies. This shows that the AGN activity is driven with the gas system.

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

The fraction of star-forming galaxies in rich clusters of galaxies increases rapidly with the redshift. This is interpreted as the result of a rapid evolution of cluster galaxies, though its mechanism is not yet clear. One hypothesis is that if galaxies run into the dense intracluster medium (ICM) regions, starbursts can be induced due to compression of molecular clouds in the galaxies by a raised "external" ICM pressure. In a merging cluster, there would be dense ICM regions and some galaxies may experience a rapid increase of external pressure. Thus, if the mechanism works, starbursts would occur in such galaxies. Such a situation is considered to be realized in Coma cluster, and galaxies showing a recent starburst have a characteristic spatial distribution; they populate an elongated region between two sub-clusters. We examined the above hypothesis by measuring the spatial distribution of blue galaxies, regarded as star-forming galaxies in Abell~168 to check whether a case like Coma is realized in this recently merging cluster. However, we could not find an enhanced blue-galaxy distribution between two sub-clusters. We discuss the causes of our result, including possibilities of surviving the hypothesis.