Source author record

Y. Taniguchi

Y. Taniguchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

156works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

156 published item(s)

preprint2022arXiv

Coherent photoproduction of the neutral pion and eta meson on the deuteron at incident energies below 1.15 GeV

Coherent photoproduction of the neutral pion and eta meson on the deuteron, $γ{d}${$\to$}$π^0η{d}$, has been experimentally studied at incident photon energies ranging from the reaction threshold to 1.15 GeV. The total cross section demonstrates a rapid rise below 1 GeV. The data are underestimated by the existing theoretical calculations based on quasi-free $π^0 η$ photoproduction on the nucleon followed by deuteron coalescence. At the same time, the data are rather well reproduced by the calculations taking into account the final-state interaction. We have also measured for the first time the differential cross sections: the $π^0 η$ invariant-mass distribution $dσ/dM_{πη}$, the $π^0 d$ invariant-mass distribution $dσ/dM_{πd}$, the $ηd$ invariant-mass distribution $dσ/dM_{ηd}$, and the distribution over the deuteron emission angle $dσ/dΩ_d$ in the overall center-of-mass frame. The measured cross section $dσ/dΩ_d$ does not exhibit strongly backward-peaking behavior predicted by the calculations. At all incident energies, an increase in $dσ/dM_{ηd}$ near the $ηd$ threshold is observed, which indicates a bound or virtual $ηd$ state resulting from a strong attraction between $η$ and a deuteron. The possibilities of using coherent $π^0η$ photoproduction on a nucleus to study the $η$-nuclear interaction are also discussed.

preprint2021arXiv

COSMOS2020: A panchromatic view of the Universe to $z\sim10$ from two complementary catalogs

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25<i<27$ reach a precision of $5\,\%$. Finally, these results are discussed in the context of previous, current, and future surveys in the COSMOS field. Compared to COSMOS2015, reaches the same photometric redshift precision at almost one magnitude deeper. Both photometric catalogs and their photometric redshift solutions and physical parameters will be made available through the usual astronomical archive systems (ESO Phase 3, IPAC IRSA, and CDS).

preprint2021arXiv

Description of isospin mixing by a generator coordinate method

Background: The isospin mixing is an interesting feature of atomic nuclei. It plays a crucial role in astrophysical nuclear reactions. However, it is not straightforward for variational nuclear structure models to describe it. Purpose: We propose a tractable method to describe the isospin mixing within a framework of the generator coordinate method and demonstrate its usability. Method: We generate the basis wave functions by applying the Fermi transition operator to the wave functions of isobars. The superposition of these basis wave functions and variationally obtained wave functions quantitatively describes the isospin mixing. Results: Using 14N as an example, we demonstrate that our method reasonably describes both T = 0 and 1 states and their mixing. Energy spectrum and E1 transition strengths are compared with the experimental data to confirm isospin mixing. Conclusion: The proposed method is effective enough to describe isospin mixing and is useful, for example, when we discuss α capture reactions of N = Z nuclei.

preprint2020arXiv

A puzzling non-detection of [O III] and [C II] from a $z \approx 7.7$ galaxy observed with ALMA

Characterizing the galaxy population in the early Universe holds the key to understanding the evolution of these objects and the role they played in cosmic reionization. However, the number of observations at the very highest redshifts are to date, few. In order to shed light on the properties of galaxies in the high-redshift Universe and their interstellar media, we observe the Lyman-$α$ emitting galaxy z7\_GSD\_3811 at $z=7.664$ with band 6 and 8 at the Atacama Large Millimeter/submillimeter Array (ALMA). We target the far-infrared [O III] 88 $μm$, [C II] 158 $μm$ emission lines and dust continuum in the star-forming galaxy z7\_GSD\_3811 with ALMA. We combine these measurements with earlier observations in the rest-frame ultraviolet (UV) in order to characterize the object, and compare results to those of earlier studies observing [O III] and [C II] emission in high-redshift galaxies. The [O III] 88 $μm$ and [C II] 158 $μm$ emission lines are undetected at the position of z7\_GSD\_3811, with $3σ$ upper limits of $1.6 \ \times \ 10^{8} \ L_{\odot}$ and $4.0 \ \times \ 10^{7} \ L_{\odot}$, respectively. We do not detect any dust continuum in band 6 nor band 8. The measured rms in the band 8 and band 6 continuum is 26 and 9.9 $μJy \ beam^{-1}$, respectively. Similar to several other high-redshift galaxies, z7\_GSD\_3811 exhibits low [C II] emission for its star formation rate compared to local galaxies. Furthermore, our upper limit on the [O III] line luminosity is lower than all the previously observed [O III] lines in high-redshift galaxies with similar ultraviolet luminosities. Our ALMA band 6 and 8 dust continuum observations imply that z7\_GSD\_3811 likely has a low dust content, and our non-detections of the [O III] and [C II] lines could indicate that z7\_GSD\_3811 has a low metallicity ($Z \lesssim 0.1 \ Z_{\odot}$).

preprint2020arXiv

alpha+28Si and 16O+16O molecular states, and their isoscalar monopole strengths

The properties of the alpha+28Si and 16O+16O molecular states which are embedded in the excited states of 32S and can have an impact on the stellar reactions are investigated using the antisymmetrized molecular dynamics. From the analysis of the cluster spectroscopic factors, the candidates of alpha+28Si and 16O+16O molecular states are identified close to and above the cluster threshold energies. The calculated properties of the alpha+28Si molecular states are consistent with those reported by the alpha+28Siresonant scattering experiments. On the other hand, the 16O+16O molecular state, which is predicted to be identical to the superdeformation of 32S, is inconsistent with the assignment proposed by an alpha inelastic scattering experiment. Our calculation suggests that the monopole transition from the ground state to the 16O+16O molecular state is rather weak and is not strongly excited by the alpha inelastic scattering.

preprint2019arXiv

$ωN$ scattering length from $ω$ photoproduction on the proton near the threshold

Photoproduction of the $ω$ meson on the proton has been experimentally studied near the threshold. The total cross sections are determined at incident energies ranging from 1.09 to 1.15 GeV. The 1/2 and 3/2 spin-averaged scattering length $a_{ωp}$ and effective range $r_{ωp}$ between the $ω$ meson and proton are estimated from the shape of the total cross section as a function of the incident photon energy: $a_{ωp} = \left(-0.97^{+0.16_{\rm stat}}_{-0.16_{\rm stat}}{}^{+0.03_{\rm syst}}_{-0.00_{\rm syst}}\right)+i \left(0.07^{+0.15_{\rm stat}}_{-0.14_{\rm stat}}{}^{+0.17_{\rm syst}}_{-0.09_{\rm syst}}\right)$ fm and $r_{ωp}=\left(+2.78^{+0.68_{\rm stat}}_{-0.54_{\rm stat}}{}^{+0.11_{\rm syst}}_{-0.13_{\rm syst}}\right)+i\left(-0.01^{+0.46_{\rm stat}}_{-0.50_{\rm stat}}{}^{+0.07_{\rm syst}}_{-0.00_{\rm syst}}\right)$ fm, resulting in a repulsive force. The real and imaginary parts for $a_{ωp}$ and $r_{ωp}$ are determined separately for the first time. A small $P$-wave contribution does not affect the obtained values.

preprint2016arXiv

ALMA observations of Lyman-alpha Blob 1: halo sub-structure illuminated from within

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

preprint2016arXiv

An event mixing technique for Bose-Einstein correlations of two pions in photoproduction around 1 GeV

We have developed an event mixing technique to observe Bose-Einstein correlations (BEC) between two identical neutral pions produced in photo-induced reactions in the non-perturbative QCD energy region. It is found that the missing-mass consistency cut and the pion-energy cut are essential for the event mixing method to effectively extract BEC observables. A Monte Carlo (MC) simulation is used to validate these constraints and confirms the efficiency of this method. Our work paves the way for similar BEC studies at lower energies where the multiplicity of emitted bosons is limited.

preprint2016arXiv

The VIMOS Ultra Deep Survey: Ly$α$ Emission and Stellar Populations of Star-Forming Galaxies at 2<z<2.5

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

preprint2015arXiv

2+1 flavor QCD simulation on a $96^4$ lattice

We generate $2+1$ flavor QCD configurations near the physical point on a $96^4$ lattice employing the 6-APE stout smeared Wilson clover action with a nonperturbative $c_{\rm SW}$ and the Iwasaki gauge action at $β=1.82$. The physical point is estimated based on the chiral perturbation theory using several data points generated by the reweighting technique from the simulation point, wherer $m_π$,$m_K$ and $m_Ω$ are used as physical inputs. The physics results include the quark masses, the hadron spectrum, the pseudoscalar meson decay constants and nucleon sigma terms, using the nonperturbative renormalization factors evaluated with the Schrodinger functional method.

preprint2015arXiv

Discovering extremely compact and metal-poor, star-forming dwarf galaxies out to z ~ 0.9 in the VIMOS Ultra-Deep Survey

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

