Deep observations have detected extended Ly$\alpha$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the sample, while CIII] is only marginally detected in a few objects. By stacking the cubes we detect emission of CIV, HeII and CIII] out to a radius of about 45 kpc. CIV and HeII show a radial decline much steeper than Ly$\alpha$, while CIII] shows a shallower profile similar to Ly$\alpha$ in the inner 45 kpc. We infer that the average metallicity of the circumgalactic gas within the central 30-50~kpc is $\sim$0.5 solar, or even higher. However, we also find evidence of a component of the Ly$\alpha$ haloes, which has much weaker metal emission lines relative to Ly$\alpha$. We suggest that the high metallicity of the circumgalactic medium within the central 30-50 kpc is associated with chemical pre-enrichment by past quasar-driven outflows and that there is a more extended component of the CGM that has much lower metallicity and likely associated with near-pristine gas accreted from the intergalactic medium. We show that our observational results are in good agreement with the expectations of the FABLE zoom-in cosmological simulations.

Recent active galactic nucleus (AGN) and quasar surveys have revealed a population showing rapid AGN luminosity variability by a factor of $\sim10$. Here we present the most drastic AGN luminosity decline by a factor of $\gtrsim 10^{3}$ constrained by a NuSTAR X-ray observation of the nearby galaxy Arp 187, which is a promising "dead" quasar whose current activity seems quiet but whose past activity of $L_\mathrm{bol} \sim 10^{46}$ erg s$^{-1}$ is still observable at a large scale by its light echo. The obtained upper bound of the X-ray luminosity is $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.2$, corresponding to $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 42.5$, indicating an inactive central engine. Even if a putative torus model with $N_\mathrm{H} \sim 1.5 \times 10^{24}$ cm$^{-2}$ is assumed, the strong upper-bound still holds with $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.8$ or $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 43.1$. Given the expected size of the narrow line region, this luminosity decrease by a factor of $\gtrsim 10^3$ must have occurred within $\lesssim 10^4$ yr. This extremely rapid luminosity/accretion shutdown is puzzling and it requires one burst-like accretion mechanism producing a clear outer boundary for an accretion disk. We raise two possible scenarios realizing such an accretion mechanism: a mass accretion 1) by the tidal disruption of a molecular cloud and/or 2) by the gas depletion as a result of vigorous nuclear starformation after rapid mass inflow to the central engine.

We explore X-ray evidence for the presence of active galactic nuclei (AGNs) in the two most actively star-forming Green Pea galaxies (GPs), SDSS J0749+3337 and SDSS J0822+2241, which have star-formation rates (SFRs) of $123~M_\odot$ yr$^{-1}$ and $78~M_\odot$ yr$^{-1}$, respectively. The GPs have red mid-infrared (MIR) spectral energy distributions and higher 22 $\mu$m luminosities than expected from a proxy of the SFR (H$\alpha$ luminosity), consistent with hosting AGNs with 2-10 keV luminosities of $\sim10^{44}$ erg s$^{-1}$. We thus obtain and analyze the first hard ($>$ 10 keV) X-ray data observed with NuSTAR and archival XMM-Newton data below 10 keV. From the NuSTAR $\approx$20 ksec data, however, we find no significant hard X-ray emission. By contrast, soft X-ray emission with 0.5--8 keV luminosities of $\approx10^{42}$ erg s$^{-1}$ is significantly detected in both targets, which can be explained only by star formation (SF). A possible reason for the lack of clear evidence is that a putative AGN torus absorbs most of the X-ray emission. Applying a smooth-density AGN torus model, we determine minimum hydrogen column densities along the equatorial plane ($N_{\rm H}^{\rm eq}$) consistent with the non-detection. The results indicate $N_{\rm H}^{\rm eq} \gtrsim 2\times10^{24}$ cm$^{-2}$ for SDSS J0749+3337 and $N_{\rm H}^{\rm eq} \gtrsim 5\times10^{24}$ cm$^{-2}$ for SDSS J0822+2241. Therefore, the GPs may host such heavily obscured AGNs. Otherwise, no AGN exists and the MIR emission is ascribed to SF. Active SF in low-mass galaxies is indeed suggested to reproduce red MIR colors. This would imply that diagnostics based on MIR photometry data alone may misidentify such galaxies as AGNs.

