This work presents the first comprehensive study of structure formation at the peak epoch of cosmic star formation over $1.4\leq z \leq 3.6$ in the COSMOS field, including the most massive high redshift galaxy proto-clusters at that era. We apply the extended COSMIC BIRTH algorithm to account for a multi-tracer and multi-survey Bayesian analysis at Lagrangian initial cosmic times. Combining the data of five different spectroscopic redshift surveys (zCOSMOS-deep, VUDS, MOSDEF, ZFIRE, and FMOS-COSMOS), we show that the corresponding unbiased primordial density fields can be inferred, if a proper survey completeness computation from the parent photometric catalogs, and a precise treatment of the non-linear and non-local evolution on the light-cone is taken into account, including (i) gravitational matter displacements, (ii) peculiar velocities, and (iii) galaxy bias. The reconstructions reveal a holistic view on the known proto-clusters in the COSMOS field and the growth of the cosmic web towards lower redshifts. The inferred distant dark matter density fields concurrently with other probes like tomographic reconstructions of the intergalactic medium will explore the interplay of gas and dark matter and are ideally suited to study structure formation at high redshifts in the light of upcoming deep surveys.

Using the Horizon-AGN hydrodynamical simulation and self-organising maps (SOMs), we show how to compress the complex data structure of a cosmological simulation into a 2-d grid which is much easier to analyse. We first verify the tight correlation between the observed 0.3$\!-\!5\mu$m broad-band colours of Horizon-AGN galaxies and their high-resolution spectra. The correlation is found to extend to physical properties such as redshift, stellar mass, and star formation rate (SFR). This direct mapping from colour to physical parameter space is shown to work also after including photometric uncertainties that mimic the COSMOS survey. We then label the SOM grid with a simulated calibration sample and estimate redshift and SFR for COSMOS-like galaxies up to $z\sim3$. In comparison to state-of-the-art techniques based on synthetic templates, our method is comparable in performance but less biased at estimating redshifts, and significantly better at predicting SFRs. In particular our "data-driven" approach, in contrast to model libraries, intrinsically allows for the complexity of galaxy formation and can handle sample biases. We advocate that obtaining the calibration for this method should be one of the goals of next-generation galaxy surveys.

We present evidence that a region of high effective Ly$\alpha$ optical depth at $z\sim5.7$ is associated with an underdense region at the tail end of cosmic reionization. We carried out a survey of Lyman-break Galaxies (LBGs) using Subaru Hyper Suprime-Cam in the field of the $z=5.98$ quasar J0148+0600, whose spectrum presents an unusually long ($\sim160 \mathrm{cMpc}$) and opaque ($\tau \gtrsim7$) Ly$\alpha$ trough at $5.5\le z\le 5.9$. LBG candidates were selected to lie within the redshift range of the trough, and the projected number densities were measured within 90~cMpc of the quasar sightline. The region within $8'$ (or $\approx 19~\mathrm{cMpc}$) of the quasar position is the most underdense of the whole field. The significance of the presence of the void is estimated to be $99\%$. This is consistent with the significant deficit of Ly$\alpha$ emitters (LAEs) at $z=5.72$ reported by Becker et al. and suggests that the paucity of LAEs is not purely due to the removal of the Ly$\alpha$ emission by the high opacity but reflects a real coherent underdensity of galaxies across the entire redshift range of the trough. These observations are consistent with scenarios in which large optical depth fluctuations arise due to fluctuations in the galaxy-dominant UV background or due to residual neutral islands that are expected from reionization that is completed at redshifts as low as $z\lesssim5.5$.

We present an observed relationship between Eddington ratio ({\lambda}Edd) and optical narrow-emission-line ratio ([N II]{\lambda}6583/H{\alpha}) of X-ray-selected broad-line active galactic nuclei (AGN) at 0.6 < z < 1.7. We use 27 near-infrared spectra from the Fiber Multi-Object Spectrograph along with 26 sources from the literature. We show that the {\lambda}Edd and [N II]{\lambda}6583/H{\alpha} ratio at 0.6 < z < 1.7 exhibits a similar anti-correlation distribution of {\lambda}Edd-[N II]{\lambda}6583/H{\alpha} as has been found for local (<z> = 0.036), ultra-hard X-ray selected AGN. The observed distribution implies that there is a consistent relationship from local to z \sim 1.7 which corresponds from the present time to 4 Gyr old. Further study of high redshift low Eddington ratio AGN (log{\lambda}Edd < -2) is necessary to determine fully whether the {\lambda}Edd-[N II]{\lambda}6583/H{\alpha} anti-correlation still holds in high-redshift AGN at low Eddington ratios.