preprint2015arXiv

Dust Heating by Low-mass Stars in Massive Galaxies at z<1

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

preprint2015arXiv

Dust properties of Lyman break galaxies at $z\sim3$

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

preprint2015arXiv

Evolution of the specific Star Formation Rate Function at z<1.4 - Dissecting the mass-SFR plane in COSMOS and GOODS

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

preprint2015arXiv

Isoscalar dipole transition as a probe for asymmetric clustering

Background: The sharp $1^-$ resonances with enhanced isoscalar dipole transition strengths are observed in many light nuclei at relatively small excitation energies, but their nature was unclear. Purpose: We show those resonances can be attributed to the cluster states with asymmetric configurations such as $α$+$^{16}{\rm O}$. We explain why asymmetric cluster states are strongly excited by the isoscalar dipole transition. We also provide a theoretical prediction of the isoscalar dipole transitions in $^{20}{\rm Ne}$ and $^{44}{\rm Ti}$. Method: The transition matrix is analytically derived to clarify the excitation mechanism. The nuclear model calculations by Brink-Bloch wave function and antisymmetrized molecular dynamics are also performed to provide a theoretical prediction for $^{20}{\rm Ne}$ and $^{44}{\rm Ti}$. Results: It is shown that the transition matrix is as large as the Weisskopf estimate even though the ground state is an ideal shell model state. Furthermore, it is considerably amplified if the ground state has cluster correlation. The nuclear model calculations predict large transition matrix to the $α$+$^{16}{\rm O}$ and $α$+$^{40}{\rm Ca}$ cluster states comparable with or larger than the Weisskopf estimate. Conclusion: We conclude that the asymmetric cluster states are strongly excited by the isoscalar dipole transition. Combined with the isoscalar monopole transition that populates the $0^+$ cluster states, the isoscalar transitions are promising probe for asymmetric clusters.

preprint2015arXiv

Mass and Axial current renormalization in the Schrödinger functional scheme for the RG-improved gauge and the stout smeared $O(a)$-improved Wilson quark actions

We present the quark mass and axial current renormalization factors for the RG-improved Iwasaki gauge action and three flavors of the stout smeared $O(a)$-improved Wilson quark action. We employ $α=0.1$ and $n_{\mathrm{step}}=6$ for the stout link smearing parameters and all links in the quark action are replaced with the smeared links. Using the Schrödinger functional scheme we evaluate the renormalization factors at $β=1.82$ where large scale simulations are being carried out.

preprint2015arXiv

Rest-frame Optical Emission Lines in Far-Infrared Selected Galaxies at z<1.7 from the FMOS-COSMOS Survey

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.

preprint2015arXiv

The composite nature of Dust-Obscured Galaxies (DOGs) at z~2-3 in the COSMOS field: I. A Far-Infrared View

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

preprint2015arXiv

The FMOS-COSMOS survey of star-forming galaxies at z~1.6 III. Survey design, performance, and sample characteristics

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.

preprint2015arXiv

The galaxy luminosity function at z ~ 6 and evidence for rapid evolution in the bright end from z ~ 7 to 5

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

preprint2015arXiv

The Lyman Continuum escape fraction of galaxies at z=3.3 in the VUDS-LBC/COSMOS field

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

preprint2015arXiv

The Quasar-LBG Two-point Angular Cross-correlation Function at z ~ 4 in the COSMOS Field

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.

preprint2015arXiv

The Subaru COSMOS 20: Subaru Optical Imaging of the HST COSMOS Field with 20 Filters

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

preprint2014arXiv

Discovery of a rich proto-cluster at z=2.9 and associated diffuse cold gas in the VIMOS Ultra-Deep Survey (VUDS)

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

preprint2014arXiv

Evolution of the Fraction of Clumpy Galaxies at 0.2<z<1.0 in the COSMOS field

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

preprint2014arXiv

The bright end of the galaxy luminosity function at z ~ 7: before the onset of mass quenching?

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

preprint2014arXiv

The evolution of clustering length, large-scale bias and host halo mass at 2<z<5 in the VIMOS Ultra Deep Survey (VUDS)

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

preprint2014arXiv

The evolving SFR-M_star relation and sSFR since z~5 from the VUDS spectroscopic survey

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

preprint2014arXiv

The FMOS-COSMOS survey of star-forming galaxies at z ~ 1.6. I. Hα-based star formation rates and dust extinction

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$α$ in emission between $1.6{\rm -}1.8 \mathrm{μm}$ with $f_{\rm Hα}\gtrsim4\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$. Here, we specifically focus on 271 sBzK-selected galaxies that yield a H$α$ 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α} \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$α$-based main sequence with a slope ($0.81\pm0.04$) and normalization similar to previous studies at these redshifts.

preprint2014arXiv

The FMOS-Cosmos Survey of Star-Forming Galaxies at z~1.6 II. The Mass-Metallicity Relation and the Dependence on Star Formation Rate and Dust Extinction

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.

preprint2014arXiv

The stellar mass - halo mass relation from galaxy clustering in VUDS: a high star formation efficiency at z~3

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

preprint2014arXiv

The VIMOS Ultra-Deep Survey (VUDS): fast increase of the fraction of strong Lyman alpha emitters from z=2 to z=6

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

preprint2014arXiv

The VIMOS Ultra-Deep Survey: ~10,000 galaxies with spectroscopic redshifts to study galaxy assembly at early epochs 2<z<~6

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

preprint2014arXiv

The zCOSMOS Redshift Survey: evolution of the light in bulges and discs since z~0.8

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

preprint2014arXiv

VIMOS Ultra-Deep Survey (VUDS): Witnessing the Assembly of a Massive Cluster at z~3.3

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

preprint2013arXiv

A Gravitational Lens Model for the Lya Emitter, LAE 221724+001716 at z=3.1 in the SSA 22 Field

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

preprint2013arXiv

A statistical relation between the X-ray spectral index and Eddington ratio of active galactic nuclei in deep surveys

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)

preprint2013arXiv

Charmed baryons at the physical point in 2+1 flavor lattice QCD

We investigate the charmed baryon mass spectrum using the relativistic heavy quark action on 2+1 flavor PACS-CS configurations previously generated on $32^3 \times 64$ lattice. The dynamical up-down and strange quark masses are tuned to their physical values, reweighted from those employed in the configuration generation. At the physical point, the inverse lattice spacing determined from the $Ω$ baryon mass gives $a^{-1}=2.194(10)$ GeV, and thus the spatial extent becomes $L = 32 a = 2.88(1)$ fm. Our results for the charmed baryon masses are consistent with experimental values, except for the mass of $Ξ_{cc}$, which has been measured by only one experimental group so far and has not been confirmed yet by others. In addition, we report values of other doubly and triply charmed baryon masses, which have never been measured experimentally.

preprint2013arXiv

Environmental Effects on the Star Formation Activity at z~0.9 in the COSMOS Field

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

preprint2013arXiv

Evolution of Galaxies and their Environments at z = 0.1 to 3 in COSMOS

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

preprint2013arXiv

Genuine Irregular Galaxies as a Relic of Building Blocks of Galaxies

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

preprint2013arXiv

Mass assembly in quiescent and star-forming galaxies since z=4 from UltraVISTA

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

preprint2013arXiv

Newly-quenched galaxies as the cause for the apparent evolution in average size of the population

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

preprint2013arXiv

The role of major mergers in shaping galaxies at 2<z<4 in the VUDS and VVDS surveys

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

preprint2012arXiv

1+1+1 flavor QCD + QED simulation at the physical point

We present the results of 1+1+1 flavor QCD+QED simulation at the physical point, in which the dynamical quark effects in QED and the up-down quark mass difference are incorporated by the reweighting technique. The physical quark masses together with the lattice spacing are determined with $m_{π^+}$, $m_{K^+}$, $m_{K^0}$ and $m_{Ω^-}$ as physical inputs. Calculations are carried out using a set of 2+1 flavor QCD configurations near the physical point generated by the non-perturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $β=1.9$ on a $32^3\times 64$ lattice. We evaluate the values of the up, down and strange quark masses individually with non-perturbative QCD renormalization.

preprint2012arXiv

Constraints on the faint end of the quasar luminosity function at z~5 in the COSMOS field

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.

preprint2012arXiv

Deep near-infrared spectroscopy of passively evolving galaxies at z>1.4

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

preprint2012arXiv

Discovery of bright z ~ 7 galaxies in the UltraVISTA survey

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

preprint2012arXiv

Ionization Source of a Minor-axis Cloud in the Outer Halo of M82

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

preprint2012arXiv

Millimeter imaging of submillimeter galaxies in the COSMOS field: Redshift distribution

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

preprint2012arXiv

Near-infrared spectroscopy of a nitrogen-loud quasar SDSS J1707+6443

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.

preprint2012arXiv

Photo-z Performance for Precision Cosmology II : Empirical Verification

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