We revisit the correlation between the mid-infrared (6 $\mu$m) and hard X-ray (2--10 keV) luminosities of active galactic nuclei (AGNs) to understand the physics behind it. We construct an X-ray flux-limited sample of 571 type 1 AGNs with $f_{0.5-2.0 \,{\rm keV}} > 2.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, drawn from the ROSAT Bright Survey catalog. Cross-matching the sample with infrared data taken from Wide-field Infrared Survey Explorer, we investigate the relation between the rest-frame 6 $\mu$m luminosity ($L_{\rm 6}$) and the rest-frame 2--10 keV luminosity ($L_{\rm X}$), where $L_{\rm 6}$ is corrected for the contamination of host galaxies by using the spectral energy distribution fitting technique. We confirm that $L_{\rm 6}$ and $L_{\rm X}$ are correlated over four orders of magnitude, in the range of $L_{\rm X} = 10^{42-46}$ erg s$^{-1}$. We investigate what kinds of physical parameters regulate this correlation. We find that $L_{\rm X}$/$L_{\rm 6}$ clearly depends on the Eddington ratio ($\lambda_{\rm Edd}$) as $\log \lambda_{\rm Edd} = -(0.56 \pm 0.10) \log \, (L_{\rm X}/L_{\rm 6}) - (1.07 \pm 0.05)$, even taking into account quasars that are undetected by ROSAT as well as those detected by XMM-Newton in the literature. We also add hyper-luminous quasars with $L_{\rm 6}$ $>$ 10$^{46}$ erg s$^{-1}$ in the literature and perform a correlation analysis. The resultant correlation coefficient is $-0.41 \pm 0.07$, indicating a moderately tight correlation between $L_{\rm X}$/$L_{\rm 6}$ and $\lambda_{\rm Edd}$. This means that AGNs with high Eddington ratios tend to have lower X-ray luminosities with respect to the mid-infrared luminosities. This dependence can be interpreted as a change in the structure of the accretion flow.

We present a near-infrared (NIR) spectrum of WISE J1042+1641, an extremely red dust-obscured galaxy (DOG), which has been observed with the LIRIS on the 4.2m William Hershel Telescope. This object was selected as a hyper-luminous DOG candidate at z ~ 2 by combining the optical and IR photometric data based on the SDSS and WISE, although its redshift had not yet been confirmed. Based on the LIRIS observation, we confirmed its redshift of 2.521 and total IR luminosity of log(L_IR/L_sun) = 14.57, which satisfies the criterion for an extremely luminous IR galaxy (ELIRG). Moreover, we indicate that this object seems to have an extremely massive black hole with M_BH = 10^10.92 M_sun based on the broad Halpha line: the host stellar mass is derived as M_star = 10^13.55 M_sun by a fit of the spectral energy distribution. Very recently, it has been reported that this object is an anomalous gravitationally lensed quasar based on near-IR high-resolution imaging data obtained with the Hubble Space Telescope. Its magnification factor has also been estimated with some uncertainty (i.e., mu = 53-122). We investigate the ratio of the black hole to galaxy mass, which is less strongly affected by a lensing magnification factor, instead of the absolute values of the luminosities and masses. We find that the M_BH/M_star ratio (i.e., 0.0140-0.0204) is significantly higher than the local relation, following a sequence of unobscured quasars instead of obscured objects (e.g., submillimeter galaxies) at the same redshift. Moreover, the LIRIS spectrum shows strongly blueshifted oxygen lines with an outflowing velocity of ~ 1100 km/s, and our Swift X-ray observation also supports that this source is an absorbed AGN with an intrinsic column density of N_H = 4.9 x 10^23 cm^-2. These results imply that WISE J1042+1641 is in a blow-out phase at the end of the buried rapid black hole growth.