Using the ligthcone from the cosmological hydrodynamical simulation Horizon-AGN, we produced a photometric catalogue over $0<z<4$ with apparent magnitudes in COSMOS, DES, LSST-like, and Euclid-like filters at depths comparable to these surveys. The virtual photometry accounts for the complex star formation history and metal enrichment of Horizon-AGN galaxies, and consistently includes magnitude errors, dust attenuation and absorption by inter-galactic medium. The COSMOS-like photometry is fitted in the same configuration as the COSMOS2015 catalogue. We then quantify random and systematic errors of photometric redshifts, stellar masses, and star-formation rates (SFR). Photometric redshifts and redshift errors capture the same dependencies on magnitude and redshift as found in COSMOS2015, excluding the impact of source extraction. COSMOS-like stellar masses are well recovered with a dispersion typically lower than 0.1 dex. The simple star formation histories and metallicities of the templates induce a systematic underestimation of stellar masses at $z<1.5$ by at most 0.12 dex. SFR estimates exhibit a dust-induced bimodality combined with a larger scatter (typically between 0.2 and 0.6 dex). We also use our mock catalogue to predict photometric redshifts and stellar masses in future imaging surveys. We stress that adding Euclid near-infrared photometry to the LSST-like baseline improves redshift accuracy especially at the faint end and decreases the outlier fraction by a factor $\sim$2. It also considerably improves stellar masses, reducing the scatter up to a factor 3. It would therefore be mutually beneficial for LSST and Euclid to work in synergy.

We present an analysis to disentangle the connection between physical quantities that characterize the conditions of ionized HII regions -- metallicity ($Z$), ionization parameter ($U$), and electron density ($n_\mathrm{e}$) -- and the global stellar mass ($M_\ast$) and specific star formation rate ($\mathrm{sSFR}=\mathrm{SFR}/M_\ast$) of the host galaxies. We construct composite spectra of galaxies at $0.027 \le z \le 0.25$ from Sloan Digital Sky Survey, separating the sample into bins of $M_\ast$ and sSFR, and estimate the nebular conditions from the emission line flux ratios. Specially, metallicity is estimated from the direct method based on the faint auroral lines [OIII]$\lambda$4363 and [OII]$\lambda\lambda$7320,7330. The metallicity estimates cover a wide range from $12+\log\mathrm{O/H}\sim7.6\textrm{--}8.9$. It is found that these three nebular parameters all are tightly correlated with the location in the $M_\ast$--sSFR plane. With simple physically-motivated ans\"atze, we derive scaling relations between these physical quantities by performing multi regression analysis. In particular, we find that $U$ is primarily controlled by sSFR, as $U \propto \mathrm{sSFR}^{0.43}$, but also depends significantly on both $Z$ and $n_\mathrm{e}$. The derived partial dependence of $U \propto Z^{-0.36}$ is weaker than the apparent correlation ($U\propto Z^{-1.52}$). The remaining negative dependence of $U$ on $n_\mathrm{e}$ is found to be $U \propto n_\mathrm{e}^{-0.29}$. The scaling relations we derived are in agreement with predictions from theoretical models and observations of each aspect of the link between these quantities. Our results provide a useful set of equations to predict the nebular conditions and emission-line fluxes of galaxies in semi-analytic models.