preprint2012arXiv

Post-starburst Tidal Tails in the Archetypical Ultra Luminous Infrared Galaxy Arp 220

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

preprint2012arXiv

Spectral Energy Distributions of Type 1 AGN in the COSMOS Survey I - The XMM-COSMOS Sample

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

preprint2012arXiv

The dominant role of mergers in the size evolution of massive early-type galaxies since z ~ 1

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

preprint2011arXiv

$ρ$ Meson Decay in 2+1 Flavor Lattice QCD

We perform a lattice QCD study of the $ρ$ meson decay from the $N_f=2+1$ full QCD configurations generated with a renormalization group improved gauge action and a non-perturbatively $O(a)$-improved Wilson fermion action. The resonance parameters, the effective $ρ\toππ$ coupling constant and the resonance mass, are estimated from the $P$-wave scattering phase shift for the isospin I=1 two-pion system. The finite size formulas are employed to calculate the phase shift from the energy on the lattice. Our calculations are carried out at two quark masses, $m_π=410\,{\rm MeV}$ ($m_π/m_ρ=0.46$) and $m_π=300\,{\rm MeV}$ ($m_π/m_ρ=0.35$), on a $32^3\times 64$ ($La=2.9\,{\rm fm}$) lattice at the lattice spacing $a=0.091\,{\rm fm}$. We compare our results at these two quark masses with those given in the previous works using $N_f=2$ full QCD configurations and the experiment.

preprint2011arXiv

Black hole accretion and host galaxies of obscured quasars in XMM-COSMOS

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

preprint2011arXiv

Charm quark system at the physical point of 2+1 flavor lattice QCD

We investigate the charm quark system using the relativistic heavy quark action on 2+1 flavor PACS-CS configurations previously generated on $32^3 \times 64$ lattice. The dynamical up-down and strange quark masses are set to the physical values by using the technique of reweighting to shift the quark hopping parameters from the values employed in the configuration generation. At the physical point, the lattice spacing equals $a^{-1}=2.194(10)$ GeV and the spatial extent $L=2.88(1)$ fm. The charm quark mass is determined by the spin-averaged mass of the 1S charmonium state, from which we obtain $m_{\rm charm}^{\msbar}(μ= m_{\rm charm}^{\msbar}) = 1.260(1)(6)(35)$ GeV, where the errors are due to our statistics, scale determination and renormalization factor. An additional systematic error from the heavy quark is of order $α_s^2 f(m_Q a)(a Λ_{QCD})$, which is estimated to be a percent level if the factor $f(m_Q a)$ analytic in $m_Q a$ is of order unity. Our results for the charmed and charmed-strange meson decay constants are $f_D=226(6)(1)(5)$ MeV, $f_{D_s}=257(2)(1)(5)$ MeV, again up to the heavy quark errors of order $α_s^2 f(m_Q a)(a Λ_{QCD})$. Combined with the CLEO values for the leptonic decay widths, these values yield $|V_{cd}| = 0.205(6)(1)(5)(9)$, $|V_{cs}| = 1.00(1)(1)(3)(3)$, where the last error is on account of the experimental uncertainty of the decay widths.

preprint2011arXiv

Dissecting Photometric redshift for Active Galactic Nuclei using XMM- and Chandra-COSMOS samples

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

preprint2011arXiv

Fe K emission from active galaxies in the COSMOS field

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

preprint2011arXiv

The redshift and nature of AzTEC/COSMOS 1: A starburst galaxy at z=4.6

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

preprint2010arXiv

A Multiwavelength Study of a Sample of 70 micron Selected Galaxies in the COSMOS Field II: The Role of Mergers in Galaxy Evolution

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

preprint2010arXiv

A Runaway Black Hole in COSMOS: Gravitational Wave or Slingshot Recoil?

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

preprint2010arXiv

A z=1.82 Analog of Local Ultra-massive Elliptical Galaxies

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

preprint2010arXiv

Bars in early- and late-type disks in COSMOS

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

preprint2010arXiv

Calculation of $ρ$ meson decay width from the PACS-CS configurations

We present preliminary results on the $ρ$ meson decay width from $N_f=2+1$ full QCD configurations generated by PACS-CS Collaboration. The decay width is estimated from the $P$-wave scattering phase shift for the isospin $I=1$ two-pion system. The finite size formula presented by Lüscher in the center of mass frame and its extension to non-zero total momentum frame by Rummukainen and Gottlieb are employed for the calculations of the phase shift. Our calculations are carried out at $m_π=410\ {\rm MeV}$ ($m_π/m_ρ=0.46$) and $a=0.091\ {\rm fm}$ on a $32^3\times 64$ ($La=2.9 {\rm fm}$) lattice.

preprint2010arXiv

Characteristic Scales of Initial Density and Velocity Fields

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

preprint2010arXiv

Decomposing Star Formation and Active Galactic Nucleus with Spitzer Mid-Infrared Spectra: Luminosity Functions and Co-Evolution

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

preprint2010arXiv

K+a galaxies in the zCOSMOS Survey: Physical properties of systems in their post-starburst phase

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

preprint2010arXiv

Low-Metallicity Star Formation in High-Redshift Galaxies at z~8

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

preprint2010arXiv

Non-perturbative renormalization of quark mass in Nf=2+1 QCD with the Schroedinger functional scheme

We present an evaluation of the quark mass renormalization factor for Nf=2+1 QCD. The Schroedinger functional scheme is employed as the intermediate scheme to carry out non-perturbative running from the low energy region, where renormalization of bare mass is performed on the lattice, to deep in the high energy perturbative region, where the conversion to the renormalization group invariant mass or the MS-bar scheme is safely carried out. For numerical simulations we adopted the Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy regions and three lattice spacings to take the continuum limit at each scale. The regularization independent step scaling function of the quark mass for the Nf=2+1 QCD is obtained in the continuum limit. Renormalization factors for the pseudo scalar density and the axial vector current are also evaluated for the same action and the bare couplings as two recent large scale Nf=2+1 simulations; previous work of the CP-PACS/JLQCD collaboration, which covered the up-down quark mass range heavier than $m_π\sim 500$ MeV and that of PACS-CS collaboration for much lighter quark masses down to $m_π=155$ MeV. The quark mass renormalization factor is used to renormalize bare PCAC masses in these simulations.

preprint2010arXiv

Physical Point Simulation in 2+1 Flavor Lattice QCD

We present the results of the physical point simulation in 2+1 flavor lattice QCD with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $β=1.9$ on a $32^3 \times 64$ lattice. The physical quark masses together with the lattice spacing is determined with $m_π$, $m_K$ and $m_Ω$ as physical inputs. There are two key algorithmic ingredients to make possible the direct simulation at the physical point: One is the mass-preconditioned domain-decomposed HMC algorithm to reduce the computational cost. The other is the reweighting technique to adjust the hopping parameters exactly to the physical point. The physics results include the hadron spectrum, the quark masses and the pseudoscalar meson decay constants. The renormalization factors are nonperturbatively evaluated with the Schr{ö}dinger functional method. The results are compared with the previous ones obtained by the chiral extrapolation method.

preprint2010arXiv

Probing the faint end of the quasar luminosity function at z ~ 4 in the COSMOS field

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.

preprint2010arXiv

The COSMOS-WIRCam near-infrared imaging survey: I: BzK selected passive and star forming galaxy candidates at z>1.4

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

preprint2010arXiv

The Opacity of Galactic Disks at z~0.7

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

preprint2010arXiv

The XMM-Newton Wide-Field Survey in the COSMOS field (XMM-COSMOS): demography and multiwavelength properties of obscured and unobscured luminous AGN

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

preprint2010arXiv

Tracking the impact of environment on the Galaxy Stellar Mass Function up to z~1 in the 10k zCOSMOS sample

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

preprint2010arXiv

Ultraluminous X-ray sources out to z~0.3 in the COSMOS field

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

preprint2009arXiv

A Multiwavelength Study of a Sample of 70 micron Selected Galaxies in the COSMOS Field I: Spectral Energy Distributions and Luminosities

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

preprint2009arXiv

Galaxy Stellar Mass Assembly between 0.2<z<2 from the S-COSMOS survey

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

preprint2009arXiv

GEMINI 3D spectroscopy of BAL+IR+Fe II QSOs: II. IRAS 04505-2958 an explosive QSO with hypershell and a new scenario for galaxy formation and galaxy end

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

preprint2009arXiv

Lyman-alpha Blobs Like Company : The Discovery of A Candidate 100 kpc Lyman-alpha Blob Near to A Radio Galaxy with A Giant Lyman-alpha halo, B3 J2330+3927 at z=3.1

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

preprint2009arXiv

On the cosmic evolution of the scaling relations between black holes and their host galaxies: Broad Line AGN in the zCOSMOS survey

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

preprint2009arXiv

Photometric properties of Ly alpha emitters at z=4.86 in the COSMOS 2 square degree field

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.

