The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

We present our ALMA Band 8 observations of a damped Ly${\alpha}$ absorption (DLA) system at $z$=3.150 observed in the spectrum of the quasar Q2233+131 at $z$=3.295. The optical counterpart of this DLA has been identified and it shows a double-peaked Ly${\alpha}$ emission line. Since one possible origin of DLAs at high redshift is an outflowing gas from star-forming galaxies, DLA2233+131 provides a good laboratory to investigate the nature of high-$z$ DLAs. Motivated by this, we have carried out ALMA band 8 observations to study the [C II] line in this system. However, we do not detect any significant emission line in the observed pass bands. Instead, we have serendipitously found three submm continuum sources in the observed sky area. One appears to be the quasar Q2233+131 itself while the other two sources are newly identified submm galaxies (SMGs), called SMG1 and SMG2 in this paper. They are located at a separation of 400.7 and 800.1 from Q2233+131, respectively. Their 646 ${\mu}$m fluxes are 6.35 mJy and 6.43 mJy, respectively, being higher than that of Q2233+131, 3.62 mJy. Since these two SMGs are not detected in the optical images obtained with the Hubble Space Telescope and the Subaru Telescope, they have a very red spectral energy distribution. It is, therefore,suggested that they are high-redshift galaxies or very dusty galaxies at intermediate redshift although we cannot rule out the possibility that they are optically very faint SMG analogs at low redshift. Follow up observations will be necessary to explore the nature of this interesting region.

We search for galaxies with a strong Balmer break (Balmer Break Galaxies; BBGs) at $z \sim 6$ over a 0.41 deg$^2$ effective area in the COSMOS field. Based on rich imaging data, including data obtained with the Atacama Large Millimeter/submillimeter Array (ALMA), three candidates are identified by their extremely red $K - [3.6]$ colors as well as by non-detection in X-ray, optical, far-infrared (FIR), and radio bands. The non-detection in the deep ALMA observations suggests that they are not dusty galaxies but BBGs at $z \sim 6$, although contamination from Active Galactic Nuclei (AGNs) at $z \sim 0$ cannot be completely ruled out for the moment. Our spectral energy distribution (SED) analyses reveal that the BBG candidates at $z \sim 6$ have stellar masses of $\approx 5 \times 10^{10} M_{\odot}$ dominated by old stellar populations with ages of $\gtrsim 700$ Myr. Assuming that all the three candidates are real BBGs at $z \sim 6$, we estimate the stellar mass density (SMD) to be $2.4^{+2.3}_{-1.3} \times 10^{4} M_{\odot}$ Mpc$^{-3}$. This is consistent with an extrapolation from the lower redshift measurements. The onset of star formation in the three BBG candidates is expected to be several hundred million years before the observed epoch of $z \sim 6$. We estimate the star-formation rate density (SFRD) contributed by progenitors of the BBGs to be 2.4 -- 12 $\times 10^{-5} M_{\odot}$ yr$^{-1} $Mpc$^{-3}$ at $z > 14$ (99.7\% confidence range). Our result suggests a smooth evolution of the SFRD beyond $z = 8$.

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

The existence of $^{12}{\rm C} + {}^{16}{\rm O}$ molecular resonances at sub-barrier energy has been a significant problem in nuclear astrophysics because they strongly affect the $^{12}{\rm C} + {}^{16}{\rm O}$ fusion reaction rate in type Ia supernovae and heavy stars. However, experimental surveys have been limited to 4~MeV and cannot access the deep sub-barrier energy due to a very small fusion cross section. Here we predict a couple of resonances with $J^\pi=0^+$, $2^+$, and $4^+$ in the deep sub-barrier energy based on the antisymmetrized molecular dynamics calculation that reproduces the known resonances and low-lying spectrum of $^{28}{\rm Si}$.

Recent supernova searches revealed a number of fast-evolving luminous transients. We perform radiation hydrodynamic simulations of light curves of several models of supernova from super--asymptotic giant branch (super-AGB, SAGB) stars with low mass envelopes ($M_{\rm env}$ $\sim$ $0.05-1.25$ $M_{\odot}$). The differences in the light curves found among the models are used to link the observed events to the most appropriate models. In particular, we propose that KSN 2015K is an electron-capture supernova. We assume "optically" thick CSM around SAGB and the circumstellar interaction powers the peak luminosity of the light curve with a short rise time. The faint tail might be influenced by the spin-down luminosity of a newborn Crab-like pulsar. Our fits indicate an ejecta mass of $0.02-0.05$ $M_{\odot}$, a circumstellar medium (CSM) mass of $0.10-0.12$ $M_{\odot}$, a radius of the CSM photosphere $\sim$ 10$^{14}$ cm, a kinetic energy of $\sim$ 3 $\times$ 10$^{50}$ erg, a photospheric velocity $v_{\rm ph} \gtrsim 10,000$ km s$^{-1}$ and a pulsar total spin energy ($2.5-4$) $\times$ 10$^{49}$ erg.

We present results of the cross-correlation Lya intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Lya emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images and the LAE samples, we further conduct extensive analysis evaluating systematics of large-scale point spread function wings, sky subtractions, and unknown errors on the basis of physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Lya intensity of the narrowband images and the LAEs. We identify very diffuse Lya emission with the $3\sigma$ ($2\sigma$) significance at $>150$ comoving kpc (ckpc) far from the LAEs at $z = 5.7$ ($6.6$), beyond a virial radius of star-forming galaxies with $M_\mathrm{h} \sim 10^{11} M_\odot$. The diffuse Lya emission possibly extends up to $1,000$ ckpc with the surface brightness of $10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We confirm that the small-scale ($<150$ ckpc) Lya radial profiles of LAEs in our Lya intensity maps are consistent with those obtained by recent MUSE observations (Leclercq et al. 2017). Comparisons with numerical simulations suggest that the large-scale ($\sim 150-1,000$ ckpc) Lya emission are not explained by unresolved faint sources of neighboring galaxies including satellites, but by a combination of Lya photons emitted from the central LAE and other unknown sources, such as a cold-gas stream and galactic outflow. We find no evolution in the Lya radial profiles of our LAEs from $z=5.7$ to $6.6$, where theoretical models predict a flattening of the profile slope made by cosmic reionization, albeit with our moderately large observational errors.

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$\alpha$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$\alpha$ emission of the LABs. The Ly$\alpha$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$\alpha$ luminosity relation of the Ly$\alpha$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass.

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

We present CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Ly{\alpha} nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 {\sigma} from one of the ALMA continuum sources associated with the Ly{\alpha} filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)x10^9 K km s^{-1} pc^2 for this CO source, which is one order of magnitude smaller than those of typical z>1 dusty star-forming galaxies (DSFGs) of similar far-infrared luminosity L(IR)~10^{12} Lsun. We derived a molecular gas mass of Mgas=(4.4^{+0.9}_{-0.6})x10^9 Msun and a star-formation rate of SFR=270\pm160 Msun yr^{-1}. We also estimated a gas depletion time of {\tau}(dep)=17\pm10 Myr, being shorter than those of typical DSFGs. It is suggested that this source is in a transition phase from DSFG to a gas-poor, early-type galaxy. From ALMA to Herschel multi-band dust continuum observations, we measured a dust emissivity index {\beta}=2.3\pm0.2, which is similar to those of local gas-poor, early-type galaxies. Such a high {\beta} can be reproduced by specific chemical compositions for interstellar dust at the submillimeter wavelengths from recent laboratory experiments. ALMA CO and multi-band dust continuum observations can constrain the evolutionary stage of high-redshift galaxies through {\tau}(dep) and {\beta}, and thus we can investigate dust chemical compositions even in the early Universe.

We present morphological properties of dusty star-forming galaxies at z=1-3 determined with high-resolution (FWHM~0"19) Atacama Large Milllimeter/submilimeter Array (ALMA) 1-mm band maps of our ASAGAO survey covering a 26-arcmin^2 area in GOODS-S. In conjunction with the ALMA archival data, the present sample consists of 42 ALMA sources with a wide rest-frame far-infrared (FIR) luminosity L_FIR range of ~10^11-10^13 Lo. To obtain an average rest-frame FIR profile, we perform individual measurements and careful stacking of the ALMA sources using the uv-visibility method that includes positional-uncertainty and smoothing-effect evaluations through Monte-Carlo simulations. We find that the dusty star-forming galaxies have the average FIR-wavelength Sersic index and effective radius of n_FIR=1.2+/-0.2 and R_e,FIR=1.0-1.3 kpc, respectively, additionally with a point source at the center, indicative of the existence of AGN. The average FIR profile agrees with a morphology of an exponential-disk clearly distinguished from a spheroidal profile (Sersic index of 4). We also examine the rest-frame optical Sersic index n_opt and effective radius R_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain n_opt=0.9+/-0.3 (~n_FIR) and R_e,opt=3.2+/-0.6 kpc (>R_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs (n_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges.

We study the heating mechanisms and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field. Dust continuum sources have been identified in 11 of the 35 LABs, all with star formation rates (SFRs) above 100 Msun/yr. Likely radio counterparts are detected in 9 out of 29 investigated LABs. The detection of submm dust emission is more linked to the physical size of the Ly{\alpha} emission than to the Ly{\alpha} luminosities of the LABs. A radio excess in the submm/radio detected LABs is common, hinting at the presence of active galactic nuclei. Most radio sources without X-ray counterparts are located at the centers of the LABs. However, all X-ray counterparts avoid the central regions. This may be explained by absorption due to exceptionally large column densities along the line-of-sight or by LAB morphologies, which are highly orientation dependent. The median Ly{\alpha} escape fraction is about 3\% among the submm-detected LABs, which is lower than a lower limit of 11\% for the submm-undetected LABs. We suspect that the large difference is due to the high dust attenuation supported by the large SFRs, the dense large-scale environment as well as large uncertainties in the extinction corrections required to apply when interpreting optical data.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.

We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.

We present rest-frame optical spectra of 12 nitrogen-loud quasars at z ~ 2.2, whose rest-frame ultraviolet (UV) spectra show strong nitrogen broad emission lines. To investigate their narrow-line region (NLR) metallicities, we measure the equivalent width (EW) of the [OIII]5007 emission line: if the NLR metallicity is remarkably high as suggested by strong UV nitrogen lines, the [OIII]5007 line flux should be very week due to the low equilibrium temperature of the ionized gas owing to significant metal cooling. In the result, we found that our spectra show moderate EW of the [OIII]5007 line similar to general quasars. This indicates nitrogen-loud quasars do not have extremely metal-rich gas clouds in NLRs. This suggests that strong nitrogen lines from broad-line regions (BLRs) are originated by exceptionally high abundances of nitrogen relative to oxygen without very high BLR metallicities. This result indicates that broad-emission lines of nitrogen are not good indicators of the BLR metallicity in some cases. On the other hand, we also investigate virial black-hole masses and Eddington ratios by using the Hbeta and CIV1549 lines for our sample. As a result, we found that black-hole masses and Eddington ratios of nitrogen-loud quasars tend to be low and high relative to normal quasars, suggesting that nitrogen-loud quasars seem to be in a rapidly-accreting phase. This can be explained in terms of a positive correlation between Eddington ratios and nitrogen abundances of quasars, that is probably caused by the connection between the mass accretion onto black holes and nuclear star formation.

We study the UV luminosity functions (LFs) at $z\sim 4$, $5$, $6,$ and $7$ based on the deep large-area optical images taken by the Hyper Suprime-Cam (HSC) Subaru strategic program (SSP). On the 100 deg$^2$ sky of the HSC SSP data available to date, we make enormous samples consisting of a total of 579,565 dropout candidates at $z\sim 4-7$ by the standard color selection technique, 358 out of which are spectroscopically confirmed by our follow-up spectroscopy and other studies. We obtain UV LFs at $z \sim 4-7$ that span a very wide UV luminosity range of $\sim 0.002 - 100 \, L_{\rm UV}^\ast$ ($-26 < M_{\rm UV} < -14$ mag) by combining LFs from our program and the ultra-deep Hubble Space Telescope legacy surveys. We derive three parameters of the best-fit Schechter function, $\phi^\ast$, $M_{\rm UV}^\ast$, and $\alpha$, of the UV LFs in the magnitude range where the AGN contribution is negligible, and find that $\alpha$ and $\phi^\ast$ decrease from $z\sim 4$ to $7$ with no significant evolution of $M_{\rm UV}^\ast$. Because our HSC SSP data bridge the LFs of galaxies and AGNs with great statistical accuracy, we carefully investigate the bright end of the galaxy UV LFs that are estimated by the subtraction of the AGN contribution either aided with spectroscopy or the best-fit AGN UV LFs. We find that the bright end of the galaxy UV LFs cannot be explained by the Schechter function fits at $> 2 \sigma$ significance, and require either double power-law functions or modified Schechter functions that consider a magnification bias due to gravitational lensing.

Deep optical imaging with both Hyper Suprime-Cam and Suprime-Cam on the 8.2 m Subaru Telescope reveals a number of outer faint structures around the archetypical Seyfert galaxy NGC 1068 (M 77). We find three ultra diffuse objects (UDOs) around NGC 1068. Since these UDOs are located within the projected distance of 45 kpc from the center of NGC 1068, they appear to be associated with NGC 1068. Hereafter, we call them UDO-SW, UDO-NE, and UDO-SE where UDO = Ultra Diffuse Object, SW = south west, NE = north west, and SE = south east; note that UDO-SE was already found in the SDSS Stripe 82 data. Among them, both UDO-NE and UDO-SW appear to show a loop or stream structure around the main body of NGC 1068, providing evidence for the physical connection to NGC 1068. We consider that UDO-SE may be a tidal dwarf galaxy. We also find another UDO-like object that is 2 magnitudes fainter and smaller by a factor of 3 to 5 than those of the three UDOs. This object may belong to a class of low surface brightness galaxy. Since this object is located along the line connecting UDO-NE and UDO-SW, it is suggested that this object is related to the past interaction event that formed the loop by UDO-NE and UDO-SW, thus implying the physical connection to NGC 1068. Another newly-discovered feature is an asymmetric outer one-arm structure emanated from the western edge of the outermost disk of NGC 1068 together with a ripple-like structure at the opposite side. These structures are expected to arise in a late phase of a minor merger according to published numerical simulations of minor mergers. All these lines of evidence show that NGC 1068 experienced a minor merger several billions years ago. We then discuss the minor-merger driven triggering of nuclear activity in the case of NGC 1068.

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

In order to investigate the physical relationship between strong Ly-alpha absorbers (logN_HI>20.0 cm^-2) such as damped Ly-alpha absorption systems (DLAs) and young star-forming galaxies at high redshift, we have conducted narrow-band observations of Ly-alpha emitters (LAEs) in a concentrated region of strong Ly-alpha absorbers at z= 2.3, the J1230+34 field. Using a catalog of Ly-alpha absorbers with logN_HI>20.0 cm^-2 based on the baryon oscillation spectroscopic survey (BOSS), we found 6 fields where 3 or more absorbers are concentrated within a (50 Mpc)^3 cubic box in the comoving scale. Among them, we focus on the J1230+34 field, where 2 DLAs and 2 sub-DLAs present. Our narrow-band imaging observations with Subaru/Suprime-Cam using a custom-made filter, NB400 (lambda_c=4003 A and FWHM=92 A) yield a sample of 149 LAEs in this field. In the large scale (~50 Mpc), we have found no differences between the obtained Ly-alpha luminosity function and those in the blank fields at similar redshifts. We also compare the frequency distribution of the Ly-alpha rest-frame equivalent width (EW_0) in the target field and other fields including both overdensity region and blank field, but find no differences. On the other hand, in the small scale (~10 Mpc), we have found a possible overdensity of LAEs around a DLA with the highest HI column density (N_HI=21.08 cm^-2) in the target field while there are no density excess around the other absorbers with a lower N_HI.

We investigate the physical conditions of ionized gas in high-z star-forming galaxies using diagnostic diagrams based on the rest-frame optical emission lines. The sample consists of 701 galaxies with an Ha detection at $1.4\lesssim z\lesssim1.7$, from the FMOS-COSMOS survey, that represent the normal star-forming population over the stellar mass range $10^{9.6} \lesssim M_\ast/M_\odot \lesssim 10^{11.6}$ with those at $M_\ast>10^{11}~M_\odot$ being well sampled. We confirm an offset of the average location of star-forming galaxies in the BPT diagram ([OIII]/Hb vs. [NII]/Ha), primarily towards higher [OIII]/Hb, compared with local galaxies. Based on the [SII] ratio, we measure an electron density ($n_e=220^{+170}_{-130}~\mathrm{cm^{-3}}$), that is higher than that of local galaxies. Based on comparisons to theoretical models, we argue that changes in emission-line ratios, including the offset in the BPT diagram, are caused by a higher ionization parameter both at fixed stellar mass and at fixed metallicity with additional contributions from a higher gas density and possibly a hardening of the ionizing radiation field. Ionization due to AGNs is ruled out as assessed with Chandra. As a consequence, we revisit the mass-metallicity relation using [NII]/Ha and a new calibration including [NII]/[SII] as recently introduced by Dopita et al. Consistent with our previous results, the most massive galaxies ($M_\ast\gtrsim10^{11}~M_\odot$) are fully enriched, while those at lower masses have metallicities lower than local galaxies. Finally, we demonstrate that the stellar masses, metallicities and star formation rates of the FMOS sample are well fit with a physically-motivated model for the chemical evolution of star-forming galaxies.

One of the unsettled, important problems in active galactic nuclei (AGNs) is the major ionization mechanism of gas clouds in their narrow-line regions (NLRs). In order to investigate this issue, we present our $J$-band spectroscopic observations for a sample of 26 nearby Seyfert galaxies. In our study, we use the flux ratio of the following two forbidden emission lines, [Fe II]1.257$\mu$m and [P II]1.188$\mu$m because it is known that this ratio is sensitive to the ionization mechanism. We obtain the [Fe II]/[P II] flux ratio or its lower limit for 19 objects. In addition to our data, we compile this flux ratio (or its lower limit) of 23 nearby Seyfert galaxies from the literature. Based on the collected data, we find that three Seyfert galaxies show very large lower limits of [Fe II]/[P II] flux ratios ($\gtrsim$10); NGC 2782, 5005, and Mrk 463. It is thus suggested that the contribution of the fast shock in the gas excitation is significantly large for them. However, more than half of Seyfert galaxies in our sample show moderate [Fe II]/[P II] flux ratios ($\sim$2), that is consistent to the pure photoionization by power-law ionizing continuum emission. We also find that the [Fe II]/[P II] flux ratio shows no clear correlation with the radio loudness, suggesting that the radio jet is not the primary origin of shocks in NLRs of Seyfert galaxies.

Recent studies have derived quasar luminosity functions (QLFs) at various redshifts. However, the faint side of the QLF at high redshifts is still too uncertain. An accurate estimate of the survey completeness is essential to derive an accurate QLF for use in studying the luminosity-dependent density evolution of the quasar population. Here we investigate how the luminosity dependence of quasar spectra (the Baldwin effect) and the attenuation model for the inter-galactic medium (IGM) affect the completeness estimates. For this purpose, we revisit the completeness of quasar surveys specifically at $z\sim4-5$, using the COSMOS images observed with Subaru/Suprime-Cam. As the result, we find that the completeness estimates are sensitive to the luminosity dependence of the quasar spectrum and difference in the IGM attenuation models. At $z\sim4$, the number density of quasars when we adopt the latest IGM model and take the luminosity dependence of spectra into account are $(3.49\pm1.62)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-24.09<M_{1450}<-23.09$ and $(5.24\pm2.13)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-23.09<M_{1450}<-22.09$, respectively, which are $\sim$$24$$\%$ lower than that estimated by the conventional method. On the other hand, at $z\sim5$, a $1\sigma$ confidence upper limit of the number density at $-24.5<M_{1450}<-22.5$ in our new estimates is $\sim$$43$$\%$ higher than that estimated previously. The results suggest that the luminosity dependence of the quasar spectrum and the IGM model are important for deriving accurate number density of high-$z$ quasars. Even taking these effects into account, the inferred luminosity-dependent density evolution of quasars is consistent with the AGN down-sizing evolutionary picture.

Using the ALMA archival data of both CO(6--5) line and 689 GHz continuum emission towards the archetypical Seyfert galaxy, NGC 1068, we identified a distinct continuum peak separated by 14 pc from the nuclear radio component S1 in projection. The continuum flux gives a gas mass of ~2x10^5 Msun and bolometric luminosity of ~10^8 Lsun, leading to a star formation rate of ~0.1 Msun/yr. Subsequent analysis on the line data suggest that the gas has a size of ~10 pc, yielding to mean H2 number density of ~10^5 cm^{-3}. We therefore refer to the gas as "massive dense gas cloud": the gas density is high enough to form a "proto starcluster" whose stellar mass of ~10^4 Msun. We found that the gas stands a unique position between galactic and extraglactic clouds in the diagrams of start formation rate (SFR) vs. gas mass proposed by Lada et al. and surface density of gas vs. SFR density by Krumholz and McKee. All the gaseous and star-formation properties may be understood in terms of the turbulence-regulated star formation scenario. Since there are two stellar populations with the ages of 300 Myr and 30 Myr in the 100 pc-scale circumnulear region, we discuss that NGC1068 has experienced at least three episodic star formation events with a tendency that the inner star-forming region is the younger. Together with several lines of evidence that the dynamics of the nuclear region is decoupled from that of the entire galactic disk, we discuss that the gas inflow towards the nuclear region of NGC 1068 may be driven by a past minor merger.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

In this paper we investigate the impact of different star formation histories (SFHs) on the relation between stellar mass M$_{*}$ and star formation rate (SFR) using a sample of galaxies with reliable spectroscopic redshift zspec>2 drawn from the VIMOS Ultra-Deep Survey (VUDS). We produce an extensive database of dusty model galaxies, calculated starting from the new library of single stellar population (SSPs) models presented in Cassara' et al. 2013 and weighted by a set of 28 different SFHs based on the Schmidt function, and characterized by different ratios of the gas infall time scale $\tau_{infall}$ to the star formation efficiency $\nu$. The treatment of dust extinction and re-emission has been carried out by means of the radiative transfer calculation. The spectral energy distribution (SED) fitting technique is performed by using GOSSIP+, a tool able to combine both photometric and spectroscopic information to extract the best value of the physical quantities of interest, and to consider the Intergalactic Medium (IGM) attenuation as a free parameter. We find that the main contribution to the scatter observed in the $SFR-M_{*}$ plane is the possibility of choosing between different families of SFHs in the SED fitting procedure, while the redshift range plays a minor role. The majority of the galaxies, at all cosmic times, are best-fit by models with SFHs characterized by a high $\tau_{\rm infall}/\nu$ ratio. We discuss the reliability of the presence of a small percentage of dusty and highly star forming galaxies, in the light of their detection in the FIR.

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 the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

We present a new method to effectively select objects which may be low-mass active black holes (BHs) at galaxy centers using high-cadence optical imaging data, and our first spectroscopic identification of an active 2.7x10^6 Msun BH at z=0.164. This active BH was originally selected due to its rapid optical variability, from a few hours to a day, based on Subaru Hyper Suprime-Cam~(HSC) g-band imaging data taken with 1-hour cadence. Broad and narrow H-alpha and many other emission lines are detected in our optical spectra taken with Subaru FOCAS, and the BH mass is measured via the broad H-alpha emission line width (1,880 km s^{-1}) and luminosity (4.2x10^{40} erg s^{-1}) after careful correction for the atmospheric absorption around 7,580-7,720A. We measure the Eddington ratio to be as low as 0.05, considerably smaller than those in a previous SDSS sample with similar BH mass and redshift, which indicates one of the strong potentials of our Subaru survey. The g-r color and morphology of the extended component indicate that the host galaxy is a star-forming galaxy. We also show effectiveness of our variability selection for low-mass active BHs.

This paper describes the first data release (DR1) of the VIMOS Ultra Deep Survey (VUDS). The DR1 includes all low-resolution spectroscopic data obtained in 276.9 arcmin2 of the CANDELS-COSMOS and CANDELS-ECFDS survey areas, including accurate spectroscopic redshifts z_spec and individual spectra obtained with VIMOS on the ESO-VLT. A total of 698 objects have a measured redshift, with 677 galaxies, two type-I AGN and a small number of 19 contaminating stars. The targets of the spectroscopic survey are selected primarily on the basis of their photometric redshifts to ensure a broad population coverage. About 500 galaxies have z_spec>2, 48 with z_spec>4, and the highest reliable redshifts reach beyond z_spec=6. This dataset approximately doubles the number of galaxies with spectroscopic redshifts at z>3 in these fields. We discuss the general properties of the sample in terms of the spectroscopic redshift distribution, the distribution of Lyman-alpha equivalent widths, and physical properties including stellar masses M_star and star formation rates (SFR) derived from spectral energy distribution fitting with the knowledge of z_spec. We highlight the properties of the most massive star-forming galaxies, noting the large range in spectral properties, with Lyman-alpha in emission or in absorption, and in imaging properties with compact, multi-component or pair morphologies. We present the catalogue database and data products. All data are publicly available and can be retrieved from a dedicated query-based database available at http://cesam.lam.fr/vuds.

We investigate morphological properties of 61 Lyman-alpha emitters (LAEs) at z = 4.86 identified in the COSMOS field, based on Hubble Space Telescope Advanced Camera for Surveys (ACS) imaging data in the F814W-band. Out of the 61 LAEs, we find the ACS counterparts for the 54 LAEs. Eight LAEs show double-component structures with a mean projected separation of 0."63 (~ 4.0 kpc at z = 4.86). Considering the faintness of these ACS sources, we carefully evaluate their morphological properties, that is, size and ellipticity. While some of them are compact and indistinguishable from the PSF half-light radius of 0."07 (~ 0.45 kpc), the others are clearly larger than the PSF size and spatially extended up to 0."3 (~ 1.9 kpc). We find that the ACS sources show a positive correlation between ellipticity and size and that the ACS sources with large size and round shape are absent. Our Monte Carlo simulation suggests that the correlation can be explained by (1) the deformation effects via PSF broadening and shot noise or (2) the source blending in which two or more sources with small separation are blended in our ACS image and detected as a single elongated source. Therefore, the 46 single-component LAEs could contain the sources which consist of double (or multiple) components with small spatial separation (i.e., < 0."3 or 1.9 kpc). Further observation with high angular resolution at longer wavelengths (e.g., rest-frame wavelengths of > 4000 A) is inevitable to decipher which interpretation is adequate for our LAE sample.

We present rapidly rising transients discovered by a high-cadence transient survey with Subaru telescope and Hyper Suprime-Cam. We discovered five transients at z=0.384-0.821 showing the rising rate faster than 1 mag per 1 day in the restframe near-ultraviolet wavelengths. The fast rising rate and brightness are the most similar to SN 2010aq and PS1-13arp, for which the ultraviolet emission within a few days after the shock breakout was detected. The lower limit of the event rate of rapidly rising transients is ~9 % of core-collapse supernova rates, assuming a duration of rapid rise to be 1 day. We show that the light curves of the three faint objects agree with the cooling envelope emission from the explosion of red supergiants. The other two luminous objects are, however, brighter and faster than the cooling envelope emission. We interpret these two objects to be the shock breakout from dense wind with the mass loss rate of ~10^{-3} Msun yr^{-1}, as also proposed for PS1-13arp. This mass loss rate is higher than that typically observed for red supergiants. The event rate of these luminous objects is >~1 % of core-collapse supernova rate, and thus, our study implies that more than ~1 % of massive stars can experience an intensive mass loss at a few years before the explosion.