We present a new data release from the Fiber Multi-Object Spectrograph (FMOS)-COSMOS survey, which contains the measurements of spectroscopic redshift and flux of rest-frame optical emission lines (H$\alpha$, [NII], [SII], H$\beta$, [OIII]) for 1931 galaxies out of a total of 5484 objects observed over the 1.7 deg$^2$ COSMOS field. We obtained $H$-band and $J$-band medium-resolution ($R\sim3000$) spectra with FMOS mounted on the Subaru telescope, which offers an in-fiber line flux sensitivity limit of $\sim 1 \times 10^{-17}~\mathrm{erg~s^{-1}~cm^{-2}}$ for an on-source exposure time of five hours. The full sample contains the main population of star-forming galaxies at $z\sim1.6$ over the stellar mass range $10^{9.5}\lesssim M_\ast/M_\odot \lesssim 10^{11.5}$, as well as other subsamples of infrared-luminous galaxies detected by Spitzer and Herschel at the same and lower ($z\sim0.9$) redshifts and X-ray emitting galaxies detected by Chandra. This paper presents an overview of our spectral analyses, a description of the sample characteristics, and a summary of the basic properties of emission-line galaxies. We use the larger sample to re-define the stellar mass--star formation rate relation based on the dust-corrected H$\alpha$ luminosity, and find that the individual galaxies are better fit with a parametrization including a bending feature at $M_\ast\approx10^{10.2}~M_\odot$, and that the intrinsic scatter increases with $M_\ast$ from 0.19 to $0.37$ dex. We also confirm with higher confidence that the massive ($M_\ast\gtrsim10^{10.5}~M_\odot$) galaxies are chemically mature as much as local galaxies with the same stellar masses, and that the massive galaxies have lower [SII]/H$\alpha$ ratios for their [OIII]/H$\beta$, as compared to local galaxies, which is indicative of enhancement in ionization parameter.

We present near-IR spectroscopy in J- and H-band for a large sample of 243 X-ray selected moderate-luminosity type-1 AGN in the COSMOS, SXDS and E-CDF-S survey fields using the multi-object spectrograph Subaru/FMOS. Our sample covers the redshift range 0.5<z<3.0 and an X-ray luminosity range of $10^{43}\lesssim L_X \lesssim 10^{45}$~erg s$^{-1}$. We provide emission-line properties and derived virial black hole mass estimates, bolometric luminosities and Eddington ratios, based on H$\alpha$ (211), H$\beta$ (63) and MgII (4). We compare line widths, luminosities and black hole mass estimates from H$\alpha$ and H$\beta$ and augment these with commensurate measurements of MgII and CIV detected in optical spectra. We demonstrate the robustness of using H$\alpha$, H$\beta$ and MgII as reliable black hole mass estimators for high-z moderate-luminosity AGN, while the use of CIV is prone to large uncertainties (>0.4 dex). We extend a recently proposed correction based on the CIV blueshift to lower luminosities and black hole masses. While our sample shows an improvement in their CIV black hole mass estimates, the deficit of high blueshift sources reduces its overall importance for moderate-luminosity AGN, compared to the most luminous quasars. In addition, we revisit luminosity correlations between $L_{\rm{bol}}$, $L_X$, $L_{\rm{[OIII]}}$, $L_{5100}$ and $L_{\rm{H}\alpha}$ and find them to be consistent with a simple empirical model, based on a small number of well-established scaling relations. We finally highlight our highest redshift AGN, CID 781, at z=4.6 which shows the lowest black hole mass ($\sim10^8$M$_\odot$) among current near-IR samples at this redshift, and is in a state of fast growth.

We present an analysis of the molecular gas properties, based on CO(2 - 1) emission, of twelve starburst galaxies at z~1.6 selected by having a boost (>~4x) in their star formation rate (SFR) above the average star-forming galaxy at an equivalent stellar mass. ALMA observations are acquired of six additional galaxies than previously reported through our effort. As a result of the larger statistical sample, we significantly detect, for the first time at high-z, a systematically lower L'_CO/L_IR ratio in galaxies lying above the star-forming `main sequence' (MS). Based on an estimate of alpha_CO (i.e., the ratio of molecular gas mass to L'_CO(1-0)), we convert the observational quantities (e.g., L'_CO/L_IR) to physical units (M_gas/SFR) that represent the gas depletion time or its inverse, the star formation efficiency. We interpret the results as indicative of the star formation efficiency increasing in a continuous fashion from the MS to the starburst regime, whereas the gas fractions remain comparable to those of MS galaxies. Although, the balance between an increase in star-formation efficiency or gas fraction depends on the adopted value of alpha_CO as discussed.