preprint2009arXiv

Precise determination of the strong coupling constant in Nf=2+1 lattice QCD with the Schrödinger functional scheme

We present an evaluation of the running coupling constant for Nf=2+1 QCD. The Schroedinger functional scheme is used as the intermediate scheme to carry out non-perturbative running from the low energy region, where physical scale is introduced, to deep in the high energy perturbative region, where conversion to the MS-bar scheme is safely performed. Possible systematic errors due to the use of perturbation theory occur only in the conversion from three-flavor to four-flavor running coupling constant near the charm mass threshold, where higher order terms beyond 5th order in the $β$ function may not be negligible. For numerical simulations we adopted Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy region and three lattice spacings to take the continuum limit at each scale. A physical scale is introduced from the previous Nf=2+1 simulation of the CP-PACS/JL-QCD collaboration, which covered the up-down quark mass range heavier than $m_π\sim 500$ MeV.

preprint2009arXiv

Properties and environment of Radio Emitting Galaxies in the VLA-zCOSMOS survey

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

preprint2009arXiv

Star formation and dust obscuration at z~2: galaxies at the dawn of downsizing

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

preprint2009arXiv

The Build-Up of the Hubble Sequence in the COSMOS Field

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

preprint2009arXiv

The optical spectra of Spitzer 24 micron galaxies in the COSMOS field: II. Faint infrared sources in the zCOSMOS-bright 10k catalogue

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

preprint2009arXiv

The zCOSMOS survey: the role of the environment in the evolution of the luminosity function of different galaxy types

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

preprint2009arXiv

Universality of the Nf=2 Running Coupling

We investigate universality of the Nf=2 running coupling in the Schödinger functional scheme, by calculating the step scaling function in lattice QCD with the renorm alization group (RG) improved gauge action at both weak(u=0.9796) and strong(u=3.3340) couplings, where u=\bar{g}^2_SF with \bar{g}_SF being the running coupling in this scheme. In our main calculations, we use the tree-level value for O(a) improvement coefficients of boundary gauge fields. In addition we employ the 1-loop value for them in order to see how scaling behaviours are affected by them. In the continuum limit, the step scaling function obtained from the RG improved gauge actions agrees with the previous result obtained from the plaquette action within errors at both couplings, though errors of our result are larger. Combined fits using all data with the RG improved action as well as the plaquette action reduce errors in the continuum limit by 2% at the weak coupling and 22% at the strong coupling.

preprint2008arXiv

2+1 flavor lattice QCD simulation with $O(a)$-improved Wilson quarks

We present simulation details and results for the light hadron spectrum in N f = 2 + 1 lattice QCD with the nonperturbatively O(a)-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at a lattice spacing of 0.09 fm on a (2.9fm)^3 box using the PACS-CS computer. We employ the Luscher's domain-decomposed HMC algorithm with several improvements to reduce the degenerate up-down quark mass toward the physical value. So far the resulting pseudoscalar meson mass is ranging from 702MeV down to 156MeV. We discuss on the stability and the efficiency of the algorithm. The light harden spectrum extrapolated at the physical point is compared with the experimental values. We also present the values of the quark masses and the pseudoscalar meson decay constants.

preprint2008arXiv

2+1 Flavor Lattice QCD toward the Physical Point

We present the first results of the PACS-CS project which aims to simulate 2+1 flavor lattice QCD on the physical point with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Numerical simulations are carried out at the lattice spacing of $a=0.0907(13)$fm on a $32^3\times 64$ lattice with the use of the DDHMC algorithm to reduce the up-down quark mass. Further algorithmic improvements make possible the simulation whose ud quark mass is as light as the physical value. The resulting PS meson masses range from 702MeV down to 156MeV, which clearly exhibit the presence of chiral logarithms. An analysis of the PS meson sector with SU(3) ChPT reveals that the NLO corrections are large at the physical strange quark mass. In order to estimate the physical ud quark mass, we employ the SU(2) chiral analysis expanding the strange quark contributions analytically around the physical strange quark mass. The SU(2) LECs ${\bar l}_3$ and ${\bar l}_4$ are comparable with the recent estimates by other lattice QCD calculations. We determine the physical point together with the lattice spacing employing $m_π$, $m_K$ and $m_Ω$ as input. The hadron spectrum extrapolated to the physical point shows an agreement with the experimental values at a few % level of statistical errors, albeit there remain possible cutoff effects. We also find that our results of $f_π=134.0(4.2)$MeV, $f_K=159.4(3.1)$MeV and $f_K/f_π=1.189(20)$ with the perturbative renormalization factors are compatible with the experimental values. For the physical quark masses we obtain $m_{\rm ud}^\msbar=2.527(47)$MeV and $m_{\rm s}^\msbar=72.72(78)$MeV extracted from the axial-vector Ward-Takahashi identity with the perturbative renormalization factors.

preprint2008arXiv

Charm quark system in 2+1 flavor lattice QCD using the PACS-CS configurations

We study heavy-heavy and heavy-light quark systems for charm with a relativistic heavy quark action in 2+1 flavor lattice QCD. Configurations are generated by the PACS-CS Collaboration at the lattice spacing is $a=0.09$ fm with the lattice size of $32^3\times 64$ employing the $O(a)$-improved Wilson quark action and the Iwasaki gauge action. We present preliminary results for the charmonium spectrum and the $D$ and $D_s$ meson decay constants evaluated at 3.5 MeV$< m_{\rm ud}<$ 12 MeV with $m_{\rm s}$ around the physical value. We investigate the dynamical quark mass dependences of the hyperfine and the orbital splittings. The decay constants are compared with the recent experimental values.

preprint2008arXiv

Chasing highly obscured QSOs in the COSMOS field

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

preprint2008arXiv

Light quark masses from unquenched lattice QCD

We calculate the light meson spectrum and the light quark masses by lattice QCD simulation, treating all light quarks dynamically and employing the Iwasaki gluon action and the nonperturbatively O(a)-improved Wilson quark action. The calculations are made at the squared lattice spacings at an equal distance a^2~0.005, 0.01 and 0.015 fm^2, and the continuum limit is taken assuming an O(a^2) discretization error. The light meson spectrum is consistent with experiment. The up, down and strange quark masses in the \bar{MS} scheme at 2 GeV are \bar{m}=(m_{u}+m_{d})/2=3.55^{+0.65}_{-0.28} MeV and m_s=90.1^{+17.2}_{-6.1} MeV where the error includes statistical and all systematic errors added in quadrature. These values contain the previous estimates obtained with the dynamical u and d quarks within the error.

preprint2008arXiv

SU(2) and SU(3) chiral perturbation theory analyses on meson and baryon masses in 2+1 flavor lattice QCD

We investigate the quark mass dependence of meson and baryon masses obtained from 2+1 flavor dynamical quark simulations performed by the PACS-CS Collaboration. With the use of SU(2) and SU(3) chiral perturbation theories up to NLO, we examine the chiral behavior of the pseudoscalar meson masses and the decay constants in terms of the degenerate up-down quark mass ranging form 3 MeV to 24 MeV and two choices of the strange quark mass around the physical value. We discuss the convergence of the SU(2) and SU(3) chiral expansions and present the results for the low energy constants. We find that the SU(3) expansion is not convergent at NLO for the physical strange quark mass. The chiral behavior of the nucleon mass is also discussed based on the SU(2) heavy baryon chiral perturbation theory up to NNLO. Our results show that the expansion is well behaved only up to m_pi^2 ~ 0.2 GeV^2.

preprint2007arXiv

2+1 flavor light hadron spectrum and quark masses with the O(a) improved Wilson-clover quark formalism

We present a summary of results of the joint CP-PACS and JLQCD project toward a 2+1 flavor full QCD simulation with the O(a)-improved Wilson quark formalism and the Iwasaki gauge action. Configurations were generated during 2002-2005 at three lattice spacings, a~0.076, 0.100 and 0.122 fm, keeping the physical volume constant at (2.0fm)^3. Up and down quark masses are taken in the range m_{PS}/m_V~0.6-0.78. We have completed the analysis for the light meson spectrum and quark masses in the continuum limit using the full configuration set. The predicted meson masses reproduce experimental values in the continuum limit at a 1% level. The average up and down, and strange quark masses turn out to be m_{ud}^{\bar{MS}}(μ=2 GeV)=3.50(14)({}^{+26}_{-15}) MeV and m_s^{\bar{MS}}(μ=2 GeV)=91.8(3.9)({}^{+6.8}_{-4.1}) MeV. We discuss our future strategy toward definitive results on hadron spectroscopy with the Wilson-clover formalism.

preprint2007arXiv

Application of chiral perturbation theory to 2+1 flavor lattice QCD with O(a)-improved Wilson quarks