Arp 187 is one of the fading active galactic nuclei (AGN), whose AGN activity is currently decreasing in luminosity. We investigate the observational signatures of AGN in Arp 187, which trace various physical scales from less than 0.1 pc to the nearly 10 kpc, to estimate the longterm luminosity change over $10^{4}$ years. The VLA 5 GHz, 8 GHz, and the ALMA 133 GHz images reveal bimodal jet lobes with $\sim$5 kpc size and the absence of the central radio-core. The 6dF optical spectrum shows that Arp 187 hosts narrow line region with the estimated size of $\sim$1 kpc, and the line strengths give the AGN luminosity of $L_{\rm bol}=1.5 \times 10^{46}$ erg s$^{-1}$. On the other hand, the current AGN activity estimated from the AGN torus emission gives the upper bound of $L_{\rm bol} < 2.2 \times 10^{43}$ erg s$^{-1}$. The absence of the radio-core gives the more strict upper bound of the current AGN luminosity of $L_{\rm bol} < 8.0 \times 10^{40}$ erg s$^{-1}$, suggesting that the central engine is already quenched. These multi-wavelength signatures indicate that Arp 187 hosts a "dying" AGN: the central engine is already dead, but the large scale AGN indicators are still observable as the remnant of the past AGN activity. The central engine has experienced the drastic luminosity decline by a factor of $\sim10^{3-5}$ fainter within $\sim10^{4}$ years, which is roughly consistent with the viscous timescale of the inner part of the accretion disk within $\sim$500 years.

Ionized outflows, revealed by broad asymmetric wings of the [OIII] line, are commonly observed in AGN but the low intrinsic spatial resolution of observations has generally prevented a detailed characterization of their properties. The MAGNUM survey aims at overcoming these limitations by focusing on the nearest AGN, including NGC 1365, a nearby Seyfert galaxy (D~17 Mpc), hosting a low-luminosity AGN (Lbol ~ 2x10^43 erg/s). We want to obtain a detailed picture of the ionized gas in the central ~5 kpc of NGC 1365 in terms of physical properties, kinematics, and ionization mechanisms. We also aim to characterize the warm ionized outflow as a function of distance from the nucleus and its relation with the nuclear X-ray wind. We employed VLT/MUSE optical integral field spectroscopic observations to investigate the warm ionized gas and Chandra ACIS-S X-ray data for the hot highly-ionized phase. We obtained flux, kinematic, and diagnostic maps of the optical emission lines, which we used to disentangle outflows from disk motions and measure the gas properties down to a spatial resolution of ~70 pc. [OIII] emission mostly traces an AGN-ionized kpc-scale biconical outflow with velocities up to ~200 km/s. H{\alpha} emission traces instead star formation in a circumnuclear ring and along the bar, where we detect non-circular motions. Soft X-rays are mostly due to thermal emission from the star-forming regions, but we could isolate the AGN photoionized component which matches the [OIII] emission. The mass outflow rate of the extended ionized outflow matches that of the nuclear X-ray wind and then decreases with radius. However, the hard X-ray emission from the circumnuclear ring suggests that star formation might contribute to the outflow. The integrated mass outflow rate, kinetic energy rate, and outflow velocity are broadly consistent with the typical relations observed in more luminous AGN.

The mass-metallicity relation (MZR) of type-2 active galactic nuclei (AGNs) at 1.2 < z < 4.0 is investigated by using high-z radio galaxies (HzRGs) and X-ray selected radio-quiet AGNs. We combine new rest-frame ultraviolet (UV) spectra of two radio-quiet type-2 AGNs obtained with FOCAS on the Subaru Telescope with existing rest-frame UV emission lines, i.e., CIV1549, HeII1640, and CIII]1909, of a sample of 16 HzRGs and 6 additional X-ray selected type-2 AGNs, whose host stellar masses have been estimated in literature. We divided our sample in three stellar mass bins and calculated averaged emission-line flux ratios of CIV1549/HeII1640 and CIII]1909/CIV1549. Comparing observed emission-line flux ratios with photoionization model predictions, we estimated narrow line region (NLR) metallicities for each mass bin. We found that there is a positive correlation between NLR metallicities and stellar masses of type-2 AGNs at z ~ 3. This is the first indication that AGN metallicities are related to their hosts, i.e., stellar mass. Since NLR metallicities and stellar masses follow a similar relation as the MZR in star-forming galaxies at similar redshifts, our results indicate that NLR metallicities are related to those of the host galaxies. This study highlights the importance of considering lower-mass X-ray selected AGNs in addition to radio galaxies to explore the metallicity properties of NLRs at high redshift.

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

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 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 present the result of our spectroscopic follow-up observation for faint quasar candidates at z~5 in a part of the Canada-France-Hawaii Telescope Legacy Survey wide field. We select nine photometric candidates and identify three z~5 faint quasars, one z~4 faint quasar, and a late-type star. Since two faint quasar spectra show Civ emission line without suffering from a heavy atmospheric absorption, we estimate the black hole mass (M$_{BH}$) and Eddington ratio (L/L$_{Edd}$) of them. The inferred log M$_{BH}$ are 9.04+/-0.14 and 8.53+/-0.20, respectively. In addition, the inferred log (L/L$_{Edd}$) are -1.00+/-0.15 and -0.42+/-0.22, respectively. If we adopt that L/L$_{Edd}$= constant or $\propto$ (1+z)^2, the seed black hole masses (M$_{seed}$) of our z~5 faint quasars are expected to be >10^5 M$_\odot$ in most cases. We also compare the observational results with a mass accretion model where angular momentum is lost due to supernova explosions (Kawakatu & Wada 2008). Accordingly, M$_{BH}$ of the z~5 faint quasars in our sample can be explained even if M$_{seed}$ is ~10^3M$_\odot$. Since z~6 luminous qusars and our z~5 faint quasars are not on the same evolutionary track, z~6 luminous quasars and our z~5 quasars are not the same populations but different populations, due to the difference of a period of the mass supply from host galaxies. Furthermore, we confirm that one can explain M$_{BH}$ of z~6 luminous quasars and our z~5 faint quasars even if their seed black holes of them are formed at z~7.