We report on the discovery of a merger-driven starburst at z = 1.52, PACS-787, based on high signal-to-noise ALMA observations. CO(5-4) and continuum emission (850um) at a spatial resolution of 0.3" reveal two compact (r_1/2 ~ 1 kpc) and interacting molecular gas disks at a separation of 8.6 kpc thus indicative of an early stage in a merger. With a SFR of 991 Msun/yr, this starburst event should occur closer to final coalescence, as usually seen in hydrodynamical simulations. From the CO size, inclination, and velocity profile for both disks, the dynamical mass is calculated through a novel method that incorporates a calibration using simulations of galaxy mergers. Based on the dynamical mass, we measure (1) the molecular gas mass, independent from the CO luminosity, (2) the ratio of the total gas mass and the CO(1 - 0) luminosity (alpha_CO = M_gas/L'_CO(1-0)), and (3) the gas-to-dust ratio, with the latter two being lower than typically assumed. We find that the high star formation, triggered in both galaxies, is caused by a set of optimal conditions: a high gas mass/fraction, a short depletion time (t_depl=85 and 67 Myrs) to convert gas into stars, and the interaction of likely counter-rotating molecular disks that may accelerate the loss of angular momentum. The state of interaction is further established by the detection of diffuse CO and continuum emission, tidal debris that bridges the two nuclei and is associated with stellar emission seen by HST/WFC3. This observation demonstrates the power of ALMA to study the dynamics of galaxy mergers at high redshift.

We investigate the evolution of the ionisation parameter of star-forming galaxies using a high-redshift ($z\sim 1.5$) sample from the FMOS-COSMOS survey and matched low-redshift samples from the Sloan Digital Sky Survey. By constructing samples of low-redshift galaxies for which the stellar mass ($\mathrm{M}_{*}$), star formation rate (SFR) and specific star formation rate (sSFR) are matched to the high-redshift sample we remove the effects of an evolution in these properties. We also account for the effect of metallicity by jointly constraining the metallicity and ionisation parameter of each sample. We find an evolution in the ionisation parameter for main-sequence, star-forming galaxies and show that this evolution is driven by the evolution of sSFR. By analysing the matched samples as well as a larger sample of $z<0.3$, star-forming galaxies we show that high ionisation parameters are directly linked to high sSFRs and are not simply the byproduct of an evolution in metallicity. Our results are physically consistent with the definition of the ionisation parameter, a measure of the hydrogen ionising photon flux relative to the number density of hydrogen atoms.

The Lyman continuum (LyC) escape fraction is a key quantity to determine the contribution of galaxies to cosmic reionization. It has been known that the escape fractions estimated by observations and numerical simulations show a large diversity. However, the origins of the diversity are still uncertain. In this work, to understand what quantities of galaxies are responsible for controlling the escape fraction, we numerically evaluate the escape fraction by performing ray-tracing calculation with simplified disc galaxy models. With a smooth disc model, we explore the dependence of the escape fraction on the disposition of ionizing sources, and find that the escape fraction varies up to $\sim 3$ orders of magnitude. It is also found that the halo mass dependence of disc scale height determines whether the escape fraction increases or decreases with halo mass. With a clumpy disc model, it turns out that the escape fraction increases as the clump mass fraction increases because the density in the inter-clump region decreases. In addition, we find that clumpiness regulates the escape fraction via two ways when the total clump mass dominates the total gas mass; the escape fraction is controlled by the covering factor of clumps if the clumps are dense sufficient to block LyC photons, otherwise the clumpiness works to reduce the escape fraction by increasing the total number of recombination events in a galaxy.

We estimate the number counts of line emitters at high redshift and their evolution with cosmic time based on a combination of photometry and spectroscopy. We predict the H$\alpha$, H$\beta$, [OII], and [OIII] line fluxes for more than $35,000$ galaxies down to stellar masses of $\sim10^9$ $M_{\odot}$ in the COSMOS and GOODS-S fields, applying standard conversions and exploiting the spectroscopic coverage of the FMOS-COSMOS survey at $z\sim1.55$ to calibrate the predictions. We calculate the number counts of H$\alpha$, [OII], and [OIII] emitters down to fluxes of $1\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the range $1.4 < z < 1.8$ covered by the FMOS-COSMOS survey. We model the time evolution of the differential and cumulative H$\alpha$ counts, steeply declining at the brightest fluxes. We expect $\sim9,300-9,700$ and $\sim2,300-2,900$ galaxies deg$^{-2}$ for fluxes $\geq1\times10^{-16}$ and $\geq2\times10^{-16}$ erg cm$^{-2}$ s$^{-1}$ over the range $0.9<z<1.8$. We show that the observed evolution of the Main Sequence of galaxies with redshift is enough to reproduce the observed counts variation at $0.2<z<2.5$. We characterize the physical properties of the H$\alpha$ emitters with fluxes $\geq2\times10^{-16}$ erg cm$^{-2}$ s$^{-1}$, including their stellar masses, UV sizes, [NII]/H$\alpha$ ratios, and H$\alpha$ equivalent widths. An aperture of $R\sim R_{\rm e}\sim0.5$" maximizes the signal-to-noise ratio for a detection, while causing a factor of $\sim2\times$ flux losses, influencing the recoverable number counts, if neglected. Our approach, based on deep and large photometric datasets, reduces the uncertainties on the number counts due to the selection and spectroscopic samplings, while exploring low fluxes. We publicly release the line flux predictions for the explored photometric samples.