We apply chiral perturbation theory to the pseudoscalar meson mass and decay constant data obtained in the PACS-CS Project toward 2+1 flavor lattice QCD simulations with the O(a)-improved Wilson quarks. We examine the existence of chiral logarithms in the quark mass range from m_{ud}=47 MeV down to 6 MeV on a (2.8 fm)^3 box with the lattice spacing a=0.09 fm. Several low energy constants are determined. We also discuss the magnitude of finite size effects based on chiral perturbation theory.

preprint2007arXiv

COSMOS : Hubble Space Telescope Observations

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

preprint2007arXiv

Deep GALEX Imaging of the HST/COSMOS Field: A First Look at the Morphology of z~0.7 Star-forming Galaxies

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

preprint2007arXiv

Deep GEMINI GMOS-IFU spectroscopy of BAL QSOs: I. Decoupling the BAL, QSO, starburst, NLR, supergiant bubbles and galactic wind in Mrk 231

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

preprint2007arXiv

Light hadron spectrum with 2+1 flavor dynamical $O(a)$-improved Wilson quarks

We present preliminary results for the light harden spectrum in $N_f=2+1$ lattice QCD obtained with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at $β=1.90$ on a $32^3 \times 64$ lattice using the PACS-CS computer. We employ Lüscher's domain-decomposed HMC algorithm to reduce the up-down quark masses toward the physical value. The pseudoscalar meson masses range from 730 MeV down to 210 MeV. We compare the light harden spectrum extrapolated to the physical point with the experimental values.

preprint2007arXiv

The First Release COSMOS Optical and Near-IR Data and Catalog

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

preprint2007arXiv

The Spitzer Legacy Survey of the HST-ACS 2 sq. deg. COSMOS Field: survey strategy and first analysis

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

preprint2006arXiv

2+1 Flavor Lattice QCD with Luescher's Domain-Decomposed HMC Algorithm

We report on a study of 2+1 flavor lattice QCD with the $O(a)$-improved Wilson quarks on a $16^3\times 32$ lattice at the lattice spacing $1/a\approx 2$GeV employing Luescher's domain-decomposed HMC(LDDHMC) algorithm. This is dedicated to a preliminary study for the PACS-CS project which plans to complete the Wilson-clover $N_f=2+1$ program lowering the up-down quark masses close to the physical values as much as possible. We focus on three issues: (i) how light quark masses we can reach with LDDHMC, (ii) efficiency of the algorithm compared with the conventional HMC, (iii) parameter choice for the production runs on PACS-CS.

preprint2006arXiv

An application of the UV-filtering preconditioner to the Polynomial Hybrid Monte Carlo algorithm

We apply the UV-filtering preconditioner, previously used to improve the Multi-Boson algorithm, to the Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm. The performance test for the algorithm is given for the plaquette gauge action and the $O(a)$-improved Wilson action at $β=5.2, c_{\mathrm{sw}}=2.02, M_π/M_ρ\sim 0.8$ and 0.7 on a $16^3\times 48$ lattice. We find that the UV-filtering reduces the magnitude of the molecular dynamics force from the pseudo fermion by a factor 3 by tuning the UV-filter parameter. Combining with the multi-time scale molecular dynamics integrator we achieve a factor 2 improvement.

preprint2006arXiv

An estimate of the eta and eta-prime meson masses in Nf=2+1 lattice QCD

Masses of the eta and eta-prime mesons are estimated in Nf=2+1 lattice QCD with the non-perturbatively O(a) improved Wilson quark action and the Iwasaki RG-improved gluon action, using CP-PACS/JLQCD configurations on a 16^3 x 32 lattice at beta=1.83 (lattice spacing is 0.122 fm). We apply a stochastic noise estimator technique combined with smearing method to evaluate correlators among flavor SU(2) singlet pseudoscalar operators and strange pseudoscalar operators for 10 combinations of up/down and strange quark masses. The correlator matrix is then diagonalized to identify signals for mass eigenstates. Masses of the ground state and the first excited state extrapolated to the physical point are m_eta= 0.545(16) GeV and m_eta-prime= 0.871(46) GeV, being close to the experimental values of the eta and eta-prime masses.

preprint2006arXiv

Kaon $B$-parameters for Generic $ΔS=2$ Four-Quark Operators in Quenched Domain Wall QCD

We present a study of $B$-parameters for generic $ΔS=2$ four-quark operators in domain wall QCD. Our calculation covers all the $B$-parameters required to study the neutral kaon mixing in the standard model (SM) and beyond it. We evaluate one-loop renormalization factors of the operators employing the plaquette and Iwasaki gauge actions. Numerical simulations are carried out in quenched QCD with both gauge actions on $16^3\times 32\times 16$ and $24^3\times 32\times 16$ at the lattice spacing $1/a\approx 2$GeV. We investigate the relative magnitudes of the non-SM $B$-parameters to the SM one, which are compared with the previous results obtained with the overlap and the clover quark actions.

preprint2006arXiv

Ly Alpha Emitters at z=5.7 in the Subaru Deep Field

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

preprint2006arXiv

Nonperturbative $O(a)$ improvement of the Wilson quark action with the RG-improved gauge action using the Schrödinger functional method

We perform a nonperturbative determination of the $O(a)$-improvement coefficient $c_{\rm SW}$ and the critical hopping parameter $κ_c$ for $N_f$=3, 2, 0 flavor QCD with the RG-improved gauge action using the Schrödinger functional method. In order to interpolate $c_{\rm SW}$ and $κ_c$ as a function of the bare coupling, a wide range of $β$ from the weak coupling region to the moderately strong coupling points used in large-scale simulations is studied. Corrections at finite lattice size of $O(a/L)$ turned out to be large for the RG-improved gauge action, and hence we make the determination at a size fixed in physical units using a modified improvement condition. This enables us to avoid $O(a)$ scaling violations which would remain in physical observables if $c_{\rm SW}$ determined for a fixed lattice size $L/a$ is used in numerical simulations.

preprint2006arXiv

The stellar content of the COSMOS field as derived from morphological and SED based gtar/galaxy Separation

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

preprint2005arXiv

Light hadron spectrum and quark masses in 2+1 flavor QCD

CP-PACS and JLQCD collaborations are carrying out a joint project of the 2+1 flavor full QCD simulation. Gauge configurations are generated for the non-perturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action using PHMC algorithm at three lattice spacings, $a\sim 0.076$, 0.010 and 0.122 fm, with a fixed physical volume $(2.0 fm)^3$. We present analysis for the light meson spectrum and quark masses in the continuum limit, which are determined using data obtained from the simulations at the two coarser lattices. Our simulations reproduce experimental values of meson masses. The ud and strange quark masses turn out to be $m_{ud}^{\bar{MS}}(μ=2 GeV)=3.34(23) MeV$ and $m_s^{\bar{MS}}(μ=2 GeV)=86.7(5.9) MeV$. We also show preliminary results at our finest lattice spacing for which simulations are still being continued.

preprint2005arXiv

Non-perturbative O(a)-improvement of Wilson quark action in three-flavor QCD with plaquette gauge action

We perform a non-perturbative determination of the O(a)-improvement coefficient c_SW for the Wilson quark action in three-flavor QCD with the plaquette gauge action. Numerical simulations are carried out in a range of β=12.0-5.2 on a single lattice size of 8^3x16 employing the Schrödinger functional setup of lattice QCD. As our main result, we obtain an interpolation formula for c_SW and the critical hopping parameter K_c as a function of the bare coupling. This enables us to remove O(a) scaling violation from physical observables in future numerical simulation in the wide range of β. Our analysis with a perturbatively modified improvement condition for c_SW suggests that finite volume effects in c_SW are not large on the 8^3x16 lattice. We investigate N_f dependence of c_SW by additional simulations for N_f=4, 2 and 0 at β=9.6. As a preparatory step for this study, we also determine c_SW in two-flavor QCD at β=5.2. At this β, several groups carried out large-scale calculations of the hadron spectrum, while no systematic determination of c_SW has been performed.

preprint2005arXiv

The PACS-CS Project

We describe our plan to develop a large-scale cluster system with a peak speed of 14.3Tflops for lattice QCD at the Center for Computational Sciences, University of Tsukuba, as a successor to the current 0.6Tflops CP-PACS computer. The system consist of 2560 nodes connected by a 16x16x10 three-dimensional hyper crossbar network. Each node has a single low-voltage 2.8GHz Xeon processor and 2GBytes of memory with 6.4GBytes/sec bandwidth, and 160 GBytes of disk in RAID1 mode. The network link in each of the three directions is made of dual Gigabit Ethernet with the peak throughput of 250MByte/sec. Hence each node has an aggregate network bandwidth of 750MByte/sec. The system will run under Linux and SCore, and an extension of the PM driver is developed for the network. The system will be developed jointly with Hitachi Limited. The installation is scheduled in the first quarter of Japanese Fiscal 2006 (April-June 2006) and the start of operation is expected in July 2006.

preprint2005arXiv

The Subaru Deep Field Project: Lyman$α$ Emitters at Redshift of 6.6

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