The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($\beta$) as well as Ly$\alpha$ equivalent widths (EW). On average, the spectroscopically measured $\beta$ (-1.36$\pm$0.02), is comparable to the photometrically measured $\beta$ (-1.32$\pm$0.02), and has smaller measurement uncertainties. The positive correlation of $\beta$ with the Spectral Energy Distribution (SED)-based measurement of dust extinction, E$_{\rm s}$(B-V), emphasizes the importance of $\beta$ as an alternative dust indicator at high redshifts. To make a proper comparison, we divide these SFGs into three subgroups based on their rest-frame Ly$\alpha$ EW: SFGs with no Ly$\alpha$ emission (SFG$_{\rm N}$; EW$\le$0\AA), SFGs with Ly$\alpha$ emission (SFG$_{\rm L}$; EW$>$0\AA), and Ly$\alpha$ emitters (LAEs; EW$\ge$20\AA). The fraction of LAEs at these redshifts is $\sim$10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SFG$_{\rm N}$, SFG$_{\rm L}$ and LAE samples. For the luminosities probed here ($\sim$L$^*$), we find that galaxies with and without Ly$\alpha$ in emission have small but significant differences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant compared to the SFR and E$_{\rm s}$(B-V). [abridged]

We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350 and 500 um) and AzTEC (1.1 mm) images. Our stacking procedure corrects the biases induced by galaxy clustering and incompleteness of our input catalogue in dense regions. We obtain the full IR spectral energy distributions (SED) of subsamples of LBGs and derive the mean IR luminosity as a function of UV luminosity, UV slope, and stellar mass. The average IRX is roughly constant over the UV luminosity range, with a mean of 7.9 (1.8 mag). However, it is correlated with UV slope, and stellar mass. We investigate using a statistically-controlled stacking analysis as a function of (stellar mass, UV slope) the dispersion of the IRX-UVslope and IRX-M* plane. Our results enable us to study the average relation between star-formation rate (SFR) and stellar mass, and we show that our LBG sample lies on the main sequence of star formation at z~3.

We present both the observations and the data reduction procedures of the Subaru COSMOS 20 project that is an optical imaging survey of the HST COSMOS field, carried out by using Suprime-Cam on the Subaru Telescope with the following 20 optical filters: 6 broad-band (B, g', V, r', i', and z'), 2 narrow-band (NB711 and NB816), and 12 intermediate-band filters (IA427, IA464, IA484, IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, and IA827). A part of this project is described in Taniguchi et al. (2007) and Capak et al. (2007) for the six broad-band and one narrow-band (NB816) filter data. In this paper, we present details of the observations and data reduction for remaining 13 filters (the 12 IA filters and NB711). In particular, we describe the accuracy of both photometry and astrometry in all the filter bands. We also present optical properties of the Suprime-Cam IA filter system in Appendix.

The Lyman continuum (LyC) flux escaping from high-z galaxies into the IGM is a fundamental quantity to understand the physical processes involved in the reionization epoch. We have investigated a sample of star-forming galaxies at z~3.3 in order to search for possible detections of LyC photons escaping from galaxy halos. UV deep imaging in the COSMOS field obtained with the prime focus camera LBC at the LBT telescope was used together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 45 galaxies at z~3.3 with L>0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 3.27<z<3.40 both in the R and deep U bands. A sub-sample of 10 galaxies apparently shows escape fractions>28% but a detailed analysis of their properties reveals that, with the exception of two marginal detections (S/N~2) in the U band, all the other 8 galaxies are most likely contaminated by the UV flux of low-z interlopers located close to the high-z targets. The average escape fraction derived from the stacking of the cleaned sample was constrained to fesc_rel<2%. The implied HI photo-ionization rate is a factor two lower than that needed to keep the IGM ionized at z~3, as observed in the Lyman forest of high-z QSO spectra or by the proximity effect. These results support a scenario where high redshift, relatively bright (L>0.5L*) star-forming galaxies alone are unable to sustain the level of ionization observed in the cosmic IGM at z~3. Star-forming galaxies at higher redshift and at fainter luminosities (L<<L*) can be the major contributors to the reionization of the Universe only if their physical properties are subject to rapid changes from z~3 to z~6-10. Alternatively, ionizing sources could be discovered looking for fainter sources among the AGN population at high-z.

We have analysed a sample of 574 Spitzer 4.5 micron-selected galaxies with [4.5]<23 and Ks_auto>24 (AB) over the UltraVISTA ultra-deep COSMOS field. Our aim is to investigate whether these mid-IR bright, near-IR faint sources contribute significantly to the overall population of massive galaxies at redshifts z>=3. By performing a spectral energy distribution (SED) analysis using up to 30 photometric bands, we have determined that the redshift distribution of our sample peaks at redshifts z~2.5-3.0, and ~32% of the galaxies lie at z>=3. We have studied the contribution of these sources to the galaxy stellar mass function (GSMF) at high redshifts. We found that the [4.5]<23, Ks_auto>24 galaxies produce a negligible change to the GSMF previously determined for Ks_auto<24 sources at 3=<z<4, but their contribution is more important at 4=<z<5, accounting for >~50% of the galaxies with stellar masses Mst>~6 x 10^10 Msun. We also constrained the GSMF at the highest-mass end (Mst>~2 x 10^11 Msun) at z>=5. From their presence at 5=<z<6, and virtual absence at higher redshifts, we can pinpoint quite precisely the moment of appearance of the first most massive galaxies as taking place in the ~0.2 Gyr of elapsed time between z~6 and z~5. Alternatively, if very massive galaxies existed earlier in cosmic time, they should have been significantly dust-obscured to lie beyond the detection limits of current, large-area, deep near-IR surveys.

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars at high redshift for the first time. In this study, we use the two-point angular cross-correlation function (CCF) for both low-luminosity quasars at $-24<M_{\rm 1450}<-22$ and Lyman-break galaxies (LBGs). Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low-luminosity quasars) in the redshift range $3.1<z<4.5$ and 835 color-selected LBGs with $z^{\prime}_{\rm LBG}<25.0$ at $z\sim4$ in the COSMOS field. We have made our analysis for the following two quasar samples; (1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample (the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity quasars at $z\sim4$ is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. We then obtain the $86\%$ upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the typical dark matter halo masses, log $(M_{\rm DM}/(h^{-1}M_{\odot}))=$$12.7$ and $13.5$, respectively. This result is not inconsistent with the predicted bias for quasars which is estimated by the major merger models.

We present a spectroscopic survey of galaxies in the COSMOS field using the Fiber Multi-Object Spectrograph (FMOS), a near-infrared instrument on the Subaru Telescope. Our survey is specifically designed to detect the Halpha emission line that falls within the H-band (1.6-1.8 um) spectroscopic window from star-forming galaxies with 1.4 < z < 1.7 and M_stellar>~10^10 Msolar. With the high multiplex capability of FMOS, it is now feasible to construct samples of over one thousand galaxies having spectroscopic redshifts at epochs that were previously challenging. The high-resolution mode (R~2600) effectively separates Halpha and [NII]6585 thus enabling studies of the gas-phase metallicity and photoionization state of the interstellar medium. The primary aim of our program is to establish how star formation depends on stellar mass and environment, both recognized as drivers of galaxy evolution at lower redshifts. In addition to the main galaxy sample, our target selection places priority on those detected in the far-infrared by Herschel/PACS to assess the level of obscured star formation and investigate, in detail, outliers from the star formation rate - stellar mass relation. Galaxies with Halpha detections are followed up with FMOS observations at shorter wavelengths using the J-long (1.11-1.35 um) grating to detect Hbeta and [OIII]5008 that provides an assessment of extinction required to measure star formation rates not hampered by dust, and an indication of embedded Active Galactic Nuclei. With 460 redshifts measured from 1153 spectra, we assess the performance of the instrument with respect to achieving our goals, discuss inherent biases in the sample, and detail the emission-line properties. Our higher-level data products, including catalogs and spectra, are available to the community.

We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$\alpha$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history.

We present the results of a search for bright (-22.7 < M_UV < -20.5) Lyman-break galaxies at z ~ 6 within a total of 1.65 square degrees of imaging in the UltraVISTA/COSMOS and UKIDSS UDS/SXDS fields. The deep near-infrared imaging available in the two independent fields, in addition to deep optical (including z'-band) data, enables the sample of z ~ 6 star-forming galaxies to be securely detected long-ward of the break (in contrast to several previous studies). We show that the expected contamination rate of our initial sample by cool galactic brown dwarfs is < 3 per cent and demonstrate that they can be effectively removed by fitting brown dwarf spectral templates to the photometry. At z ~ 6 the galaxy surface density in the UltraVISTA field exceeds that in the UDS by a factor of ~ 1.8, indicating strong cosmic variance even between degree-scale fields at z > 5. We calculate the bright end of the rest-frame Ultra-Violet (UV) luminosity function (LF) at z ~ 6. The galaxy number counts are a factor of ~1.7 lower than predicted by the recent LF determination by Bouwens et al.. In comparison to other smaller area studies, we find an evolution in the characteristic magnitude between z ~ 5 and z ~ 7 of dM* ~ 0.4 mag, and show that a double power-law or a Schechter function can equally well describe the LF at z = 6. Furthermore, the bright-end of the LF appears to steepen from z ~ 7 to z ~ 5, which could indicate the onset of mass quenching or the rise of dust obscuration, a conclusion supported by comparing the observed LFs to a range of theoretical model predictions.

Dust-Obscured galaxies (DOGs) are bright 24 um-selected sources with extreme obscuration at optical wavelengths. They are typically characterized by a rising power-law continuum of hot dust (T_D ~ 200-1000K) in the near-IR indicating that their mid-IR luminosity is dominated by an an active galactic nucleus (AGN). DOGs with a fainter 24 um flux display a stellar bump in the near-IR and their mid-IR luminosity appears to be mainly powered by dusty star formation. Alternatively, it may be that the mid-IR emission arising from AGN activity is dominant but the torus is sufficiently opaque to make the near-IR emission from the AGN negligible with respect to the emission from the host component. In an effort to characterize the astrophysical nature of the processes responsible for the IR emission in DOGs, this paper exploits Herschel data (PACS + SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colors; dust temperatures and masses) based on SED fitting. Of particular interest are the 24 um-bright DOGs (F_24um > 1mJy). They present bluer far-IR/mid-IR colors than the rest of the sample, unveiling the potential presence of an AGN. The AGN contribution to the total 8-1000 um flux increases as a function of the rest-frame 8 um-luminosity irrespective of the redshift. This confirms that faint DOGs (L_8 um< 10^12 L_sun) are dominated by star-formation while brighter DOGs show a larger contribution from an AGN.

We have used FMOS on Subaru to obtain near-infrared spectroscopy of 123 far-infrared selected galaxies in COSMOS and obtain the key rest-frame optical emission lines. This is the largest sample of infrared galaxies with near-infrared spectroscopy at these redshifts. The far-infrared selection results in a sample of galaxies that are massive systems that span a range of metallicities in comparison with previous optically selected surveys, and thus has a higher AGN fraction and better samples the AGN branch. We establish the presence of AGN and starbursts in this sample of (U)LIRGs selected as Herschel-PACS and Spitzer-MIPS detections in two redshift bins (z~0.7 and z~1.5) and test the redshift dependence of diagnostics used to separate AGN from star-formation dominated galaxies. In addition, we construct a low redshift (z~0.1) comparison sample of infrared selected galaxies and find that the evolution from z~1.5 to today is consistent with an evolving AGN selection line and a range of ISM conditions and metallicities from the models of Kewley et al. (2013b). We find that a large fraction of (U)LIRGs are BPT-selected AGN using their new, redshift-dependent classification line. We compare the position of known X-ray detected AGN (67 in total) with the BPT selection and find that the new classification line accurately selects most of these objects (> 70%). Furthermore, we identify 35 new (likely obscured) AGN not selected as such by their X-ray emission. Our results have direct implications for AGN selection at higher redshift with either current (MOSFIRE, KMOS) or future (PFS, MOONS) spectroscopic efforts with near-infrared spectral coverage.

We report the discovery of 31 low-luminosity (-14.5 > M_{AB}(B) > -18.8), extreme emission line galaxies (EELGs) at 0.2 < z < 0.9 identified by their unusually high rest-frame equivalent widths (100 < EW[OIII] < 1700 A) as part of the VIMOS Ultra Deep Survey (VUDS). VIMOS optical spectra of unprecedented sensitivity ($I_{AB}$ ~ 25 mag) along with multiwavelength photometry and HST imaging are used to investigate spectrophotometric properties of this unique sample and explore, for the first time, the very low stellar mass end (M* < 10^8 M$_{\odot}$) of the luminosity-metallicity (LZR) and mass-metallicity (MZR) relations at z < 1. Characterized by their extreme compactness (R50 < 1 kpc), low stellar mass and enhanced specific star formation rates (SFR/M* ~ 10^{-9} - 10^{-7} yr^{-1}), the VUDS EELGs are blue dwarf galaxies likely experiencing the first stages of a vigorous galaxy-wide starburst. Using T_e-sensitive direct and strong-line methods, we find that VUDS EELGs are low-metallicity (7.5 < 12+log(O/H) < 8.3) galaxies with high ionization conditions, including at least three EELGs showing HeII 4686A emission and four EELGs of extremely metal-poor (<10% solar) galaxies. The LZR and MZR followed by EELGs show relatively large scatter, being broadly consistent with the extrapolation toward low luminosity and mass from previous studies at similar redshift. However, we find evidences that galaxies with younger and more vigorous star formation -- as characterized by their larger EWs, ionization and sSFR -- tend to be more metal-poor at a given stellar mass.

The relation between the stellar mass and the star formation rate characterizes how the instantaneous star formation is determined by the galaxy past star formation history and by the growth of the dark matter structures. We deconstruct the M-SFR plane by measuring the specific SFR functions in several stellar mass bins from z=0.2 out to z=1.4. Our analysis is primary based on a MIPS 24$\mu m$ selected catalogue combining the COSMOS and GOODS surveys. We estimate the SFR by combining mid- and far-infrared data for 20500 galaxies. The sSFR functions are derived in four stellar mass bins within the range 9.5<log(M/Msun)<11.5. First, we demonstrate the importance of taking into account selection effects when studying the M-SFR relation. Secondly, we find a mass-dependent evolution of the median sSFR with redshift varying as $sSFR \propto (1+z)^{b}$, with $b$ increasing from $b=2.88$ to $b=3.78$ between $M=10^{9.75}Msun$ and $M=10^{11.1}Msun$, respectively. At low masses, this evolution is consistent with the cosmological accretion rate and predictions from semi-analytical models (SAM). This agreement breaks down for more massive galaxies showing the need for a more comprehensive description of the star-formation history in massive galaxies. Third, we obtain that the shape of the sSFR function is invariant with time at z<1.4 but depends on the mass. We observe a broadening of the sSFR function ranging from 0.28 dex at $M=10^{9.75}Msun$ to 0.46 dex at $M=10^{11.1}Msun$. Such increase in the scatter of the M-SFR relation suggests an increasing diversity of SFHs as the stellar mass increases. Finally, we find a gradual decline of the sSFR with mass as $log(sSFR) \propto -0.17M$. We discuss the numerous physical processes, as gas exhaustion in hot gas halos or secular evolution, which can gradually reduce the sSFR and increase the SFH diversity.

Using the Hubble Space Telescope/Wide Field Camera 3 imaging data and multi-wavelength photometric catalog, we investigated the dust temperature of passively evolving and star-forming galaxies at 0.2<z<1.0 in the CANDELS fields. We estimated the stellar radiation field by low-mass stars from the stellar mass and surface brightness profile of these galaxies and then calculated their steady-state dust temperature. At first, we tested our method using nearby early-type galaxies with the deep FIR data by the Herschel Virgo cluster survey and confirmed that the estimated dust temperatures are consistent with the observed temperatures within the uncertainty. We then applied the method to galaxies at 0.2<z<1.0, and found that most of passively evolving galaxies with Mstar > 10^{10} Msun have a relatively high dust temperature of Tdust > 20 K, for which the formation efficiency of molecular hydrogen on the surface of dust grains in the diffuse ISM is expected to be very low from the laboratory experiments. The fraction of passively evolving galaxies strongly depends on the expected dust temperature at all redshifts and increases rapidly with increasing the temperature around Tdust ~ 20 K. These results suggest that the dust heating by low-mass stars in massive galaxies plays an important role for the continuation of their passive evolution, because the lack of the shielding effect of the molecular hydrogen on the UV radiation can prevent the gas cooling and formation of new stars.

The relation between the galaxy stellar mass M_star and the dark matter halo mass M_h gives important information on the efficiency in forming stars and assembling stellar mass in galaxies. We present the stellar mass to halo mass ratio (SMHR) measurements at redshifts 2<z<5, obtained from the VIMOS Ultra Deep Survey. We use halo occupation distribution (HOD) modelling of clustering measurements on ~3000 galaxies with spectroscopic redshifts to derive the dark matter halo mass M_h, and SED fitting over a large set of multi-wavelength data to derive the stellar mass M_star and compute the SMHR=M_star/M_h. We find that the SMHR ranges from 1% to 2.5% for galaxies with M_star=1.3x10^9 M_sun to M_star=7.4x10^9 M_sun in DM halos with M_h=1.3x10^{11} M_sun} to M_h=3x10^{11} M_sun. We derive the integrated star formation efficiency (ISFE) of these galaxies and find that the star formation efficiency is a moderate 6-9% for lower mass galaxies while it is relatively high at 16% for galaxies with the median stellar mass of the sample ~7x10^9 M_sun. The lower ISFE at lower masses may indicate that some efficient means of suppressing star formation is at work (like SNe feedback), while the high ISFE for the average galaxy at z~3 is indicating that these galaxies are efficiently building-up their stellar mass at a key epoch in the mass assembly process. We further infer that the average mass galaxy at z~3 will start experiencing star formation quenching within a few hundred millions years.

We present a new chemical evolution model for dwarf spheroidal galaxies (dSphs) in the Local Universe. Our main aim is to explain both their observed star formation histories and metallicity distribution functions simultaneously. Applying our new model for the four local dSphs, that is, Fornax, Sculptor, Leo II, and Sextans, we find that our new model reproduces the observed chemical properties of the dSphs consistently. Our results show that the dSphs have evolved with both a low star formation efficiency and a large gas outflow efficiency compared with the Milky Way, as suggested by previous works. Comparing the observed [alpha/Fe]--[Fe/H] relation of the dSphs with the model predictions, we find that our model favors a longer onset time of Type Ia supernovae (i.e., 0.5 Gyr) than that suggested in previous studies (i.e., 0.1 Gyr). We discuss the origin of this discrepancy in detail.

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$\alpha$ halo at z=4.11. We used a sample of Ly$\alpha$ emitters (LAEs) down to $\log(L_{\rm Ly\alpha} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $\delta_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$\alpha$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Ly\alpha}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$\alpha$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings.

Using new spectroscopic observations obtained as part of the VIMOS Ultra-Deep Survey (VUDS), we perform a systematic search for overdense environments in the early universe ($z>2$) and report here on the discovery of Cl J0227-0421, a massive protocluster at $z=3.29$. This protocluster is characterized by both the large overdensity of spectroscopically confirmed members, $\delta_{gal}=10.5\pm2.8$, and a significant overdensity in photometric redshift members. The halo mass of this protocluster is estimated, by a variety of methods, to be roughly $3\times10^{14}$ $\mathcal{M}_{\odot}$ at $z\sim3.3$, which, evolved to $z=0$ results in a halo mass rivaling or exceeding that of the Coma cluster. The properties of 19 spectroscopically confirmed member galaxies are compared with a large sample of VUDS/VVDS galaxies in lower density field environments at similar redshifts. We find tentative evidence for an excess of redder, brighter, and more massive galaxies within the confines of the protocluster relative to the field population, which suggests that we may be observing the beginning of environmentally-induced quenching. The properties of these galaxies are investigated, including a discussion of the brightest protocluster galaxy which appears to be undergoing vigorous coeval nuclear and starburst activity. The remaining member galaxies appear to have characteristics which are largely similar to the field population. Though we find weaker evidence of the suppression of the median star formation rates amongst and differences in stacked spectra of member galaxies with respect to the field, we defer any conclusions of these trends to future work with the ensemble of protostructures that are found in the full VUDS sample.

We investigate the evolution of galaxy clustering for galaxies in the redshift range 2.0<$z$<5.0 using the VIMOS Ultra Deep Survey (VUDS). We present the projected (real-space) two-point correlation function $w_p(r_p)$ measured by using 3022 galaxies with robust spectroscopic redshifts in two independent fields (COSMOS and VVDS-02h) covering in total 0.8 deg$^2$. We quantify how the scale dependent clustering amplitude $r_0$ changes with redshift making use of mock samples to evaluate and correct the survey selection function. Using a power-law model $\xi(r) = (r/r_0)^{-\gamma}$ we find that the correlation function for the general population is best fit by a model with a clustering length $r_0$=3.95$^{+0.48}_{-0.54}$ h$^{-1}$Mpc and slope $\gamma$=1.8$^{+0.02}_{-0.06}$ at $z$~2.5, $r_0$=4.35$\pm$0.60 h$^{-1}$Mpc and $\gamma$=1.6$^{+0.12}_{-0.13}$ at $z$~3.5. We use these clustering parameters to derive the large-scale linear galaxy bias $b_L^{PL}$, between galaxies and dark matter. We find $b_L^{PL}$ = 2.68$\pm$0.22 at redshift $z$~3 (assuming $\sigma_8$ = 0.8), significantly higher than found at intermediate and low redshifts. We fit an HOD model to the data and we obtain that the average halo mass at redshift $z$~3 is $M_h$=10$^{11.75\pm0.23}$ h$^{-1}$M$_{\odot}$. From this fit we confirm that the large-scale linear galaxy bias is relatively high at $b_L^{HOD}$ = 2.82$\pm$0.27. Comparing these measurements with similar measurements at lower redshifts we infer that the star-forming population of galaxies at $z$~3 should evolve into the massive and bright ($M_r$<-21.5) galaxy population which typically occupy haloes of mass $\langle M_h\rangle$ = 10$^{13.9}$ h$^{-1}$ $M_{\odot}$ at redshift $z$=0.

We study the evolution of the star formation rate (SFR) - stellar mass (M_star) relation and specific star formation rate (sSFR) of star forming galaxies (SFGs) since a redshift z~5.5 using 2435 (4531) galaxies with highly reliable (reliable) spectroscopic redshifts in the VIMOS Ultra-Deep Survey (VUDS). It is the first time that these relations can be followed over such a large redshift range from a single homogeneously selected sample of galaxies with spectroscopic redshifts. The log(SFR) - log(M_star) relation for SFGs remains roughly linear all the way up to z=5 but the SFR steadily increases at fixed mass with increasing redshift. We find that for stellar masses M_star>3.2 x 10^9 M_sun the SFR increases by a factor ~13 between z=0.4 and z=2.3. We extend this relation up to z=5, finding an additional increase in SFR by a factor 1.7 from z=2.3 to z=4.8 for masses M_star > 10^10 M_sun. We observe a turn-off in the SFR-M_star relation at the highest mass end up to a redshift z~3.5. We interpret this turn-off as the signature of a strong on-going quenching mechanism and rapid mass growth. The sSFR increases strongly up to z~2 but it grows much less rapidly in 2<z<5. We find that the shape of the sSFR evolution is not well reproduced by cold gas accretion-driven models or the latest hydrodynamical models. Below z~2 these models have a flatter evolution (1+z)^{Phi} with Phi=2-2.25 compared to the data which evolves more rapidly with Phi=2.8+-0.2. Above z~2, the reverse is happening with the data evolving more slowly with Phi=1.2+-0.1. The observed sSFR evolution over a large redshift range 0<z<5 and our finding of a non linear main sequence at high mass both indicate that the evolution of SFR and M_star is not solely driven by gas accretion. The results presented in this paper emphasize the need to invoke a more complex mix of physical processes {abridge}

(arXiv abridged abstract) The observed UV rest-frame spectra of distant galaxies are the result of their intrinsic emission combined with absorption along the line of sight produced by the inter-galactic medium (IGM). Here we analyse the evolution of the mean IGM transmission Tr(Ly_alpha) and its dispersion along the line of sight for 2127 galaxies with 2.5<z<5.5 in the VIMOS Ultra Deep Survey (VUDS). We fit model spectra combined with a range of IGM transmission to the galaxy spectra using the spectral fitting algorithm GOSSIP+. We use these fits to derive the mean IGM transmission towards each galaxy for several redshift slices from z=2.5 to z=5.5. We find that the mean IGM transmission defined as Tr(Ly_alpha)=e^{-tau} (with tau the HI optical depth) is 79%, 69%, 59%, 55% and 46% at redshifts 2.75, 3,22, 3.70, 4.23, 4.77, respectively. We compare these results to measurements obtained from quasars lines of sight and find that the IGM transmission towards galaxies is in excellent agreement with quasar values up to redshift z~4. We find tentative evidence for a higher IGM transmission at z>= 4 compared to results from QSOs, but a degeneracy between dust extinction and IGM prevents to draw firm conclusions if the internal dust extinction for star-forming galaxies at z>4 takes a mean value significantly in excess of E(B-V)>0.15. Most importantly, we find a large dispersion of IGM transmission along the lines of sight towards distant galaxies with 68% of the distribution within 10 to 17% of the median value in delta z=0.5 bins, similar to what is found on the LOS towards QSOs. We demonstrate the importance of taking into account this large range of IGM transmission when selecting high redshift galaxies based on their colour properties (e.g. LBG or photometric redshift selection) or otherwise face a significant incompleteness in selecting high redshift galaxy populations.

Aims. The aim of this work is to constrain the evolution of the fraction of Lya emitters among UV selected star forming galaxies at 2<z<6, and to measure the stellar escape fraction of Lya photons over the same redshift range. Methods. We exploit the ultradeep spectroscopic observations collected by the VIMOS Ultra Deep Survey (VUDS) to build an unique, complete and unbiased sample of 4000 spectroscopically confirmed star forming galaxies at 2<z<6. Our galaxy sample UV luminosities brighter than M* at 2<z<6, and luminosities down to one magnitude fainter than M* at 2<z<3.5. Results. We find that 80% of the star forming galaxies in our sample have EW0(Lya)<10A, and correspondingly fesc(Lya)<1%. By comparing these results with literature, we conclude that the bulk of the Lya luminosity at 2<z<6 comes from galaxies that are fainter in the UV than those we sample in this work. The strong Lya emitters constitute, at each redshift, the tail of the distribution of the galaxies with extreme EW0(Lya) and fesc(Lya) . This tail of large EW0 and fesc(Lya) becomes more important as the redshift increases, and causes the fraction of Lya with EW0> 25A to increase from 5% at z=2 to 30% at z=6, with the increase being relatively stronger beyond z=4. We observe no difference, for the narrow range of UV luminosities explored in this work, between the fraction of strong Lya emitters among galaxies fainter or brighter than M*, although the fraction for the FUV faint galaxies evolves faster, at 2<z<3.5, than for the bright ones. We do observe an anticorrelation between E(B-V) and fesc(Lya): generally galaxies with high fesc(Lya) have also small amounts of dust (and viceversa). However, when the dust content is low (E(B-V)<0.05) we observe a very broad range of fesc(Lya), ranging from 10^-3 to 1. This implies that the dust alone is not the only regulator of the amount of escaping Lya photons.

We present the first observations of foreground Lyman-$\alpha$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$\alpha$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$\alpha$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

[Abridged] We characterise a massive proto-cluster at z=2.895 that we found in the COSMOS field using the spectroscopic sample of the VIMOS Ultra-Deep Survey (VUDS). This is one of the rare structures at z~3 not identified around AGNs or radio galaxies, so it is an ideal laboratory to study galaxy formation in dense environments. The structure comprises 12 galaxies with secure spectroscopic redshift in an area of 7'x8', in a z bin of Dz=0.016. The measured galaxy number overdensity is delta_g=12+/-2. This overdensity has total mass of M~8.1x10^(14)M_sun in a volume of 13x15x17 Mpc^3. Simulations indicate that such an overdensity at z~2.9 is a proto-cluster that will collapse in a cluster of total mass M~2.5x10^(15)M_sun at z=0. We compare the properties of the galaxies within the overdensity with a control sample at the same z but outside the overdensity. We did not find any statistically significant difference between the properties (stellar mass, SFR, sSFR, NUV-r, r-K) of the galaxies inside and outside the overdensity. The stacked spectrum of galaxies in the overdensity background shows a significant absorption feature at the wavelength of Lya redshifted at z=2.895 (lambda=4736 A), with a rest frame EW = 4+/- 1.4 A. Stacking only background galaxies without intervening sources at z~2.9 along their line of sight, we find that this absorption feature has a rest frame EW of 10.8+/-3.7 A, with a detection S/N of ~4. These EW values imply a high column density (N(HI)~3-20x10^(19)cm^(-2)), consistent with a scenario where such absorption is due to intervening cold gas streams, falling into the halo potential wells of the proto-cluster galaxies. However, we cannot exclude the hypothesis that this absorption is due to the diffuse gas within the overdensity.

Using the first 50% of data collected for the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) observations on the 1.8 deg$^2$ Cosmological Evolution Survey (COSMOS) we estimate the masses and star formation rates of 3398 $M_*>10^{10}M_\odot $ star-forming galaxies at $4 < z < 6$ with a substantial population up to $M_* \gtrsim 10^{11.5} M_\odot$. We find that the strong correlation between stellar mass and star formation rate seen at lower redshift (the "main sequence" of star-forming galaxies) extends to $z\sim6$. The observed relation and scatter is consistent with a continued increase in star formation rate at fixed mass in line with extrapolations from lower-redshift observations. It is difficult to explain this continued correlation, especially for the most massive systems, unless the most massive galaxies are forming stars near their Eddington-limited rate from their first collapse. Furthermore, we find no evidence for moderate quenching at higher masses, indicating quenching either has not occurred prior to $z \sim 6$ or else occurs rapidly, so that few galaxies are visible in transition between star-forming and quenched.

We investigate the relationships between stellar mass, gas-phase oxygen abundance (metallicity), star formation rate, and dust content of star-forming galaxies at z$\sim$1.6 using Subaru/FMOS spectroscopy in the COSMOS field. The mass-metallicity relation at $z\sim1.6$ is steeper than the relation observed in the local Universe. The steeper MZ relation at $z\sim1.6$ is mainly due to evolution in the stellar mass where the MZ relation begins to turnover and flatten. This turnover mass is 1.2 dex larger at $z\sim1.6$. The most massive galaxies at $z\sim1.6$ ($\sim 10^{11}M_\odot$) are enriched to the level observed in massive galaxies in the local Universe. The mass-metallicity relation we measure at $z\sim1.6$ supports the suggestion of an empirical upper metallicity limit that does not significantly evolve with redshift. We find an anti-correlation between metallicity and star formation rate for galaxies at a fixed stellar mass at $z\sim1.6$ which is similar to trends observed in the local Universe. We do not find a relation between stellar mass, metallicity and star formation rate that is independent of redshift; our data suggest that there is redshift evolution in this relation. We examine the relation between stellar mass, metallicity and dust extinction. We find that at a fixed stellar mass dustier galaxies tend to be more metal rich. From examination of the stellar masses, metallicities, SFRs and dust extinctions we conclude that stellar mass is most closely related to dust extinction.

We studied the chronology of galactic bulge and disc formation by analysing the relative contributions of these components to the B-band rest-frame luminosity density at different epochs. We present the first estimate of the evolution of the fraction of rest-frame B-band light in galactic bulges and discs since redshift z~0.8. We performed a bulge-to-disc decomposition of HST/ACS images of 3266 galaxies in the zCOSMOS-bright survey with spectroscopic redshifts in the range 0.7 < z < 0.9. We find that the fraction of B-band light in bulges and discs is $(26 \pm 4)%$ and $(74 \pm 4)%$, respectively. When compared with rest-frame B-band measurements of galaxies in the local Universe in the same mass range ($10^{9} M_{\odot}\lessapprox M \lessapprox 10^{11.5} M_{\odot}$), we find that the B-band light in discs decreases by ~30% from z~0.7-0.9 to z~0, while the light from the bulge increases by ~30% over the same period of time. We interpret this evolution as the consequence of star formation and mass assembly processes, as well as morphological transformation, which gradually shift stars formed at half the age of the Universe from star-forming late-type/irregular galaxies toearlier types and ultimately into spheroids.

We present the first results from a near-IR spectroscopic survey of the COSMOS field, using the Fiber Multi-Object Spectrograph on the Subaru telescope, designed to characterize the star-forming galaxy population at $1.4<z<1.7$. The high-resolution mode is implemented to detect H$\alpha$ in emission between $1.6{\rm -}1.8 \mathrm{\mu m}$ with $f_{\rm H\alpha}\gtrsim4\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$. Here, we specifically focus on 271 sBzK-selected galaxies that yield a H$\alpha$ detection thus providing a redshift and emission line luminosity to establish the relation between star formation rate and stellar mass. With further $J$-band spectroscopy for 89 of these, the level of dust extinction is assessed by measuring the Balmer decrement using co-added spectra. We find that the extinction ($0.6\lesssim A_\mathrm{H\alpha} \lesssim 2.5$) rises with stellar mass and is elevated at high masses compared to low-redshift galaxies. Using this subset of the spectroscopic sample, we further find that the differential extinction between stellar and nebular emission \hbox{$E_\mathrm{star}(B-V)/E_\mathrm{neb}(B-V)$} is 0.7--0.8, dissimilar to that typically seen at low redshift. After correcting for extinction, we derive an H$\alpha$-based main sequence with a slope ($0.81\pm0.04$) and normalization similar to previous studies at these redshifts.

We present the VIMOS Ultra Deep Survey (VUDS), a spectroscopic redshift survey of ~10.000 very faint galaxies to study the major phase of galaxy assembly 2<z<~6. The survey covers 1 deg^2 in 3 separate fields: COSMOS, ECDFS and VVDS-02h, with targets selection based on an inclusive combination of photometric redshifts and color properties. Spectra covering 3650<lambda<9350 A are obtained with VIMOS on the ESO-VLT with integration times of 14h. Here we present the survey strategy, the target selection, the data processing, as well as the redshift measurement process, emphasizing the specific methods adapted to this high redshift range. The spectra quality and redshift reliability are discussed, and we derive a completeness in redshift measurement of 91%, or 74% for the most reliable measurements, down to i_AB=25, and measurements are performed all the way down to i_AB=27. The redshift distribution of the main sample peaks at z=3-4 and extends over a large redshift range mainly in 2 < z < 6. At 3<z<5, the galaxies cover a large range of luminosities -23< M_U < -20.5, stellar mass 10^9 M_sun< M_star < 10^{11} M_sun, and star formation rates 1 M_sun/yr< SFR < 10^3 M_sun/yr. We discuss the spectral properties of galaxies using individual as well as stacked spectra. The comparison between spectroscopic and photometric redshifts as well as color selection demonstrate the effectiveness of our selection scheme. With ~6000 galaxies with reliable spectroscopic redshifts in 2<z<6 expected when complete, this survey is the largest at these redshifts and offers the opportunity for unprecedented studies of the star-forming galaxy population and its distribution in large scale structures during the major phase of galaxy assembly.

Using the Hubble Space Telescope/Advanced Camera for Surveys data in the COSMOS field, we systematically searched clumpy galaxies at 0.2<z<1.0 and investigated the fraction of clumpy galaxies and its evolution as a function of stellar mass, star formation rate (SFR), and specific SFR (SSFR). The fraction of clumpy galaxies in star-forming galaxies with Mstar > 10^9.5 Msun decreases with time from ~0.35 at 0.8<z<1.0 to ~0.05 at 0.2<z<0.4 irrespective of the stellar mass, although the fraction tends to be slightly lower for massive galaxies with Mstar > 10^10.5 Msun at each redshift. On the other hand, the fraction of clumpy galaxies increases with increasing both SFR and SSFR in all the redshift ranges we investigated. In particular, we found that the SSFR dependences of the fractions are similar among galaxies with different stellar masses, and the fraction at a given SSFR does not depend on the stellar mass in each redshift bin. The evolution of the fraction of clumpy galaxies from z~0.9 to z~0.3 seems to be explained by such SSFR dependence of the fraction and the evolution of SSFRs of star-forming galaxies. The fraction at a given SSFR also appears to decrease with time, but this can be due to the effect of the morphological K-correction. We suggest that these results are understood by the gravitational fragmentation model for the formation of giant clumps in disk galaxies, where the gas mass fraction is a crucial parameter.

We present the results of a new search for bright star-forming galaxies at z ~ 7 within the UltraVISTA DR2 and UKIDSS UDS DR10 data, which together provide 1.65 sq deg of near-infrared imaging with overlapping optical and Spitzer data. Using a full photo-z analysis to identify high-z galaxies and reject contaminants, we have selected a sample of 34 luminous (-22.7 < M_UV < -21.2) galaxies with 6.5 < z < 7.5. Crucially, the deeper imaging provided by UltraVISTA DR2 confirms all of the robust objects previously uncovered by Bowler et al. (2012), validating our selection technique. Our sample includes the most massive galaxies known at z ~ 7, with M_* ~ 10^{10} M_sun, and the majority are resolved, consistent with larger sizes (r_{1/2} ~ 1 - 1.5 kpc) than displayed by less massive galaxies. From our final sample, we determine the form of the bright end of the rest-frame UV galaxy luminosity function (LF) at z ~ 7, providing strong evidence that the bright end of the z = 7 LF does not decline as steeply as predicted by the Schechter function fitted to fainter data. We consider carefully, and exclude the possibility that this is due to either gravitational lensing, or significant contamination of our galaxy sample by AGN. Rather, our results favour a double power-law form for the galaxy LF at high z, or, more interestingly, a LF which simply follows the form of the dark-matter halo mass function at bright magnitudes. This suggests that the physical mechanism which inhibits star-formation activity in massive galaxies (i.e. AGN feedback or some other form of `mass quenching') has yet to impact on the observable galaxy LF at z ~ 7, a conclusion supported by the estimated masses of our brightest galaxies which have only just reached a mass comparable to the critical `quenching mass' of M_* = 10 ^{10.2} M_sun derived from studies of the mass function of star-forming galaxies at lower z.

A provocative summary of the conference on "Galaxy Mergers in an Evolving Universe" is presented. The conference includes the following issues; 1) Mergers at low redshift, 2) Merger dynamics and numerical simulations, 3) Mergers in galaxy evolution, 4) Mergers and active galactic nuclei (AGNs), 5) Mergers at high redshift and submm galaxies, 6) Feedback from starbursts and AGNs, 6) Chemical evolution of galaxies, and 7) Results from new surveys. Although there are many important issues related to galaxy mergers and nuclear activities in galaxies, in this summary, l et us focus on following the question; {\it Is there any unified model for triggering nuclear activities ?} Given key observational properties of galaxies with AGNs, a possible unified model is discussed as an evolutionary model from starbursts to AGNs based on minor and major mergers between/among galaxies.