Low-metallicity active galactic nuclei (AGNs) are interesting to study the early phase of the AGN evolution. However most AGNs are chemically matured and accordingly low-metallicity AGNs are extremely rare. One approach to search for low-metallicity AGNs systematically is utilizing the so-called BPT diagram that consists of the [O III]$\lambda$5007/H$\beta$$\lambda4861$ and [N II]$\lambda6584$/H$\alpha$$\lambda6563$ flux ratios. Specifically, photoionization models predict that low-metallicity AGNs show a high [O III]$\lambda$5007/H$\beta$$\lambda$4861 ratio and a relatively low [N II]$\lambda$6584/H$\alpha$$\lambda$6563 ratio, that corresponds to the location between the sequence of star-forming galaxies and that of usual AGNs on the BPT diagram (hereafter "the BPT valley"). However, other populations of galaxies such as star-forming galaxies and AGNs with a high electron density or a high ionization parameter could be also located in the BPT valley, not only low-metallicity AGNs. In this paper, we examine whether most of emission-line galaxies at the BPT valley are low-metallicity AGNs or not. We select 70 BPT-valley objects from 212,866 emission line galaxies obtained by the Sloan Digital Sky Survey. Among the 70 BPT-valley objects, 43 objects show firm evidence of the AGN activity; i.e., the He II$\lambda$4686 emission and/or weak but significant broad H$\alpha$ emission. Our analysis shows that those 43 BPT-valley AGNs are not characterized by a very high gas density nor ionization parameter, inferring that at least 43 among 70 BPT-valley objects (i.e., $>60$ %) are low-metallicity AGNs. This suggests that the BPT diagram is an efficient tool to search for low-metallicity AGNs.

We investigate the nature of far-infrared (70 um) and hard X-ray (3-24 keV) selected galaxies in the COSMOS field detected with both Spitzer and Nuclear Spectroscopic Telescope Array (NuSTAR). By matching the Spitzer-COSMOS catalog against the NuSTAR-COSMOS catalog, we obtain a sample consisting of a hyperluminous infrared galaxy with log(L_IR/L_sun) > 13, 12 ultraluminous infrared galaxies with 12 < log(L_IR/L_sun) < 13, and 10 luminous infrared galaxies with 11 < log(L_IR/L_sun) < 12, i.e., 23 Hy/U/LIRGs in total. Using their X-ray hardness ratios, we find that 12 sources are obscured active galactic nuclei (AGNs) with absorption column densities of N_H > 10^22 cm^-2, including several Compton-thick (N_H ~ 10^24 cm^-2) AGN candidates. On the basis of the infrared (60 um) and intrinsic X-ray luminosities, we examine the relation between star-formation (SF) and AGN luminosities of the 23 Hy/U/LIRGs. We find that the correlation is similar to that of optically-selected AGNs reported by Netzer (2009), whereas local, far-infrared selected U/LIRGs show higher SF-to-AGN luminosity ratios than the average of our sample. This result suggests that our Hy/U/LIRGs detected both with Spitzer and NuSTAR are likely situated in a transition epoch between AGN-rising and cold-gas diminishing phases in SF-AGN evolutional sequences. The nature of a Compton-thick AGN candidate newly detected above 8 keV with NuSTAR (ID 245 in Civano et al. 2015) is briefly discussed.