We study the properties of dark matter halos that contain star-forming galaxies at $1.43 \le z \le 1.74$ using the FMOS-COSMOS survey. The sample consists of 516 objects with a detection of the H$\alpha$ emission line, that represent the star-forming population at this epoch having a stellar mass range of $10^{9.57}\le M_\ast/M_\odot \lesssim 10^{11.4}$ and a star formation rate range of $15\lesssim \mathrm{SFR}/(M_\odot \mathrm{yr^{-1}}) \lesssim 600$. We measure the projected two-point correlation function while carefully taking into account observational biases, and find a significant clustering amplitude at scales of $0.04$-$10~h^{-1}~\mathrm{cMpc}$, with a correlation length $r_0 = 5.21^{+0.70}_{-0.67}~h^{-1}~\mathrm{cMpc}$ and a bias $b=2.59^{+0.41}_{-0.34}$. We interpret our clustering measurement using a halo occupation distribution model. The sample galaxies appear to reside in halos with mass $M_\mathrm{h} = 4.6^{+1.1}_{-1.6}\times10^{12}~h^{-1}M_\odot$ on average that will likely become present-day halos of mass $M_\mathrm{h} (z=0) \sim2\times10^{13}~h^{-1}M_\odot$, equivalent to the typical halo mass scale of galaxy groups. We then confirm the decline of the stellar-to-halo mass ratio at $M_\mathrm{h}<10^{12}~M_\odot$, finding $M_\ast/M_\mathrm{h} \approx 5\times10^{-3}$ at $M_\mathrm{h}=10^{11.86}~M_\odot$, which is lower by a factor of 2-4 than those measured at higher masses. Finally, we use our results to illustrate the future capabilities of Subaru's Prime-Focus Spectrograph, a next-generation instrument that will provide strong constraints on the galaxy-formation scenario by obtaining precise measurements of galaxy clustering at $z>1$.

We present rest-frame optical spectra from the FMOS-COSMOS survey of twelve $z \sim 1.6$ \textit{Herschel} starburst galaxies, with Star Formation Rate (SFR) elevated by $\times$8, on average, above the star-forming Main Sequence (MS). Comparing the H$\alpha$ to IR luminosity ratio and the Balmer Decrement we find that the optically-thin regions of the sources contain on average only $\sim 10$ percent of the total SFR whereas $\sim90$ percent comes from an extremely obscured component which is revealed only by far-IR observations and is optically-thick even in H$\alpha$. We measure the [NII]$_{6583}$/H$\alpha$ ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity-SFR anticorrelation observed at fixed stellar mass for the MS population. The [SII]$_{6732}$/[SII]$_{6717}$ ratio from the average spectrum indicates an electron density $n_{\rm e} \sim 1,100\ \mathrm{cm}^{-3}$, larger than what estimated for MS galaxies but only at the 1.5$\sigma$ level. Our results provide supporting evidence that high-$z$ MS outliers are the analogous of local ULIRGs, and are consistent with a major merger origin for the starburst event.

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.