preprint2004arXiv

A scaling study of the step scaling function in SU(3) gauge theory with improved gauge actions

We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the renormalization-group improved Iwasaki gauge action and the perturbatively improved Lüscher-Weisz gauge action. We confirm that the step scaling functions from the improved gauge actions agree with that previously obtained from the plaquette action within errors in the continuum limit at both weak and strong coupling regions. We also investigate how different choices of boundary counter terms for the improved gauge actions affect the scaling behavior. In the extrapolation to the continuum limit, we observe that the cut off dependence becomes moderate for the Iwasaki action, if a perturbative reduction of scaling violations is applied to the simulation results. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.

preprint2004arXiv

A scaling study of the step scaling function of quenched QCD with improved gauge actions

We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the Iwasaki gauge action and the Luescher-Weisz gauge action. In particular, we test the choice of boundary counter terms and apply a perturbative procedure for removal of lattice artifacts for the simulation results in the extrapolation procedure. We confirm the universality of the step scaling functions at both weak and strong coupling regions. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.

preprint2004arXiv

Lattice QCD calculation of the proton decay matrix element in the continuum limit

We present a quenched lattice QCD calculation of the αand βparameters of the proton decay matrix element. The simulation is carried out using the Wilson quark action at three values of the lattice spacing in the range a\approx 0.1-0.064 fm to study the scaling violation effect. We find only mild scaling violation when the lattice scale is determined by the nucleon mass. We obtain in the continuum limit, |α(NDR,2GeV)|=0.0090(09)(^{+5}_{-19})GeV^3 and |β(NDR,2GeV)|=0.0096(09)(^{+6}_{-20})GeV^3 with αand βin a relatively opposite sign, where the first error is statistical and the second is due to the uncertainty in the determination of the physical scale.

preprint2004arXiv

Non-perturbative renormalization of meson decay constants in quenched QCD for a renormalization group improved gauge action

Renormalization constants ($Z$-factors) of vector and axial-vector currents are determined non-perturbatively in quenched QCD for a renormalization group improved gauge action and a tadpole improved clover quark action using the Schrödinger functional method. Non-perturbative values of $Z$-factors turn out to be smaller than one-loop perturbative values by $O(15%)$ at lattice spacing of $a^{-1}\approx$ 1 GeV. The pseudoscalar and vector meson decay constants calculated with the non-perturbative $Z$-factors show a much better scaling behavior compared to previous results obtained with tadpole improved one-loop $Z$-factors. In particular, the non-perturbative $Z$-factors normalized at infinite physical volume show that scaling violation of the decay constants are within about 10% up to the lattice spacing $a^{-1}\sim 1$ GeV. The continuum estimates obtained from data in the range $a^{-1}=$ 1 -- 2 GeV agree with those determined from finer lattices ($a^{-1}\sim 2-4$ GeV) with the standard action.

preprint2004arXiv

The Subaru Deep Field: The Optical Imaging Data

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy survey over a blank field as large as 34'x27'. The program consists of very deep multi-band optical imaging, near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific goals of the project are to construct large samples of Lyman-break galaxies at z~4-5 and Lyman alpha emitters at z~5.7 and 6.6, and to make detailed studies these very high-redshift galaxy populations. In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked images and object catalogs in seven bandpasses.The optical imaging was made through five broad-band filters, B, V, R, i', z', and two narrow-band filters, NB816 (lambda_c=8150A) and NB921 (lambda_c=9196A) with almost 10 hours long integrations for each band. The limiting magnitudes measured at 3-sigma on a 2" aperture are B=28.45, V=27.74, R=27.80, i'=27.43, z'=26.62, NB816=26.63, and NB921=26.54 in the AB system. The object catalog constructed for each of the seven bands contains more than 10^5 objects. The galaxy number counts corrected for detection incompleteness and star count contribution are found to be consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/.

preprint2003arXiv

Calculation of Non-Leptonic Kaon Decay Amplitudes from $K\toπ$ Matrix Elements in Quenched Domain-Wall QCD

We explore application of the domain wall fermion formalism of lattice QCD to calculate the $K\toππ$ decay amplitudes in terms of the $K\toπ$ and $K\to 0$ hadronic matrix elements through relations derived in chiral perturbation theory. Numerical simulations are carried out in quenched QCD using domain-wall fermion action for quarks and an RG-improved gauge action for gluons on a $16^3\times 32\times 16$ and $24^3\times 32\times 16$ lattice at $β=2.6$ corresponding to the lattice spacing $1/a\approx 2$GeV. Quark loop contractions which appear in Penguin diagrams are calculated by the random noise method, and the $ΔI=1/2$ matrix elements which require subtractions with the quark loop contractions are obtained with a statistical accuracy of about 10%. We confirm the chiral properties required of the $K\toπ$ matrix elements. Matching the lattice matrix elements to those in the continuum at $μ=1/a$ using the perturbative renormalization factor to one loop order, and running to the scale $μ=m_c=1.3$ GeV with the renormalization group for $N_f=3$ flavors, we calculate all the matrix elements needed for the decay amplitudes. With these matrix elements, the $ΔI=3/2$ decay amplitude shows a good agreement with experiment in the chiral limit. The $ΔI=1/2$ amplitude, on the other hand, is about 50--60% of the experimental one even after chiral extrapolation. In view ofthe insufficient enhancement of the $ΔI=1/2$ contribution, we employ the experimental values for the real parts of the decay amplitudes in our calculation of $ε'/ε$. We find that the $ΔI=3/2$ contribution is larger than the $ΔI=1/2$ contribution so that $ε'/ε$ is negative and has a magnitude of order $10^{-4}$. Possible reasons for these unsatisfactory results are discussed.

preprint2003arXiv

Flavor Singlet Meson Mass in the Continuum Limit in Two-Flavor Lattice QCD

We present results for the mass of the eta-prime meson in the continuum limit for two-flavor lattice QCD, calculated on the CP-PACS computer, using a renormalization-group improved gauge action, and Sheikholeslami and Wohlert's fermion action with tadpole-improved csw. Correlation functions are measured at three values of the coupling constant beta corresponding to the lattice spacing a approx. 0.22, 0.16, 0.11 fm and for four values of the quark mass parameter kappa corresponding to mpi over mrho approx. 0.8, 0.75, 0.7 and 0.6. For each beta, kappa pair, 400-800 gauge configurations are used. The two-loop diagrams are evaluated using a noisy source method. We calculate eta-prime propagators using local sources, and find that excited state contributions are much reduced by smearing. A full analysis for the smeared propagators gives metaprime=0.960(87)+0.036-0.248 GeV, in the continuum limit, where the second error represents the systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations.

preprint2003arXiv

Non-perturbative calculation of $Z_V$ and $Z_A$ in domain-wall QCD on a finite box

We report on a non-perturbative evaluation of the renormalization factors for the vector and axial-vector currents, $Z_V$ and $Z_A$, in the quenched domain-wall QCD (DWQCD) with plaquette and renormalization group improved gauge actions. We take the Dirichlet boundary condition for both gauge and domain-wall fermion fields on the finite box, and introduce the flavor-chiral Ward-Takahashi identities to calculate the renormalization factors. As a test of the method, we numerically confirm the expected relation that $Z_V \simeq Z_A$ in DWQCD. Employing two different box sizes for the numerical simulations at several values of the gauge coupling constant $g^2$ and the domain-wall height $M$, we extrapolate $Z_V$ to the infinite volume to remove $a/L$ errors. We finally give the interpolation formula of $Z_V$ in the infinite volume as a function of $g^2$ and $M$.

preprint2003arXiv

Non-perturbative renormalization of vector and axial vector currents in quenched QCD for a renormalization group improved gauge action

Renormalization constants of vector ($Z_V$) and axial-vector ($Z_A$) currents are determined non-perturbatively in quenched QCD for an RG-improved gauge action and a tadpole-improved clover quark action using the Schrödinger functional method. Meson decay constants $f_ρ$ and $f_π$ show much better scaling when $Z_V$ and $Z_A$ estimated for infinite physical volume are used instead of $Z$-factors from tadpole-improved one-loop perturbation theory.

preprint2003arXiv

The Discovery of Two Lyman$α$ Emitters Beyond Redshift 6 in the Subaru Deep Field

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

preprint2003arXiv

Two flavors of dynamical quarks on anisotropic lattices

We report on our study of two-flavor full QCD on anisotropic lattices using $O(a)$-improved Wilson quarks coupled with an RG-improved glue. The bare gauge and quark anisotropies corresponding to the renormalized anisotropy $ξ=a_s/a_t = 2$ are determined as functions of $β$ and $κ$, which covers the region of spatial lattice spacings $a_s\approx 0.28$--0.16 fm and $m_{PS}/m_V\approx 0.6$--0.9. The calibrations of the bare anisotropies are performed with the Wilson loop and the meson dispersion relation at 4 lattice cutoffs and 5--6 quark masses. Using the calibration results we calculate the meson mass spectrum and the Sommer scale $r_0$. We confirm that the values of $r_0$ calculated for the calibration using pseudo scalar and vector meson energy momentum dispersion relation coincide in the continuum limit within errors. This work serves to lay ground toward studies of heavy quark systems and thermodynamics of QCD including the extraction of the equation of state in the continuum limit using Wilson-type quark actions.