We systematically investigate the near- (NIR) to far-infrared (FIR) photometric properties of a nearly complete sample of local active galactic nuclei (AGN) detected in the Swift/Burst Alert Telescope (BAT) all-sky ultra hard X-ray (14-195 keV) survey. Out of 606 non-blazar AGN in the Swift/BAT 70-month catalog at high galactic latitude of $|b|>10^{\circ}$, we obtain IR photometric data of 604 objects by cross-matching the AGN positions with catalogs from the WISE, AKARI, IRAS, and Herschel infrared observatories. We find a good correlation between the ultra-hard X-ray and mid-IR (MIR) luminosities over five orders of magnitude ($41 < \log (L_{14-195}/{\rm erg}~{\rm s}^{-1})< 46$). Informed by previous measures of the intrinsic spectral energy distribution of AGN, we find FIR pure-AGN candidates whose FIR emission is thought to be AGN-dominated with low starformation activity. We demonstrate that the dust covering factor decreases with the bolometric AGN luminosity, confirming the luminosity-dependent unified scheme. We also show that the completeness of the WISE color-color cut in selecting Swift/BAT AGN increases strongly with 14-195 keV luminosity.

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.

The ability to accurately derive black hole (BH) masses at progressively higher redshifts and over a wide range of continuum luminosities has become indispensable in the era of large-area extragalactic spectroscopic surveys. In this paper we present an extension of existing comparisons between rest-frame UV and optical virial BH mass estimators to intermediate redshifts and luminosities comparable to the local H$\beta$ reverberation mapped active galactic nuclei (AGN). We focus on the MgII, CIV, and CIII] broad emission lines and compare them to both H$\alpha$ and H$\beta$. We use newly acquired near-infrared spectra from the FMOS instrument on the Subaru telescope for 89 broad-lined AGN at redshifts between 0.3 and 3.5, complemented by data from the AGES survey. We employ two different prescriptions for measuring the emission line widths and compare the results. We confirm that MgII shows a tight correlation with H$\alpha$ and H$\beta$, with a scatter of ~0.25 dex. The CIV and CIII] estimators, while showing larger scatter, are viable virial mass estimators after accounting for a trend with the UV-to-optical luminosity ratio. We find an intrinsic scatter of ~0.37 dex between Balmer and carbon virial estimators by combining our dataset with previous high redshift measurements. This updated comparison spans a total of 3 decades in BH mass. We calculate a virial factor for CIV/CIII] logf(CIV/CIII])=0.87 with an estimated systematic uncertainty of ~0.4 dex and find excellent agreement between the local reverberation mapped AGN sample and our high-z sample.

We investigate optical, infrared, and radio active galactic nucleus (AGN) signs in the merger remnant Arp 187, which hosts luminous jets launched in the order of $10^5$ yr ago but whose present-day AGN activity is still unknown. We find AGN signs from the optical BPT diagram and infrared [OIV]25.89 $\mu$m line, originating from the narrow line regions of AGN. On the other hand, Spitzer/IRS show the host galaxy dominated spectra, suggesting that the thermal emission from the AGN torus is considerably small or already diminished. Combining the black hole mass, the upper limit of radio luminosity of the core, and the fundamental plane of the black hole enable us to estimate X-ray luminosity, which gives $<10^{40}$ erg s$^{-1}$. Those results suggest that the AGN activity of Arp 187 has already been quenched, but the narrow line region is still alive owing to the time delay of emission from the past AGN activity.