In order to understand nature of building blocks of galaxies in the early universe, we investigate "genuine irregular galaxies (GIGs)" in the nearby universe. Here, GIGs are defined as isolated galaxies without regular structures (spheroid, bulge, disk, bar, spiral arm, and nucleus). Using the results of two excellent studies on galaxy morphology based on the Sloan Digital Sky Survey (SDSS), we obtain a sample of 66 irregular galaxies. We carry out new classification of them into GIGs and non-GIGs which have regular structure or show evidence for galaxy interaction, by using the SDSS Data Release 10 images. We then find that a half of these irregular galaxies (33/66) are GIGs and obtain an unambiguous sample of 33 GIGs for the first time. We discuss their observational properties by comparing them with those of elliptical, S0, spiral galaxies, and irregular galaxies without the GIGs. We find that our GIGs have smaller sizes, lower optical luminosities, bluer rest-frame optical colors, lower surface stellar mass densities, and lower gas metallicity than normal galaxies. All these properties suggest that they are in chemically and dynamically younger phases even in the nearby universe.

Abridged. We use COSMOS to study in a self-consistent way the change in the number densities of quenched early-type galaxies (Q-ETGs) of a given size over the interval 0.2 < z < 1.0 to study the claimed size evolution of these galaxies. At 10^10.5<Mgalaxy<10^11 Msun, we see no change in the number density of compact Q-ETGs, while at >10^11 Msun we find a decrease by 30%. In both mass bins, the increase of the median sizes of Q-ETGs with time is primarily caused by the addition to the size function of larger and more diffuse Q-ETGs. At all masses, compact Q-ETGs become systematically redder towards later epochs, with a (U-V) difference consistent with passive evolution of their stellar populations, indicating that they are a population that does not appreciably evolve in size. At all epochs, the larger Q-ETGs (at least in the lower mass bin) have average rest-frame colors systematically bluer than those of the more compact Q-ETGs, suggesting that the former are younger than the latter. The idea that new, large, Q-ETGs are responsible for the observed growth in the median size of the population at a given mass is supported by the sizes and number of the star-forming galaxies that are expected to be progenitors of the new Q-ETGs over the same period. In the low mass bin, the new Q-ETG have 30% smaller sizes than their star-forming progenitors. This is likely due to the fading of their disks after they cease star-formation. Comparison with higher z shows that the median size of newly-quenched galaxies roughly scales, at constant mass, as (1+z)^-1. The dominant cause of the size evolution seen in the Q-ETG population is thus that the average sizes of individual Q-ETGs scale with the average density of the Universe at the time when they were quenched, with subsequent size changes in individual objects through eg merging of secondary importance, especially at masses <10^11 Msun.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

The mass assembly of galaxies can proceed through different physical processes. Here we report on the spectroscopic identification of close physical pairs of galaxies at redshifts 2 <z < 4 and discuss the impact of major mergers in shaping galaxies at these early cosmological times. We aim to identify and characterize close physical pairs of galaxies destined to merge and use their properties to infer the contribution of merging processes to the early mass assembly of galaxies. We search for galaxy pairs with a transverse separation r_p < 25 h^{-1} kpc and a velocity difference Delta v < 500 km s^{-1} using the VIMOS VLT Deep Survey (VVDS) and early data from the VIMOS Ultra Deep Survey (VUDS) that comprise a sample of 1111 galaxies with spectroscopic redshifts measurements at redshifts 1.8 < z < 4 in the COSMOS, ECDFS, and VVDS-02h fields. We analyse their spectra and associated visible and near-infrared photometry to assess the main properties of merging galaxies that have an average stellar mass M_{star}=2.3 x 10^{10} M_{sun} at these redshifts. Using the 12 physical pairs found in our sample we obtain a first estimate of the merger fraction at these rdshifts, fmerg ~(15-20)%. The pair separations indicate that these pairs will merge within 1 Gyr, on average, with each producing a more massive galaxy by the time the cosmic star formation peaks at z~1-2. From the average mass ratio between galaxies in the pairs, the stellar mass of the resulting galaxy after merging will be ~60% larger than the most massive galaxy in the pair before merging. We conclude that major merging of galaxy pairs is on-going at 2<z<4 and significantly contributing to the major assembly phase of galaxies at this early epoch.

We investigated the fraction of [OII] emitters in galaxies at z~0.9 as a function of the local galaxy density in the Hubble Space Telescope (HST) COSMOS 2 square degree field. [OII] emitters are selected by the narrow-band excess technique with the NB711-band imaging data taken with Suprime-Cam on the Subaru telescope. We carefully selected 614 photo-z selected galaxies with M_U3500 < -19.31 at z=0.901-0.920, which includes 195 [OII] emitters, to directly compare results with our previous study at z~1.2. We found that the fraction is almost constant at 0.3 Mpc^-2 < Sigma_10th < 10 Mpc^-2. We also checked the fraction of galaxies with blue rest-frame colors of NUV-R < 2 in our photo-z selected sample, and found that the fraction of blue galaxies does not significantly depend on the local density. On the other hand, the semi-analytic model of galaxy formation predicted that the fraction of star-forming galaxies at z~0.9 decreases with increasing the projected galaxy density even if the effects of the projection and the photo-z error in our analysis were taken into account. The fraction of [OII] emitters decreases from ~60% at z~1.2 to ~30% at z~0.9 independent of the galaxy environment. The decrease of the [OII] emitter fraction could be explained mainly by the rapid decrease of the star formation activity in the universe from z~1.2 to z~0.9.