preprint2002arXiv

A Multi-band Photometric Study of Tidal Debris in A Compact Group of Galaxies: Seyfert's Sextet

In order to investigate the properties of the prominent tidal debris feature extending to the northeast of the compact group of galaxies Seyfert's Sextet, we analyzed multi-band (U, B, V, VR, R, I, J, H and K') photometric imaging data and obtained the following results: 1) The radial surface brightness distribution of this tidal debris in Seyfert's Sextet (TDSS) in each band appears to be well approximated by an exponential profile. 2) The observed B-V color of TDSS is similar to those of dwarf elliptical galaxies in nearby clusters. 3) Comparing the spectral energy distribution (SED) of TDSS with theoretical photometric evolution models and with the SED of the stars in the outer part of HCG 79b, we find that its SED is comparable to that of a $\sim$ 10 Gyr-old stellar population with solar metallicity, similar to the stellar population in the outer part of HCG 79b. This suggests that TDSS consists of stars that may have been liberated from HCG 79b by strong galaxy interactions, not a pre-existing dwarf galaxy previously thought.

preprint2002arXiv

Charmonium Spectrum from Quenched Anisotropic Lattice QCD

We present a detailed study of the charmonium spectrum using anisotropic lattice QCD. We first derive a tree-level improved clover quark action on the anisotropic lattice for arbitrary quark mass. The heavy quark mass dependences of the improvement coefficients, i.e. the ratio of the hopping parameters $ζ=K_t/K_s$ and the clover coefficients $c_{s,t}$, are examined at the tree level. We then compute the charmonium spectrum in the quenched approximation employing $ξ= a_s/a_t = 3$ anisotropic lattices. Simulations are made with the standard anisotropic gauge action and the anisotropic clover quark action at four lattice spacings in the range $a_s$=0.07-0.2 fm. The clover coefficients $c_{s,t}$ are estimated from tree-level tadpole improvement. On the other hand, for the ratio of the hopping parameters $ζ$, we adopt both the tree-level tadpole-improved value and a non-perturbative one. We calculate the spectrum of S- and P-states and their excitations. The results largely depend on the scale input even in the continuum limit, showing a quenching effect. When the lattice spacing is determined from the $1P-1S$ splitting, the deviation from the experimental value is estimated to be $\sim$30% for the S-state hyperfine splitting and $\sim$20% for the P-state fine structure. Our results are consistent with previous results at $ξ= 2$ obtained by Chen when the lattice spacing is determined from the Sommer scale $r_0$. We also address the problem with the hyperfine splitting that different choices of the clover coefficients lead to disagreeing results in the continuum limit.

preprint2002arXiv

Exploring QCD at small sea quark masses with improved Wilson-type quarks

We explore the region of small sea quark masses below $m_{PS}/m_V=0.5$ in two-flavor QCD using a mean-field improved clover quark action and an RG-improved gauge action at $a \simeq 0.2$ fm on $12^3 \times 24$ and $16^3 \times 24$ lattices. We find that instability of the standard BiCGStab algorithm at small quark masses can be mostly removed by the BiCGStab(DS-$L$) algorithm, which employs $L$-th minimal residual polynomials with a dynamical selection of $L$. We also find singular spikes of $ΔH$ in the HMC algorithm at moderate values of $Δτ$. Nature of the spike is studied. We also study finite-size effects and chiral properties of meson masses.

preprint2002arXiv

I=2 Pion Scattering Phase Shift with Wilson Fermions

We present a lattice QCD calculation of the scattering phase shift for the I=2 $S$-wave two-pion system using the finite size method proposed by Lüscher. We work in the quenched approximation employing the standard plaquette action at $β=5.9$ for gluons and the Wilson fermion action for quarks. The phase shift is extracted from the energy eigenvalues of the two-pion system, which are obtained by a diagonalization of the pion 4-point function evaluated for a set of relative spatial momenta. In order to change momentum of the two-pion system, calculations are carried out on $24^3\times 60$, $32^3\times 60$, and $48^3\times 60$ lattices. The phase shift is successfully calculated over the momentum range $0 < p^2 < 0.3 {\rm GeV}^2$.

preprint2002arXiv

Runaway Merging of Black Holes: Analytical Constraint on the Timescale

Following the discovery of a black hole (BH) with a mass of 10^3-10^6 M(sun) in a starburst galaxy M82, we study formation of such a BH via successive merging of stellar-mass BHs within a star cluster. The merging has a runaway characteristic. This is because massive BHs sink into the cluster core and have a high number density, and because the merging probability is higher for more massive BHs. We use the Smoluchowski equation to study analytically the evolution of the BH mass distribution. Under favorable conditions, which are expected for some star clusters in starburst galaxies, the timescale of the runaway merging is at most of order 10^7 yr. This is short enough to account for the presence of a BH heavier than 10^3 M(sun) in an ongoing starburst region.

preprint2001arXiv

Charmonium spectrum from quenched QCD on anisotropic lattices

We present our final results of the charmonium spectrum in quenched QCD on anisotropic lattices. Simulations are made with the plaquette gauge action and a tadpole improved clover quark action employing $ξ= a_s/a_t = 3$. We calculate the spectrum of S- and P-states and their excitation, and study the scaling behavior of mass splittings. Comparison is made with the experiment and previous lattice results. The issue of hyperfine splitting for different choices of the clover coefficients obtained by Klassen is discussed.

preprint2001arXiv

Chiral property of domain-wall fermion from eigenvalues of 4D Wilson-Dirac Operator

We investigate a chiral property of the domain-wall fermion (DWF) system using the four-dimensional hermitian Wilson-Dirac operator $H_W$. A formula expressing the Ward-Takahashi identity quark mass $m_{5q}$ with eigenvalues of this operator is derived, which well explains the $N_5$ dependence of $m_{5q}$ observed in previous numerical simulations. We further discuss the chiral property of DWF in the large volume in terms of the spectra of $H_W$.

preprint2001arXiv

Equation of state for pure SU(3) gauge theory on anisotropic lattices

We present results for the equation of state for pure SU(3) gauge theory obtained on anisotropic lattices with the anisotropy $ξ\equiv a_s/a_t = 2$. The pressure and energy density are calculated on $N_t / ξ= 4, 5$ and 6 lattices with the integral method. They are found to satisfy the leading $1/N_t^2$ scaling from our coarsest lattice $N_t/ξ=4$. This enables us to carry out well controlled continuum extrapolations. We find that the pressure and energy density agree with those obtained using the isotropic plaquette action, but have smaller and more reliable errors.

preprint2001arXiv

Equation of state in finite-temperature QCD with two flavors of improved Wilson quarks

We present results of a first study of equation of state in finite-temperature QCD with two flavors of Wilson-type quarks. Simulations are made on lattices with temporal size $N_t=4$ and 6, using an RG-improved action for the gluon sector and a meanfield-improved clover action for the quark sector. The lines of constant physics corresponding to fixed values of the ratio $m_{\rm PS}/m_{\rm V}$ of the pseudo-scalar to vector meson masses at zero temperature are determined, and the beta functions which describe the renormalization-group flow along these lines are calculated. Using these results, the energy density and the pressure are calculated as functions of temperature along the lines of constant physics in the range $m_{\rm PS}/m_{\rm V} = 0.65$--0.95. The quark mass dependence in the equation of state is found to be small for $m_{\rm PS}/m_{\rm V} \simlt 0.8$. Comparison of results for $N_t=4$ and $N_t=6$ lattices show significant scaling violation present in the $N_t=4$ results. At high temperatures the results for $N_t=6$ are quite close to the continuum Stefan-Boltzmann limit, suggesting the possibility of a precise continuum extrapolation of thermodynamic quantities from simulations at $N_t\simgt 6$.

preprint2001arXiv

Far-Infrared Census of Starburst-Seyfert Connection

Far-infrared flux densities are newly extracted from the IRAS database for the RSA and CfA complete samples of Seyfert galaxies. These data are used to classify the Seyfert galaxies into those where the far-infrared continuum emission is dominated by the active galactic nucleus (AGN), circumnuclear starburst, or host galaxy. While AGN-dominant objects consist of comparable numbers of Seyfert 1 and 2 galaxies, starburst- and host-dominant objects consist preferentially of Seyfert 2 galaxies. Thus, in addition to the dusty torus, the circumnuclear starburst region and host galaxy are important in hiding the broad-line region. Morphologically, starburst-dominant Seyfert galaxies are of later types and more strongly interacting than AGN-dominant Seyfert galaxies. In a later-type galaxy, the AGN central engine has a lower Eddington luminosity, and the gaseous content is higher. The gas is efficiently supplied to the starburst via a galaxy-galaxy interaction. Morphologies of host-dominant Seyfert galaxies are of various types. Since starbursts in Seyfert galaxies are older than those in classical starburst galaxies, we propose an evolution from starburst to starburst-dominant Seyfert to host-dominant Seyfert for a late-type galaxy. An evolution from AGN-dominant Seyfert to host-dominant Seyfert is proposed for an early-type galaxy. These sequences have durations of a few x 10^8 yr and occur repeatedly within a galaxy during its evolution from a late type to an early type.