We present the results of new infrared spectroscopic observations of 37 quasars at z~3, selected based on the optical r'-band magnitude and the availability of nearby bright stars for future imaging follow-up with Adaptive Optics system. The supermassive black hole (SMBH) masses (M_BH) were successfully estimated in 28 out of 37 observed objects from the combination of the H_beta emission linewidth and continuum luminosity at rest-frame 5100A. Comparing these results with those from previous studies of quasars with similar redshift, our sample exhibited slightly lower (~ -0.11 dex in median) Eddington ratios; and, the SMBH masses are slightly (~ 0.38 dex in median) higher. The SMBH growth time, t_grow, was calculated by dividing the estimated SMBH mass by the mass accretion rate measured using optical luminosity. We found, given reasonable assumptions, that t_grow was smaller than the age of the universe at the redshift of individual quasars for a large fraction of observed sources, suggesting that the SMBHs in many of our observed quasars are in growing phase with high accretion rates. A comparison of the SMBH masses derived from our H_beta data and archived CIV data indicated considerable scattering, as indicated in previous studies. All quasars with measured SMBH masses have at least one nearby bright star, such that they are suitable targets for adaptive optics observations to study the mass relationship between SMBHs and host galaxies' stellar component at high redshift.

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 investigate the relation between star formation (SF) and black hole accretion luminosities, using a sample of 492 type-2 active galactic nuclei (AGNs) at z < 0.22, which are detected in the far-infrared (FIR) surveys with AKARI and Herschel. We adopt FIR luminosities at 90 and 100 um as SF luminosities, assuming the proposed linear proportionality of star formation rate with FIR luminosities. By estimating AGN luminosities from [OIII]5007 and [OI]6300 emission lines, we find a positive linear trend between FIR and AGN luminosities over a wide dynamical range. This result appears to be inconsistent with the recent reports that low-luminosity AGNs show essentially no correlation between FIR and X-ray luminosities, while the discrepancy is likely due to the Malmquist and sample selection biases. By analyzing the spectral energy distribution, we find that pure-AGN candidates, of which FIR radiation is thought to be AGN-dominated, show significantly low-SF activities. These AGNs hosted by low-SF galaxies are rare in our sample (~ 1%). However, the low fraction of low-SF AGN is possibly due to observational limitations since the recent FIR surveys are insufficient to examine the population of high-luminosity AGNs hosted by low-SF galaxies.

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 use Chandra data to study the incidence and properties of Active Galactic Nuclei (AGN) in 16 intermediate redshift ($0.5 < z < 1.1$) X-ray-selected galaxy groups in the Chandra Deep Field-South. We measure an AGN fraction of $f(L_{X,H} > 10^{42}; M_R<-20) = 8.0_{-2.3}^{+3.0}\%$ at $\bar{z} \sim 0.74$, approximately a factor of two higher than the AGN fraction found for rich clusters at comparable redshift. This extends the trend found at low redshift for groups to have higher AGN fractions than clusters. Our estimate of the AGN fraction is also more than a factor of 3 higher than that of low redshift X-ray-selected groups. Using optical spectra from various surveys, we also constrain the properties of emission-line selected AGN in these groups. Contrary to the large population of X-ray AGN ($N(L_{X,H} > 10^{41}$ erg/s) = 25), we find only 4 emission-line AGN, 3 of which are also X-ray bright. Furthermore, most of the X-ray AGN in our groups are optically-dull (i.e. lack strong emission-lines) similar to those found in low redshift X-ray groups and clusters of galaxies. This contrasts with the AGN population found in low redshift optically-selected groups which are dominated by emission-line AGN. The differences between the optically and X-ray-selected AGN populations in groups are consistent with a scenario where most AGN in the densest environments are currently in a low accretion state.

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.

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)

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.

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 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 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 present J-band long-slit spectroscopic observation of NGC 1068 classified as a Seyfert 2 galaxy. J-band observations with OAO/ISLE provide clear detection of spatially extended [Fe II]1.257um and [P II]1.188um lines. We found that [Fe II]1.257um/[P II]1.188um increases with distance from a central continuum peak. Observed line ratios around the nucleus (continuum peak) are consistent with a typical value expected from photoionization models, while the ratios at 3" - 4" (210-280 pc) east and west of the nucleus are slightly higher than this. In the off nucleus region of NGC 1068 we also found a possible association between [Fe II]1.257um/[P II]1.188um and the radio continuum. This suggests a mild contribution of shock ionization induced by a radio jet outside nucleus while photoionization by the central energy source is dominant near the nucleus.

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.

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.

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