During the course of our Lyman continuum imaging survey, we found that the spectroscopically confirmed Lya emitter LAE 221724+001716 at z = 3.10 in the SSA 22 field shows strong Lyman continuum emission (lambda_rest ~ 900 A) that escapes from this galaxy. However, another recent spectroscopic survey revealed that the supposed Lyman continuum emission could arise from a foreground galaxy at z = 1.76 if the emission line newly detected from the galaxy at lambda_obs ~ 3360 A is Lya. If this is the case, as the angular separation between these two galaxies is very small (~ 0.6"), LAE 221724+001716 at z = 3.10 could be amplified by the gravitational lensing caused by this intervening galaxy. Here we present a possible gravitational lens model for the system of LAE 221724+001716. First, we estimate the stellar mass of the intervening galaxy as Mstar ~ 3.5x10^9 Msun from its UV luminosity and ~ 3.0x10^7--2.4x10^9 Msun through the SED fitting. Then, we find that the gravitational magnification factor ranges from 1.01 to 1.16 using the so-called singular isothermal sphere model for strong lensing. While LAE 221724+001716 is the first system of an LAE-LAE lensing reported so far, the estimated magnification factor is not so significant because the stellar mass of the intervening galaxy is small.

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

In this paper we measure the merger fraction and rate, both minor and major, of massive early-type galaxies (M_star >= 10^11 M_Sun) in the COSMOS field, and study their role in mass and size evolution. We use the 30-band photometric catalogue in COSMOS, complemented with the spectroscopy of the zCOSMOS survey, to define close pairs with a separation 10h^-1 kpc <= r_p <= 30h-1 kpc and a relative velocity Delta v <= 500 km s^-1. We measure both major (stellar mass ratio mu = M_star,2/M_star,1 >= 1/4) and minor (1/10 <= mu < 1/4) merger fractions of massive galaxies, and study their dependence on redshift and on morphology. The merger fraction and rate of massive galaxies evolves as a power-law (1+z)^n, with major mergers increasing with redshift, n_MM = 1.4, and minor mergers showing little evolution, n_mm ~ 0. When split by their morphology, the minor merger fraction for early types is higher by a factor of three than that for spirals, and both are nearly constant with redshift. Our results show that massive early-type galaxies have undergone 0.89 mergers (0.43 major and 0.46 minor) since z ~ 1, leading to a mass growth of ~30%. We find that mu >= 1/10 mergers can explain ~55% of the observed size evolution of these galaxies since z ~ 1. Another ~20% is due to the progenitor bias (younger galaxies are more extended) and we estimate that very minor mergers (mu < 1/10) could contribute with an extra ~20%. The remaining ~5% should come from other processes (e.g., adiabatic expansion or observational effects). This picture also reproduces the mass growth and velocity dispersion evolution of these galaxies. We conclude from these results that merging is the main contributor to the size evolution of massive ETGs at z <= 1, accounting for ~50-75% of that evolution in the last 8 Gyr. Nearly half of the evolution due to mergers is related to minor (mu < 1/4) events.

The M82 `cap' is a gas cloud at a projected radius of 11.6 kpc along the minor axis of this well known superwind source. The cap has been detected in optical line emission and X-ray emission and therefore provides an important probe of the wind energetics. In order to investigate the ionization source of the cap, we observed it with the Kyoto3DII Fabry-Perot instrument mounted on the Subaru Telescope. Deep continuum, Ha, [NII]6583/Ha, and [SII]6716,6731/Ha maps were obtained with sub-arcsecond resolution. The superior spatial resolution compared to earlier studies reveals a number of bright Ha emitting clouds within the cap. The emission line widths (< 100 km s^-1 FWHM) and line ratios in the newly identified knots are most reasonably explained by slow to moderate shocks velocities (v_shock = 40--80 km s^-1) driven by a fast wind into dense clouds. The momentum input from the M82 nuclear starburst region is enough to produce the observed shock. Consequently, earlier claims of photoionization by the central starburst are ruled out because they cannot explain the observed fluxes of the densest knots unless the UV escape fraction is very high (f_esc > 60%), i.e., an order of magnitude higher than observed in dwarf galaxies to date. Using these results, we discuss the evolutionary history of the M82 superwind. Future UV/X-ray surveys are expected to confirm that the temperature of the gas is consistent with our moderate shock model.

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

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

We have exploited the new, deep, near-infrared UltraVISTA imaging of the COSMOS field, in tandem with deep optical and mid-infrared imaging, to conduct a new search for luminous galaxies at redshifts z ~ 7. The unique multi-wavelength dataset provided by VISTA, CFHT, Subaru, HST and Spitzer over a common area of 1 deg^2 has allowed us to select galaxy candidates at z > 6.5 by searching first for Y+J-detected (< 25 AB mag) objects which are undetected in the CFHT+HST optical data. This sample was then refined using a photometric redshift fitting code, enabling the rejection of lower-redshift galaxy contaminants and cool galactic M,L,T dwarf stars.The final result of this process is a small sample of (at most) ten credible galaxy candidates at z > 6.5 which we present in this paper. The first four of these appear to be robust galaxies at z > 6.5, and fitting to their stacked SED yields z = 6.98+-0.05 with a stellar mass M* = 5x10^9 Msun, and rest-frame UV spectral slope beta = -2.0+-0.2. The next three are also good candidates for z > 6.5 galaxies, but the possibility that they are low-redshift galaxies or dwarf stars cannot be excluded. Our final subset of three additional candidates is afflicted not only by potential dwarf-star contamination, but also contains objects likely to lie at redshifts just below z = 6.5. We show that the three even-brighter z > 7 galaxy candidates reported in the COSMOS field by Capak et al. (2011) in fact all lie at z ~ 1.5-3.5. Consequently the new z ~ 7 galaxies reported here are the first credible z ~ 7 Lyman-break galaxies discovered in the COSMOS field and, as the most UV-luminous discovered to date at these redshifts, are prime targets for deep follow-up spectroscopy. We explore their physical properties, and briefly consider the implications of their inferred number density for the form of the galaxy luminosity function at z = 7.

We present the result of our low-luminosity quasar survey in the redshift range of 4.5 < z < 5.5 in the COSMOS field. Using the COSMOS photometric catalog, we selected 15 quasar candidates with 22 < i' < 24 at z~5, that are ~ 3 mag fainter than the SDSS quasars in the same redshift range. We obtained optical spectra for 14 of the 15 candidates using FOCAS on the Subaru Telescope and did not identify any low-luminosity type-1 quasars at z~5 while a low-luminosity type-2 quasar at z~5.07 was discovered. In order to constrain the faint end of the quasar luminosity function at z~5, we calculated the 1sigma confidence upper limits of the space density of type-1 quasars. As a result, the 1sigma confidence upper limits on the quasar space density are Phi< 1.33*10^{-7} Mpc^{-3} mag^{-1} for -24.52 < M_{1450} < -23.52 and Phi< 2.88*10^{-7} Mpc^{-3} mag^{-1} for -23.52 < M_{1450} < -22.52. The inferred 1sigma confidence upper limits of the space density are then used to provide constrains on the faint-end slope and the break absolute magnitude of the quasar luminosity function at z~5. We find that the quasar space density decreases gradually as a function of redshift at low luminosity (M_{1450} ~ -23), being similar to the trend found for quasars with high luminosity (M_{1450}<-26). This result is consistent with the so-called downsizing evolution of quasars seen at lower redshifts.

We present new IRAM PdBI 1.3mm continuum observations at ~1.5" resolution of 28 SMGs previously discovered with the 870um bolometer LABOCA at APEX within the central 0.7deg2 of the COSMOS field. 19 out of the 28 LABOCA sources were detected with the PdBI at a >~3sigma level of ~1.4mJy/b. A combined analysis of this new sample with existing interferometrically identified SMGs in the COSMOS field yields the following results: 1) >~15%, and possibly up to ~40% of single-dish detected SMGs consist of multiple sources, 2) statistical identifications of multi-wavelength counterparts to the single-dish SMGs yield that only ~50% of these single-dish SMGs have real radio or IR counterparts, 3) ~18% of interferometric SMGs have only radio or even no multi-wavelength counterpart at all, and 4) ~50-70% of z>~3 SMGs have no radio counterparts down to an rms of 7-12uJy at 1.4GHz. Using the exact interferometric positions to identify proper multi-wavelength counterparts allows us to determine accurate photometric redshifts for these sources. The redshift distributions of the combined and the individual 1.1mm and 870um selected samples have a higher mean and broader width than the redshift distributions derived in previous studies. Our sample supports the previous tentative trend that on average brighter and/or mm-selected SMGs are located at higher redshifts. There is a tentative offset between the mean redshift for the 1.1mm (<z>=3.1+/-0.4) and 870um (<z>=2.6+/-0.4) selected samples, with the 1.1mm sources lying on average at higher redshifts. Based on our nearly complete sample of AzTEC 1.1mm SMGs within a uniform 0.15deg2 area we infer a higher surface density of z>~4 SMGs than predicted by current cosmological models. In summary, our findings imply that (sub-)millimeter interferometric identifications are crucial to build statistically complete and unbiased samples of SMGs.

[Abridged] We present the results of new near-IR spectroscopic observations of passive galaxies at z>1.4 in a concentration of BzK-selected galaxies in the COSMOS field. The observations have been conducted with Subaru/MOIRCS, and have resulted in absorption lines and/or continuum detection for 18 out of 34 objects. This allows us to measure spectroscopic redshifts for a sample almost complete to K(AB)=21. COSMOS photometric redshifts are found in fair agreement overall with the spectroscopic redshifts, with a standard deviation of ~0.05; however, ~30% of objects have photometric redshifts systematically underestimated by up to ~25%. We show that these systematic offsets in photometric redshifts can be removed by using these objects as a training set. All galaxies fall in four distinct redshift spikes at z=1.43, 1.53, 1.67 and 1.82, with this latter one including 7 galaxies. SED fits to broad-band fluxes indicate stellar masses in the range of ~4-40x10^10Msun and that star formation was quenched ~1 Gyr before the cosmic epoch at which they are observed. The spectra of several individual galaxies have allowed us to measure their Hdelta_F and Dn4000 indices, which confirms their identification as passive galaxies, as does a composite spectrum resulting from the coaddition of 17 individual spectra. The effective radii of the galaxies have been measured on the HST/ACS F814W image, confirming the coexistence at these redshifts of passive galaxies which are substantially more compact than their local counterparts with others that follow the local size-stellar mass relation. For the galaxy with best S/N spectrum we were able to measure a velocity dispersion of 270+/-105 km/s, indicating that this galaxy lies closely on the virial relation given its stellar mass and effective radius.

We present near-infrared spectroscopy of the z=3.2 quasar SDSS J1707+6443, obtained with MOIRCS on the Subaru Telescope. This quasar is classified as a "nitrogen-loud" quasar because of the fairly strong NIII] and NIV] semi-forbidden emission lines from the broad-line region (BLR) observed in its rest-frame UV spectrum. However, our rest-frame optical spectrum from MOIRCS shows strong [OIII] emission from the narrow-line region (NLR) suggesting that, at variance with the BLR, NLR gas is not metal-rich. In order to reconcile these contradictory results, there may be two alternative possibilities; (1) the strong nitrogen lines from the BLR are simply due to a very high relative abundance of nitrogen rather than to a very high BLR metallicity, or (2) the BLR metallicity is not representative of the metallicity of the host galaxy, better traced by the NLR. In either case, the strong broad nitrogen lines in the UV spectrum are not indication of a chemically enriched host galaxy. We estimated the black hole mass and Eddington ratio of this quasar from the velocity width of both CIV and H_beta, that results in log(M_BH/M_sun) = 9.50 and log(L_bol/L_Edd) = -0.34. The relatively high Eddington ratio is consistent with our earlier result that strong nitrogen emission from BLRs is associated with high Eddington ratios. Finally, we detected significant [NeIII] emission from the NLR, implying a quite high gas density of n~10^6 cm^-3 and suggesting a strong coupling between quasar activity and dense interstellar clouds in the host galaxy.

We present our new deep optical imaging and long-slit spectroscopy for Arp 220 that is the archetypical ULIRG in the local universe. Our sensitive Ha imaging has newly revealed large-scale, Ha absorption, i.e., post-starburst regions in this merger; one is found in the eastern superbubble and the other is in the two tidal tails that are clearly reveled in our deep optical imaging. The size of Ha absorption region in the eastern bubble is 5 kpc x 7.5 kpc and the observed Ha equivalent widths are ~2 A +- 0.2 A. The sizes of the northern and southern Ha-absorption tidal tails are ~5 kpc x 10 kpc and ~6 kpc x 20 kpc, respectively. The observed Ha equivalent widths range from 4 A to 7 A. In order to explain the presence of the two post-starburst tails, we suggest a possible multiple-merger scenario for Arp 220 in which two post-starburst disk-like structures merged into one, and then caused the two tails. This favors that Arp 220 is a multiple merging system composed of four or more galaxies, arising from a compact group of galaxies. Taking our new results into account, we discuss a star formation history in the last 1 Gyr in Arp 220.

The success of future large scale weak lensing surveys will critically depend on the accurate estimation of photometric redshifts of very large samples of galaxies. This in turn depends on both the quality of the photometric data and the photo-z estimators. In a previous study, (Bordoloi et al. 2010) we focussed primarily on the impact of photometric quality on photo-z estimates and on the development of novel techniques to construct the N(z) of tomographic bins at the high level of precision required for precision cosmology, as well as the correction of issues such as imprecise corrections for Galactic reddening. We used the same set of templates to generate the simulated photometry as were then used in the photo-z code, thereby removing any effects of "template error". In this work we now include the effects of "template error" by generating simulated photometric data set from actual COSMOS photometry. We use the trick of simulating redder photometry of galaxies at higher redshifts by using a bluer set of passbands on low z galaxies with known redshifts. We find that "template error" is a rather small factor in photo-z performance, at the photometric precision and filter complement expected for all-sky surveys. With only a small sub-set of training galaxies with spectroscopic redshifts, it is in principle possible to construct tomographic redshift bins whose mean redshift is known, from photo-z alone, to the required accuracy of 0.002(1+z).

In order to understand environmental effects on star formation in high-redshift galaxies, we investigate the physical relationships between the star formation activity, stellar mass, and environment for z ~1.2 galaxies in the 2 deg^2 COSMOS field. We estimate star formation using the [OII] emission line and environment from the local galaxy density. Our analysis shows that for massive galaxies M_*>10^10 M_sun, the fraction of [OII] emitters in high-density environments is 1.7 times higher than in low-density environments, while the [OII] emitter fraction does not depend on environment for low-mass M_* < 10^10 M_sun galaxies. In order to understand what drives these trends, we investigate the role of companion galaxies in our sample. We find that the fraction of [OII] emitters in galaxies with companions is 2.4 times as high as that in galaxies without companions at M_* > 10^10 M_sun. In addition, massive galaxies are more likely to have companions in high-density environments. However, although the "number" of star forming galaxies increases for massive galaxies with close companions and in dense environments, the "average" star formation rate of star forming galaxies at a given mass is independent of environment and the presence/absence of a close companion. These results suggest that interactions and/or mergers in high-density environment could induce star formation in massive galaxies at z~1.2, increasing the fraction of star-forming galaxies with M_* > 10^10 M_sun.

We present a rest-frame spectral stacking analysis of ~1000 X-ray sources detected in the XMM-COSMOS field in order to investigate the iron K line properties of active galaxies beyond redshift z~1. In Type I AGN that have a typical X-ray luminosity of Lx~1.5e44 erg/s and z~1.6, the cold Fe K at 6.4 keV is weak (EW~0.05keV), in agreement with the known trend. In contrast, high-ionization lines of Fe XXV and Fe XXVI are pronounced. These high-ionization Fe K lines appear to have a connection with high accretion rates. While no broad Fe emission is detected in the total spectrum, it might be present, albeit at low significance, when the X-ray luminosity is restricted to the range below 3e44 erg/s, or when an intermediate range of Eddington ratio around 0.1 is selected. In Type II AGN, both cold and high-ionzation lines become weak with increasing X-ray luminosity. However, strong high-ionization Fe K (EW~0.3 keV) is detected in the spectrum of objects at z>2, while no 6.4 keV line is found. It is then found that the primary source of the high-ionization Fe K emission is those objects detected with Spitzer-MIPS at 24 micron. Given their median redshift of z=2.5, their bolometric luminosity is likely to reach 10^13 Lsun and the MIPS-detected emission most likely originates from hot dust heated by embedded AGN, probably accreting at high Eddington ratio. These properties match those of rapidly growing black holes in ultra-luminous infrared galaxies at the interesting epoch (z=2-3) of galaxy formation.

We explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS survey and the physical properties of their host galaxies. We study a bolometric regime (<Lbol > 8 x 10^45 erg/s) where several theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. We confirm that obscured quasars mainly reside in massive galaxies (Mstar>10^10 Msun) and that the fraction of galaxies hosting such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical properties of the hosts. We find that at z ~1, ~62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ~71% at z ~2, and 100% at z ~3. We also find that the the evolution from z ~1 to z ~3 of the specific star-formation rate of the Type-2 QSO hosts is in excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated galaxies tend to host Type-2 QSOs with low Eddington ratios (lambda<0.1), while disk-dominated or merging galaxies have at their centers BHs accreting at high Eddington ratios (lambda > 0.1).

With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

We present a deep optical spectrum of TN J0924-2201, the most distant radio galaxy at z = 5.19, obtained with FOCAS on the Subaru Telescope. We successfully detect, for the first time, the CIV1549 emission line from the narrow-line region (NLR). In addition to the emission-line fluxes of Ly alpha and CIV, we set upper limits on the NV and HeII emissions. We use these line detections and upper limits to constrain the chemical properties of TN J0924-2201. By comparing the observed emission-line flux ratios with photoionization models, we infer that the carbon-to-oxygen relative abundance is already [C/O] > -0.5 at a cosmic age of ~ 1.1 Gyr. This lower limit on [C/O] is higher than the ratio expected at the earliest phases of the galaxy chemical evolution, indicating that TN J0924-2201 has already experienced significant chemical evolution at z = 5.19.

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

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

We present Subaru/FOCAS spectropolarimetry of two active galaxies in the Cosmic Evolution Survey. These objects were selected to be optically dull, with the bright X-ray emission of an AGN but missing optical emission lines in our previous spectroscopy. Our new observations show that one target has very weak emission lines consistent with an optically dull AGN, while the other object has strong emission lines typical of a host-diluted Type 2 Seyfert galaxy. In neither source do we observe polarized emission lines, with 3-sigma upper limits of P_BLR < 2%. This means that the missing broad emission lines (and weaker narrow emission lines) are not due to simple anisotropic obscuration, e.g., by the canonical AGN torus. The weak-lined optically dull AGN exhibits a blue polarized continuum with P = 0.78 +/- 0.07% at 4400 A < lambda_rest < 7200 A (P = 1.37 +/- 0.16% at 4400 A < lambda_rest < 5050 A). The wavelength dependence of this polarized flux is similar to that of an unobscured AGN continuum and represents the intrinsic AGN emission, either as synchrotron emission or the outer part of an accretion disk reflected by a clumpy dust scatterer. Because this intrinsic AGN emission lacks emission lines, this source is likely to have a radiatively inefficient accretion flow.

We have identified a complete, flux-limited, (S_160>120 mJy), sample of 160 micron-selected sources from Spitzer observations of the 1-deg^2 ISO Deep Field region in the Lockman Hole. Ground-based UV, optical and near-infrared (NIR) photometry and optical spectroscopy have been used to determine colors, redshifts and masses for the complete sample of 40 galaxies. Spitzer-IRAC+MIPS photometry, supplemented by ISOPHOT data at 90 micron and 170 micron, has been used to calculate accurate total infrared luminosities, LIR(8-1000 micron), and to determine the IR luminosity function (LF) of luminous infrared galaxies (LIRGs). The maximum observed redshift is z~0.80 and the maximum total infrared luminosity is log(L_IR/L_Sun)=12.74. Over the luminosity range log(L_IR/L_Sun)=10-12, the LF for LIRGs in the Lockman Hole Deep Field is similar to that found previously for local sources at similar infrared luminosities. The mean host galaxy mass, log(M/M_Sun)=10.7, and dominance of HII-region spectral types, is also similar to what has been found for local LIRGs, suggesting that intense starbursts likely power the bulk of the infrared luminosity for sources in this range of LIR. However for the most luminous sources, log(L_IR/L_Sun)>12.0, we find evidence for strong evolution in the LF \propto (1+z)6\pm1, assuming pure number density evolution. These ultraluminous infrared galaxies (ULIRGs) have a larger mean host mass, log(M/M_Sun)=11.0, and exhibit disturbed morphologies consistent with strong-interactions/mergers, and they are also more likely to be characterized by starburst-AGN composite or AGN spectral types.

We searched for quasars that are ~ 3 mag fainter than the SDSS quasars in the redshift range 3.7 < z < 4.7 in the COSMOS field to constrain the faint end of the quasar luminosity function. Using optical photometric data, we selected 31 quasar candidates with 22 < i' < 24 at z ~ 4. We obtained optical spectra for most of these candidates using FOCAS on the Subaru telescope, and identified 8 low-luminosity quasars at z ~ 4. In order to derive the quasar luminosity function (QLF) based on our spectroscopic follow-up campaign, we estimated the photometric completeness of our quasar survey through detailed Monte Carlo simulations. Our QLF at z ~ 4 has a much shallower faint-end slope beta = -1.67^{+0.11}_{-0.17} than that obtained by other recent surveys in the same redshift. Our result is consistent with the scenario of downsizing evolution of active galactic nuclei inferred by recent optical and X-ray quasar surveys at lower redshifts.

We report the discovery of a significant excess of candidate Halpha emitters (HAEs) in the field of the radio galaxy 4C 23.56 at z=2.483. Using the MOIRCS near-infrared imager on the Subaru Telescope we found 11 candidate emission-line galaxies to a flux limit of ~7.5 10^-17 erg s-1 cm-2, which is about 5 times excess from the expected field counts with ~3-sigma significance. Three of these are spectroscopically confirmed as redshifted Halpha at z=2.49. The distribution of candidate emitters on the sky is tightly confined to a 1.2-Mpc-radius area at z=2.49, locating 4C 23.56 at the western edge of the distribution. Analysis of the deep Spitzer MIPS 24 mu m imaging shows that there is also an excess of faint MIPS sources. All but two of the 11 HAEs are also found in the MIPS data. The inferred star-formation rate (SFR) of the HAEs based on the extinction-corrected Halpha luminosity (median SFR >~100 M_solar yr-1) is similar to those of HAEs in random fields at z~2. On the other hand, the MIPS-based SFR for the HAEs is on average 3.6 times larger, suggesting the existence of the star-formation significanly obscured by dust. The comparison of the Halpha-based star-formation activities of the HAEs in the 4C 23.56 field to those in another proto-cluster around PKS 1138-262 at z=2.16 reveals that the latter tend to have fainter Halpha emission despite similar K-band magnitudes. This suggests that star-formation may be suppressed in the PKS 1138-262 protocluster relative to the 4C 23.56 protocluster. This difference among the HAEs in the two proto-clusters at z > 2 may imply that some massive cluster galaxies are just forming at these epochs with some variation among clusters.

Active galactic nuclei (AGNs) are characterized by a clear correlation between luminosity and metallicity (L_AGN-Z_AGN relation). The origin of this correlation is not clear. It may result from a relation between the black hole mass (M_BH) and metallicity, or from a relation between the accretion rate (L/L_Edd) and metallicity. To investigate the origin of the L_AGN-Z_AGN relation, we use optical spectra of 2383 quasars at 2.3 < z < 3.0 from the Sloan Digital Sky Survey. By using this data set, we have constructed composite spectra of 33 subsamples in intervals of both M_BH and L/L_Edd. From these composite spectra we measure emission-line flux ratios that are sensitive to the metallicity of the broad line region (BLR); specifically, NV/CIV, NV/HeII, (SiIV+OIV])/CIV, and AlIII/CIV. We find that there is a significant correlation between M_BH and Z_BLR as inferred from all four metallicity-sensitive emission-line flux ratios. This result strongly suggests that the observed L_AGN-Z_AGN relation is mostly a consequence of the M_BH-Z_AGN relation. The relation between M_BH and Z_BLR is likely a consequence of both the M_BH-M_bul relation and of the mass-metallicity relation in the host galaxy. We also find that L/L_Edd correlates with the emission line flux ratios involving NV (more specifically, NV/CIV and NV/HeII), while it does not correlate with the other two metallicity sensitive emission line flux ratios, i.e., (SiIV+OIV])/CIV and AlIII/CIV. These correlations indicate that the emission-line flux ratios involving NV depend on both metallicity and relative abundance of nitrogen. We suggest that the relation between L/L_Edd and those line ratios involving nitrogen, is caused by a delay of the black hole accretion rate relative to the onset of nuclear star formation of about 10^8 years, which is the timescale required for the nitrogen enrichment.

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)<11.7 with high-resolution HST/ACS imaging from the COSMOS field. The visual analysis of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted morphologies in the AGN and control samples, i.e. quantifying the signature of recent mergers which might potentially be responsible for fueling/triggering the AGN. We find that (1) the vast majority (>85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

Based on the recent very deep near-infrared imaging of the Hubble Ultra Deep Field with WFC3 on the Hubble Space Telescope, five groups published most probable samples of galaxies at z~8, selected by the so-called dropout method or photometric redshift; e.g., Y_105-dropouts (Y_105-J_125 > 0.8). These studies are highly useful for investigating both the early star formation history of galaxies and the sources of cosmic re-ionization. In order to better understand these issues, we carefully examine if there are low-$z$ interlopers in the samples of z~8 galaxy candidates. We focus on the strong emission-line galaxies at z~2 in this paper. Such galaxies may be selected as Y_105-dropouts since the [OIII] lambda 5007 emission line is redshifted into the J_125-band. We have found that the contamination from such low-$z$ interlopers is negligibly small. Therefore, all objects found by the five groups are free from this type of contamination. However, it remains difficult to extract real z~8 galaxies because all the sources are very faint and the different groups have found different candidates. With this in mind, we construct a robust sample of eight galaxies at z~8 from the objects found by the five groups: each of these eight objects has been selected by at least two groups. Using this sample, we discuss their UV continuum slope. We also discuss the escape fraction of ionizing photons adopting various metallicities. Our analysis suggests that massive stars forming in low-metallicity gas (Z~5 \times 10^-4 Z_sun) can be responsible for the completion of cosmic re-ionization if the escape fraction of ionizing continuum from galaxies is as large as 0.5, and this is consistent with the observed blue UV continua.

We study the impact of the environment on the evolution of galaxies in the zCOSMOS 10k sample in the redshift range 0.1<z<1.0 over an area of ~1.5 deg2. The considered sample of secure spectroscopic redshifts contains about 8500 galaxies, with their stellar masses estimated by SED fitting of the multiwavelength optical to NIR photometry. The evolution of the galaxy stellar mass function (GSMF) in high and low density regions provides a tool to study the mass assembly evolution in different environments; moreover, the contributions to the GSMF from different galaxy types, as defined by their SEDs and their morphologies, can be quantified. At redshift z~1, the GSMF is only slightly dependent on environment, but at lower redshifts the shapes of the GSMFs in high- and low-density environments become extremely different, with high density regions exhibiting a marked bimodality. As a result, we infer that galaxy evolution depends on both the stellar mass and the environment, the latter setting the probability of a galaxy to have a given mass: all the galaxy properties related to the stellar mass show a dependence on environment, reflecting the difference observed in the mass functions. The shapes of the GSMFs of early- and late-type galaxies are almost identical for the extremes of the density contrast we consider. The evolution toward z=0 of the mass at which the early- and late-type GSMFs match each other is more rapid in high density environments. The comparison of the observed GSMFs to the same quantities derived from a set of mock catalogues shows that blue galaxies in sparse environments are overproduced in the semi-analytical models at intermediate and high masses, because of a deficit of star formation suppression, while at z<0.5 an excess of red galaxies is present in dense environments at intermediate and low masses, because of the overquenching of satellites. ABRIDGED

We present deep optical 18-medium-band photometry from the Subaru telescope over the ~30' x 30' Extended Chandra Deep Field-South (ECDF-S), as part of the Multiwavelength Survey by Yale-Chile (MUSYC). This field has a wealth of ground- and space-based ancillary data, and contains the GOODS-South field and the Hubble Ultra Deep Field. We combine the Subaru imaging with existing UBVRIzJHK and Spitzer IRAC images to create a uniform catalog. Detecting sources in the MUSYC BVR image we find ~40,000 galaxies with R_AB<25.3, the median 5 sigma limit of the 18 medium bands. Photometric redshifts are determined using the EAZY code and compared to ~2000 spectroscopic redshifts in this field. The medium band filters provide very accurate redshifts for the (bright) subset of galaxies with spectroscopic redshifts, particularly at 0.1 < z < 1.2 and at z > 3.5. For 0.1 < z < 1.2, we find a 1 sigma scatter in \Delta z/(1+z) of 0.007, similar to results obtained with a similar filter set in the COSMOS field. As a demonstration of the data quality, we show that the red sequence and blue cloud can be cleanly identified in rest-frame color-magnitude diagrams at 0.1 < z < 1.2. We find that ~20% of the red-sequence-galaxies show evidence of dust-emission at longer rest-frame wavelengths. The reduced images, photometric catalog, and photometric redshifts are provided through the public MUSYC website.

We present Spitzer 7-38um spectra for a 24um flux limited sample of galaxies at z~0.7 in the COSMOS field. The detailed high-quality spectra allow us to cleanly separate star formation (SF) and active galactic nucleus (AGN) in individual galaxies. We first decompose mid-infrared Luminosity Functions (LFs). We find that the SF 8um and 15um LFs are well described by Schechter functions. AGNs dominate the space density at high luminosities, which leads to the shallow bright-end slope of the overall mid-infrared LFs. The total infrared (8-1000um) LF from 70um selected galaxies shows a shallower bright-end slope than the bolometrically corrected SF 15um LF, owing to the intrinsic dispersion in the mid-to-far-infrared spectral energy distributions. We then study the contemporary growth of galaxies and their supermassive black holes (BHs). Seven of the 31 Luminous Infrared Galaxies with Spitzer spectra host luminous AGNs, implying an AGN duty cycle of 23+/-9%. The time-averaged ratio of BH accretion rate and SF rate matches the local M_BH-M_bulge relation and the M_BH-M_host relation at z ~ 1. These results favor co-evolution scenarios in which BH growth and intense SF happen in the same event but the former spans a shorter lifetime than the latter. Finally, we compare our mid-infrared spectroscopic selection with other AGN identification methods and discuss candidate Compton-thick AGNs in the sample. While only half of the mid-infrared spectroscopically selected AGNs are detected in X-ray, ~90% of them can be identified with their near-infrared spectral indices.

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts).

We analyze the morphological properties of a large sample of 1503 70 micron selected galaxies in the COSMOS field spanning the redshift range 0.01<z< 3.5 with a median redshift of 0.5 and an infrared luminosity range of 10^8<L_IR<10^14L_sun with a median luminosity of 10^11.4 L_sun. In general these galaxies are massive, with a stellar mass range of 10^10-10^12 M_sun, and luminous, with -25<M_K<-20. We find a strong correlation between the fraction of major mergers and L_IR, with the fraction at the highest luminosity being up to 50%. We also find that the fraction of spirals drops dramatically with L_IR. Minor mergers likely play a role in boosting the infrared luminosity for sources with low luminosities. The precise fraction of mergers in any given L_IR bin varies by redshift due to sources at z>1 being difficult to classify and subject to the effects of band pass shifting, therefore, these numbers can only be considered lower limits. At z<1, where the morphological classifications are most robust, major mergers clearly dominate the ULIRG population (50-80%) and are important for the LIRG population (25-40%). At z>1 the fraction of major mergers is at least 30-40% for ULIRGs. Although the general morphological trends agree with what has been observed for local (U)LIRGs, the fraction of major mergers is slightly lower than seen locally. This is in part due to the difficulty of identifying merger signatures at high redshift. We argue that given the number of major gas-rich mergers observed and the relatively short timescale that they would be observable in the (U)LIRG phase that it is plausible for the observed red sequence of massive ellipticals (<10^12 M_sun) to have been formed entirely by gas-rich major mergers.

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

We present observations of a very massive galaxy at z=1.82 which show that its morphology, size, velocity dispersion and stellar population properties that are fully consistent with those expected for passively evolving progenitors of today's giant ellipticals. These findings are based on a deep optical rest-frame spectrum obtained with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru telescope of a high-z passive galaxy candidate (pBzK) from the COSMOS field, for which we accurately measure its redshift of z=1.8230 and obtain an upper limit on its velocity dispersion sigma_star<326 km/s. By detailed stellar population modeling of both the galaxy broad-band SED and the rest-frame optical spectrum we derive a star-formation-weighted age and formation redshift of t_sf~1-2 Gyr and z_form~2.5-4, and a stellar mass of M_star~(3-4)x10^{11} M_sun. This is in agreement with a virial mass limit of M_vir<7x10^{11}M_sun, derived from the measured sigma_star value and stellar half-light radius, as well as with the dynamical mass limit based on the Jeans equations. In contrast with previously reported super-dense passive galaxies at z~2, the present galaxy at z=1.82 appears to have both size and velocity dispersion similar to early-type galaxies in the local Universe with similar stellar mass. This suggests that z~2 massive and passive galaxies may exhibit a wide range of properties, then possibly following quite different evolutionary histories from z~2 to z=0.

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust.

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.

Using Chandra observations we have identified a sample of seven off-nuclear X-ray sources, in the redshift range z=0.072-0.283, located within optically bright galaxies in the COSMOS Survey. Using the multi-wavelength coverage available in the COSMOS field, we study the properties of the host galaxies of these ULXs. In detail, we derived their star formation rate from H_alpha measurements and their stellar masses using SED fitting techniques with the aim to compute the probability to have an off-nuclear source based on the host galaxy properties. We divide the host galaxies in different morphological classes using the available ACS/HST imaging. We find that our ULXs candidates are located in regions of the SFR versus M$_star$ plane where one or more off-nuclear detectable sources are expected. From a morphological analysis of the ACS imaging and the use of rest-frame colours, we find that our ULXs are hosted both in late and early type galaxies. Finally, we find that the fraction of galaxies hosting a ULX ranges from ~0.5% to ~0.2% going from L[0.5-2 keV]=3 x 10^39 erg s^-1 to L[0.5-2 keV]= 2 x 10^40 erg s^-1.

For the initial fields of the density contrast and peculiar velocity, we theoretically calculate the differential and integral length scales, i.e., statistical measures that respectively characterize the small- and large-scale fluctuations of a random field. These length scales and the associated mass scales explain the length and mass scales observed for (1) halos of young galaxies at z > 5, (2) halos of galaxies at z = 0, and (3) the largest structures in the galaxy distribution at z = 0. We thereby discuss that such observed scales are fossil imprints of the characteristic scales of the initial fields.

The identities of the main processes triggering and quenching star-formation in galaxies remain unclear. A key stage in evolution, however, appears to be represented by post-starburst galaxies. To investigate their impact on galaxy evolution, we initiated a multiwavelength study of galaxies with k+a spectral features in the COSMOS field. We examine a mass-selected sample of k+a galaxies at z=0.48-1.2 using the spectroscopic zCOSMOS sample. K+a galaxies occupy the brightest tail of the luminosity distribution. They are as massive as quiescent galaxies and populate the green valley in the colour versus luminosity (or stellar mass) distribution. A small percentage (<8%) of these galaxies have radio and/or X-ray counterparts (implying an upper limit to the SFR of ~8Msun/yr). Over the entire redshift range explored, the class of k+a galaxies is morphologically a heterogeneous population with a similar incidence of bulge-dominated and disky galaxies. This distribution does not vary with the strength of the Hdelta absorption line but instead with stellar mass in a way reminiscent of the well-known mass-morphology relation. Although k+a galaxies are also found in underdense regions, they appear to reside typically in a similarly rich environment as quiescent galaxies on a physical scale of ~2-8Mpc, and in groups they show a morphological early-to-late type ratio similar to the quiescent galaxy class. With the current data set, we do not find evidence of statistical significant evolution in either the number/mass density of k+a galaxies at intermediate redshift with respect to the local values, or the spectral properties. Those galaxies, which are affected by a sudden quenching of their star-formation activity, may increase the stellar mass of the red-sequence by up to a non-negligible level of ~10%.

(abridged) We present a new near-infrared survey covering the 2 deg sq COSMOS field. Combining our survey with Subaru B and z images we construct a deep, wide-field optical-infrared catalogue. At Ks<23 (AB magnitudes) our survey completeness is greater than 90% and 70% for stars and galaxies respectively and contains 143,466 galaxies and 13,254 stars. At z~2 our catalogues contain 3931 quiescent and 25,757 star-forming BzK-selected galaxies representing the largest and most secure sample of these objects to date. Our counts of quiescent galaxies turns over at Ks~22 an effect which we demonstrate cannot be due to sample incompleteness. In our survey both the number of faint and bright quiescent objects exceeds the predictions of a semi-analytic galaxy formation model, indicating potentially the need for further refinements in the amount of merging and AGN feedback at z~2 in these models. We measure the angular correlation function for each sample and find that at small scales the correlation function for passive BzK galaxies exceeds the clustering of dark matter. We use 30-band photometric redshifts to derive the spatial correlation length and the redshift distributions for each object class. At Ks<22 we find r_0^{\gamma/1.8}=7.0 +/-0.5h^{-1} Mpc for the passive BzK candidates and 4.7+/-0.8h^{-1} Mpc for the star-forming BzK galaxies. Our pBzK galaxies have an average photometric redshift of z_p~1.4, in approximate agreement with the limited spectroscopic information currently available. The stacked Ks image will be made publicly available from IRSA.

We follow the galaxy stellar mass assembly by morphological and spectral type in the COSMOS 2-deg^2 field. We derive the stellar mass functions and stellar mass densities from z=2 to z=0.2 using 196,000 galaxies selected at F(3.6 micron) > 1 microJy with accurate photometric redshifts (sigma_((zp-zs)/(1+zs))=0.008 at i<22.5). Using a spectral classification, we find that z~1 is an epoch of transition in the stellar mass assembly of quiescent galaxies. Their stellar mass density increases by 1.1 dex between z=1.5-2 and z=0.8-1 (Delta t ~2.5 Gyr), but only by 0.3 dex between z=0.8-1 and z~0.1 (Delta t ~ 6 Gyr). Then, we add the morphological information and find that 80-90% of the massive quiescent galaxies (log(M)~11) have an elliptical morphology at z<0.8. Therefore, a dominant mechanism links the shutdown of star formation and the acquisition of an elliptical morphology in massive galaxies. Still, a significant fraction of quiescent galaxies present a Spi/Irr morphology at low mass (40-60% at log(M)~9.5), but this fraction is smaller than predicted by semi-analytical models using a ``halo quenching'' recipe. We also analyze the evolution of star-forming galaxies and split them into ``intermediate activity'' and ``high activity'' galaxies. We find that the most massive ``high activity'' galaxies end their high star formation rate phase first. Finally, the space density of massive star-forming galaxies becomes lower than the space density of massive elliptical galaxies at z<1. As a consequence, the rate of ``wet mergers'' involved in the formation of the most massive ellipticals must decline very rapidly at z<1, which could explain the observed slow down in the assembly of these quiescent and massive sources.

We present a large robust sample of 1503 reliable and unconfused 70microm selected sources from the multiwavelength data set of the Cosmic Evolution Survey (COSMOS). Using the Spitzer IRAC and MIPS photometry, we estimate the total infrared luminosity, L_IR (8--1000 microns), by finding the best fit template from several different template libraries. The long wavelength 70 and 160 micron data allow us to obtain a reliable estimate of L_IR, accurate to within 0.2 and 0.05 dex, respectively. The 70 micron data point enables a significant improvement over the luminosity estimates possible with only a 24 micron detection. The full sample spans a wide range in L_IR, L_IR ~ 10^8-10^14 L_sun, with a median luminosity of 10^11.4 L_sun. We identify a total of 687 luminous, 303 ultraluminous, and 31 hyperluminous infrared galaxies (LIRGs, ULIRGs, and HyLIRGs) over the redshift range 0.01<z<3.5 with a median redshift of 0.5. Presented here are the full spectral energy distributions for each of the sources compiled from the extensive multiwavelength data set from the ultraviolet (UV) to the far-infrared (FIR). Using SED fits we find possible evidence for a subset of cooler ultraluminous objects than observed locally. However, until direct observations at longer wavelengths are obtained, the peak of emission and the dust temperature cannot be well constrained. We use these SEDs, along with the deep radio and X-ray coverage of the field, to identify a large sample of candidate active galactic nuclei (AGN). We find that the fraction of AGN increases strongly with L_IR, as it does in the local universe, and that nearly 70% of ULIRGs and all HyLIRGs likely host a powerful AGN.

We constrain the ratio of black hole (BH) mass to total stellar mass of type-1 AGN in the COSMOS survey at 1<z<2. For 10 AGN at mean redshift z~1.4 with both HST/ACS and HST/NICMOS imaging data we are able to compute total stellar mass M_(*,total), based on restframe UV-to-optical host galaxy colors which constrain mass-to-light ratios. All objects have virial BH mass-estimates available from the COSMOS Magellan/IMACS and zCOSMOS surveys. We find zero difference between the M_BH--M_(*,total)-relation at z~1.4 and the M_BH--M_(*,bulge)-relation in the local Universe. Our interpretation is: (a) If our objects were purely bulge-dominated, the M_BH--M_(*,bulge)-relation has not evolved since z~1.4. However, (b) since we have evidence for substantial disk components, the bulges of massive galaxies (logM_(*,total)=11.1+-0.25 or logM_BH~8.3+-0.2) must have grown over the last 9 Gyrs predominantly by redistribution of disk- into bulge-mass. Since all necessary stellar mass exists in the galaxy at z=1.4, no star-formation or addition of external stellar material is required, only a redistribution e.g. induced by minor and major merging or through disk instabilities. Merging, in addition to redistributing mass in the galaxy, will add both BH and stellar/bulge mass, but does not change the overall final M_BH/M_(*,bulge) ratio. Since the overall cosmic stellar and BH mass buildup trace each other tightly over time, our scenario of bulge-formation in massive galaxies is independent of any strong BH-feedback and means that the mechanism coupling BH and bulge mass until the present is very indirect.

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate density determination. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of the 8.0-4.5 $\mu$m infrared color--specific star formation plane, from which we extract the corresponding control samples. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive AGN increases from ~2% in underdense regions to ~15% for overdensities (1+delta) greater than 10. This trend is also present as a function of richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, while in other cases it is determined by the structure (group or cluster) in which the galaxy resides.

We have used the zCOSMOS-bright 10k sample to identify 3244 Spitzer/MIPS 24-micron-selected galaxies with 0.06< S(24um)< 0.50 mJy and I(AB)<22.5, over 1.5 deg^2 of the COSMOS field, and studied different spectral properties, depending on redshift. At 0.2<z<0.3, we found that different reddening laws of common use in the literature explain the dust extinction properties of around 80% of our infrared (IR) sources, within the error bars. For up to 16% of objects, instead, the Halpha/Hbeta ratios are too high for their IR/UV attenuations, which is probably a consequence of inhomogenous dust distributions. In only a few of our galaxies at 0.2<z<0.3 the IR emission could be mainly produced by dust heated by old rather than young stars. Besides, the line ratios of ~22% of our galaxies suggest that they might be star-formation/nuclear-activity composite systems. At 0.5<z<0.7, we estimated galaxy metallicities for 301 galaxies: at least 12% of them are securely below the upper-branch mass-metallicity trend, which is consistent with the local relation. Finally, we performed a combined analysis of the Hdelta equivalent-width versus Dn(4000) diagram for 1722 faint and bright 24um galaxies at 0.6<z<1.0, spanning two decades in mid-IR luminosity. We found that, while secondary bursts of star formation are necessary to explain the position of the most luminous IR galaxies in that diagram, quiescent, exponentially-declining star formation histories can well reproduce the spectral properties of ~40% of the less luminous sources. Our results suggest a transition in the possible modes of star formation at total IR luminosities L(TIR)=(3 +/-2)x10^11 Lsun.

(Abriged) We report on the measurement of the rest frame K-band luminosity and total stellar mass of the hosts of 89 broad line Active Galactic Nuclei detected in the zCOSMOS survey in the redshift range 1<z<2.2. The unprecedented multiwavelength coverage of the survey field allows us to disentangle the emission of the host galaxy from that of the nuclear black hole in their Spectral Energy Distributions. We derive an estimate of black hole masses through the analysis of the broad Mg II emission lines observed in the medium-resolution spectra taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared to the local value, the average black hole to host galaxy mass ratio appears to evolve positively with redshift, with a best fit evolution of the form (1+z)^{0.68 \pm0.12 +0.6 -0.3}, where the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, in the calibration of the virial relation used to estimate BH masses and in the mean QSO SED adopted. A thorough analysis of observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion that an evolution of the MBH-M* relation must ensue for actively growing black holes at early times: either its overall normalization, or its intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as signature of either a more rapid growth of supermassive black holes at high redshift, a change of structural properties of AGN hosts at earlier times, or a significant mismatch between the typical growth times of nuclear black holes and host galaxies.

We present infrared, optical, and X-ray data of 48 X-ray bright, optically dull AGNs in the COSMOS field. These objects exhibit the X-ray luminosity of an active galactic nucleus (AGN) but lack broad and narrow emission lines in their optical spectrum. We show that despite the lack of optical emission lines, most of these optically dull AGNs are not well-described by a typical passive red galaxy spectrum: instead they exhibit weak but significant blue emission like an unobscured AGN. Photometric observations over several years additionally show significant variability in the blue emission of four optically dull AGNs. The nature of the blue and infrared emission suggest that the optically inactive appearance of these AGNs cannot be caused by obscuration intrinsic to the AGNs. Instead, up to ~70% of optically dull AGNs are diluted by their hosts, with bright or simply edge-on hosts lying preferentially within the spectroscopic aperture. The remaining ~30% of optically dull AGNs have anomalously high f_x/f_o ratios and are intrinsically weak, not obscured, in the optical. These optically dull AGNs are best described as a weakly accreting AGN with a truncated accretion disk from a radiatively inefficient accretion flow.

(Abridged) We studied the evolution in the B band luminosity function to z~1 in the zCOSMOS 10k sample, for which both accurate galaxy classifications and a detailed description of the local density field are available. The global LF exhibits a brightening of ~0.7 mag in M* from z~0.2 to z~0.9. At low z, late types dominate at faint magnitudes, while the bright end is populated mainly by early types. At higher z, late-type galaxies evolve significantly and, at z~1, the contribution from the various types to the bright end of the LF is comparable. The evolution for early types is in both luminosity and normalization. A similar behaviour is exhibited by late types, but with an opposite trend for the normalization. Studying the role of the environment, we find that the global LF of galaxies in overdense regions has always a brighter M* and a flatter slope. In low density environments, the main contribution to the LF is from blue galaxies, while for high density environments there is an important contribution from red galaxies to the bright end. The differences between the global LF in the two environments are not due to only a difference in the relative numbers of red and blue galaxies, but also to their relative luminosity distributions: the value of M* for both types in underdense regions is always fainter than in overdense environments. The "specular" evolution of late- and early-type galaxies is consistent with a scenario where a part of blue galaxies is transformed in red galaxies with increasing cosmic time, without significant changes in the fraction of intermediate-type galaxies. The bulk of this tranformation in overdense regions probably happened before z~1, while it is still ongoing at lower z in underdense environments.

We present the discovery of a candidate of giant radio-quiet Lyman-alpha (Lya) blob (RQLAB) in a large-scale structure around a high-redshift radio galaxy (HzRG) lying in a giant Lya halo, B3 J2330+3927 at redshift z=3.087. We obtained Lya imaging around B3 J2330+3927 with Subaru/Suprime-Cam to search for Lya emitters (LAEs) and absorbers (LAAs) at redshift z=3.09+-0.03. We detected candidate 127 LAEs and 26 LAAs in the field of view of 31' x 24'. We found that B3 J2330+3927 is surrounded by a 130 kpc Lya halo and a large-scale (60 x 20 comoving Mpc) filamentary structure. The large-scale structure contains one prominent local density peak with an overdensity of greater than 5, which is 8' (15 comoving Mpc) away from B3 J2330+3927. In this peak, we discovered a candidate 100 kpc RQLAB. The existence of both types of Lya nebulae in the same large-scale structure suggests that giant Lya nebulae need special large-scale environments to form. On smaller scales, however, the location of B3 J2330+3927 is not a significant local density peak in this structure, in contrast to the RQLAB. There are two possible interpretations of the difference of the local environments of these two Lya nebulae. Firstly, RQLAB may need a prominent (delta ~ 5) density peak of galaxies to form through intense star-bursts due to frequent galaxy interactions/mergers and/or continuous gas accretion in an overdense environment. On the other hand, Lya halo around HzRG may not always need a prominent density peak to form if the surrounding Lya halo is mainly powered by its radio and AGN activities. Alternatively, both RQLAB and Lya halo around HzRG may need prominent density peaks to form but we could not completely trace the density of galaxies because we missed evolved and dusty galaxies in this survey.

We investigate if the discrepancy between estimates of the total baryon mass fraction obtained from observations of the cosmic microwave background (CMB) and of galaxy groups/clusters persists when a large sample of groups is considered. To this purpose, 91 candidate X-ray groups/poor clusters at redshift 0.1 < z < 1 are selected from the COSMOS 2 deg^2 survey, based only on their X-ray luminosity and extent. This sample is complemented by 27 nearby clusters with a robust, analogous determination of the total and stellar mass inside R_500. The total sample of 118 groups and clusters with z < 1 spans a range in M_500 of ~10^13--10^15 M_sun. We find that the stellar mass fraction associated with galaxies at R_500 decreases with increasing total mass as (M_500)^-0.37 \pm 0.04, independent of redshift. Estimating the total gas mass fraction from a recently derived, high quality scaling relation, the total baryon mass fraction (f_500^stars+gas=f_500^stars+f_500^gas) is found to increase by ~ 25% when M_500 increases from <M>=5 X 10^13 M_sun to <M> = 7 X 10^14 M_sun. After consideration of a plausible contribution due to intra--cluster light (11--22% of the total stellar mass), and gas depletion through the hierarchical assembly process (10% of the gas mass), the estimated values of the total baryon mass fraction are still lower than the latest CMB measure of the same quantity (WMAP5), at a significance level of 3.3\sigma for groups of <M>=5 X 10^13~M_sun. The discrepancy decreases towards higher total masses, such that it is 1\sigma at <M>= 7 X 10^14~M_sun. We discuss this result in terms of non--gravitational processes such as feedback and filamentary heating.

We present detailed morphological properties of Lyman alpha emitters (LAEs) at z~ 5.7 in the COSMOS field, based on {\it Hubble Space Telescope} Advanced Camera for Surveys (ACS) data. The ACS imaging in the F814W filter covered 85 LAEs of the 119 LAEs identified in the full two square degree field, and 47 LAEs of them are detected in the ACS images. Nearly half of them are spatially extended with a size larger than 0.15 arcsec (~0.88 kpc at z=5.7) up to 0.4 arcsec (~2.5 kpc at z=5.7). The others are nearly unresolved compact objects. Two LAEs show double-component structures, indicating interaction or merging of building components to form more massive galaxies. By stacking the ACS images of all the detected sources, we obtain a Sersic parameter of n~0.7 with a half-light radius of 0.13 arcsec (0.76 kpc), suggesting that the majority of ACS detected LAEs have not spheroidal-like but disk-like or irregular light profiles. Comparing ACS F814W magnitudes (I_814) with Subaru/Suprime-Cam magnitudes in the NB816, i', and z' bands, we find that the ACS imaging in the F814W band mainly probes UV continuum rather than Lyman alpha line emission. UV continuum sizes tend to be larger for LAEs with larger Ly\alpha emission regions as traced by the NB816 imaging. The non-detection of 38 LAEs in the ACS images is likely due to the fact that their surface brightness is even too low both in the UV continuum and Lyalpha emission. Estimating I_814 for the LAEs with ACS non-detection from the z' and NB816 magnitudes, we find that 16 of these are probably LAEs with a size larger than 0.15 arcsec in UV continuum. All these results suggest that our LAE sample contains systematically larger LAEs in UV continuum size than those previously studied at z~6.

We investigate the relation between the star-formation activity in galaxies and environment at z=1.2 in the COSMOS field, using the fraction of [OII] emitters and the local galaxy density. The fraction of [OII] emitters appears to be almost constant over the surface density of galaxies between 0.2 and 10 Mpc^-2. This trend is different from that seen in the local universe where the star-formation activity is weaker in higher density regions. To understand this difference between z~1 and z~0, we study the fraction of non-isolated galaxies as a function of local galaxy density. We find that the fraction of non-isolated galaxies increases with increasing density. Our results suggest that the star formation in galaxies at z~1 is triggered by galaxy interaction and/or mergers.

We present new deep optical spectra of 9 high-z radio galaxies (HzRGs) at z > 2.7 obtained with FORS2 on VLT. These rest-frame ultraviolet spectra are used to infer the metallicity of the narrow-line regions (NLRs) in order to investigate the chemical evolution of galaxies in high-z universe. We focus mainly on the CIV/HeII and CIII]/CIV flux ratios that are sensitive to gas metallicity and ionization parameter. Although the NV emission has been widely used to infer the gas metallicity, it is often too weak to be measured accurately for NLRs. By combining our new spectra with data from the literature, we examine the possible redshift evolution of the NLR metallicity for 57 HzRGs at 1 < z < 4. Based on the comparison between the observed emission-line flux ratios and the results of our photoionization model calculations, we find no significant metallicity evolution in NLRs of HzRGs, up to z ~ 4. Our results imply that massive galaxies had almost completed their chemical evolution at much higher redshift (z > 5). Finally, although we detect strong NV emission lines in 5 HzRGs at z > 2.7, we point out that high NV/HeII ratios are not indicative of high metallicities but correspond to high ionization parameters of gas clouds in NLRs.

We present first results of a study aimed to constrain the star formation rate and dust content of galaxies at z~2. We use a sample of BzK-selected star-forming galaxies, drawn from the COSMOS survey, to perform a stacking analysis of their 1.4 GHz radio continuum as a function of different stellar population properties, after removing AGN contaminants from the sample. Dust unbiased star formation rates are derived from radio fluxes assuming the local radio-IR correlation. The main results of this work are: i) specific star formation rates are constant over about 1 dex in stellar mass and up to the highest stellar mass probed; ii) the dust attenuation is a strong function of galaxy stellar mass with more massive galaxies being more obscured than lower mass objects; iii) a single value of the UV extinction applied to all galaxies would lead to grossly underestimate the SFR in massive galaxies; iv) correcting the observed UV luminosities for dust attenuation based on the Calzetti recipe provide results in very good agreement with the radio derived ones; v) the mean specific star formation rate of our sample steadily decreases by a factor of ~4 with decreasing redshift from z=2.3 to 1.4 and a factor of ~40 down the local Universe. These empirical SFRs would cause galaxies to dramatically overgrow in mass if maintained all the way to low redshifts, we suggest that this does not happen because star formation is progressively quenched, likely starting from the most massive galaxies.

From a study of BAL + IR + Fe II QSOs (using deep Gemini GMOS-IFU spectroscopy) new results are presented: for IRAS 04505-2958. Specifically, we have studied in detail the out flow (OF) process and their associated structures, mainly at two large galactic scales: (i) two blobs/shells (S1, S2) at radius r = 1.1 and 2.2 kpc; and (ii) an external hypergiant shell (S3) at r = 11 kpc. In addition, the presence of two very extended hypergiant shells (S4, S5) at r = 80 kpc is discussed. From this GMOS study the following main results were obtained: (i) For the external hypergiant shell S3 the kinematics GMOS maps of the ionized gas show very similar features to those observed for the prototype of exploding external supergiant shell: in NGC 5514. (ii) The main knots K1, K2 and K3 -of this hypergiant shell S3- show a stellar population and emission line ratios associated with the presence of a starburst + OF/shocks. (iii) The internal shells S1 and S2 show structures, OF components and properties very similar to those detected in the nuclear shells of Mrk 231. (iv) The shells S1+S2 and S3 are aligned at PA = 131: i.e. suggesting that the OF process is in the blow-out phase with bipolar structure. In addition, the shells S4 and S5 (at 80-100 kpc scale) are aligned at PA = 40, i.e.: a bipolar OF perpendicular to the internal OF. Finally, the generation of UHE cosmic rays and neutrino/ dark-matter -associated with HyNe in BAL + IR + Fe II QSOs- is discussed.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra} program that has imaged the central 0.5 sq.deg of the COSMOS field (centered at 10h, +02deg) with an effective exposure of ~160ksec, and an outer 0.4sq.deg. area with an effective exposure of ~80ksec. The limiting source detection depths are 1.9e-16 erg cm(-2) s(-1) in the Soft (0.5-2 keV) band, 7.3e(-16) erg cm^-2 s^-1 in the Hard (2-10 keV) band, and 5.7e(-16) erg cm(-2) s(-1) in the Full (0.5-10 keV) band. Here we describe the strategy, design and execution of the C-COSMOS survey, and present the catalog of 1761 point sources detected at a probability of being spurious of <2e(-5) (1655 in the Full, 1340 in the Soft, and 1017 in the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing positions with the ACIS-I imager, a remarkably uniform (to 12%) exposure across the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower flux limit. The widely different PSFs obtained in each exposure at each point in the field required a novel source detection method, because of the overlapping tiling strategy, which is described in a companion paper. (Puccetti et al. Paper II). This method produced reliable sources down to a 7-12 counts, as verified by the resulting logN-logS curve, with sub-arcsecond positions, enabling optical and infrared identifications of virtually all sources, as reported in a second companion paper (Civano et al. Paper III). The full catalog is described here in detail, and is available on-line.

When, and how, did the first galaxies and supermassive black holes (SMBH) form, and how did they reionization the Universe? First galaxy formation and cosmic reionization are among the last frontiers in studies of cosmic structure formation. We delineate the detailed astrophysical probes of early galaxy and SMBH formation afforded by observations at centimeter through submillimeter wavelengths. These observations include studies of the molecular gas (= the fuel for star formation in galaxies), atomic fine structure lines (= the dominant ISM gas coolant), thermal dust continuum emission (= an ideal star formation rate estimator), and radio continuum emission from star formation and relativistic jets. High resolution spectroscopic imaging can be used to study galaxy dynamics and star formation on sub-kpc scales. These cm and mm observations are the necessary compliment to near-IR observations, which probe the stars and ionized gas, and X-ray observations, which reveal the AGN. Together, a suite of revolutionary observatories planned for the next decade from centimeter to X-ray wavelengths will provide the requisite panchromatic view of the complex processes involved in the formation of the first generation of galaxies and SMBHs, and cosmic reionization.

We present results of a survey for Ly alpha emitters at z=4.86 based on optical narrowband (lambda_c=7126 angstrom, Delta lambda=73 angstrom) and broadband (B, V, r', i', and z') observations of the Cosmic Evolution Survey (COSMOS) field using Suprime-Cam on the Subaru Telescope. We find 79 LAE candidates at z=4.86 over a contiguous survey area of 1.83 deg^2, down to the Ly alpha line flux of 1.47 x 10^-17 ergs s^-1 cm^-2. We obtain the Ly alpha luminosity function with a best-fit Schechter parameters of log L^*=42.9^+0.5_-0.3 ergs s^-1 and phi^* = 1.2^+8.0_-1.1 x 10^-4 Mpc^-3 for alpha=-1.5 (fixed). The two-point correlation function for our LAE sample is xi(r) = (r/4.4^+5.7_-2.9 Mpc)^-1.90+/-0.22. In order to investigate the field-to-field variations of the properties of Ly alpha emitters, we divide the survey area into nine tiles of 0.5^circ x 0.5^circ each. We find that the number density varies with a factor of ~ 2 from field to field with high statistical significance. However, we find no significant field-to-field variance when we divide the field into four tiles with 0.7^circ x 0.7^circ each. We conclude that at least 0.5 deg^2 survey area is required to derive averaged properties of LAEs at z~5, and our survey field is wide enough to overcome the cosmic variance.

We study the spatial clustering of 538 X-ray selected AGN in the 2 deg^2 XMM-COSMOS field that are spectroscopically identified to I_{AB}<23 and span the redshift range z=0.2-3.0. The median redshift and luminosity of the sample are z = 0.98 and L_{0.5-10}=6.3 x 10^{43} erg/s, respectively. A strong clustering signal is detected at ~18sigma level, which is the most significant measurement obtained to date for clustering of X-ray selected AGN. By fitting the projected correlation function w(r_p) with a power law on scales of r_p=0.3-40 Mpc/h, we derive a best fit comoving correlation length of r_0 = 8.6 +- 0.5 Mpc/h and slope of gamma=1.88 +- 0.07 (Poissonian errors; bootstrap errors are about a factor of 2 larger). An excess signal is observed in the range r_p~5-15 Mpc/h, which is due to a large scale structure at z ~ 0.36 containing about 40 AGN. When removing the z ~ 0.36 structure, or computing w(r_p) in a narrower range around the peak of the redshift distribution (e.g. z=0.4-1.6), the correlation length decreases to r_0 ~ 5-6 Mpc/h, which is consistent with that observed for bright optical QSOs at the same redshift. We investigate the clustering properties of obscured and unobscured AGN separately. Within the statistical uncertainties, we do not find evidence that AGN with broad optical lines (BLAGN) cluster differently from AGN without broad optical lines (non-BLAGN). The correlation length measured for XMM-COSMOS AGN at z~1 is similar to that of massive galaxies (stellar mass M_*> 3 x 10^{10} M_sun) at the same redshift. This suggests that AGN at z~1 are preferentially hosted by massive galaxies, as observed both in the local and in the distant (z~2) Universe. (shortened)

We report the detection of CO molecular line emission in the z=4.5 millimeter-detected galaxy COSMOS_J100054+023436 (hereafter: J100+0234) using the IRAM Plateau de Bure interferometer (PdBI) and NRAO's Very Large Array (VLA). The CO(4-3) line as observed with PdBI has a full line width of ~1000 km/s, an integrated line flux of 0.66 Jy km/s, and a CO luminosity of 3.2e10 L_sun. Comparison to the 3.3sigma detection of the CO(2-1) line emission with the VLA suggests that the molecular gas is likely thermalized to the J=4-3 transition level. The corresponding molecular gas mass is 2.6e10 M_sun assuming an ULIRG-like conversion factor. From the spatial offset of the red- and blue-shifted line peaks and the line width a dynamical mass of 1.1e11 M_sun is estimated assuming a merging scenario. The molecular gas distribution coincides with the rest-frame optical and radio position of the object while being offset by 0.5'' from the previously detected Ly$\alpha$ emission. J1000+0234 exhibits very typical properties for lower redshift (z~2) sub-millimeter galaxies (SMGs) and thus is very likely one of the long sought after high redshift (z>4) objects of this population. The large CO(4-3) line width taken together with its highly disturbed rest-frame UV geometry suggest an ongoing major merger about a billion years after the Big Bang. Given its large star formation rate (SFR) of >1000 M_sun/yr and molecular gas content this object could be the precursor of a 'red-and-dead' elliptical observed at a redshift of z=2.

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.

We present photometric redshifts and spectral energy distribution (SED) classifications for a sample of 1542 optically identified sources detected with XMM in the COSMOS field. Our template fitting classifies 46 sources as stars and 464 as non-active galaxies, while the remaining 1032 require templates with an AGN contribution. High accuracy in the derived photometric redshifts was accomplished as the result of 1) photometry in up to 30 bands with high significance detections, 2) a new set of SED templates including 18 hybrids covering the far-UV to mid-infrared, which have been constructed by the combination of AGN and non-active galaxies templates, and 3) multi-epoch observations that have been used to correct for variability (most important for type 1 AGN). The reliability of the photometric redshifts is evaluated using the sub-sample of 442 sources with measured spectroscopic redshifts. We achieved an accuracy of $\sigma_{\Delta z/(1+z_{spec})} = 0.014$ for i$_{AB}^*<$22.5 ($\sigma_{\Delta z/(1+z_{spec})} \sim0.015$ for i$_{AB}^*<$24.5). The high accuracies were accomplished for both type 2 (where the SED is often dominated by the host galaxy) and type 1 AGN and QSOs out to $z=4.5$. The number of outliers is a large improvement over previous photometric redshift estimates for X-ray selected sources (4.0% and 4.8% outliers for i$_{AB}^*<$22.5 and i$_{AB}^*<$24.5, respectively). We show that the intermediate band photometry is vital to achieving accurate photometric redshifts for AGN, whereas the broad SED coverage provided by mid infrared (Spitzer/IRAC) bands is important to reduce the number of outliers for normal galaxies.

We present accurate photometric redshifts in the 2-deg2 COSMOS field. The redshifts are computed with 30 broad, intermediate, and narrow bands covering the UV (GALEX), Visible-NIR (Subaru, CFHT, UKIRT and NOAO) and mid-IR (Spitzer/IRAC). A chi2 template-fitting method (Le Phare) was used and calibrated with large spectroscopic samples from VLT-VIMOS and Keck-DEIMOS. We develop and implement a new method which accounts for the contributions from emission lines (OII, Hbeta, Halpha and Ly) to the spectral energy distributions (SEDs). The treatment of emission lines improves the photo-z accuracy by a factor of 2.5. Comparison of the derived photo-z with 4148 spectroscopic redshifts (i.e. Delta z = zs - zp) indicates a dispersion of sigma_{Delta z/(1+zs)}=0.007 at i<22.5, a factor of 2-6 times more accurate than earlier photo-z in the COSMOS, CFHTLS and COMBO-17 survey fields. At fainter magnitudes i<24 and z<1.25, the accuracy is sigma_{Delta z/(1+zs)}=0.012. The deep NIR and IRAC coverage enables the photo-z to be extended to z~2 albeit with a lower accuracy (sigma_{Delta z/(1+zs)}=0.06 at i~24). The redshift distribution of large magnitude-selected samples is derived and the median redshift is found to range from z=0.66 at 22<i<22.5 to z=1.06 at 24.5<i<25. At i<26.0, the multi-wavelength COSMOS catalog includes approximately 607,617 objects. The COSMOS-30 photo-z enable the full exploitation of this survey for studies of galaxy and large scale structure evolution at high redshift.

We present an analysis of the radio properties of large samples of Lyman Break Galaxies (LBGs) at $z \sim 3$, 4, and 5 from the COSMOS field. The median stacking analysis yields a statistical detection of the $z \sim 3$ LBGs (U-band drop-outs), with a 1.4 GHz flux density of $0.90 \pm 0.21 \mu$Jy. The stacked emission is unresolved, with a size $< 1"$, or a physical size $< 8$kpc. The total star formation rate implied by this radio luminosity is $31\pm 7$ $M_\odot$ year$^{-1}$, based on the radio-FIR correlation in low redshift star forming galaxies. The star formation rate derived from a similar analysis of the UV luminosities is 17 $M_\odot$ year$^{-1}$, without any correction for UV dust attenuation. The simplest conclusion is that the dust attenuation factor is 1.8 at UV wavelengths. However, this factor is considerably smaller than the standard attenuation factor $\sim 5$, normally assumed for LBGs. We discuss potential reasons for this discrepancy, including the possibility that the dust attenuation factor at $z \ge 3$ is smaller than at lower redshifts. Conversely, the radio luminosity for a given star formation rate may be systematically lower at very high redshift. Two possible causes for a suppressed radio luminosity are: (i) increased inverse Compton cooling of the relativistic electron population due to scattering off the increasing CMB at high redshift, or (ii) cosmic ray diffusion from systematically smaller galaxies. The radio detections of individual sources are consistent with a radio-loud AGN fraction of 0.3%. One source is identified as a very dusty, extreme starburst galaxy (a 'submm galaxy').

We present our analysis on the clustering properties of star-forming galaxies selected by narrow-band excesses in the Subaru Deep Field. Specifically we focus on Halpha emitting galaxies at z = 0.24 and z = 0.40 in the same field, to investigate possible evolutionary signatures of clustering properties of star-forming galaxies. Based on the analysis on 228 Halpha emitting galaxies with 39.8 < log L(Halpha) < 40.8 at z = 0.40, we find that their two-point correlation function is estimated as xi = (r/1.62^{+0.64}_{-0.50} Mpc)^{-1.84 +/- 0.08}. This is similar to that of Halpha emitting galaxies in the same Halpha luminosity range at z = 0.24, xi = (r/1.88^{+0.60}_{-0.49} Mpc)^{-1.89 +/- 0.07}. These correlation lengths are smaller than those for the brighter galaxy sample studied by Meneux et al. (2006) in the same redshift range. The evolution of correlation length between z = 0.24 and z = 0.40 is interpreted by the gravitational growth of the dark matter halos.

We make a search for Halpha emitting galaxies at z=0.24 in the Subaru Deep Field (SDF) using the archival data set obtained with the Subaru Telescope. We carefully select Halpha emitters in the narrowband filter NB816, using B, V, Rc, i', and z' broad-band colors. We obtain a sample of 258 emitting galaxies with observed equivalent widths of (Halpha+[NII]6548,6584) greater than 12 angstrom. We also analyze a sample of Halpha emitters taken from the Sloan Digital Sky Survey (SDSS) to constrain the luminous end of Halpha luminosity function. Using the same selection criteria as for the SDF, and after excluding AGNs, we obtain 317 Halpha emitting star-forming galaxies. Combining these two samples of Halpha emitters found in both SDF and SDSS, we derive a Halpha luminosity function with best-fit Schechter function parameters of alpha = -1.31^+0.17_-0.17, log phi^* = -2.46^+0.34_-0.40 Mpc^-3, log L^* = 41.99^+0.08_-0.07 ergs s^-1. An extinction-corrected Halpha luminosity density is 4.45^+2.96_-1.75 x 10^39 ergs s^-1 Mpc^-3. Using the Kennicutt relation between the Halpha luminosity and star formation rate, the star formation rate density in the survey volume is estimated as 0.035^+0.024_-0.014 M_sun yr^-1 Mpc^-3. The angular two-point correlation function of Halpha emitters over 875 arcmin^2 at z = 0.24 is well fitted by a power-law form with w(theta) = 0.047^+0.017_-0.013 theta^-0.66 +- 0.08, corresponding to the correlation function of xi(r) = (r/2.6^+1.0_-0.8 Mpc)^(-1.66 +- 0.08). The small correlation length of Halpha emitters may imply the weak clustering of active star-forming galaxies.

New results of a Programme of study of BAL + IR + Fe II QSOs (at low and high redshift) are presented. Which are based mainly on deep Gemini GMOS integral field unit (IFU/3D) spectroscopy. We have performed a detailed study of the kinematics, morphological, and physical conditions, in the BAL + IR + Fe II QSO: IRAS 04505-2958. From this study, some selected results are presented, mainly for the 3 expanding giant shells (observed with Gemini). In particular, the GMOS data suggest that the outflow (OF) process -in this IR QSO- generated multiple expanding hypergiant shells (from 10, to 100 kpc), in several extreme explosive events. These new Gemini GMOS data are in good agreement with our evolutionary, explosive and composite Model: where part of the ISM of the host galaxy is ejected in the form of multiple giant shells, mainly by HyN explosions. This process could generate satellite/companion galaxies, and even could expel a high fraction -or all- the host galaxy. In addition, this Model for AGN could give important clue about the physical processes that could explain the origin -in AGNs- of very energetic cosmic rays, detected by the P. Auger Observatory.

We present recent progress in searching for galaxies at redshift from z = 5 to z = 10. Wide-field and senstive surveys with 8m class telescopes have been providing more than several hundreds of star forming galaxies at z =5 - 7 that are probed in the optical window. These galaxies are used to study the early cosmic star formation activity as well as the early structure formation in the universe. Moreover, near infrared deep imaging and spectropscopic surveys have found probable candidates of galaxies from z = 7 to z = 10. Although these candidates are too faint to be identified unambiguously, we human being are now going to the universe beyond 13 billion light years, close to the epoch of first-generations stars; i.e., Population III stars. We also mention about challanges to find Population III-dominated galaxies in the early universe.

We present a new photometric search for high-z galaxies hosting Population III (PopIII) stars based on deep intermediate-band imaging observations obtained in the Subaru Deep Field (SDF), by using Suprime-Cam on the Subaru Telescope. By combining our new data with the existing broad-band and narrow-band data, we searched for galaxies which emit strongly both in Ly_alpha and in HeII 1640 (``dual emitters'') that are promising candidates for PopIII-hosting galaxies, at 3.93<z<4.01 and 4.57<z<4.65. Although we found 10 ``dual emitters'', most of them turn out to be [OII]-[OIII] dual emitters or H_beta-(H_alpha+[NII]) dual emitters at z<1, as inferred from their broad-band colors and from the ratio of the equivalent widths. No convincing candidate of Ly_alpha-HeII dual emitter of SFR_PopIII > 2 Msun/yr was found by our photometric search in 4.03 x 10^5 Mpc^3 in the SDF. This result disfavors low feedback models for PopIII star clusters, and implies an upper-limit of the PopIII SFR density of SFRD_PopIII < 5 x 10^-6 Msun/yr/Mpc^3. This new selection method to search for PopIII-hosting galaxies should be useful in future narrow-band surveys to achieve the first observational detection of PopIII-hosting galaxies at high redshifts.

To derive a new H$\alpha$ luminosity function and to understand the clustering properties of star-forming galaxies at $z \approx 0.24$, we have made a narrow-band imaging survey for H$\alpha$ emitting galaxies in the HST COSMOS 2 square degree field. We used the narrow-band filter NB816 ($\lambda_c = 8150$ \AA, $\Delta \lambda = 120$ \AA) and sampled H$\alpha$ emitters with $EW_{\rm obs}(\rm H\alpha + [N\textsc{ii}]) > 12$ \AA in a redshift range between $z=0.233$ and $z=0.251$ corresponding to a depth of 70 Mpc. We obtained 980 H$\alpha$ emitting galaxies in a sky area of 5540 arcmin$^2$, corresponding to a survey volume of $3.1 \times 10^4 {\rm Mpc^3}$. We derive a H$\alpha$ luminosity function with a best-fit Schechter function parameter set of $\alpha = -1.35^{+0.11}_{-0.13}$, $\log\phi_* = -2.65^{+0.27}_{-0.38}$, and $\log L_* ({\rm erg s^{-1}}) = 41.94^{+0.38}_{-0.23}$. The H$\alpha$ luminosity density is $2.7^{+0.7}_{-0.6} \times 10^{39}$ ergs s$^{-1}$ Mpc$^{-3}$. After subtracting the AGN contribution (15 %) to the H$\alpha$ luminosity density, the star formation rate density is evaluated as $1.8^{+0.7}_{-0.4} \times 10^{-2}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. The angular two-point correlation function of H$\alpha$ emitting galaxies of $\log L({\rm H\alpha}) > 39.8$ is well fit by a power law form of $w(\theta) = 0.013^{+0.002}_{-0.001} \theta^{-0.88 \pm 0.03}$, corresponding to the correlation function of $\xi(r) = (r/1.9{\rm Mpc})^{-1.88}$. We also find that the H$\alpha$ emitters with higher H$\alpha$ luminosity are more strongly clustered than those with lower luminosity.

In this work we present the first results of a study of BAL QSOs (at low and high redshift), based on very deep Gemini GMOS integral field spectroscopy. In particular, the results obtained for the nearest BAL IR QSO Mrk 231 are presented. Very deep three-dimensional spectra and maps clearly show that the BAL systems I and II are extended (reaching 1.4-1.6" = 1.2-1.3 kpc, from the nucleus) and clearly elongated at the position angle close to the radio jet PA. Which suggest that the BAL systems I and II are both associated with the radio jet, and supporting the bipolar jet-wind model for some BALs. For the nuclear region of Mrk 231, the QSO and host-galaxy components were modelled, using a new technique of decoupling 3D spectra. From this study, the following main results were found: (i) in the pure host galaxy spectrum an extreme nuclear starburst component was clearly observed, mainly as a very strong increase in the flux, at the blue wavelengths; (ii) the BAL system I is observed in the spectrum of the host galaxy; (iii) in the clean/pure QSO emission spectrum, only broad lines were detected. 3D GMOS individual spectra (specially in the IR Ca II triplet) and maps confirm the presence of an extreme and young nuclear starburst (8 < age < 15 Myr), which was detected mainly in a ring or toroid with a radius r = 0.3" = 200 pc, around the very nucleus. The physical properties of the four expanding nuclear bubbles were analysed, using the GMOS 3D spectra and maps. These results suggest that an important part of the nuclear NLR is generated by the OF process and the associated low velocity ionizing shocks.

In this work we analise the role and evidence of exploding BAL + IR + Fe II QSOs, and their relation with new -and previous- explosive models for evolution, formation and end of galaxies.

We measure the fraction of luminous galaxies in pairs at projected separations of 5-20 kpc out to z=1.2 in the COSMOS field using ACS images and photometric redshifts derived from an extensive multiwavelength dataset. Analysis of a complete sample of 106,188 galaxies more luminous than M_V=-19.8 (~ L_V*) in the redshift range 0.1 < z < 1.2 yields 1,749 galaxy pairs. These data are supplemented by a local z=0-0.1 value for the galaxy pair fraction derived from the Sloan Digital Sky Survey (SDSS). After statistically correcting the COSMOS pair sample for chance line-of-sight superpositions, the evolution in the pair fraction is fit by a power law \propto (1+z)^{n=3.1 \pm 0.1}. If this strongly evolving pair fraction continues out to higher redshift, ~ 50% of all luminous galaxies at z ~ 2 are in close pairs. This clearly signifies that galaxy mergers are a very significant and possibly dominant mechanism for galaxy evolution during the epoch of galaxy formation at z=1 to 3.

Using the initial peculiar velocity field, we analytically study the hierarchical formation of gravitationally bound objects. The field is smoothed over a scale that corresponds to the mass of a given class of objects. Through the Zel'dovich approximation, the smoothed field determines how the objects cluster together to form a new class of more massive objects. The standard cosmological parameters lead to the evolution of primordial clouds with 10^6 M(sun) -> galaxies with 10^12 M(sun) -> clusters of galaxies with 10^15 M(sun) -> superclusters of galaxies with 10^16 M(sun). The epochs obtained for the formation of these classes of objects are consistent with observations.

We present the results of six-epoch Very Long Baseline Array (VLBA) observations of 3C~66A. The high-resolution Very Long Baseline Interferometer (VLBI) maps obtained at multi-frequency (2.3, 8.4, and 22.2 GHz) simultaneously enabled us to identify the brightest compact component with the core. We find that the spectrum of the core can be reasonably fitted by the synchrotron self-absorption model. Our VLBA maps show that the jet of 3C~66A has two bendings at about 1.2 and 4 mas from the core. We also give possible identifications of our jet components with the components in previous VLBA observations by analysing their proper motions. We find consistent differences of the position from the core in one component between different frequencies at six epochs.

We present measurements of the two-point galaxy angular correlation function w(\theta) in the COSMOS field. Independent determinations of w(\theta) as a function of magnitude limit are presented for both the HST ACS catalog and also for the ground-based data from Subaru and the CFHT. Despite having significantly different masks, these three determinations agree well. At bright magnitudes (IAB<22), our data generally match very well with existing measurements and with mock catalogs based on semi-analytic galaxy formation calculations of Kitzbichler and White from the Millennium Simulation. The exception is that our result is at the upper end of the expected cosmic variance scatter for \theta > 10 arcmin, which we attribute to a particularly rich structure known to exist at z~0.8. For fainter samples, however, the level of clustering is somewhat higher than reported by some previous studies: in all three catalogues we find w(\theta=1')~0.014 at a median IAB magnitude of 24. At these very faintest magnitudes, our measurements agree well with the latest determinations from the Canada-France Legacy Survey. This level of clustering is approximately double what is predicted by the semi-analytic catalogs (at all angles). The semi-analytic results allow an estimate of cosmic variance, which is too small to account for the discrepancy. We therefore conclude that the mean amplitude of clustering at this level is higher than previously estimated.

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

We have carried out a wide-field imaging survey for [OII]3727 emitting galaxies at z~1.2 in the HST COSMOS 2 square degree field using the Suprime-Cam on the Subaru Telescope. The survey covers a sky area of 6700 arcmin^2 in the COSMOS field, and a redshift range between 1.17 and 1.20 (Delta_z = 0.03), corresponding to a survey volume of 5.56*10^5 Mpc^3. We obtain a sample of 3176 [OII] emitting galaxies with observed emission-line equivalent widths greater than 26 AA. Since our survey tends to sample brighter [OII]3727 emitting galaxies, we also analyze a sample of fainter [OII]3727 emitting galaxies found in the Subaru Deep Field (SDF). We find an extinction-corrected [OII] luminosity density of 10^{40.35^+0.08_-0.06} ergs s^-1 Mpc-3, corresponding to star formation rate density of 0.32^+0.06_-0.04 M_sun yr-1 Mpc^-3 in the COSMOS field at z~1.2. This is the largest survey for [OII]3727 emitters beyond z=1 currently available.

We present the first identification of large-scale structures (LSS) at z $< 1.1$ in the Cosmic Evolution Survey (COSMOS). The structures are identified from adaptive smoothing of galaxy counts in the pseudo-3d space ($\alpha,\delta$,z) using the COSMOS photometric redshift catalog. The technique is tested on a simulation including galaxies distributed in model clusters and a field galaxy population -- recovering structures on all scales from 1 to 20\arcmin without {\it a priori} assumptions for the structure size or density profile. Our procedure makes {\bf no} {\it a priori} selection on galaxy spectral energy distribution (SED, for example the Red Sequence), enabling an unbiased investigation of environmental effects on galaxy evolution. The COSMOS photometric redshift catalog yields a sample of $1.5\times10^5$ galaxies with redshift accuracy, $\Delta z_{FWHM}/(1+z) \leq 0.1$ at z $< 1.1$ down to I$_{AB} \leq 25$ mag. Using this sample of galaxies, we identify 42 large-scale structures and clusters. abstract truncated for astroph 25 line limit -- see preprint

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

We present results from a narrow-band optical survey of a contiguous area of 1.95 deg^2, covered by the Cosmic Evolution Survey (COSMOS). Both optical narrow-band (lambda_c = 8150 AA and Delta_lambda = 120 AA) and broad-band (B, V, g', r', i', and z') imaging observations were performed with the Subaru prime-focus camera, Suprime-Cam on the Subaru Telescope. We provide the largest contiguous narrow-band survey, targetting Ly alpha emitters (LAEs) at z~5.7. We find a total of 119 LAE candidates at z~5.7. Over the wide-area covered by this survey, we find no strong evidence for large scale clustering of LAEs. We estimate a star formation rate (SFR) density of ~7*10^-4 M_sun yr^-1 Mpc^-3 for LAEs at z~5.7, and compare it with previous measurements.

We report new follow-up spectroscopy of i-dropout galaxies with an NB921-band depression found in the Subaru Deep Field. The NB921-depressed i-dropout selection method is expected to select galaxies with large equivalent width Ly alpha emission over a wide redshift range, 6.0<z<6.5. Two of four observed targets show a strong emission line with a clear asymmetric profile, identified as Ly alpha emitters at z=6.11 and 6.00. Their rest-frame equivalent widths are 153A and 114A, which are lower limits on the intrinsic equivalent widths. Through our spectroscopic observations (including previous ones) of NB921-depressed i-dropout galaxies, we identified 5 galaxies in total with a rest-frame equivalent width larger than 100A at 6.0<z<6.5 out of 8 photometric candidates, which suggests that the NB921-depressed i-dropout selection method is possibly an efficient way to search for Ly alpha emitters with a large Ly alpha equivalent width, in a wider redshift range than usual narrow-band excess techniques. By combining these findings with our previous observational results, we infer that the fraction of broad-band selected galaxies having a rest-frame equivalent width larger than 100A is significantly higher at z~6 (the cosmic age of ~1 Gyr) than that at z~3 (~2 Gyr), being consistent with the idea that the typical stellar population of galaxies is significantly younger at z~6 than that at z~3. The NB921-depressed i-dropout galaxies may be interesting candidates for hosts of massive, zero-metallicity Population III stars.

ABRIDGED: We study the evolution since z~1 of the rest-frame B luminosity function of the early-type galaxies (ETGs) in ~0.7 deg^2 in the COSMOS field. In order to identify ALL progenitors of local ETGs we construct the sample of high-z galaxies using two complementary criteria: (i) A morphological selection based on the Zurich Estimator of Structural Types, and (ii) A photometric selection based on the galaxy properties in the (U-V)-M_V color-magnitude diagram. We furthermore constrain both samples so as to ensure that the selected progenitors of ETGs are compatible with evolving into systems which obey the mu_B-r_{hl} Kormendy relation. Assuming the luminosity evolution derived from studies of the fundamental plane for high-z ETGs, our analysis shows no evidence for a decrease in the number density of the most massive ETGs out to z~ 0.7: Both the morphologically- and the photometrically-selected sub-samples show no evolution in the number density of bright (~L>2.5L*) ETGs. Allowing for different star formation histories, and cosmic variance, we estimate a maximum decrease in the number density of massive galaxies at that redshift of ~30%. We observe, however, in both the photometrical and morphological samples, a deficit of up to ~2-3 of fainter ETGs over the same cosmic period. Our results argue against a significant contribution of recent dissipationless ``dry'' mergers to the formation of the most massive ETGs. We suggest that the mass growth in low luminosity ETGs can be explained with a conversion from z~0.7 to z=0 of blue, irregular and disk galaxies into low- and intermediate-mass ``red'' ETGs, possibly also through gas rich mergers.

We use observations from the Giant Metrewave Radio Telescope (GMRT) to measure the atomic hydrogen gas content of star-forming galaxies at z = 0.24 (i.e. a look-backtime of ~3 Gyr). The sample of galaxies studied were selected from Halpha-emitting field galaxies detected in a narrow-band imaging survey with the Subaru Telescope. The Anglo-Australian Telescope was used to obtain precise optical redshifts for these galaxies. We then coadded the HI 21 cm emission signal for all the galaxies within the GMRT spectral line data cube. From the coadded signal of 121 galaxies, we measure an average atomic hydrogen gas mass of (2.26 +- 0.90)*10^9 solar masses. We translate this HI signal into a cosmic density of neutral gas at z = 0.24 of Omega_gas = (0.91 +- 0.42)*10^-3. This is the current highest redshift at which Omega_gas has been constrained from 21 cm emission and our value is consistent with that estimated from damped Lyman-alpha systems around this redshift. We also find that the correlations between the Halpha luminosity and the radio continuum luminosity and between the star formation rate and the HI gas content in star-forming galaxies at z = 0.24 are consistent with the correlations found at z = 0. These two results suggest that the star formation mechanisms in field galaxies ~3 Gyr ago were not substantially different from the present, even though the star formation rate is 3 times higher.

Ordinary baryonic particles (such as protons and neutrons) account for only one-sixth of the total matter in the Universe. The remainder is a mysterious "dark matter" component, which does not interact via electromagnetism and thus neither emits nor reflects light. As dark matter cannot be seen directly using traditional observations, very little is currently known about its properties. It does interact via gravity, and is most effectively probed through gravitational lensing: the deflection of light from distant galaxies by the gravitational attraction of foreground mass concentrations. This is a purely geometrical effect that is free of astrophysical assumptions and sensitive to all matter -- whether baryonic or dark. Here we show high fidelity maps of the large-scale distribution of dark matter, resolved in both angle and depth. We find a loose network of filaments, growing over time, which intersect in massive structures at the locations of clusters of galaxies. Our results are consistent with predictions of gravitationally induced structure formation, in which the initial, smooth distribution of dark matter collapses into filaments then into clusters, forming a gravitational scaffold into which gas can accumulate, and stars can be built.

We study the enviromental dependence and the morphological composition of the galaxy color-magnitude diagram at z~0.7, using a pilot sub-sample of ~2000 galaxies from the COSMOS surve, with I_AB<24, photometric redshift within 0.61<z<0.85 and galaxy morphologies based on the HST-ACS data. The C-M diagram shows a clear red-sequence dominated by early-type galaxies and also a remarkably well-defined "blue sequence'" described by late-type objects. While the percentage of objects populating the two sequences is a function of environment, with a clear morphology/color-density relation, their normalization and slope are independent of local density. We identify and study a number of objects with "anomalous'' colors, given their morphology. Red late-type galaxies are found to be mostly highly-inclined spiral galaxies, with colors dominated by internal reddening (33% contamination with respect to truly passive spheroidals). Conversely, the population of blue early-type galaxies is composed by objects of moderate luminosity and mass, concurring to only ~5% of the mass in spheroidal galaxies. The majority of them (~70%) occupy a position in the surface-brightness/effective-radius plane not consistent with them being precursors of current epoch elliptical galaxies. In a color-mass diagram, color sequences are even better defined, with red galaxies covering in general a wider range of masses at nearly constant color, and blue galaxies showing a more pronounced dependence of color on mass. While the red sequence is adequately reproduced by models of passive evolution, the blue sequence is better interpreted as a specific star-formation sequence. Its substantial invariance with respect to local density suggests that the overall, "secular'' star formation is driven more by galaxy mass than by environment. [...]

We present a study of the morphological nature of redshift z~0.7 star-forming galaxies using a combination of HST/ACS, GALEX and ground-based images of the COSMOS field. Our sample consists of 8,146 galaxies, 5,777 of which are detected in the GALEX near-ultraviolet band down to a limiting magnitude of 25.5 (AB). We make use of the UV to estimate star formation rates, correcting for the effect of dust using the UV-slope, and compute, from the ACS F814W images, the C,A,S,G,M20 morphological parameters for all objects in our sample. We observe a morphological bimodality in the galaxy population and show that it has a strong correspondence with the FUV - g color bimodality. We conclude that UV-optical color predominantly evolves concurrently with morphology. We observe many of the most star-forming galaxies to have morphologies approaching that of early-type galaxies, and interpret this as evidence that strong starburst events are linked to bulge growth and constitute a process through which galaxies can be brought from the blue to the red sequence while simultaneously modifying their morphology accordingly. We conclude that the red sequence has continued growing at z~<0.7. We also observe z~0.7 galaxies to have physical properties similar to that of local galaxies, except for higher star formation rates. Whence we infer that the dimming of star-forming galaxies is responsible for most of the evolution in the star formation rate density of the Universe since that redshift, although our data are also consistent with a mild number evolution. [abridged]

We have identified a large-scale structure at z~0.73 in the COSMOS field, coherently described by the distribution of galaxy photometric redshifts, an ACS weak-lensing convergence map and the distribution of extended X-ray sources in a mosaic of XMM observations. The main peak seen in these maps corresponds to a rich cluster with Tx= 3.51+0.60/-0.46 keV and Lx=(1.56+/-0.04) x 10^{44} erg/s ([0.1-2.4] keV band). We estimate an X-ray mass within $r500$ corresponding to M500~1.6 x 10^{14} Msun and a total lensing mass (extrapolated by fitting a NFW profile) M(NFW)=(6+/-3) x 10^15 Msun. We use an automated morphological classification of all galaxies brighter than I_AB=24 over the structure area to measure the fraction of early-type objects as a function of local projected density Sigma_10, based on photometric redshifts derived from ground-based deep multi-band photometry. We recover a robust morphology-density relation at this redshift, indicating, for comparable local densities, a smaller fraction of early-type galaxies than today. Interestingly, this difference is less strong at the highest densities and becomes more severe in intermediate environments. We also find, however, local "inversions'' of the observed global relation, possibly driven by the large-scale environment. In particular, we find direct correspondence of a large concentration of disk galaxies to (the colder side of) a possible shock region detected in the X-ray temperature map and surface brightness distribution of the dominant cluster. We interpret this as potential evidence of shock-induced star formation in existing galaxy disks, during the ongoing merger between two sub-clusters.

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

We present observations at 1.4 and 250 GHz of the $z\sim 5.7$ Ly$\alpha$ emitters (LAE) in the COSMOS field found by Murayama et al.. At 1.4 GHz there are 99 LAEs in the lower noise regions of the radio field. We do not detect any individual source down to 3$\sigma$ limits of $\sim 30\mu$Jy beam$^{-1}$ at 1.4 GHz, nor do we detect a source in a stacking analysis, to a 2$\sigma$ limit of $2.5\mu$Jy beam$^{-1}$. At 250 GHz we do not detect any of the 10 LAEs that are located within the central regions of the COSMOS field covered by MAMBO ($20' \times 20'$) to a typical 2$\sigma$ limit of $S_{250} < 2$mJy. The radio data imply that there are no low luminosity radio AGN with $L_{1.4} > 6\times 10^{24}$ W Hz$^{-1}$ in the LAE sample. The radio and millimeter observations also rule out any highly obscured, extreme starbursts in the sample, ie. any galaxies with massive star formation rates $> 1500$ M$_\odot$ year$^{-1}$ in the full sample (based on the radio data), or 500 M$_\odot$ year$^{-1}$ for the 10% of the LAE sample that fall in the central MAMBO field. The stacking analysis implies an upper limit to the mean massive star formation rate of $\sim 100$ M$_\odot$ year$^{-1}$.

We present a detailed spectral analysis of point-like X-ray sources in the XMM-COSMOS field. Our sample of 135 sources only includes those that have more than 100 net counts in the 0.3-10 keV energy band and have been identified through optical spectroscopy. The majority of the sources are well described by a simple power-law model with either no absorption (76%) or a significant intrinsic, absorbing column (20%).As expected, the distribution of intrinsic absorbing column densities is markedly different between AGN with or without broad optical emission lines. We find within our sample four Type-2 QSOs candidates (L_X > 10^44 erg/s, N_H > 10^22 cm^-2), with a spectral energy distribution well reproduced by a composite Seyfert-2 spectrum, that demonstrates the strength of the wide field XMM/COSMOS survey to detect these rare and underrepresented sources.

We report on the stellar content of the COSMOS two degree field, as derived from a rigorous star/galaxy separation approach developed for using stellar sources to define the point spread function variation map used in a study of weak galaxy lensing. The catalog obtained in one filter from the ACS (Advanced Camera for Survey on the Hubble Space Telescope) is cross-identified with ground based multi-wavelength catalogs. The classification is reliable to magnitude $F_{814W}=24$ and the sample is complete even fainter. We construct a color-magnitude diagram and color histograms and compare them with predictions of a standard model of population synthesis. We find features corresponding to the halo subdwarf main sequence turnoff, the thick disk, and the thin disk. This data set provides constraints on the thick disk and spheroid density laws and on the IMF at low mass. We find no evidence of a sharp spheroid edge out to this distance. We identify a blue population of white dwarfs with counts that agree with model predictions. We find a hint for a possible slight stellar overdensity at about 22-34 kpc but the data are not strong enough at present to claim detection of a stream feature in the halo (abridged).

We present the results of a search for galaxy clusters in the first 36 XMM-Newton pointings on the COSMOS field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 3x10-15 ergs cm-2 s-1, having one of the widest XMM-Newton contiguous raster surveys, covering an area of 2.1 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. Verification of the cluster candidates is done based on a galaxy concentration analysis in redshift slices of thickness of 0.1-0.2 in redshift, using the multi-band photometric catalog of the COSMOS field and restricting the search to z<1.3 and i_AB < 25. We identify 72 clusters and derive their properties based on the X-ray cluster scaling relations. A statistical description of the survey in terms of the cumulative log(N>S)-lg(S) distribution compares well with previous results, although yielding a somewhat higher number of clusters at similar fluxes. The X-ray luminosity function of COSMOS clusters matches well the results of nearby surveys, providing a comparably tight constraint on the faint end slope of alpha=1.93+/-0.04. For the probed luminosity range of 8x10+42 - 2x10+44 ergs s-1, our survey is in agreement with and adds significantly to the existing data on the cluster luminosity function at high redshifts and implies no substantial evolution at these luminosities to z=1.3.

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift <1.5, while the redshift distribution of the full BL AGN population peaks at z~1.5. Our analysis also suggests that the Type 2/Type 1 ratio decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral analysis and the most recent modeling of the X-ray luminosity function evolution.

We measure photometric redshifts and spectral types for galaxies in the COSMOS survey. We use template fitting technique combined with luminosity function priors and with the option to simultaneously estimate dust extinction (i.e. E(B-V)) for each galaxy.Our estimated redshifts are accurate to i<25 and z~1.2. Using simulations with sampling and noise characteristics similar to those in COSMOS, the accuracy and reliability is estimated for the photometric redshifts as a function of the magnitude limits of the sample, S/N ratios and the number of bands used. From the simulations we find that the ratio of derived 95% confidence interval in the redshift probability distribution to the estimated photometric redshift (D95) can be used to identify and exclude the catastrophic failures in the photometric redshift estimates. We compare the derived redshifts with high-reliability spectroscopic redshifts for a sample of 868 normal galaxies with z < 1.2 from zCOSMOS. Considering different scenarios, depending on using prior, no prior and/or extinction, we compare the photometric and spectroscopic redshifts for this sample. This corresponds to an rms scatter of 0.031, with a small number of outliers (<2.5%). We also find good agreement (rms=0.10) between photometric and spectroscopic redshifts for Type II AGNs. We compare results from our photometric redshift procedure with three other independent codes and find them in excellent agreement. We show preliminary results, based on photometric redshifts for the entire COSMOS sample (to i < 25 mag.).

We report on the discovery of a new potential galaxy threshing system in the COSMOS 2 square degree field using the prime-focus camera, Suprime-Cam, on the 8.2 m Subaru Telescope. This system consists of a giant elliptical galaxy with $M_V \approx -21.6$ and a tidally disrupted satellite galaxy with $M_V \approx -17.7$ at a photometric redshift of $z \approx 0.08$. This redshift is consistent with the spectroscopic redshift of 0.079 for the giant elliptical galaxy obtained from the Sloan Digital Sky Survey (SDSS) archive. The luminosity masses of the two galaxies are $3.7 \times 10^{12} \cal{M}_{\odot}$ and $3.1 \times 10^{9} \cal{M}_{\odot}$, respectively. The distance between the two galaxies is greater than 100 kpc. The two tidal tails emanating from the satellite galaxy extend over 150 kpc. This system would be the second well-defined galaxy threshing system found so far.

The Cosmic Evolution Survey (COSMOS) is designed to probe the correlated evolution of galaxies, star formation, active galactic nuclei (AGN) and dark matter (DM) with large-scale structure (LSS) over the redshift range z $> 0.5 $ to 6. The survey includes multi-wavelength imaging and spectroscopy from X-ray to radio wavelengths covering a 2 $\sq$\deg area, including HST imaging. Given the very high sensitivity and resolution of these datasets, COSMOS also provides unprecedented samples of objects at high redshift with greatly reduced cosmic variance, compared to earlier surveys. Here we provide a brief overview of the survey strategy, the characteristics of the major COSMOS datasets, and summarize the science goals.

We present deep optical imaging observations of 2 square degree area, covered by the Cosmic Evolution Survey (COSMOS), made by the prime-focus Camera (Supreme-Cam) on the 8.2m Subaru Telescope. Observations were done in six broad-band [B (4459.7 AA), g' (4723.1 AA), V (5483.8 AA), r' (6213.0 AA), i' (7640.8 AA), z' (8855.0 AA)], and one narrow-band (NB816) filters. A total of 10^6 galaxies were detected to i'~26.5 mag. These data, combined with observations at u* and K-band are used to construct the photometric catalogs for the COSMOS and to measure their photometric redshifts, multi-band spectral energy distributions, stellar masses and identification of high redshift candidates. This catalog provides multi-waveband data for scientific analysis of the COSMOS survey.

zCOSMOS is a large redshift survey that is being undertaken in the COSMOS field using 600 hours of observation with the VIMOS spectrograph on the 8-m VLT. The survey is designed to characterise the environments of COSMOS galaxies from the 100 kpc scales of galaxy groups up to the 100 Mpc scale of the cosmic web and to produce diagnostic information on galaxies and active galactic nuclei. The zCOSMOS survey consists of two parts: (a) zCOSMOS-bright, a magnitude-limited I-band IAB < 22.5 sample of about 20,000 galaxies with 0.1 < z < 1.2 covering the whole 1.7 deg2 COSMOS ACS field and designed to mimic the parameters of the 2dfGRS; and (b) zCOSMOS-deep, a survey of approximately 10,000 galaxies selected through colour-selection criteria to have 1.4 < z < 3.0, within the central 1 deg2. This paper describes the survey design and the construction of the target catalogues, and briefly outlines the observational program and the data pipeline. In the first observing season, spectra of 1303 zCOSMOS-bright targets and of 977 zCOSMOS-deep targets have been obtained. These are briefly analysed to demonstrate the characteristics that may be expected from zCOSMOS, and particularly zCOSMOS-bright, when it is finally completed between 2008-2009. The power of combining spectroscopic and photometric redshifts is demonstrated, especially in correctly identifying the emission line in single-line spectra and in determining which of the less reliable spectroscopic redshifts are correct and which are incorrect. Our zCOSMOS-deep spectra demonstrate the effectiveness of our selection techniques to isolate high redshift galaxies at 1.4 < z < 3.0 and of VIMOS to measure their redshifts using ultraviolet absorption lines.

(ABRIDGED) Motivated by the desire to reliably and automatically classify structure of thousands of COSMOS galaxies, we present ZEST, the Zurich Estimator of Structural Types. To classify galaxy structure, ZEST uses: (i) Five non-parametric diagnostics: asymmetry, concentration, Gini coefficient, 2nd-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (ii) The exponent n of single--Sersic fits to the 2D surface brightness distributions. To fully exploit the wealth of information while reducing the redundancy present in these diagnostics, ZEST performs a principal component (PC) Analysis. We use a sample of ~56,000 I<24 COSMOS galaxies to show that the first three PCs fully describe the key aspects of the galaxy structure, i.e., to calibrate a three-dimensional classification grid of axis PC_1, PC_2, and PC_3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. (1996). The ZEST classification breaks most of the degeneracy between different galaxy populations that affects morphological classifications based on only some of the diagnostics included in ZEST. As a first application, we present the evolution since z~1 of the Luminosity Functions of COSMOS galaxies of early, disk and irregular galaxies and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find that the LF up to a redshift z=1 is consistent with a pure-luminosity evolution (of about 0.95 magnitudes at z \~0.7). We highlight however two trends, that are in general agreement with a down-sizing scenario for galaxy formation: (1.) A deficit of a factor of about two at z~0.7 of MB>-20.5 structurally--classified early--type galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of irregular galaxies.

We have identified a complex galaxy cluster system in the COSMOS field via a wide angle tail (WAT) radio galaxy consistent with the idea that WAT galaxies can be used as tracers of clusters. The WAT galaxy, CWAT-01, is coincident with an elliptical galaxy resolved in the HST-ACS image. Using the COSMOS multiwavelength data set, we derive the radio properties of CWAT-01 and use the optical and X-ray data to investigate its host environment. The cluster hosting CWAT-01 is part of a larger assembly consisting of a minimum of four X-ray luminous clusters within ~2 Mpc distance. We apply hydrodynamical models that combine ram pressure and buoyancy forces on CWAT-01. These models explain the shape of the radio jets only if the galaxy's velocity relative to the intra-cluster medium (ICM) is in the range of about 300-550 km/s which is higher than expected for brightest cluster galaxies (BCGs) in relaxed systems. This indicates that the CWAT-01 host cluster is not relaxed, but is possibly dynamically young. We argue that such a velocity could have been induced through subcluster merger within the CWAT-01 parent cluster and/or cluster-cluster interactions. Our results strongly indicate that we are witnessing the formation of a large cluster from an assembly of multiple clusters, consistent with the hierarchical scenario of structure formation. We estimate the total mass of the final cluster to be approximately 20% of the mass of the Coma cluster.

Simulations of nearby (0.015 < z < 0.025) SDSS galaxies have been used to reproduce as accurately as possible the appearance that they would have on COSMOS ACS images if they had been observed at z ~ 0.7 and z ~ 1.2. By adding the SDSS galaxies to random locations in the COSMOS images, we simulate the effects of chance superpositions of high redshift galaxies with unrelated foreground or background objects. We have used these simulated images, together with those of real COSMOS galaxies at these same redshifts, to undertake a "blind" morphological classification of galaxies to identify those that appear to be undergoing mergers and thus to estimate the change in merger fraction with redshift. We find that real mergers are harder to recognize at high redshift, and also that the chance superposition of unrelated galaxies often produces the appearance of mergers where in reality none exists. In particular, we estimate that 1.5 - 2.0% of objects randomly added to ACS images are misclassified as mergers due to projection with unrelated objects, and as a result, that 40% of the apparent mergers in COSMOS at z=0.7 are likely to be spurious. We find that the fraction of galaxies undergoing mergers increases as (1+z)^3.8+/-1.2 to z ~ 0.7 and that this trend appears to continue to z = 1.2. Merger candidates at z ~ 0.7 are bluer than the parent population, especially when the statistical effects of the chance projections are accounted for. Merger candidates are more asymmetric than the population as a whole, and are often associated with irregular morphology. Nevertheless, the majority (~60%) of the merger candidates appear to be associated with spiral galaxies although in this case we cannot correct for the effects of chance projections.

We present ZEBRA, the Zurich Extragalactic Bayesian Redshift Analyzer. The current version of ZEBRA combines and extends several of the classical approaches to produce accurate photometric redshifts down to faint magnitudes. In particular, ZEBRA uses the template-fitting approach to produce Maximum Likelihood and Bayesian redshift estimates based on: (1.) An automatic iterative technique to correct the original set of galaxy templates to best represent the SEDs of real galaxies at different redshifts; (2.) A training set of spectroscopic redshifts for a small fraction of the photometric sample; and (3.) An iterative technique for Bayesian redshift estimates, which extracts the full two-dimensional redshift and template probability function for each galaxy. We demonstrate the performance of ZEBRA by applying it to a sample of 866 I_AB <= 22.5 COSMOS galaxies with available u*, B, V, g', r', i', z' and K_s photometry and zCOSMOS spectroscopic redshifts in the range 0 < z < 1.3. Adopting a 5-sigma-clipping that excludes less than 10 galaxies, both the Maximum Likelihood and Bayesian estimates for this sample have an accuracy dz/1+z better than 0.03. Similar accuracies are recovered using mock galaxies.

We study a sample of approximately 16,500 galaxies with I_AB <= 22.5 in the COSMOS field. Structural information on the galaxies is derived by fitting single Sersic models to their two-dimensional surface brightness distributions. We investigate the evolution of the number density of disk galaxies larger than 5 kpc between redshift z~1 and the present epoch. To this end, we use the measurements of the half-light radii to construct, as a function of redshift, the size function of both the total disk galaxy population and of disk galaxies split in four bins of bulge-to-disk ratio. Furthermore, we use a selected sample of roughly 1800 SDSS galaxies to calibrate our results with respect to the local universe. We find that: (i) The number density of disk galaxies with intermediate sizes (r_{1/2}~5-7 kpc) remains nearly constant from z~1 to today. (ii) The number density of the largest disks (r_{1/2}>7 kpc) decreases by a factor of about two out to z~1. (iii) There is a constancy in the number density of large bulgeless disks out to z~1; the deficit of large disks at early epochs seems to arise from a smaller number of bulged disks. Our results indicate that the bulk of the large disk galaxy population has completed its growth by z~1, and support the hypothesis that secular evolution processes produce - or at least add stellar mass to - the bulge components of disk galaxies.

There is observed a trend that a lower mass galaxy forms stars at a later epoch. This downsizing of star-forming galaxies has been attributed to hydrodynamical or radiative feedback processes that regulate star formation. However, here we explain the downsizing by gravitational processes alone, in the bottom-up scenario where galaxies evolve from subgalactic-scale objects. Within a region of the initial density field that is to evolve into a lower mass galaxy, subgalactic-scale fluctuation is of a smaller amplitude. The formation of subgalactic-scale objects, i.e., gravitational collapse of the subgalactic-scale fluctuation, and the subsequent onset of star formation accordingly occur at a later epoch for a lower mass galaxy. As a function of galaxy mass, we calculate the peak epoch of formation of subgalactic-scale objects. The peak epoch is consistent with the peak epoch of star formation derived from observations. [abridged]

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.

Recently, Zheng et al. (2005) found evidence for an overdensity of galaxies around a radio-loud quasar, SDSS J0836+0054, at z=5.8 (a five arcmin$^2$ region). We have examined our deep optical imaging data (B, V, r', i', z', and NB816) taken with the Suprime-Cam on the Subaru Telescope. The NB816 narrow-band filter (lambda_c = 815 nm and $\Delta\lambda = 12$ nm) is suitable for searching for Ly$\alpha$ emitters at $z\approx 5.7$. We have found a new strong Ly$\alpha$ emitter at $z \approx 5.7$ close to object B identified by Zheng et al. Further, the non detection of the nine objects selected by Zheng et al. (2005) in our B, V, and r' images provides supporting evidence that they are high-z objects.

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

We discuss observational properties of strong emission-line galaxies at low redshift found by our deep imaging survey for high-redshift Ly alpha emitters. In our surveys, we used the narrowband filter, NB816 (lambda_center=8150A with FWHM = 120A), and the intermediate-band filter, IA827 (lambda_center = 8270A with FWHM = 340A). In this survey, 62 NB816-excess (> 0.9 mag) and 21 IA827-excess (> 0.8 mag) objects were found. Among them, we found 20 NB816-excess and 4 IA827-excess Ly alpha emitter candidates. Therefore, it turns out that 42 NB816-excess and 17 IA827-excess objects are strong emission-line objects at lower redshift. Since 4 objects in the two low-z samples are common, the total number of strong low-z emitters is 55. Applying our photometric redshift technique, we identify 7 H alpha emitters at z~0.24, 20 H beta-[OIII] ones at z~0.65, and 11 [OII] ones at z~1.19. However, we cannot determine reliable photometric redshifts of the remaining 17 emitters. The distributions of their rest frame equivalent widths are consistently understood with recent studies of galaxy evolution from z~1 to z~0.

We present results of our intermediate-band optical imaging survey for high-$z$ Ly$\alpha$ emitters (LAEs) using the prime focus camera, Suprime-Cam, on the 8.2m Subaru Telescope. In our survey, we use eleven filters; four broad-band filters ($B$, $R_{\rm c}$, $i^\prime$, and $z^\prime$) and seven intermediate-band filters covering from 500 nm to 720 nm; we call this imaging program as the Mahoroba-11. The seven intermediate-band filters are selected from the IA filter series that is the Suprime-Cam intermediate-band filter system whose spectral resolution is $R = 23$. Our survey has been made in a $34^\prime \times 27^\prime$ sky area in the Subaru XMM Newton Deep Survey field. We have found 409 IA-excess objects that provide us a large photometric sample of strong emission-line objects. Applying the photometric redshift method to this sample, we obtained a new sample of 198 LAE candidates at $3 < z < 5$. We found that there is no evidence for evolution of the number density and the star formation rate density for LAEs with $\log L({\rm Ly}\alpha) ({\rm erg s^{-1}}) > 42.67$ between $z \sim 3$ and 5.

We report the discovery of a large-scale coherent filamentary structure of Lyman alpha emitters in a redshift space at z=3.1. We carried out spectroscopic observations to map the three dimensional structure of the belt-like feature of the Lyman alpha emitters discovered by our previous narrow-band imaging observations centered on the protocluster at z=3.1. The feature was found to consist of at least three physical filaments connecting with each other. The result is in qualitative agreement with the prediction of the 'biased' galaxy-formation theories that galaxies preferentially formed in large-scale filamentary or sheet-like mass overdensities in the early Universe. We also found that the two known giant Lyman alpha emission-line nebulae showing high star-formation activities are located near the intersection of these filaments, which presumably evolves into a massive cluster of galaxies in the local Universe. This may suggest that massive galaxy formation occurs at the characteristic place in the surrounding large-scale structure at high redshift.

We present results from optical narrow-band lambda_c = 8150A ~ and Delta lambda = 120A) observations of the Great Observatories Origins Deep Survey (GOODS) fields, using Suprime-Cam on the Subaru Telescope. Using these narrow-band data, we then perform a survey of Lyman alpha Emitters (LAEs) at z~5.7. The LAE survey covers an area of approx 320 arcmin^2 and a co-moving volume of ~8.0 x 10^4 Mpc^3. We found a total of 10 (GOODS-N) and 4 (GOODS-S) LAE candidates at z~5.7. We perform a study of the spatial distribution, space density, and star formation properties of the LAEs at z~5.7.

We discuss on the early stage of galaxy formation based on recent deep surveys for very high-redshift galaxies, mostly beyond redshift of 6. These galaxies are observed to be strong Lyman$\alpha$ emitters, indicating bursts of massive star formation in them. The fraction of such star-forming system appears to increase with increasing redshift. On the other hand, the star formation rate density derived from Lyman$\alpha$ emitters tends to decrease with increasing redshift. It is thus suggested that the major epoch of initial starbursts may occur around $z \sim$ 6 -- 7. In order to understand the early stage of galaxy formation, new surveys for galaxies beyond redshift of 7 will be important in near future.

We present a very deep near-infrared spectroscopic observation of a strong Ly_alpha emitter at z=6.33, SDF J132440.6+273607, which we used to search for HeII 1640. This emission line is expected if the target hosts a significant number of population III stars. Even after 42 ksec of integration with the Subaru/OHS spectrograph, no emission-line features are detected in the JH band, which confirms that SDF J132440.6+273607 is neither an active galactic nucleus nor a low-$z$ emission-line object. We obtained a 2sigma upper-limit of 9.06e-18 ergs/s/cm^2 on the HeII 1640 emission line flux, which corresponds to a luminosity of 4.11e42 ergs/s. This upper-limit on the HeII 1640 luminosity implies that the upper limit on population III star-formation rate is in the range 4.9--41.2 M_sun/yr if population III stars suffer no mass loss, and in the range 1.8--13.2 M_sun/yr if strong mass loss is present. The non-detection of HeII in SDF J132440.6+273607 at z=6.33 may thus disfavor weak feedback models for population III stars.

We present the results from a submm survey of a sample of 23 giant Lya emitting nebulae in the overdensity at z=3.09 in the SA22 field. These objects, which have become known as Lya Blobs (LABs) have a diverse range of morphology and surface brightness, but the nature of their power source is unclear - with cooling flows and/or AGN/starburst ionised winds being possibilities. Using the SCUBA submm camera we measure the 850um flux of a sample of LABs, detecting four LABs at >3.5sigma individually, and a modest statistical detection of the full sample at about 3mJy. These fluxes correspond to bolometric luminosities in the ultraluminous regime, with star-formation rates of about 1e3 Msun/yr. We show there is a trend between Lya luminosity and bolometric output, which suggests that a galactic scale superwind generated from starbursts of age 10-100Myr may be responsible for the Lya emission. We estimate the star-formation rate density in SA22 to be >3 Msun/yr/Mpc^3 - greater than the field at this epoch, and note that there are now 7 submm galaxies in the SA22 structure, making this region the richest association of these intensely active galaxies. Finally we suggest that Lya haloes may be a common feature of the submm population in general, and have an important role in the heating and enrichment of the intergalactic medium.

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.

In order to search for high-redshift galaxies beyond $z = 6.6$ in the Subaru Deep Field, we have investigated NB921-dropout galaxies where NB921 is the narrowband filter centered at 919.6 nm with FWHM of 13.2 nm for the Suprime-Cam on the Subaru Telescope. There are no secure NB921-dropout candidates brighter than $z^\prime = 25.5$. Based on this result, we discuss the UV luminosity function of star-forming galaxies at $z > 6.6$.

By exploiting the far-infrared(FIR) and radio correlation, we have performed the Likelihood-Ratio analysis to identify optical counterparts to the far-infrared sources in the Lockman Hole. Using the likelihood ratio analysis and the associated reliability, 44 FIR sources have been identified with radio sources. Redshifts have been obtained for 29 out of 44 identified sources. One hyper-luminous infrared galaxy (HyLIRG) with and four ultraluminous infrared galaxies (ULIRGs) are identified in our sample. The space density of the FIR sources at z = 0.3-0.6 is 4.6\times 10^{-5}Mpc^{-3}, implying a rapid evolution of the ULIRG population. Most of \ISO FIR sources have their FIR-radio ratios similar to star-forming galaxies ARP 220 and M82. At least seven of our FIR sources show evidence for the presence of an active galactic nucleus (AGN) in optical emission lines, radio continuum excess, or X-ray activity. Three out of five (60%) of the ULIRG/HyLIRGs are AGN galaxies. Five of the seven AGN galaxies are within the ROSAT X-ray survey field, and two are within the XMM-Newton survey fields. X-ray emission has been detected in only one source, 1EX030, which is optically classified as a quasar. The non-detection in the XMM-Newton 2-10 keV band suggests a very thick absorption obscuring the central source of the two AGN galaxies. Several sources have an extreme FIR luminosity relative to the optical R-band, L(90\mu\mathrm{m})/L(R) > 500, which is rare even among the local ULIRG population. While source confusion or blending might offer an explanation in some cases, they may represent a new population of galaxies with an extreme activity of star formation in an undeveloped stellar system -- i.e., formation of bulges or young ellipticals.

In order to search for high-redshift galaxies around $z \sim 6$ in the Subaru Deep Field, we have investigated NB816-dropout galaxies where NB816 is the narrowband filter centered at 815nm with FWHM of 12.5 nm for the Suprime-Cam on the Subaru Telescope. Since the NB816 imaging is so deep we can detect 10 well-defined NB816-dropout galaxies that we identify to lie at $z$ = 5.8 -- 6.5. We discuss their observational properties.

The Cosmic Evolution Survey (COSMOS) is a Hubble Space Telescope (HST) treasury project.The COSMOS aims to perform a 2 square degree imaging survey of an equatorial field in $I$(F814W) band, using the Advanced Camera for Surveys (ACS). Such a wide field survey, combined with ground-based photometric and spectroscopic data, is essential to understand the interplay between large scale structure, evolution and formation of galaxies and dark matter. In 2004, we have obtained high-quality, broad band images of the COSMOS field ($B, V, r^\prime, i^\prime,$ and $ z^\prime$) using Suprime-Cam on the Subaru Telescope, and we have started our new optical multi-band program, COSMOS-21 in 2005. Here, we present a brief summary of the current status of the COSMOS project together with contributions from the Subaru Telescope. Our future Subaru program, COSMOS-21, is also discussed briefly.

We propose a new method to search for high-redshift galaxies that is based on an intermediate-band dropout technique rather than the usual broad-band dropout one. In this method, we use an intermediate-band filter whose central wavelength is longer than 7000 \AA. This new method makes it possible to distinguish both very late-type stars such as L and T dwarfs and dusty galaxies at intermediate redshift from real high-$z$ Lyman break galaxies. The reason for this is that such interlopers do not show strong intermediate-band depression although they have very red broad-band colors that are indicative of Lyman break galaxies. Applying our new method to imaging data sets obtained with the Suprime-Cam on the Subaru Telescope, we find a new sample of Lyman break galaxies at $z \simeq 5$.

We present new results of a deep optical imaging survey using a narrowband filter ($NB921$) centered at $\lambda =$ 9196 \AA ~ together with $B$, $V$, $R$, $i^\prime$, and $z^\prime$ broadband filters in the sky area of the Subaru Deep Field which has been promoted as one of legacy programs of the 8.2m Subaru Telescope. We obtained a photometric sample of 58 Ly$\alpha$ emitter candidates at $z \approx$ 6.5 -- 6.6 among $\sim 180$ strong $NB921$-excess ($z^\prime - NB921 > 1.0$) objects together with a color criterion of $i^\prime - z^\prime > 1.3$. We then obtained optical spectra of 20 objects in our $NB921$-excess sample and identified at least nine Ly$\alpha$ emitters at $z \sim 6.5$ -- 6.6 including the two emitters reported by Kodaira et al. (2003). Since our Ly$\alpha$ emitter candidates are free from strong amplification of gravitational lensing, we are able to discuss their observational properties from a statistical point of view. Based on these new results, we obtain a lower limit of the star formation rate density of $\rho_{\rm SFR} \simeq 5.5 \times 10^{-4}$ $h_{0.7}$ $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ at $z \approx 6.6$, being consistent with our previous estimate. We discuss the nature of star-formation activity in galaxies beyond $z=6$.

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

We present deep spatially-resolved optical spectroscopy of the NW companion galaxy of the quasar BR 1202-0725 at $z=4.7$. Its rest-frame UV spectrum shows star-forming activity in the nuclear region. The Ly$\alpha$ emission profile is symmetric with wavelength while its line width is unusually wide (FWHM $\simeq 1100$ km s$^{-1}$) for such a high-$z$ star-forming galaxy. Spectrum taken along the Ly$\alpha$ nebula elongation, which is almost along the minor axis of the companion host galaxy, reveals that off-nuclear Ly$\alpha$ nebulae have either flat-topped or multi-peaked profiles along the extension. All these properties can be understood in terms of superwind activity in the companion galaxy. We also find a diffuse continuum component around the companion, which shows similar morphology to that of Ly$\alpha$ nebula, and is most likely due to scattering of the quasar light at dusty halo around the companion. We argue that the superwind could expel dusty material out to the halo region, making a dusty halo for scattering.

A study of morphology, kinematics and ionization structure of the IR merger NGC 5514, is presented. This study is based mainly on Integral 2D spectroscopy obtained on 4.2 m William Herschel Telescope (La Palma, Spain) and data from CASLEO (Argentina). Clear evidence of two extra-nuclear starbursts with young outflows and LINER activity are reported. One of these outflow generate a supergiant bubble, where the emission lines and kinematics maps show 4 extended ejections. These results suggest that the bubble is in the rupture phase.

We report on the result of our high quality spectropolarimetric observation for the narrow-line region in a nearby Seyfert 2 galaxy, Mrk 573, by the Subaru Telescope. The polarized flux spectrum of Mrk 573 shows not only prominent scattered broad H$\alpha$ emission but also various narrow forbidden emission lines. We find that the measured polarization degree of the observed forbidden emission lines is positively correlated with the ionization potential of the corresponding ions and the critical density of the corresponding transitions. We discuss some possible origins of these correlations, and then we point out that the correlations are caused due to the obscuration of the stratified narrow-line region in Mrk 573 by the optically and geometrically thick dusty torus, just similar to the previous study on NGC 4258.

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.

In this Work we report, for the extreme starburst in the IR merger NGC 3256: (i) The detection of 4 galactic bubbles, associated with SN explosions.(ii) The first analysis of the spatial distribution of young star clusters (YSC) candidates.(iii) The kinematic study of the ionized gas in the core of of the main optical nucleus, performed with HST STIS spectra. The shape of the rotation curve and the emission line profile could be explained by the presence in the core of YSC with outflow.

We present the results of a survey for Ly-alpha emitters at z ~ 5.8 using a new intermediate-band filter centered at lambda_c = 8275 AA with Delta-lambda_FWHM = 340 AA (i.e., the spectroscopic resolution is R ~ 23) with a combination with a traditional narrow-band centered at lambda_c = 8150 AA with Delta-lambda_FWHM = 120 AA (R ~ 68). Our observations were made with use of the Subaru Prime Focus Camera, Suprime-Cam, on the 8.2 m Subaru telescope in a sky area surrounding the high redshift quasar, SDSSp J104433.04-012522.2 at z=5.74, covering an effective sky area with ~ 720 arcmin^2. In this survey, we have found four Ly-alpha-emitter candidates from the intermediate-band image (z \~ 5.8 with Delta z ~ 0.3). Combined with our previous results based on the NB816 imaging, we discuss the star formation activity in galaxies between z ~ 5.7 and z ~ 5.9.

We report the properties of the 35 robust candidates of Ly-alpha blobs (LABs), which are larger than 16 arcsec^2 in isophotal area and brighter than 0.7 x 10^-16 ergs s^-1 cm^-2, searched in and around the proto-cluster region at redshift z=3.1 discovered by Steidel et al. in the SSA22 field, based on wide-field (31'x23') and deep narrow-band (NB497; 4977/77) and broad-band (B,V, and R) images taken with the prime-focus camera on the Subaru telescope. The two previously known giant LABs are the most luminous and the largest ones in our survey volume of 1.3 x 10^5 Mpc^3. We revealed the internal structures of the two giant LABs and discovered some bubble-like features, which suggest that intensive starburst and galactic superwind phenomena occurred in these objects in the past. The rest 33 LABs have isophotal area of about 16-78 arcsec^2 and flux of 0.7-7 x 10^-16 ergs s^-1 cm^-2. These 35 LABs show a continuous distribution of isophotal area and emission line flux. The distributions of average surface brightness and morphology are widespread from relatively compact high surface brightness objects to very diffuse low surface brightness ones. The physical origins of these LABs may be (i) photo-ionization by massive stars, or active galactic nuclei, or (ii) cooling radiation from gravitationally heated gas, or (iii) shock heating by starburst driven galactic superwind. One third of them are apparently not associated with ultra-violet continuum sources that are bright enough to produce Ly-alpha emission, assuming a Salpeter initial mass function. The 90% of these LABs are located inside the high surface density region of the 283 relatively compact and strong Ly-alpha emitters selected in our previous study. This suggests that these LABs may be the phenomena related to dense environment at high redshift.

Powerful radio galaxies at high redshift are highly useful in studies of early evolution of AGN-hosting galaxies because their observed optical and near infrared light are dominated by their stellar population rather than the nonthermal continuum emitted by the central engine of AGNs. In addition, the presence of AGN activity in them implies that a supermassive black hole has been already made in their nuclei. These properties allow us to investigate a possible starburst-AGN connection in early universe and then provide some crucial hints for the formation mechanism of supermassive black holes. Taking observational properties of high-z powerful radio galaxies into account, we discuss a possible formation mechanism of supermassive black holes in their nuclei.

We report the discovery of the scattered emission from a hidden broad-line region (BLR) in a Seyfert 2 galaxy, Mrk 573, based on our recent spectropolarimetric observation performed at the Subaru Telescope. This object has been regarded as a type 2 AGN without a hidden BLR by the previous observations. However, our high quality spectrum of the polarized flux of Mrk 573 shows prominent broad (~3000 km/s) H_alpha emission, broad weak H_beta emission, and subtle Fe II multiplet emission. Our new detection of these indications for the presence of the hidden BLR in the nucleus of Mrk 573 is thought to be owing to the high signal-to-noise ratio of our data, but the possibility of a time variation of the scattered BLR emission is also mentioned. Some diagnostic quantities such as the IRAS color, the radio power, and the line ratio of the emission from the narrow-line region of Mrk 573 are consistent with the distributions of such quantities of type 2 AGNs with a hidden BLR. Mrk 573 is thought to be an object whose level of the AGN activity is the weakest among the type 2 AGNs with a hidden BLR. In terms of the systematic differences between the type 2 AGNs with and without a hidden BLR, we briefly comment on an interesting Seyfert 2 galaxy, Mrk 266SW, which may possess a hidden BLR but has been treated as a type 2 AGNs without a hidden BLR.

Galaxy interactions are studied in terms of the starburst-Seyfert connection. The starburst requires a high rate of gas supply. Since the efficiency for supplying the gas is high in a galaxy interaction, although the companion is not necessarily discernible, Seyfert galaxies with circumnuclear starbursts are expected to be interacting. Since the large amounts of circumnuclear gas and dust obscure the broad-line region, they are expected to be observed as Seyfert 2. The active galactic nucleus itself does not require a high rate of gas supply. Seyfert galaxies without circumnuclear starbursts are not necessarily expected to be interacting even at the highest luminosities. They are not necessarily expected to evolve from Seyfert galaxies with circumnuclear starbursts. We derive these and other theoretical expectations and confirm them with statistics on observational data of magnitude-limited samples of Seyfert galaxies.

We report results of our new spatially-resolved, optical spectroscopy of the giant Ly$\alpha$ nebula around a powerful radio galaxy 1243+036 (4C+03.24) at $z=3.57$. The nebula is extended over $\sim 30$ kpc from the nucleus, and forms a pair of cones or elongated bubbles. The high-velocity ($\sim -1000$ km s$^{-1}$; blueshifted with respect to the systemic velocity) Ly$\alpha$-emitting components are detected at both sides of the nucleus along its major axis. The northwestern nebula is more spectacular in its velocity shift (blueshifted by -1000 km s$^{-1}$ to -1400 km s$^{-1}$) and in its width ($\simeq 1900$ km s$^{-1}$ FWHM) over $\simeq 30$ kpc scale. We discuss possible origin of the nebula; 1) the shock-heated expanding bubble or outflowing cone associated with the superwind activity of the host galaxy, 2) halo gas photoionized by the anisotropic radiation from the active galactic nuclei (AGN), and 3) the jet-induced star-formation or shock. The last possibility may not be likely because Ly$\alpha$ emission is distributed out of the narrow channel of the radio jet. We show that the superwind model is most plausible since it can explain both the characteristics of the morphology (size and shape) and the kinematical structures (velocity shift and line width) of the nebula although the photoionization by AGN may contribute to the excitation to some extent.

The 8m class telescopes in the ground-based optical astronomy together with help from the ultra-sharp eye of the Hubble Space Telescope have enabled us to observe forming galaxies beyond redshift $z=5$. In particular, more than twenty Ly$\alpha$-emitting galaxies have already been found at $z > 5$. These findings provide us with useful hints to investigate how galaxies formed and then evolved in the early universe. Further, detailed analysis of Ly$\alpha$ emission line profiles are useful in exploring the nature of the intergalactic medium because the trailing edge of cosmic reionization could be close to $z \sim 6$ -- 7, at which forming galaxies have been found recently. We also discuss the importance of superwinds from forming galaxies at high redshift, which has an intimate relationship between galaxies and the intergalactic medium. We then give a review of early cosmic star formation history based on recent progress in searching for Ly$\alpha$-emitting young galaxies beyond redshift 5.

We present a summary of optical/NIR deep surveys for very high-$z$ galaxies using the 8.2m Subaru Telescope operated by National Astronomical Observatory of Japan. The prime focus mosaic CCD camera, Suprime-Cam, with a very wide field of view, 34$^\prime \times 27^\prime$, allows us to carry out efficient optical deep surveys. In particular, the Subaru Deep Field project has provided us a number of Lyman$\alpha$ emitters beyond $z=6$. We discuss the star formation history in the early universe based on this project.

We discuss the nature of faint 6.7um galaxies detected with the mid-infrared camera ISOCAM on board the Infrared Space Observatory (ISO). The 23 hour integration on the Hawaii Deep Field SSA13 has provided a sample of 65 sources down to 6uJy at 6.7um. For 57 sources, optical or near-infrared counterparts were found with a statistical method. All four Chandra sources, three SCUBA sources, and one VLA/FIRST source in this field were detected at 6.7um with high significance. Using their optical to mid-infrared colors, we divided the 6.7um sample into three categories: low redshift galaxies with past histories of rapid star formation, high redshift ancestors of these, and other star forming galaxies. Rapidly star forming systems at high redshifts dominate the faintest end. Spectroscopically calibrated photometric redshifts were derived from fits to a limited set of template SEDs. They show a high redshift tail in their distribution with faint (<30uJy) galaxies at z>1. The 6.7um galaxies tend to have brighter K magnitudes and redder I-K colors than the blue dwarf population at intermediate redshifts. Stellar masses of the 6.7um galaxies were estimated from their rest-frame near-infrared luminosities. Massive galaxies (M_star~10e11M_sun) were found in the redshift range of z=0.2-3. Epoch dependent stellar mass functions indicate a decline of massive galaxies' comoving space densities with redshift. Even with such a decrease, the contributions of the 6.7um galaxies to the stellar mass density in the universe are found to be comparable to those expected from UV bright galaxies detected in deep optical surveys.

We present an optical multicolor-imaging study of the galaxy cluster Abell 521 at $z = 0.25$, using Suprime-Cam on the Subaru Telescope, covering an area of $32 \times 20$ arcmin$^2$ ($9.4 \times 5.8 h_{50}^{-2}$ Mpc$^2$ at $z = 0.25$). Our imaging data taken with both a narrow-band filter, $NB816$ ($\lambda_0 = 8150$\AA and $\Delta \lambda = 120$\AA), and broad-band filters, $B,V,R_{\rm C}, i^\prime$, and $z^\prime$ allow us to find 165 H$\alpha$ emitters. We obtain the H$\alpha$ luminosity function (LF) for the cluster galaxies within 2 Mpc; the Schechter parameters are $\alpha = -0.75 \pm 0.23$, $\phi^\star = 10^{-0.25 \pm 0.20}$ Mpc$^{-3}$, and $L^\star = 10^{42.03 \pm 0.17}$ erg s$^{-1}$. Although the faint end slope, $\alpha$, is consistent with that of the local cluster H$\alpha$ LFs, the characteristic luminosity, $L^\star$, is about 6 times (or $\approx 2$ mag) brighter. This strong evolution implies that Abell 521 contains more active star-forming galaxies than the local clusters, being consistent with the observed Butcher-Oemler effect. However, the bright $L^\star$ of Abell 521 may be, at least in part, due to the dynamical condition of this cluster.

We present the results of a survey for Ly$\alpha$ emitters at $z\approx 5.7$ based on optical narrow-band ($\lambda_{\rm c} = 8150$ \AA ~ and $\Delta\lambda = 120$ \AA), and broad-band ($B$, $R_{\rm C}$, $I_{\rm C}$, and $z^\prime$) observations of the field surrounding the high redshift quasar, SDSSp J104433.04$-$012522.2, on the 8.2 m Subaru Telescope with the Subaru Prime Focus Camera, Suprime-Cam. This survey covers a sky area of $\approx 720$ arcmin$^2$ and a co-moving volume of $\simeq 2 \times 10^5 ~ h_{0.7}^{-3}$ Mpc$^3$. We have found 20 candidates of Ly$\alpha$ emitters at $z \approx$ 5.7 with $\Delta z \approx 0.1$. Two of them have been confirmed star-forming galaxies at $z=5.655$ and $z=5.687$ from our follow-up optical spectroscopy. We discuss star-formation properties of the 20 objects from a statistical point of view. Our survey leads to a new estimate of the star formation rate density at $z \approx 5.7$, $\sim 1.2 \times 10^{-3} h_{0.7} M_\odot$ yr$^{-1}$ Mpc$^{-3}$.

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

We present results of a deep mid-infrared survey in the SSA13 field with the Infrared Space Observatory (ISO). In order to probe the near-infrared light at high redshifts, we surveyed the field with the broad band LW2 (5-8.5um) filter of the mid-infrared camera ISOCAM. Adopting a highly redundant imaging strategy for the 23 hour observation and carefully treating gradual changes in the detector responsivity caused by a very high rate of cosmic ray impacts, we succeeded in reaching an 80% completeness limit of 16uJy in the central 7 arcmin2 region. Utilizing the signal-to-noise ratio map, we detected 65 sources down to 6uJy in the 16 arcmin2 field. Integral galaxy number counts at 6.7um are then derived, reaching 1.3e4/deg2 at the faint limit with a slope of -1.6 between 13uJy and 130uJy. Integrating individual sources in this flux range, the resolved fraction of the extragalactic background light at 6.7um is estimated to be 0.56nW/m2/sr. These results, which reach a flux limit three times fainter than those in the Hubble Deep Fields, are in fairly good agreement with a model prediction by Franceschini et al. (1997). Finally, we discuss the relation of distant massive E/S0 galaxies to the faint 6.7um galaxy population.

We report on a study of forbidden emission-line spectrum of nearby Seyfert 1.5 galaxy NGC 4151 based on the high-resolution (R ~ 45,000) optical spectrum obtained by using the High Dispersion Spectrograph boarded on the Subaru Telescope. The profile parameters such as the emission-line widths, the velocity shifts from the recession velocity of the host galaxy, and the asymmetry indices, for emission lines including very faint ones such as [Ar IV]4712,4740 and [Fe VI]5631,5677 are investigated. Statistically significant correlations between the measured profile parameters and the critical densities of transitions are found while there are no meaningful correlations between the profile parameters and the ionization potentials of ions. By comparing the results with photoionization model calculations, we remark that a simple power-law distribution of the gas density which is independent of the radius from the nucleus cannot explain the observed correlation between the emission-line widths and the critical densities of the transitions. Taking the additional dense gas component expected to exist at the innermost of the narrow-line region into account, the observed correlations between the emission-line width and the critical density of the transitions can be understood since high-critical-density emission lines can arise at such relatively inner regions even if their ionization potentials are low. The observed correlation between the blueshift amounts of emission lines and the critical densities of the ions is also explained if such dense gas clouds located closer to the nucleus have larger outflowing velocities.

Supermassive black hole binaries may exist in the centers of active galactic nuclei like quasars and radio galaxies and mergers between galaxies may result in the formation of supermassive binaries during the course of galactic evolution. Using the very-long-baseline interferometer, we imaged the radio galaxy 3C 66B at radio frequencies and found that the unresolved radio core of 3C 66B shows well-defined elliptical motions with a period of 1.05 +/- 0.03 years, which provides a direct detection of a supermassive black hole binary.

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 have performed a deep optical imaging survey using a narrowband filter ($NB921$) centered at $\lambda =$ 9196 \AA ~ together with $i^\prime$ and $z^\prime$ broadband filters covering an 814 arcmin$^2$ area of the Subaru Deep Field. We obtained a sample of 73 strong $NB921$-excess objects based on the following two color criteria; $z^\prime - NB921 > 1$ and $i^\prime - z^\prime > 1.3$. We then obtained optical spectroscopy of nine objects in our $NB921$-excess sample, and identified at least two Ly$\alpha$ emitters at$z=6.541 \pm 0.002$ and $z=6.578 \pm 0.002$, each of which shows the characteristic sharp cutoff together with the continuum depression at wavelengths shortward of the line peak. The latter object is more distant than HCM-6A at $z=6.56$ and thus this is the most distant known object found so far. These new data allow us to estimate the first meaningful lower limit of the star formation rate density beyond redshift 6; $\rho_{\rm SFR} \sim 5.2 \times 10^{-4} M_\odot$ yr$^{-1}$ Mpc$^{-3}$. Since it is expected that the actual density is higher by a factor of several than this value, our new observation reveals that a moderately high level of star formation activity already occurred at $z \sim$ 6.6.

We have carried out a deep imaging survey for H$\alpha$ emitting galaxies at $z\approx$0.24 using a narrowband filter tuned with the redshifted line. The total sky area covered is 706 arcmin$^2$ within a redshift range from 0.234 to 0.252 ($\delta z$=0.018). This corresponds to a volume of 3.9$\times10^3$ Mpc$^3$ when $\Omega_{\rm matter}=0.3$, $\Omega_\Lambda=0.7$, and $H_{\mathrm{0}}$=70 km s$^{-1}$ Mpc$^{-1}$ are adopted. We obtain a sample of 348 H$\alpha$ emitting galaxies whose observed emission-line equivalent widths are greater than 12 \AA. We find an extinction-corrected H$\alpha$ luminosity density of $10^{39.65^{+0.08}_{-0.12}}$ ergs s$^{-1}$ Mpc$^{-3}$. Using the Kennicutt relation between the H$\alpha$ luminosity and star formation rate, the star formation rate density in the covered volume is estimated as $0.036^{+0.006}_{-0.012}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. This value is higher by a factor of 3 than the local SFR density.

We report the discovery of a large-scale structure of Lyman alpha emitters (LAEs) at z=4.86 based on wide-field imaging with the prime-focus camera (Suprime-Cam) on the Subaru telescope. We observed a 25' x 45' area of the Subaru Deep Field in a narrow band (NB, lambdaC=7126 A and FWHM=73 A) together with R and i'. We isolate from these data 43 LAE candidates down to NB=25.5 mag using color criteria. Follow-up spectroscopy of five candidates suggests the contamination by low-z objects to be ~ 20%. We find that the LAE candidates are clustered in an elongated region on the sky of 20 Mpc in width and 50 Mpc in length at z=4.86, which is comparable in size to present-day large-scale structures (we adopt H0=70 km/s/Mpc, Omega0=0.3, lambda0=0.7). This elongated region includes a circular region of 12 Mpc radius of higher surface overdensity (delta=2), which may be the progenitor of a cluster of galaxies. Assuming this circular region to be a sphere with a spatial overdensity of 2, we compare our observation with predictions by Cold Dark Matter models. We find that an Omega0=0.3 flat model with sigma8=0.9 predicts the number of such spheres consistent with the observed number (one sphere in our survey volume) if the bias parameter of LAEs is b \simeq 6. This value suggests that the typical mass of dark haloes hosting LAEs at z \simeq 5 is of the order of 10^12 Msolar. Such a large mass poses an interesting question about the nature of LAEs.

We report on the discovery of a very narrow-line star forming object beyond redshift of 5. Using the prime-focus camera, Suprime-Cam, on the 8.2 m Subaru telescope together with a narrow-passband filter centered at $\lambda_{\rm c}$ = 8150 \AA with passband of $\Delta\lambda$ = 120 \AA, we have obtained a very deep image of the field surrounding the quasar SDSSp J104433.04$-$012502.2 at a redshift of 5.74. Comparing this image with optical broad-band images, we have found an object with a very strong emission line. Our follow-up optical spectroscopy has revealed that this source is at a redshift of $z=5.655\pm0.002$, forming stars at a rate $\sim 13 ~ h_{0.7}^{-2} ~ M_\odot$ yr$^{-1}$. Remarkably, the velocity dispersion of Ly$\alpha$-emitting gas is only 22 km s$^{-1}$. Since a blue half of the Ly$\alpha$ emission could be absorbed by neutral hydrogen gas, perhaps in the system, a modest estimate of the velocity dispersion may be $\gtrsim$ 44 km s$^{-1}$. Together with a linear size of 7.7 $h_{0.7}^{-1}$ kpc, we estimate a lower limit of the dynamical mass of this object to be $\sim 2 \times 10^9 M_\odot$. It is thus suggested that LAE J1044$-$0123 is a star-forming dwarf galaxy (i.e., a subgalactic object or a building block) beyond redshift 5 although we cannot exclude a possibility that most Ly$\alpha$ emission is absorbed by the red damping wing of neutral intergalactic matter.

We present a summary of the new optical intermediate-band filter system for the prime-focus camera, Suprime-Cam, on the Subaru telescope at Mauna Kea Observatories. We also discuss a future plan to promote a new deep survey with this filter system (the MAHOROBA project).

In order to examine whether or not high-ionization nuclear emission-line regions (HINERs) in narrow-line regions of active galactic nuclei are dusty, we focus on two high-ionization forbidden emission lines, [Fe VII]6087 and [Ne V]3426. We perform photoionization model calculations to investigate possible dependences of the flux ratio of [Fe VII]6087/[Ne V]3426 on various gas properties, in order to investigate how useful this flux ratio to explore the dust abundances in HINERs. Based on our photoionization model calculations, we show that the observed range of the flux ratio of [Fe VII]6087/[Ne V]3426 is consistent with the dust-free models while that is hard to be explained by the dusty models. This suggests that iron is not depleted at HINERs, which implies that the HINERs are not dusty. This results is consistent with the idea that the HINERs are located closer than the dust-sublimation radius (i.e., inner radius of dusty tori) and thus can be hidden by dusty tori when seen from a edge-on view toward the tori, which has been also suggested by the AGN-type dependence of the visibility of high-ionization emission lines.

The tight correlation between the masses of central black holes and their host spheroids in nearby galaxies and active galactic nuclei suggests that black hole growth is closely related to their spheroid formation. Based on our previous work regarding such a joint evolutionary scheme and the consequential black hole to bulge mass correlation, we use the X-ray luminosity function of AGN and the cosmological evolution rate which are from ROSAT X-ray surveys to estimate the cosmic star formation history associated with the black hole growth. By the Basic assumption that the major black hole growth occurs during the luminous AGN phase, the luminosity function of AGNs as a function of redshift traces not only the accretion history of the black holes but also the cosmic star formation history of the spheroids. Although the space density of the especially luminous QSOs is very low, we show that the total amount of star formation associated with the massive black hole growth is almost the same as that of Lyman Break Galaxies detected by the current optical deep surveys. We thus argue that the optical deep surveys may miss about half of the net star formation in our Universe, This is probably due to in part significant dust extinction as well as the small field of view of previous optical surveys which cannot sample such rare events with relatively short time scale. However, the far infrared emission from the dust heated by star fomation on-going during the black hole growth could sufficiently account for the observed SCUBA number counts, and would be the probable dominating energy source of the SUBA population.

During the course of our optical deep survey program on L$\alpha$ emitters at $z \approx 5.7$ in the sky area surrounding the quasar SDSSp J104433.04$-$012502.2 at $z=5.74$, we found that a faint galaxy with $m_B$(AB) $\approx 25$ is located at \timeform{1".9} southwest of the quasar. Its broad-band color properties from $B$ to $z^\prime$ suggest that the galaxy is located at a redshift of $z \sim 1.5$ -- 2.5. This is consistent with no strong emission line in our optical spectroscopy. Since the counter image of the quasar cannot be seen in our deep optical images, the magnification factor seems not to be very high. Our modest estimate is that this quasar is gravitationally magnified by a factor of 2.

We give a review on possible starburst-AGN (active galactic nuclei) connections from high redshift to the present day. First, we give a historical review on some basic ideas related to the starburst-AGN connection published in the literature. Second, we focus our attention to the so-called Magorrian relation which is the close relationship between the nuclear black hole mass and the bulge mass, established in nearby galaxies. If the Magorrian relation is universal, we obtain an important implication that any supermassive black holes were made through successive merging processes of starburst remnants (i.e., neutron stars and stellar-sized black holes), providing a channel of the starburst-AGN connection. Third, we briefly discuss a possible new scenario for the formation of quasar nuclei at very high redshift based on an idea that the successive mergers of starburst remnants formed in subgalactic gaseous clouds.

The orbital decay and tidal disruption of a star cluster in a galaxy is studied in an analytical manner. Owing to dynamical friction, the star cluster spirals in toward the center of the galaxy. Simultaneously, the galactic tidal field strips stars from the outskirts of the star cluster. Under an assumption that the star cluster undergoes a self-similar evolution, we obtain the condition and timescale for the star cluster to reach the galaxy center before its disruption. The result is used to discuss the fate of so-called intermediate-mass black holes with >10^3 M(sun) found recently in young star clusters of starburst galaxies and also the mass function of globular clusters in galaxies.

We present the results of a survey for emission-line objects based on optical intermediate-band ($\lambda_{\rm c}$ = 5736 \AA ~ and $\Delta\lambda$ = 280 \AA) and broad-band ($B$, $V$, $R$, and $i^\prime$) observations of the Subaru/XMM-Newton Deep Field on the 8.2 m Subaru telescope with the Subaru Prime Focus Camera, Suprime-Cam. All the data were obtained during the guaranteed time observations of the Suprime-Cam instrument. The intermediate-band image covered a sky area with 10\minpoint62 $\times 12\minpoint40 \approx 132$ arcmin$^2$ in the Subaru/XMM-Newton Deep Field (Ouchi et al.). Using this image, we have found 23 emission-line sources whose observed emission-line equivalent widths are greater than 250 \AA. Their optical multicolor properties indicate that six emission-line sources are Ly$\alpha$ emitters at $z \approx$ 3.7 ($\Delta z \approx 0.22$). They are either intense starburst galaxies or active galactic nuclei like quasars at $z \approx$ 3.7. Two more emission-line sources may also be Ly$\alpha$ emitters at $z \approx$ 3.7 although their multicolor properties are marginal. Among the remaining 15 emission-line objects, eight objects appear strong emission-line galaxies at lower redshift; e.g., [O {\sc ii}] $\lambda$3727 emitters at $z \approx 0.54$, H$\beta$ at $z \approx 0.18$, or [O {\sc iii}]$\lambda$5007 emitters at $z \approx 0.15$. The remaining seven objects are unclassified because they are too faint to be detected in broad-band images. We discuss observational properties of these strong emission-line sources. In particular, our data allow us to estimate the star formation density at $z \approx 3.7$ for the first time.

Faint submillimeter sources detected with the Submillimeter Common-User Bolometer Array on the James Clerk Maxwell Telescope have faced an identification problem due to the telescope's broad beam profile. Here we propose a new method to identify such submillimeter sources with a mid-infrared image having a finer point spread function. The Infrared Space Observatory has provided a very deep 6.7 um image of the Hawaii Deep Field SSA13. All three faint 850 um sources in this field have their 6.7 um counterparts. They have been identified with interacting galaxy pairs in optical images. These pairs are also detected in the radio. Two of them are optically faint and very red (I>24, I-K>4), one of which has a hard X-ray detection with the Chandra satellite. As these observing properties are similar to those of local ultraluminous infrared galaxies, their photometric redshifts are derived based on submillimeter to mid-infrared flux ratios assuming a spectral energy distribution (SED) of Arp220. Other photometric redshifts are obtained via chi^2 minimization between the available photometry data and template SEDs. Both estimates are in the range z=1-2, in good agreement with a spectroscopic redshift and a millimetric one. The reconstructed Arp220 SEDs with these redshift estimates are consistent with all the photometry data except Chandra's hard X-ray detection. The sources would be a few times more luminous than Arp220. With an assumption that AGN contributions are negligible, it appears that extremely high star formation rates are occurring in galaxies at high redshifts with massive stellar contents already in place.

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.