preprint2001arXiv

Kaon B parameter from quenched domain-wall QCD

We report on a calculation of $B_K$ with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at $β=2.6$ and 2.9 corresponding to $a^{-1}\approx 2$GeV and 3GeV. Effects due to finite fifth dimensional size $N_5$ and finite spatial size $N_σ$ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain $B_K(μ= 2 GeV)= 0.5746(61)(191)$, where the first error is statistical and the second error represents an estimate of scaling violation and ${\cal O}(α^2)$ errors in the renormalization factor added in quadrature, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is consistent, albeit somewhat small, with $B_K(μ= 2 {GeV})= 0.628(42)$ obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for light quark masses are also reported.

preprint2001arXiv

Maximum entropy analysis of hadron spectral functions and excited states in quenched lattice QCD

Employing the maximum entropy method we extract the spectral functions from meson correlators at four lattice spacings in quenched QCD with the Wilson quark action. We confirm that the masses and decay constants, obtained from the position and the area of peaks, agree well with the results from the conventional exponential fit. For the first excited state, we obtain $m_{π_1} = 660(590)$ MeV, $m_{ρ_1} = 1540(570)$ MeV, and $f_{ρ_1} = 0.085(36)$ in the continuum limit.

preprint2001arXiv

Spectral function and excited states in lattice QCD with maximum entropy method

We apply the maximum entropy method to extract the spectral functions for pseudoscalar and vector mesons from hadron correlators previously calculated at four different lattice spacings in quenched QCD with the Wilson quark action. We determine masses and decay constants for the ground and excited states of the pseudoscalar and vector channels from position and area of peaks in the spectral functions. We obtain the results, $m_{π_1} = 660(590)$ MeV and $m_{ρ_1} = 1540(570)$ MeV for the masses of the first excited state masses, in the continuum limit of quenched QCD. We also find unphysical states which have infinite mass in the continuum limit, and argue that they are bound states of two doublers of the Wilson quark action. If the interpretation is correct, this is the first time that the state of doublers is identified in lattice QCD numerical simulations.

preprint2000arXiv

Calculation of $K\to π$ matrix elements in quenched domain-wall QCD

We explore the possibility for an evaluation of non-leptonic $ΔS=1$ $K$ decay amplitudes through the calculation of $K\toπ$ matrix elements using domain-wall QCD. The relation between the physical $K\toππ$ matrix elements and $K\toπ$ matrix elements deduced from chiral perturbation theory is recapitulated. Quenched numerical simulations are performed on an $16^3\times 32\times 16$ lattice at a lattice spacing $a^{-1}\approx 2$GeV for both the standard plaquette gauge action and a renormalization-group improved gauge action, and reasonable signals for $K\toπ$ matrix elements are obtained. Preliminary results are reported on the $K\toππ$ matrix elements, and results from two actions are compared.

preprint2000arXiv

Eigenvalues of the hermitian Wilson-Dirac operator and chiral properties of the domain-wall fermion

Chiral properties of QCD formulated with the domain-wall fermion (DWQCD) are studied using the anomalous quark mass m_{5q} and the spectrum of the 4-dimensional Wilson-Dirac operator. Numerical simulations are made with the standard plaquette gauge action and a renormalization-group improved gauge action. Results are reported on the density of zero eigenvalue obtained with the accumulation method, and a comparison is made with the results for m_{5q}.

preprint2000arXiv

Erratum: Excitation Mechanism of Near-Infrared [Fe II] Emission in Seyfert and Starburst Galaxies (ApJ, 528, 186 [2000])

In the paper "Excitation Mechanism of Near-Infrared [Fe II] Emission in Seyfert and Starburst Galaxies" by Hideaki Mouri, Kimiaki Kawara, and Yoshiaki Taniguchi (ApJ, 528, 186 [2000]), the two panels of Figure 6 were printed into different sizes, as the result of an error in the printing process. The correct version of Figure 6 appears here.

preprint2000arXiv

Kaon B parameter from quenched domain-wall QCD

We report on a calculation of $B_K$ with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at $β=2.6$ and 2.9 corresponding to $a^{-1}\approx 2$GeV and 3GeV. Effects due to finite fifth dimensional size $N_5$ and finite spatial size $N_σ$ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain $B_K(μ= 2 {GeV})= 0.5746(61)$, where the error is statistical only, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is smaller but consistent with $B_K(μ= 2 {GeV})= 0.628(42)$ obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for strange quark mass are also reported.

preprint2000arXiv

Quenched charmonium spectrum on anisotropic lattices

We present the results of quenched charmonium spectrum for S- and P-states, obtained by a relativistic heavy quark method on anisotropic lattices. Simulations are carried out using the standard plaquette gauge action and a meanfield-improved clover quark action at $a_t^{-1} = 3$--6 GeV with the renormalized anisotropy fixed to $ξ\equiv a_s/a_t =3$. We study the scaling of our fine and hyperfine mass splittings, and compare with previous results.

preprint2000arXiv

Sub-parsec-scale Accleration of the Radio Jet in the Powerful Radio Galaxy NGC 6251

In order to investigate the genesis of powerful radio jet, we have mapped the central 10 pc region of the nearby radio galaxy NGC 6251 with a 0.2 pc resolution using Very Long Baseline Interferometer (VLBI) at two radio frequencies, 5 GHz and 15 GHz, we have found the sub-parsec-scale counterjet for the first time in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic beaming model.

preprint1999arXiv

Lattice QCD with Domain-Wall Fermions

We study the quenched lattice QCD using domain-wall fermions at $β=6.0$. Behaviors of both pion mass and the explicit breaking term in the axial Ward-Takahashi identity support the existence of the chiral zero modes. We observe a good agreement between the pion decay constants $f_π$ from both the conserved axial current and the local current perturbatively renormalized at 1-loop. Finally the possible existence of the parity broken phase is also examined in this model.

preprint1999arXiv

Quenched QCD with domain-wall fermions on coarse lattices

We investigate the existence of chiral zero modes at a^{-1} \simeq 1 GeV in quenched domain-wall QCD. Simulations are carried out for the plaquette and an RG-improved gauge actions on a 12^3x24xN_s lattice with N_s=10-50. We find that the pion mass in the chiral limit remains non-vanishing as N_s\to\infty for both gauge actions. Possible origins of this non-vanishing pion mass are discussed.

preprint1998arXiv

A Geometrical Relationship between Broad-Line Clouds and an Accretion Disk around Active Galactic Nuclei

Recent hard X-ray spectroscopy of active galactic nuclei has strongly suggested that double-peaked, very broad Fe K emission arises from an accretion disk around the central engine. Model fitting of the observed Fe K emission line profile makes it possible to estimate a probable inclination angle of the accretion disk. In order to study the geometrical relationship between the accretion disk and broad emission-line regions (BLRs), we investigate the correlation between the inclination angle of the accretion disk and the velocity width of BLRs for 18 type-1 Seyfert galaxies. We found that there may be a negative correlation between them, i.e., Seyfert nuclei with a more face-on accretion disk tend to have larger BLR velocity widths, suggesting that the BLRs are not coplanar with respect to the accretion disk. The most probable interpretation may be that the BLRs arise from outer parts ({\it r} $\sim$ 0.01 pc) of a warped accretion disk illuminated by the central engine.

preprint1998arXiv

Dusty star forming galaxies at high redshift

The global star formation rate in high redshift galaxies, based on optical surveys, shows a strong peak at a redshift of z=1.5, which implies that we have already seen most of the formation. High redshift galaxies may, however, emit most of their energy at submillimeter wavelengths if they contain substantial amounts of dust. The dust would absorb the starlight and reradiate it as far-infrared light, which would be redshifted to the submillimeter range. Here we report a deep survey of two blank regions of sky performed at submillimeter wavelengths (450 and 850-micron). If the sources we detect in the 850-micron band are powered by star formation, then each must be converting more than 100 solar masses of gas per year into stars, which is larger than the maximum star formation rates inferred for most optically-selected galaxies. The total amount of high redshift star formation is essentially fixed by the level of background light, but where the peak occurs in redshift for the submillimeter is not yet established. However, the background light contribution from only the sources detected at 850-micron is already comparable to that from the optically-selected sources. Establishing the main epoch of star formation will therefore require a combination of optical and submillimeter studies.