We use $\sim$83,000 star-forming galaxies at $0.04<z<0.3$ from the Sloan Digital Sky Survey to study the so-called fundamental metallicity relation (FMR) and report on the disappearance of its anti-correlation between metallicity and star formation rate (SFR) when using the new metallicity indicator recently proposed by Dopita et al. In this calibration, metallicity is primarily sensitive to the emission line ratio [NII]$\lambda$6584 / [SII]$\lambda\lambda$6717, 6731 that is insensitive to dilution by pristine infalling gas that may drive the FMR anti-correlation with SFR. Therefore, we conclude that the apparent disappearance of the FMR (using this new metallicity indicator) does not rule out its existence.

Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star-forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure the molecular gas content of seven high-redshift (z ~ 1.6) starburst galaxies with the Atacama Large (sub-)Millimeter Array and IRAM/Plateau de Bure Interferometer. Our targets are selected from the sample of Herschel far-infrared detected galaxies having star formation rates (~300-800 Msolar/yr) elevated (>4x) above the star-forming `main sequence', and included in the FMOS-COSMOS near-infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. We detect CO emission in all cases at high levels of significance, indicative of high gas fractions (~30-50%). Even more compelling, we firmly establish with a clean and systematic selection that starbursts, identified as main-sequence outliers, at high redshift generally have a lower ratio of CO to total infrared luminosity as compared to typical main-sequence star-forming galaxies, although with a smaller offset than expected based on past studies of local starbursts. We put forward a hypothesis that there exists a continuous increase in star formation efficiency with elevation from the main sequence with galaxy mergers as a possible physical driver. Along with a heightened star formation efficiency, our high-redshift sample is similar in other respects to local starbursts such as being metal rich and having a higher ionization state of the interstellar medium.

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 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 examine the mass-metallicity relation for $z\lesssim 1.6$. The mass-metallicity relation follows a steep slope with a turnover or `knee' at stellar masses around $10^{10} M_\odot$. At stellar masses higher than the characteristic turnover mass, the mass-metallicity relation flattens as metallicities begin to saturate. We show that the redshift evolution of the mass-metallicity relation depends only on evolution of the characteristic turnover mass. The relationship between metallicity and the stellar mass normalized to the characteristic turnover mass is independent of redshift. We find that the redshift independent slope of the mass-metallicity relation is set by the slope of the relationship between gas mass and stellar mass. The turnover in the mass-metallicity relation occurs when the gas-phase oxygen abundance is high enough that the amount of oxygen locked up in low mass stars is an appreciable fraction of the amount of oxygen produced by massive stars. The characteristic turnover mass is the stellar mass where the stellar-to-gas mass ratio is unity. Numerical modeling suggests that the relationship between metallicity and stellar-to-gas mass ratio is a redshift independent, universal relationship followed by all galaxies as they evolve. The mass-metallicity relation originates from this more fundamental universal relationship between metallicity and stellar-to-gas mass ratio. We test the validity of this universal metallicity relation in local galaxies where stellar mass, metallicity and gas mass measurements are available. The data are consistent with a universal metallicity relation. We derive an equation for estimating the hydrogen gas mass from measurements of stellar mass and metallicity valid for $z\lesssim1.6$ and predict the cosmological evolution of galactic gas masses.

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 compare various star formation rate (SFR) indicators for star-forming galaxies at $1.4<z<2.5$ in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR-stellar mass ($M_*$) relation, compared to that of the dominant main sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR-{$M_*$} relation is $\sim0.8-0.9$, in agreement with that derived from the dust-corrected UV-luminosity. Exploiting deeper 24$\mu$m-Spitzer data we have characterized a sub-sample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the $B$ band. The combination of Herschel with Spitzer data have allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the H$\alpha$ luminosity of a sample of star-forming galaxies at $1.4<z<1.7$. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.

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.

The standard models of inflation predict statistically homogeneous and isotropic primordial fluc- tuations, which should be tested by observations. In this paper we test the statistical isotropy of the mean of the CMB temperature fluctuations measured by the WMAP in the spherical harmonic space. A classical method to test a mean, like the simple student's-t test, is not appropriate for this purpose because the WMAP data contain anisotropic instrumental noise and suffer from the effect of the mask for the foreground emissions which breaks the statistical independence. Here we perform a band-power analysis with Monte-Carlo simulations in which we take into account the anisotropic noise and the mask. We find evidence of a non-zero mean at 99.93 % confidence level in a particular range of multipoles. The evidence against the zero-mean assumption as a whole is still significant at the 99 % confidence level even if the fact is taken into account that we have tested multiple ranges.