Source author record

D. B. Sanders

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.HE astro-ph.IM

Catalog footprint

What is connected

92works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1).

preprint2022arXiv

COSMOS2020: Manifold Learning to Estimate Physical Parameters in Large Galaxy Surveys

We present a novel method to estimate galaxy physical properties from spectral energy distributions (SEDs), alternate to template fitting techniques and based on self-organizing maps (SOM) to learn the high-dimensional manifold of a photometric galaxy catalog. The method has been previously tested with hydrodynamical simulations in Davidzon et al. (2019) while here is applied to real data for the first time. It is crucial for its implementation to build the SOM with a high quality, panchromatic data set, which we elect to be the "COSMOS2020" galaxy catalog. After the training and calibration steps with COSMOS2020, other galaxies can be processed through SOM to obtain an estimate of their stellar mass and star formation rate (SFR). Both quantities result to be in good agreement with independent measurements derived from more extended photometric baseline, and also their combination (i.e., the SFR vs. stellar mass diagram) shows a main sequence of star forming galaxies consistent with previous studies. We discuss the advantages of this method compared to traditional SED fitting, highlighting the impact of having, instead of the usual synthetic templates, a collection of empirical SEDs built by the SOM in a "data-driven" way. Such an approach also allows, even for extremely large data sets, an efficient visual inspection to identify photometric errors or peculiar galaxy types. Considering in addition the computational speed of this new estimator, we argue that it will play a valuable role in the analysis of oncoming large-area surveys like Euclid or the Legacy Survey of Space and Time at the Vera Cooper Rubin Telescope.

preprint2022arXiv

COSMOS2020: The cosmic evolution of the stellar-to-halo mass relation for central and satellite galaxies up to z~5

We use the COSMOS2020 catalogue to measure the stellar-to-halo mass relation (SHMR) divided by central and satellite galaxies from $z=0.2$ to $z = 5.5$. Starting from accurate photometric redshifts we measure the near-infrared selected two-point angular correlation and stellar mass functions in ten redshift bins and fit them with an HOD-based model. At each redshift, we measure the ratio of stellar mass to halo mass, $M_*/M_h$, which shows the characteristic strong dependence of halo mass with a peak at $M_h^{\rm peak} \sim 10^{12}\, M_{\odot}$. Our results are in accordance with the scenario in which the peak of star-formation efficiency moves towards more massive halos at higher redshifts. We also measure the fraction of satellites as a function of stellar mass and redshift. For all stellar mass thresholds the satellite fraction decreases at higher redshifts. At a given redshift there is a higher fraction of low-mass satellites. The satellite contribution to the total stellar mass budget in halos becomes more important than centrals at halo masses of about $M_h > 10^{13} \, M_{\odot}$ and always stays below by peak, indicating that quenching mechanisms are present in massive halos that keep the star-formation efficiency low. Finally, we compare our results with three hydrodynamical simulations Horizon-AGN, Illustris-TNG-100 and EAGLE. We find that the most significant discrepancy is at the high mass end, where the simulations generally show that satellites have a higher contribution to the total stellar mass budget than the observations. This, together with the finding that the fraction of satellites is higher in the simulations, indicates that the feedback mechanisms acting in group-and cluster-scale halos appear to be less efficient in quenching the mass assembly of satellites, and/or that quenching occurs much later in the simulations.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

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

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2020arXiv

The Molecular Gas in the NGC 6240 Merging Galaxy System at the Highest Spatial Resolution

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

preprint2019arXiv

Star-forming Clumps in Local Luminous Infrared Galaxies

We present $HST$ narrow-band near-infrared imaging of Pa$α$ and Pa$β$ emission of 48 local Luminous Infrared Galaxies (LIRGs) from the Great Observatories All-Sky LIRG Survey (GOALS). These data allow us to measure the properties of 810 spatially resolved star-forming regions (59 nuclei and 751 extra-nuclear clumps), and directly compare their properties to those found in both local and high-redshift star-forming galaxies. We find that in LIRGs, the star-forming clumps have radii ranging from $\sim90-900$ pc and star formation rates (SFRs) of $\sim1\times10^{-3}$ to 10 M$_\odot$yr$^{-1}$, with median values for extra-nuclear clumps of 170 pc and 0.03 M$_\odot$yr$^{-1}$. The detected star-forming clumps are young, with a median stellar age of $8.7$ Myrs, and a median stellar mass of $5\times10^{5}$ M$_\odot$. The SFRs span the range of those found in normal local star-forming galaxies to those found in high-redshift star-forming galaxies at $\rm{z}=1-3$. The luminosity function of the LIRG clumps has a flatter slope than found in lower-luminosity, star-forming galaxies, indicating a relative excess of luminous star-forming clumps. In order to predict the possible range of star-forming histories and gas fractions, we compare the star-forming clumps to those measured in the MassiveFIRE high-resolution cosmological simulation. The star-forming clumps in MassiveFIRE cover the same range of SFRs and sizes found in the local LIRGs and have total gas fractions that extend from 10 to 90%. If local LIRGs are similar to these simulated galaxies, we expect future observations with ALMA will find a large range of gas fractions, and corresponding star formation efficiencies, among the star-forming clumps in LIRGs.

preprint2016arXiv

NuSTAR unveils a heavily obscured low-luminosity Active Galactic Nucleus in the Luminous Infrared Galaxy NGC 6286

We report the detection of a heavily obscured Active Galactic Nucleus (AGN) in the luminous infrared galaxy (LIRG) NGC 6286, identified in a 17.5 ks NuSTAR observation. The source is in an early merging stage, and was targeted as part of our ongoing NuSTAR campaign observing local luminous and ultra-luminous infrared galaxies in different merger stages. NGC 6286 is clearly detected above 10 keV and, by including the quasi-simultaneous Swift/XRT and archival XMM-Newton and Chandra data, we find that the source is heavily obscured [$N_{\rm\,H}\simeq (0.95-1.32)\times 10^{24}\rm\,cm^{-2}$], with a column density consistent with being Compton-thick [CT, $\log (N_{\rm\,H}/\rm cm^{-2})\geq 24$]. The AGN in NGC 6286 has a low absorption-corrected luminosity ($L_{2-10\rm\,keV}\sim 3-20\times 10^{41}\rm\,erg\,s^{-1}$) and contributes $\lesssim$1\% to the energetics of the system. Because of its low-luminosity, previous observations carried out in the soft X-ray band ($<10$ keV) and in the infrared did not notice the presence of a buried AGN. NGC 6286 has multi-wavelength characteristics typical of objects with the same infrared luminosity and in the same merger stage, which might imply that there is a significant population of obscured low-luminosity AGN in LIRGs that can only be detected by sensitive hard X-ray observations.

preprint2015arXiv

A massive, distant proto-cluster at z=2.47 caught in a phase of rapid formation?

Numerical simulations of cosmological structure formation show that the Universe's most massive clusters, and the galaxies living in those clusters, assemble rapidly at early times (2.5 < z < 4). While more than twenty proto-clusters have been observed at z > 2 based on associations of 5-40 galaxies around rare sources, the observational evidence for rapid cluster formation is weak. Here we report observations of an asymmetric, filamentary structure at z = 2.47 containing seven starbursting, submillimeter-luminous galaxies and five additional AGN within a comoving volume of 15000 Mpc$^{3}$. As the expected lifetime of both the luminous AGN and starburst phase of a galaxy is ~100 Myr, we conclude that these sources were likely triggered in rapid succession by environmental factors, or, alternatively, the duration of these cosmologically rare phenomena is much longer than prior direct measurements suggest. The stellar mass already built up in the structure is $\sim10^{12}M_{\odot}$ and we estimate that the cluster mass will exceed that of the Coma supercluster at $z \sim 0$. The filamentary structure is in line with hierarchical growth simulations which predict that the peak of cluster activity occurs rapidly at z > 2.

preprint2015arXiv

A Turnover in the Galaxy Main Sequence of Star Formation at $M_{*} \sim 10^{10} M_{\odot}$ for Redshifts $z < 1.3$

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $μ$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

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

Kinematic classifications of local interacting galaxies: implications for the merger/disk classifications at high-z

The classification of galaxy mergers and isolated disks is key for understanding the relative importance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic properties of galaxies as traced by emission lines have been used to suggest the existence of a significant population of high-z star-forming galaxies consistent with isolated rotating disks. However, recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. To further investigate the robustness of merger/disk classifications based on kinematic properties, we carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially redshift the WiFeS observations of these local (U)LIRGs to z=1.5 to make a realistic comparison with observations at high-z, and also to ensure that all galaxies have the same spatial sampling of ~900 pc. Using both kinemetry-based and visual classifications, we find that the reliability of kinematic classification shows a strong trend with the interaction stage of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties compared to isolated disks, whereas a significant population of the interacting disks and merger remnants are indistinguishable from isolated disks. The high fraction of late-stage mergers showing disk-like kinematics reflects the complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified disks and mergers depend on the definition of kinematic asymmetries and the classification threshold when using kinemetry-based classifications. Our results suggest that additional indicators such as morphologies traced by stars or molecular gas are required to further constrain the merger/disk classifications at high-z.

preprint2015arXiv

Merger Signatures in the Dynamics of Star-forming Gas

Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion.

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 Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

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

A comparison of the morphological properties between local and z~1 infrared luminous galaxies. Are local and high-z (U)LIRGs different?

Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.

preprint2014arXiv

ALMA Observations of Warm Dense Gas in NGC 1614 --- Breaking of Star Formation Law in the Central kpc

We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)% of the total CO(6-5) flux measured by Herschel. The molecular ring, located between 100pc < r < 350pc from the nucleus, looks clumpy and includes seven unresolved (or marginally resolved) knots with median velocity dispersion of 40 km/s. These knots are associated with strong star formation regions with Σ_{SFR} 100 M_\sun/yr/kpc^{2} and Σ_{Gas} 1.0E4 M_\sun/pc^{2}. The non-detections of the nucleus in both the CO (6-5) line emission and the 435um continuum rule out, with relatively high confidence, a Compton-thick AGN in NGC 1614. Comparisons with radio continuum emission show a strong deviation from an expected local correlation between Σ_{Gas} and Σ_{SFR}, indicating a breakdown of the Kennicutt-Schmidt law on the linear scale of 100 pc.

preprint2014arXiv

ALMA Observations of Warm Molecular Gas and Cold Dust in NGC 34

We present ALMA Cycle-0 observations of the CO (6-5) line emission (rest-frame frequency = 691.473 GHz) and of the 435$μm$ dust continuum emission in the nuclear region of NGC 34, a local luminous infrared galaxy (LIRG) at a distance of 84 Mpc (1" = 407 pc) which contains a Seyfert 2 active galactic nucleus (AGN) and a nuclear starburst. The CO emission is well resolved by the ALMA beam ($\rm 0."26\times 0."23$), with an integrated flux of $\rm f_{CO~(6-5)} = 1004\; (\pm 151) \; Jy\; km\; s^{-1}$. Both the morphology and kinematics of the CO (6-5) emission are rather regular, consistent with a compact rotating disk with a size of 200 pc. A significant emission feature is detected on the red-shifted wing of the line profile at the frequency of the $\rm H^{13}CN\; (8-7)$ line, with an integrated flux of $\rm 17.7 \pm 2.1 (random) \pm 2.7 (sysmatic)\; Jy\;km\; s^{-1}$. However, it cannot be ruled out that the feature is due to an outflow of warm dense gas with a mean velocity of $\rm 400\; km\; s^{-1}$. The continuum is resolved into an elongated configuration, and the observed flux corresponds to a dust mass of $\rm M_{dust} = 10^{6.97\pm 0.13}\; M_{sun}$. An unresolved central core ($\rm radius \simeq 50\; pc$) contributes $28\%$ of the continuum flux and $19\%$ of the CO (6-5) flux, consistent with insignificant contributions of the AGN to both emissions. Both the CO (6-5) and continuum spatial distributions suggest a very high gas column density ($\rm >= 10^4\; M_{sun}\; pc^{-2}$) in the nuclear region at $\rm radius <= 100\; pc$.

preprint2014arXiv

Are Dusty Galaxies Blue? Insights on UV Attenuation from Dust-Selected Galaxies

Galaxies' rest-frame ultraviolet (UV) properties are often used to directly infer the degree to which dust obscuration affects the measurement of star formation rates. While much recent work has focused on calibrating dust attenuation in galaxies selected at rest-frame ultraviolet wavelengths, locally and at high-$z$, here we investigate attenuation in dusty, star-forming galaxies (DSFGs) selected at far-infrared wavelengths. By combining multiwavelength coverage across 0.15--500\,$μ$m in the COSMOS field, in particular making use of {\it Herschel} imaging, and a rich dataset on local galaxies, we find a empirical variation in the relationship between rest-frame UV slope ($β$) and ratio of infrared-to-ultraviolet emission ($L_{\rm IR}/L_{\rm UV}\equiv\,IRX$) as a function of infrared luminosity, or total star formation rate, SFR. Both locally and at high-$z$, galaxies above SFR$\gt$50\,M$_\odot$\,yr$^{-1}$ deviate from the nominal $IRX-β$ relation towards bluer colors by a factor proportional to their increasing IR luminosity. We also estimate contamination rates of DSFGs on high-$z$ dropout searches of $\ll1$\%\ at $z\lt4-10$, providing independent verification that contamination from very dusty foreground galaxies is low in LBG searches. Overall, our results are consistent with the physical interpretation that DSFGs, e.g. galaxies with $>50$\,M$_\odot$\,yr$^{-1}$, are dominated at all epochs by short-lived, extreme burst events, producing many young O and B stars that are primarily, yet not entirely, enshrouded in thick dust cocoons. The blue rest-frame UV slopes of DSFGs are inconsistent with the suggestion that most DSFGs at $z\sim2$ exhibit steady-state star formation in secular disks.

preprint2014arXiv

Extended [CII] Emission in Local Luminous Infrared Galaxies

We present Herschel/PACS observations of extended [CII]157.7μm line emission detected on ~ 1 - 10 kpc scales in 60 local luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS). We find that most of the extra-nuclear emission show [CII]/FIR ratios >~ 4 x 10^-3, larger than the mean ratio seen in the nuclei, and similar to those found in the extended disks of normal star-forming galaxies and the diffuse inter-stellar medium (ISM) of our Galaxy. The [CII] "deficits" found in the most luminous local LIRGs are therefore restricted to their nuclei. There is a trend for LIRGs with warmer nuclei to show larger differences between their nuclear and extra-nuclear [CII]/FIR ratios. We find an anti-correlation between [CII]/FIR and the luminosity surface density, Σ_IR, for the extended emission in the spatially-resolved galaxies. However, there is an offset between this trend and that found for the LIRG nuclei. We use this offset to derive a beam filling-factor for the star-forming regions within the LIRG disks of ~ 6 % relative to their nuclei. We confront the observed trend to photo-dissociation region (PDR) models and find that the slope of the correlation is much shallower than the model predictions. Finally, we compare the correlation found between [CII]/FIR and Σ_IR with measurements of high-redshift starbursting IR-luminous galaxies.

preprint2014arXiv

Is There a Maximum Star Formation Rate in High-Redshift Galaxies?

We use the James Clerk Maxwell Telescope's SCUBA-2 camera to image a 400 arcmin^2 area surrounding the GOODS-N field. The 850 micron rms noise ranges from a value of 0.49 mJy in the central region to 3.5 mJy at the outside edge. From these data, we construct an 850 micron source catalog to 2 mJy containing 49 sources detected above the 4-sigma level. We use an ultradeep (11.5 uJy at 5-sigma) 1.4 GHz image obtained with the Karl G. Jansky Very Large Array together with observations made with the Submillimeter Array to identify counterparts to the submillimeter galaxies. For most cases of multiple radio counterparts, we can identify the correct counterpart from new and existing Submillimeter Array data. We have spectroscopic redshifts for 62% of the radio sources in the 9 arcmin radius highest sensitivity region (556/894) and 67% of the radio sources in the GOODS-N region (367/543). We supplement these with a modest number of additional photometric redshifts in the GOODS-N region (30). We measure millimetric redshifts from the radio to submillimeter flux ratios for the unidentified submillimeter sample, assuming an Arp 220 spectral energy distribution. We find a radio flux dependent K-z relation for the radio sources, which we use to estimate redshifts for the remaining radio sources. We determine the star formation rates (SFRs) of the submillimeter sources based on their radio powers and their submillimeter and find that they agree well. The radio data are deep enough to detect star-forming galaxies with SFRs >2000 solar masses per year to z~6. We find galaxies with SFRs up to ~6,000 solar masses per year over the redshift range z=1.5-6, but we see evidence for a turn-down in the SFR distribution function above 2000 solar masses per year.

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

Warm Molecular Gas in Luminous Infrared Galaxies

We present our initial results on the CO rotational spectral line energy distribution (SLED) of the $J$ to $J$$-$1 transitions from $J=4$ up to $13$ from Herschel SPIRE spectroscopic observations of 65 luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). The observed SLEDs change on average from one peaking at $J \le 4$ to a broad distribution peaking around $J \sim\,$6$-$7 as the IRAS 60-to-100 um color, $C(60/100)$, increases. However, the ratios of a CO line luminosity to the total infrared luminosity, $L_{\rm IR}$, show the smallest variation for $J$ around 6 or 7. This suggests that, for most LIRGs, ongoing star formation (SF) is also responsible for a warm gas component that emits CO lines primarily in the mid-$J$ regime ($5 \lesssim J \lesssim 10$). As a result, the logarithmic ratios of the CO line luminosity summed over CO (5$-$4), (6$-$5), (7$-$6), (8$-$7) and (10$-$9) transitions to $L_{\rm IR}$, $\log R_{\rm midCO}$, remain largely independent of $C(60/100)$, and show a mean value of $-4.13$ ($\equiv \log R^{\rm SF}_{\rm midCO}$) and a sample standard deviation of only 0.10 for the SF-dominated galaxies. Including additional galaxies from the literature, we show, albeit with small number of cases, the possibility that galaxies, which bear powerful interstellar shocks unrelated to the current SF, and galaxies, in which an energetic active galactic nucleus contributes significantly to the bolometric luminosity, have their $R_{\rm midCO}$ higher and lower than $R^{\rm SF}_{\rm midCO}$, respectively.

preprint2013arXiv

A Herschel Survey of the [N II] 205 micron Line in Local Luminous Infrared Galaxies --- The [N II] 205 micron Emission as a Star Formation Rate Indicator

We present, for the first time, a statistical study of [N II] 205 mciron line emission for a large sample of local luminous infrared galaxies using Herschel Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE FTS) data. For our sample of galaxies, we investigate the correlation between the [N II] luminosity (LNII) and the total infrared luminosity (LIR), as well as the dependence of LNII/LIR ratio on LIR, far infrared colors (IRAS $f_{60}/f_{100}$) and the [O III] 88 micron to [N II] luminosity ratio. We find that LNII correlates almost linearly with LIR for non AGN galaxies (all having $L_{IR} < 10^{12} L_solar$) in our sample, which implies that LNII can serve as a SFR tracer which is particularly useful for high redshift galaxies which will be observed with forthcoming submm spectroscopic facilities such as the Atacama Large Millimeter/submillimeter Array. Our analysis shows that the deviation from the mean LNII-LIR relation correlates with tracers of the ionization parameter, which suggests the scatter in this relation is mainly due to the variations in the hardness, and/or ionization parameter, of the ambient galactic UV field among the sources in our sample.

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

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

Explaining the [CII]158um Deficit in Luminous Infrared Galaxies - First Results from a Herschel/PACS Study of the GOALS Sample

We present the first results of a survey of the [CII]158um emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board Herschel. The [CII] luminosities of the LIRGs in GOALS range from ~10^7 to 2x10^9 Lsun. We find that LIRGs show a tight correlation of [CII]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ~10^-2 to 10^-4, as the average temperature of dust increases. We find correlations between the [CII]/FIR ratio and the strength of the 9.7um silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (Sigma_MIR), suggesting that warmer, more compact starbursts have substantially smaller [CII]/FIR ratios. Pure star-forming (SF) LIRGs have a mean [CII]/FIR ~ 4x10^-3, while galaxies with low 6.2um PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [CII]/FIR. However, we show that even when only pure SF galaxies are considered, the [CII]/FIR ratio drops by an order of magnitude, from 10^-2 to 10^-3, with Sigma_MIR and Sigma_IR, implying that the [CII] luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [CII]/FIR >= 10^-3 typical of high 6.2um PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [CII]/FIR and the specific SFR (SSFR) for SF LIRGs. Finally, we present predictions for the starburst size based on the observed [CII] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.

preprint2013arXiv

Hubble Space Telescope ACS Imaging of the GOALS Sample: Quantitative Structural Properties of Nearby Luminous Infrared Galaxies with L_IR > 10^{11.4} L_sun

A {\it Hubble Space Telescope} ({\it HST}) / Advanced Camera for Surveys (ACS) study of the structural properties of 85 luminous and ultraluminous ($L_{\rm IR} > 10^{11.4}$ L$_\odot$) infrared galaxies (LIRGs and ULIRGs) in the Great Observatories All-sky LIRG Survey (GOALS) sample is presented. Two-dimensional GALFIT analysis has been performed on F814W "{\it I}-band" images to decompose each galaxy, as appropriate, into bulge, disk, central PSF and stellar bar components. The fraction of bulge-less disk systems is observed to be higher in LIRGs (35%) than in ULIRGs (20%), with the disk+bulge systems making up the dominant fraction of both LIRGs (55%) and ULIRGs (45%). Further, bulge+disk systems are the dominant late-stage merger galaxy type and are the dominant type for LIRGs and ULIRGs at almost every stage of galaxy-galaxy nuclear separation. The mean {\it I}-band host absolute magnitude of the GOALS galaxies is $-22.64\pm$0.62 mag (1.8$^{+1.4}_{-0.4}$ L$^*_I$), and the mean bulge absolute magnitude in GOALS galaxies is about 1.1 magnitude fainter than the mean host magnitude. Almost all ULIRGs have bulge magnitudes at the high end (-20.6 to - 23.5 mag) of the GOALS bulge magnitude range. Mass ratios in the GOALS binary systems are consistent with most of the galaxies being the result of major mergers, and an examination of the residual-to-host intensity ratios in GOALS binary systems suggests that smaller companions suffer more tidal distortion than the larger companions. We find approximately twice as many bars in GOALS disk+bulge systems (32.8%) than in pure-disk mergers (15.9%) but most of the disk+bulge systems that contain bars are disk-dominated with small bulges. The bar-to-host intensity ratio, bar half-light radius, and bar ellipticity in GOALS galaxies are similar to those found in nearby spiral galaxies.(abridged)

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

Mid-Infrared Properties of Nearby Luminous Infrared Galaxies I: Spitzer IRS Spectra for the GOALS Sample

The Great Observatories All-Sky LIRG Survey (GOALS) is a multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here we present low resolution Spitzer spectra covering 5-38um and provide a basic analysis of the mid-IR spectral properties for nearby LIRGs. In a companion paper, we discuss detailed fits to the spectra. The GOALS sample of 244 nuclei in 180 luminous and 22 ultraluminous IR galaxies represents a complete subset of the IRAS RBGS and covers a range of merger stages, morphologies and spectral types. The majority (>60%) of GOALS LIRGs have high 6.2um PAH equivalent widths (EQW > 0.4um) and low levels of silicate absorption (s_9.7um >-1.0). There is a general trend among the U/LIRGs for silicate depth and MIR slope to increase with LIR. U/LIRGs in the late stages of a merger also have on average steeper MIR slopes and higher levels of dust obscuration. Together these trends suggest that as gas & dust is funneled towards the center of a coalescing merger, the nuclei become more compact and obscured. The sources that depart from these correlations have very low PAH EQW (EQW < 0.1um) consistent with their MIR emission being dominated by an AGN. The most heavily dust obscured sources are the most compact in their MIR emission, suggesting that the obscuring (cool) dust is associated with the outer regions of the starburst. As the merger progresses a marked decline is seen for the fraction of high EQW (star formation dominated) sources while the fraction of composite sources increases but the fraction of AGN-dominated sources remains low. When compared to the MIR spectra of submillimeter galaxies (SMGs) at z~2, the average GOALS LIRG is more absorbed at 9.7um and has more PAH emission. However, when the AGN contributions to both the local LIRGs and the high-z SMGs are removed, the average local starbursting LIRG closely resembles the starbursting SMGs.

preprint2013arXiv

Multi-Wavelength SEDs of Herschel Selected Galaxies in the COSMOS Field

We combine Herschel PACS and SPIRE maps of the full 2 deg^2 COSMOS field with existing multi-wavelength data to obtain template and model-independent optical-to-far-infrared spectral energy distributions (SEDs) for 4,218 Herschel-selected sources with log(L_IR/L_sun) = 9.4-13.6 and z = 0.02-3.54. Median SEDs are created by binning the optical to far-infrared (FIR) bands available in COSMOS as a function of infrared luminosity. Herschel probes rest-frame wavelengths where the bulk of the infrared radiation is emitted, allowing us to more accurately determine fundamental dust properties of our sample of infrared luminous galaxies. We find that the SED peak wavelength (lambda_peak) decreases and the dust mass (M_dust) increases with increasing total infrared luminosity (L_IR). In the lowest infrared luminosity galaxies (log(L_IR/L_sun) = 10.0-11.5), we see evidence of Polycyclic Aromatic Hydrocarbons (PAH) features (lambda=7-9 um), while in the highest infrared luminosity galaxies (L_IR > 10^12 L_sun) we see an increasing contribution of hot dust and/or power-law emission, consistent with the presence of heating from an active galactic nucleus (AGN). We study the relationship between stellar mass and star formation rate of our sample of infrared luminous galaxies and find no evidence that Herschel-selected galaxies follow the SFR/M_* "main sequence" as previously determined from studies of optically selected, star-forming galaxies. Finally, we compare the mid-infrared (MIR) to FIR properties of our infrared luminous galaxies using the previously defined diagnostic, IR8 = L_IR / L_8, and find that galaxies with L_IR > 10^11.3 L_sun tend to systematically lie above (~3-5) the IR8 "infrared main sequence", suggesting either suppressed PAH emission or an increasing contribution from AGN heating.

preprint2013arXiv

The Build-Up of Nuclear Stellar Cusps in Extreme Starburst Galaxies and Major Mergers

Nuclear stellar cusps are defined as central excess light component in the stellar light profiles of galaxies and are suggested to be stellar relics of intense compact starbursts in the central ~100-500 pc region of gas-rich major mergers. Here we probe the build-up of nuclear cusps during the actual starburst phase for a complete sample of Luminous Infrared Galaxy systems (85 LIRGs, with 11.4<log[LIR/L_sun]<12.5) in the GOALS sample. Cusp properties are derived via 2-dimensional fitting of the nuclear stellar light imaged in the near-infrared by the Hubble Space Telescope and have been combined with mid-IR diagnostics for AGN/starburst characterization. We find that nuclear stellar cusps are resolved in 76% of LIRGs (merger and non-interacting galaxies). The cusp strength and luminosity increases with far-IR luminosity (excluding AGN) and merger stage, confirming theoretical models that starburst activity is associated with the build-up of nuclear stellar cusps. Evidence for ultra compact nuclear starbursts is found in ~13% of LIRGs, which have a strong unresolved central NIR light component but no significant contribution of an AGN. The nuclear near-IR surface density (measured within 1 kpc radius) increases by a factor of ~5 towards late merger stages. A careful comparison to local early-type galaxies with comparable masses reveals (a) that local (U)LIRGs have a significantly larger cusp fraction and (b) that the majority of the cusp LIRGs have host galaxy luminosities (H-band) similar to core ellipticals which is roughly one order in magnitude larger than for cusp ellipticals.

preprint2013arXiv

Unveiling The Sigma-Discrepancy II: Revisiting the Evolution of ULIRGs & The Origin of Quasars

We present the first central velocity dispersions (sigma_o) measured from the 0.85 micron Calcium II Triplet (CaT) for 8 advanced (i.e. single nuclei) local (z < 0.15) Ultraluminous Infrared Galaxies (ULIRGs). First, these measurements are used to test the prediction that the "sigma-Discrepancy," in which the CaT sigma_o is systematically larger than the sigma_o obtained from the 1.6 or 2.3 micron stellar CO band-heads, extends to ULIRG luminosities. Next, we combine the CaT data with rest-frame I-band photometry obtained from archival Hubble Space Telescope data and the Sloan Digital Sky Survey (SDSS) to derive dynamical properties for the 8 ULIRGs. These are then compared to the dynamical properties of 9,255 elliptical galaxies from the SDSS within the same redshift volume and of a relatively nearby (z < 0.4) sample of 53 QSO host galaxies. A comparison is also made between the I-band and H-band dynamical properties of the ULIRGs. We find four key results: 1) the sigma-Discrepancy extends to ULIRG luminosities; 2) at I-band ULIRGs lie on the Fundamental Plane (FP) in a region consistent with the most massive elliptical galaxies and not low-intermediate mass ellipticals as previously reported in the near-infrared; 3) the I-band M/L of ULIRGs are consistent with an old stellar population, while at H-band ULIRGs appear significantly younger and less massive; and 4) we derive an I-band Kormendy Relation from the SDSS ellipticals and demonstrate that ULIRGs and QSO host galaxies are dynamically similar.

preprint2012arXiv

A Population of z>2 Far-Infrared Herschel-SPIRE selected Starbursts

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

preprint2012arXiv

A Redshift Survey of Herschel Far-Infrared Selected Starbursts and Implications for Obscured Star Formation

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

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

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

preprint2012arXiv

Identifying Luminous AGN in Deep Surveys: Revised IRAC Selection Criteria

Spitzer IRAC selection is a powerful tool for identifying luminous AGN. For deep IRAC data, however, the AGN selection wedges currently in use are heavily contaminated by star-forming galaxies, especially at high redshift. Using the large samples of luminous AGN and high-redshift star-forming galaxies in COSMOS, we redefine the AGN selection criteria for use in deep IRAC surveys. The new IRAC criteria are designed to be both highly complete and reliable, and incorporate the best aspects of the current AGN selection wedges and of infrared power-law selection while excluding high redshift star-forming galaxies selected via the BzK, DRG, LBG, and SMG criteria. At QSO-luminosities of log L(2-10 keV) (ergs/s) > 44, the new IRAC criteria recover 75% of the hard X-ray and IRAC-detected XMM-COSMOS sample, yet only 38% of the IRAC AGN candidates have X-ray counterparts, a fraction that rises to 52% in regions with Chandra exposures of 50-160 ks. X-ray stacking of the individually X-ray non-detected AGN candidates leads to a hard X-ray signal indicative of heavily obscured to mildly Compton-thick obscuration (log N_H (cm^-2) = 23.5 +/- 0.4). While IRAC selection recovers a substantial fraction of luminous unobscured and obscured AGN, it is incomplete to low-luminosity and host-dominated AGN.

preprint2012arXiv

Investigation of Dual Active Nuclei, Outflows, Shock-Heated Gas, and Young Star Clusters in Markarian 266

Results of observations with the Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space telescopes are presented for the Luminous Infrared Galaxy (LIRG) merger Mrk 266. The SW (Seyfert 2) and NE (LINER) nuclei reside in galaxies with Hubble types SBb (pec) and S0/a (pec), respectively. Both galaxies have L > L*, and they are inferred to each contain a ~2.5x10^8 M_sun black hole. Mrk 266 SW is likely the primary source of a bright Fe K-alpha line detected from the system, consistent with the reflection-dominated X-ray spectrum of a heavily obscured AGN. Optical knots embedded in an arc with aligned radio continuum radiation, combined with luminous H_2 line emission, provide evidence for a radiative bow shock in an AGN-driven outflow surrounding the NE nucleus. Soft X-ray emission modeled as shock-heated plasma is co-spatial with radio continuum emission between the galaxies. Mid-IR diagnostics indicate roughly equal contributions of AGN and starburst radiation powering the bolometric luminosity. Approximately 120 star clusters have been detected, with most having estimated ages < 50 Myr. Detection of 24 micron emission extending ~34 arcsec (20 kpc) north of the galaxies is interpreted as ~2x10^7 M_sun of dust entrained in an outflowing superwind; at optical wavelengths this region is resolved into a fragmented morphology indicative of Rayleigh-Taylor instabilities in an expanding shell of ionized gas. Mrk 266 demonstrates that the dust "blow-out" phase can begin in a LIRG well before the galaxies fully coalesce during a subsequent ULIRG phase, and rapid gas consumption in luminous dual AGNs (kpc scale separations) early in the merger process may explain the paucity of detected binary QSOs (sub-pc scale orbits) in large surveys. An evolutionary sequence is proposed representing a progression from dual to binary AGNs, accompanied by an increase in observed L_x/L_ir ratios by 10^4 or more.

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 WIRCam Deep Survey I: Counts, colours and mass-functions derived from near-infrared imaging in the CFHTLS deep fields

We present a new near-infrared imaging survey in the four CFHTLS deep fields: the WIRCam Deep Survey (WIRDS). WIRDS comprises extremely deep, high quality (FWHM ~0.6") J, H and K imaging covering a total effective area of 2.1 deg^2 and reaching AB 50% completeness limits of ~24.5. We combine our images with the CFHTLS to create a unique eight-band ugrizJHK photometric catalogues in the CFHTLS deep fields; these four separate fields allow us to make a robust estimate of the effect of cosmic variance for all our measurements. We use these catalogues to estimate precise photometric redshifts, galaxy types and stellar masses for a unique sample of ~1.8 million galaxies. Our JHK number counts are consistent with previous studies. We apply the BzK selection to our gzK filter set and find that the star forming BzK selection successfully selects 76% of star-forming galaxies in the redshift range 1.4<z<2.5 in our photometric catalogue. The passive BzK selection returns 52% of the passive 1.4<z<2.5 population identified in the photometric catalogue. We present the galaxy stellar mass function as a function of redshift up to z=2 and present fits using double Schechter functions. A mass-dependent evolution of the mass function is seen with the numbers of galaxies with masses of log(M)<10.75 still evolving at z<1, but galaxies of higher mass reaching their present day numbers by z~0.8-1. This is consistent with the present picture of downsizing in galaxy evolution. We compare our results with the predictions of the GALFORM semi-analytical galaxy formation model and find that the simulations provide a relatively successful fit to the observed mass functions at intermediate masses (i.e. 10<log(M)<11). However, the GALFORM results under-predict the mass function at low masses, whilst the fit as a whole degrades beyond redshifts of z~1.2.

preprint2011arXiv

C-GOALS: Chandra observations of a complete sample of luminous infrared galaxies from the IRAS Revised Bright Galaxy Survey

We present X-ray data for a complete sample of 44 luminous infrared galaxies (LIRGs), obtained with the Chandra X-ray Observatory. These are the X-ray observations of the high luminosity portion of the Great Observatory All-sky LIRG Survey (GOALS), which includes the most luminous infrared selected galaxies, log (Lir/Lsun) > 11.73, in the local universe, z < 0.088. X-rays were detected from 43 out of 44 objects, and their arcsec-resolution images, spectra, and radial brightness distributions are presented. With a selection by hard X-ray colour and the 6.4 keV iron line, AGN are found in 37% of the objects, with higher luminosity sources more likely to contain an AGN. These AGN also tend to be found in late-stage mergers. The AGN fraction would increase to 48% if objects with mid-IR [Ne V] detection are included. Double AGN are clearly detected only in NGC 6240 among 24 double/triple systems. Other AGN are found either in single nucleus objects or in one of the double nuclei at similar rates. Objects without conventional X-ray signatures of AGN appear to be hard X-ray quiet, relative to the X-ray to far-IR correlation for starburst galaxies, as discussed elsewhere. Most objects also show extended soft X-ray emission, which is likely related to an outflow from the nuclear region, with a metal abundance pattern suggesting enrichment by core collapse supernovae, as expected for a starburst.

preprint2011arXiv

Chandra Discovery of a Binary Active Galactic Nucleus in Mrk 739

We have discovered a binary AGN in the galaxy Mrk 739 using Chandra and Swift BAT. We find two luminous L_2-10 keV=1.1*10^43 and 1.0*10^42 erg/s, unresolved nuclei with a projected separation of 3.4 kpc (5.8\pm0.1) coincident with two bulge components in the optical image. The western X-ray source (Mrk 739W) is highly variable (2.5x) during the 4-hour Chandra observation and has a very hard spectrum consistent with an AGN. While the eastern component was already known to be an AGN based on the presence of broad optical recombination lines, Mrk 739W shows no evidence of being an AGN in optical, UV, and radio observations, suggesting the critical importance of high spatial resolution hard X-ray observations (>2 keV) in finding these binary AGN. A high level of star formation combined with a very low L_[O III]/L_2-10 keV ratio cause the AGN to be missed in optical observations. CO observations of the (3-2) and (2-1) lines indicate large amounts of molecular gas in the system that could be driven towards the black holes during the violent galaxy collision and be key to fueling the binary AGN. Mrk 739E has a high Eddington ratio of 0.71 and a small black hole (log M_BH=7.05\pm0.3) consistent with an efficiently accreting AGN. Other than NGC 6240, this stands as the nearest case of a binary AGN discovered to date.

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

Multi-wavelength GOALS Observations of Star Formation and Active Galactic Nucleus Activity in the Luminous Infrared Galaxy IC 883

New optical HST, Spitzer, GALEX, and Chandra observations of the single-nucleus, luminous infrared galaxy (LIRG) merger IC 883 are presented. The galaxy is a member of the Great Observatories All-sky LIRG Survey (GOALS), and is of particular interest for a detailed examination of a luminous late-stage merger due to the richness of the optically-visible star clusters and the extended nature of the nuclear X-ray, mid-IR, CO and radio emission. In the HST ACS images, the galaxy is shown to contain 156 optically visible star clusters distributed throughout the nuclear regions and tidal tails of the merger, with a majority of visible clusters residing in an arc ~ 3-7 kpc from the position of the mid-infrared core of the galaxy. The luminosity functions of the clusters have an alpha_F435W ~ -2.17+/-0.22 and alpha_F814W ~ -2.01+/-0.21. Further, the colors and absolute magnitudes of the majority of the clusters are consistent with instantaneous burst population synthesis model ages in the range of a few x10^7 - 10^8 yrs (for 10^5 M_sun clusters), but may be as low as few x10^6 yrs with extinction factored in. The X-ray and mid-IR spectroscopy are indicative of predominantly starburst-produced nuclear emission, and the star formation rate is ~ 80 M_sun / yr. The kinematics of the CO emission and the morphology of both the CO and radio emission are consistent with the nuclear starburst being situated in a highly inclined disk 2 kpc in diameter with an infrared surface brightness mu_IR ~ 2x10^11 L_sun kpc^-2, a factor of 10 less than that of the Orion star-forming region. Finally, the detection of the [Ne V] 14.32 um emission line is evidence that an AGN is present. The faintness of the line (i.e., [Ne V] / [Ne II] um ~ 0.01) and the small equivalent width of the 6.2 um PAH feature ($= 0.39μ$m) are both indicative of a relatively weak AGN. (abridged)

preprint2011arXiv

Spectroscopy of luminous z>7 galaxy candidates and sources of contamination in z>7 galaxy searches

We present three bright z+ dropout candidates selected from deep Near-Infrared (NIR) imaging of the COSMOS 2 square degree field. All three objects match the 0.8-8um colors of other published z>7 candidates but are three magnitudes brighter, facilitating further study. Deep spectroscopy of two of the candidates covering 0.64-1.02um with Keck-DEIMOS and all three covering 0.94-1.10um and 1.52-1.80um with Keck-NIRSPEC detects weak spectral features tentatively identified as Ly-alpha at z=6.95 and z=7.69 in two of the objects. The third object is placed at z~1.6 based on a 24um and weak optical detection. A comparison with the spectral energy distributions of known z<7 galaxies, including objects with strong spectral lines, large extinction, and large systematic uncertainties in the photometry yields no objects with similar colors. However, the lambda>1um properties of all three objects can be matched to optically detected sources with photometric redshifts at z~1.8, so the non-detection in the i+ and z+ bands are the primary factors which favors a z>7 solution. If any of these objects are at z~7 the bright end of the luminosity function is significantly higher at z>7 than suggested by previous studies, but consistent within the statistical uncertainty and the dark matter halo distribution. If these objects are at low redshift, the Lyman-Break selection must be contaminated by a previously unknown population of low redshift objects with very strong breaks in their broad band spectral energy distributions and blue NIR colors. The implications of this result on luminosity function evolution at high redshift is discussed. We show that the primary limitation of z>7 galaxy searches with broad filters is the depth of the available optical data.

preprint2011arXiv

The location of an active nucleus and the soft X-ray shadowing by a tidal tail in the ULIRG Mrk 273

Analysis of data from the Chandra X-ray Observatory for the double nucleus ULIRG Mrk 273 reveals an absorbed hard X-ray source coincident with the southwest nucleus, implying that this unresolved near infrared source is where an active nucleus resides while the northern nuclear region contains a powerful starburst which dominates the far infrared luminosity. There is evidence of a slight image extension in the 6-7 keV band, where a Fe K line is present, towards the northern nucleus. A large-scale, diffuse emission nebula detected in soft X-rays contains a dark lane that spatially coincides with a high surface-brightness tidal tail extending ~50 arcsec (40 kpc) to the south. The soft X-ray source is likely located behind the tidal tail which absorbs X-ray photons along the line of sight. The estimated column density of cold gas in the tidal tail responsible for shadowing the soft X-rays is nH >= 6e+21 cm-2, consistent with the tidal tail having an edge-on orientation.

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 Far-IR Characterization of 24 micron Selected Galaxies at 0<z<2.5 using Stacking at 70 microns and 160 microns in the COSMOS Field

We present a study of the average properties of luminous infrared galaxies detected directly at 24 $μ$m in the COSMOS field using a median stacking analysis at 70$μ$m and 160 $μ$m. Over 35000 sources spanning 0<z<3 and 0.06 mJy<S_{24}<3.0 mJy are stacked, divided into bins of both photometric redshift and 24 $μ$m flux. We find no correlation of $S_{70}/S_{24}$ flux density ratio with $S_{24}$, but find that galaxies with higher $S_{24}$ have a lower $S_{160}/S_{24}$ flux density ratio. These observed ratios suggest that 24 $μ$m selected galaxies have warmer SEDs at higher mid-IR fluxes, and therefore have a possible higher fraction of AGN. Comparisons of the average $S_{70}/S_{24}$ and $S_{160}/S_{24}$ colors with various empirical templates and theoretical models show that the galaxies detected at 24 $μ$m are consistent with "normal" star-forming galaxies and warm mid-IR galaxies such as Mrk 231, but inconsistent with heavily obscured galaxies such as Arp 220. We perform a $χ^{2}$ analysis to determine best fit galactic model SEDs and total IR luminosities for each of our bins. We compare our results to previous methods of estimating $L_{\rm{IR}}$ and find that previous methods show considerable agreement over the full redshift range, except for the brightest $S_{24}$ sources, where previous methods overpredict the bolometric IR luminosity at high redshift, most likely due to their warmer dust SED. We present a table that can be used as a more accurate and robust method for estimating bolometric infrared luminosity from 24 $μ$m flux densities.

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

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

preprint2010arXiv

Herschel observations of water vapour in Markarian 231

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

preprint2010arXiv

No Evolution in the IR-Radio Relation for IR-Luminous Galaxies at z<2 in the COSMOS Field

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution.

preprint2010arXiv

Radio galaxy feedback in X-ray selected groups from COSMOS: the effect on the ICM

We quantify the importance of the mechanical energy released by radio-galaxies inside galaxy groups. We use scaling relations to estimate the mechanical energy released by 16 radio-AGN located inside X-ray detected galaxy groups in the COSMOS field. By comparing this energy output to the host groups' gravitational binding energy, we find that radio galaxies produce sufficient energy to unbind a significant fraction of the intra-group medium. This unbinding effect is negligible in massive galaxy clusters with deeper potential wells. Our results correctly reproduce the breaking of self-similarity observed in the scaling relation between entropy and temperature for galaxy groups.

preprint2010arXiv

The Buried Starburst in the Interacting Galaxy II Zw 096 as Revealed by the Spitzer Space Telescope

An analysis of data from the Spitzer Space Telescope, Hubble Space Telescope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is presented for the z=0.036 merging galaxy system II Zw 096 (CGCG 448-020). Because II Zw 096 has an infrared luminosity of log(L_IR/L_sun) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories All-sky LIRG Survey (GOALS). The Spitzer data suggest that 80% of the total infrared luminosity comes from an extremely compact, red source not associated with the nuclei of the merging galaxies. The Spitzer mid-infrared spectra indicate no high-ionization lines from a buried active galactic nucleus in this source. The strong detection of the 3.3 micron and 6.2 micron PAH emission features in the AKARI and Spitzer spectra also implies that the energy source of II Zw 096 is a starburst. Based on Spitzer infrared imaging and AKARI near-infrared spectroscopy, the star formation rate is estimated to be 120 M_sun/yr and > 45 M_sun/yr, respectively. Finally, the high-resolution B, I, and H-band images show many star clusters in the interacting system. The colors of these clusters suggest at least two populations - one with an age of 1-5 Myr and one with an age of 20-500 Myr, reddened by 0-2 magnitudes of visual extinction. The masses of these clusters span a range between 10^6-10^8 M_sun. This starburst source is reminiscent of the extra-nuclear starburst seen in NGC 4038/9 (the Antennae Galaxies) and Arp 299 but approximately an order of magnitude more luminous than the Antennae. The source is remarkable in that the off-nuclear infrared luminosity dominates the enitre system.

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 Nuclear Structure in Nearby Luminous Infrared Galaxies: HST NICMOS Imaging of the GOALS Sample

We present results of Hubble Space Telescope NICMOS H-band imaging of 73 of most luminous (i.e., log[L_IR/L_0]>11.4) Infrared Galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). This dataset combines multi-wavelength imaging and spectroscopic data from space (Spitzer, HST, GALEX, and Chandra) and ground-based telescopes. In this paper we use the high-resolution near-infrared data to recover nuclear structure that is obscured by dust at optical wavelengths and measure the evolution in this structure along the merger sequence. A large fraction of all galaxies in our sample possess double nuclei (~63%) or show evidence for triple nuclei (~6%). Half of these double nuclei are not visible in the HST B-band images due to dust obscuration. The majority of interacting LIRGs have remaining merger timescales of 0.3 to 1.3 Gyrs, based on the projected nuclear separations and the mass ratio of nuclei. We find that the bulge luminosity surface density increases significantly along the merger sequence (primarily due to a decrease of the bulge radius), while the bulge luminosity shows a small increase towards late merger stages. No significant increase of the bulge Sersic index is found. LIRGs that show no interaction features have on average a significantly larger bulge luminosity, suggesting that non merging LIRGs have larger bulge masses than merging LIRGs. This may be related to the flux limited nature of the sample and the fact that mergers can significantly boost the IR luminosity of otherwise low luminosity galaxies. We find that the projected nuclear separation is significantly smaller for ULIRGs (median value of 1.2 kpc) than for LIRGs (mean value of 6.7 kpc), suggesting that the LIRG phase appears earlier in mergers than the ULIRG phase.

preprint2010arXiv

The Spatial Extent of (U)LIRGs in the mid-Infrared I: The Continuum Emission

We present an analysis of the extended mid-infrared (MIR) emission of the Great Observatories All-Sky LIRG Survey (GOALS) sample based on 5-15um low resolution spectra obtained with the IRS on Spitzer. We calculate the fraction of extended emission as a function of wavelength for the galaxies in the sample, FEE_lambda. We can identify 3 general types of FEE_lambda: one where it is constant, one where features due to emission lines and PAHs appear more extended than the continuum, and a third which is characteristic of sources with deep silicate absorption at 9.7um. More than 30% of the galaxies have a median FEE_lambda larger than 0.5 implying that at least half of their MIR emission is extended. Luminous Infrared Galaxies (LIRGs) display a wide range of FEE in their warm dust continuum (0<=FEE_13.2um<=0.85). The large values of FEE_13.2um that we find in many LIRGs suggest that their extended MIR continuum emission originates in scales up to 10kpc. The mean size of the LIRG cores at 13.2um is 2.6kpc. However, once the LIR of the systems reaches the threshold of ~10^11.8Lsun, all sources become clearly more compact, with FEE_13.2um<=0.2, and their cores are unresolved. Our estimated upper limit for the core size of ULIRGs is less than 1.5kpc. The analysis indicates that the compactness of systems with LIR>~10^11.25Lsun strongly increases in those classified as mergers in their final stage of interaction. The FEE_13.2um is also related to the contribution of an active galactic nucleus (AGN) to the MIR. Galaxies which are more AGN-dominated are less extended, independently of their LIR. We finally find that the extent of the MIR continuum emission is correlated with the far-IR IRAS log(f_60um/f_100um) color. This enables us to place a lower limit to the area in a galaxy from where the cold dust emission may originate, a prediction which can be tested soon with the Herschel Space Telescope.

preprint2010arXiv

The VLA-COSMOS Perspective on the IR-Radio Relation. I. New Constraints on Selection Biases and the Non-Evolution of the IR/Radio Properties of Star Forming and AGN Galaxies at Intermediate and High Redshift

VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5<z<5. At redshift z<1.4 we observe that radio-quiet AGN populate the locus of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- selected populations and show that it can reconcile apparently discrepant measurements presented in the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity.

preprint2009arXiv

A Deep HST H-Band Imaging Survey of Massive Gas-Rich Mergers. II. The QUEST PG QSOs

We report the results from a deep HST NICMOS H-band imaging survey of 28 z < 0.3 QSOs from the Palomar-Green (PG) sample. This program is part of QUEST (Quasar / ULIRG Evolution STudy) and complements a similar set of data on 26 highly-nucleated ULIRGs presented in Paper I. Our analysis indicates that the fraction of QSOs with elliptical hosts is higher among QSOs with undetected far-infrared (FIR) emission, small infrared excess, and luminous hosts. The hosts of FIR-faint QSOs show a tendency to have less pronounced merger-induced morphological anomalies and larger QSO-to-host luminosity ratios on average than the hosts of FIR-bright QSOs, consistent with late-merger evolution from FIR-bright to FIR-faint QSOs. The spheroid sizes and total host luminosities of the radio-quiet PG QSOs in our sample are statistically indistinguishable from the ULIRG hosts presented in Paper I, while those of radio-loud PG QSOs are systematically larger and more luminous. ULIRGs and PG QSOs with elliptical hosts fall near, but not exactly on, the fundamental plane of inactive spheroids. We confirm the systematic trend noted in Paper I for objects with small (< 2 kpc) spheroids to be up to ~1 mag. brighter than inactive spheroids. The host colors and wavelength dependence of their sizes support the idea that these deviations are due at least in part to non-nuclear star formation. However, the amplitudes of these deviations does not depend on host R-H colors. Taken at face value (i.e., no correction for extinction or the presence of a young stellar population), the H-band spheroid-host luminosities imply BH masses ~5 -- 200 x 10^7 M_sun and sub-Eddington mass accretion rates for both QSOs and ULIRGs. These results are compared with published BH mass estimates derived from other methods. (abridged)

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

GOALS: The Great Observatories All-Sky LIRG Survey

The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full range of nuclear spectral types defined via traditional optical line-ratio diagrams as well as interaction stages. They provide an unbiased picture of the processes responsible for enhanced infrared emission in galaxies in the local Universe. As an example of the analytic power of the multi-wavelength GOALS dataset, we present data for the interacting system VV 340 (IRAS F14547+2449). Between 80-95% of the total far-infrared emission (or about 5E11 solar luminosities) originates in VV 340 North. While the IRAC colors of VV 340 North and South are consistent with star-forming galaxies, both the Spitzer IRS and Chandra ACIS data indicate the presence of a buried AGN in VV 340 North. The GALEX far and near-UV fluxes imply a extremely large infrared "excess" (IRX) for the system (IR/FUV = 81) which is well above the correlation seen in starburst galaxies. Most of this excess is driven by VV 340 N, which alone has an IR excess of nearly 400. The VV 340 system seems to be comprised of two very different galaxies - an infrared luminous edge-on galaxy (VV 340 North) that dominates the long-wavelength emission from the system and which hosts a buried AGN, and a face-on starburst (VV 340 South) that dominates the short-wavelength emission.

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

Spitzer 70 and 160-micron Observations of the COSMOS Field

We present Spitzer 70 and 160 micron observations of the COSMOS Spitzer survey (S-COSMOS). The data processing techniques are discussed for the publicly released products consisting of images and source catalogs. We present accurate 70 and 160 micron source counts of the COSMOS field and find reasonable agreement with measurements in other fields and with model predictions. The previously reported counts for GOODS-North and the extragalactic First Look Survey are updated with the latest calibration, and counts are measured based on the large area SWIRE survey to constrain the bright source counts. We measure an extragalactic confusion noise level of sigma_c = 9.4+/-3.3 mJy (q=5) for the MIPS 160-micron band based on the deep S-COSMOS data and report an updated confusion noise level of sigma_c = 0.35+/-0.15 mJy (q=5) for the MIPS 70-micron band.

preprint2009arXiv

Spitzer Quasar and ULIRG Evolution Study (QUEST). IV. Comparison of 1-Jy Ultraluminous Infrared Galaxies with Palomar-Green Quasars

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)

preprint2009arXiv

The Bimodal Galaxy Stellar Mass Function in the COSMOS Survey to z~1: A Steep Faint End and a New Galaxy Dichotomy

We present a new analysis of stellar mass functions (MF) in the COSMOS field to fainter limits than has been previously probed to z~1. Neither the total nor the passive or star-forming MF can be well fit with a single Schechter function once one probes below 3e9 Msun. We observe a dip or plateau at masses ~1e10 Msun, and an upturn towards a steep faint-end slope of -1.7 at lower mass at any z<1. This bimodal nature of the MF is not solely a result of the blue/red dichotomy. The blue MF is by itself bimodal at z~1. This suggests a new dichotomy in galaxy formation that predates the appearance of the red sequence. We propose two interpretations for this bimodality. If the gas fraction increases towards lower mass, galaxies with M_baryon~1e10 Msun would shift to lower stellar masses, creating the observed dip. This would indicate a change in star formation efficiency, perhaps linked to supernovae feedback becoming much more efficient. Therefore, we investigate whether the dip is present in the baryonic (stars+gas) MF. Alternatively, the dip could be created by an enhancement of the galaxy assembly rate at ~1e11 Msun, a phenomenon that naturally arises if the baryon fraction peaks at M_halo ~1e12 Msun. In this scenario, galaxies occupying the bump around M* would be identified with central galaxies and the second fainter component having a steep faint-end slope with satellites. While the dip is apparent in the total MF at any z, it appears to shift from the blue to red population, likely as a result of transforming high-mass blue galaxies into red ones. At the same time, we detect a drastic upturn in the number of low-mass red galaxies. Their increase with time reflects a decrease in the number of blue systems and so we tentatively associate them with satellite dwarf galaxies that have undergone quenching.

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 Role of Starburst-AGN composites in Luminous Infrared Galaxy Mergers: Insights from the New Optical Classification Scheme

We investigate the fraction of starbursts, starburst-AGN composites, Seyferts, and LINERs as a function of infrared luminosity (L_IR) and merger progress for ~500 infrared-selected galaxies. Using the new optical classifications afforded by the extremely large data set of the Sloan Digital Sky Survey, we find that the fraction of LINERs in IR-selected samples is rare (< 5%) compared with other spectral types. The lack of strong infrared emission in LINERs is consistent with recent optical studies suggesting that LINERs contain AGN with lower accretion rates than in Seyfert galaxies. Most previously classified infrared-luminous LINERs are classified as starburst-AGN composite galaxies in the new scheme. Starburst-AGN composites appear to "bridge" the spectral evolution from starburst to AGN in ULIRGs. The relative strength of the AGN versus starburst activity shows a significant increase at high infrared luminosity. In ULIRGs (L_IR >10^12 L_odot), starburst-AGN composite galaxies dominate at early--intermediate stages of the merger, and AGN galaxies dominate during the final merger stages. Our results are consistent with models for IR-luminous galaxies where mergers of gas-rich spirals fuel both starburst and AGN, and where the AGN becomes increasingly dominant during the final merger stages of the most luminous infrared objects.

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.

preprint2007arXiv

An unusually brilliant transient in the galaxy Messier 85

Historically, variable and transient sources have both surprised astronomers and provided new views of the heavens. Here we report the discovery of an optical transient in the outskirts of the lenticular galaxy Messier 85 in the Virgo Cluster. With a peak absolute R magnitude of -12 this event is distinctly brighter than novae, but fainter than type Ia supernovae (expected from a population of old stars in lenticular galaxies). Archival images of the field do not show a luminous star at that position with an upper limit of ~-4.1, so it is unlikely to be a giant eruption from a luminous blue variable star. Over a two month period the transient emitted radiation energy of almost 10E47 erg and subsequently faded in the opical sky. It is similar to, but more luminous at peak by a factor of 6 than, an enigmatic transient in the galaxy M31. A possible origin of M85 OT2006-1 is a stellar merger. If so, searches for similar events in nearby galaxies will not only allow to study the physics of hyper-Eddington sources, but also probe an important phase in the evolution of stellar binary systems.

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

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.

preprint2005arXiv

The HST Cosmos Project: Contribution from the Subaru Telescope

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

preprint2003arXiv

Imaging Large Scale Structure in the X-ray Sky

We present the first results from a wide solid angle, moderately deep {\it Chandra} survey of the Lockman Hole North-West region. Our 9 ACIS-I fields cover an effective solid angle of 0.4 deg$^{2}$ and reach a depth of $3 \times 10^{-16}$ \ergpcmsqps in the 0.4--2 keV band and $3 \times 10^{-15}$ \ergpcmsqps in the 2--8 keV band. The best fit logN-logS for the entire field, the largest contiguous {\it Chandra} field yet observed, matches well onto that of the {\it Chandra} Deep Field North. We show that the full range of the `cosmic variance' previously seen in different {\it Chandra} fields is reproduced in this small region of the sky. Counts-in-cells analysis shows that the hard band sources are more strongly correlated than the soft band sources.

preprint2003arXiv

The IRAS Revised Bright Galaxy Sample (RBGS)

IRAS flux densities, redshifts, and infrared luminosities are reported for all sources identified in the IRAS Revised Bright Galaxy Sample (RBGS), a complete flux-limited survey of all extragalactic objects with total 60 micron flux density greater than 5.24 Jy, covering the entire sky surveyed by IRAS at Galactic latitude |b| > 5 degrees. The RBGS includes 629 objects, with a median (mean) sample redshift of 0.0082 (0.0126) and a maximum redshift of 0.0876. The RBGS supersedes the previous two-part IRAS Bright Galaxy Samples, which were compiled before the final ("Pass 3") calibration of the IRAS Level 1 Archive in May 1990. The RBGS also makes use of more accurate and consistent automated methods to measure the flux of objects with extended emission. Basic properties of the RBGS sources are summarized, including estimated total infrared luminosities, as well as updates to cross-identifications with sources from optical galaxy catalogs established using the NASA/IPAC Extragalactic Database (NED). In addition, an atlas of images from the Digitized Sky Survey with overlays of the IRAS position uncertainty ellipse and annotated scale bars is provided for ease in visualizing the optical morphology in context with the angular and metric size of each object. The revised bolometric infrared luminosity function, phi(L_ir), for infrared bright galaxies in the local Universe remains best fit by a double power law, phi(L_ir) ~ L_ir^alpha, with alpha = -0.6 (+/- 0.1), and alpha = -2.2 (+/- 0.1) below and above the "characteristic" infrared luminosity L_ir ~ 10^{10.5} L_solar, respectively. (Abridged)

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.

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.

D. B. Sanders

What is connected

Connect this record

See the researcher in context

Building this map preview

92 published item(s)

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

COSMOS2020: Manifold Learning to Estimate Physical Parameters in Large Galaxy Surveys

COSMOS2020: The cosmic evolution of the stellar-to-halo mass relation for central and satellite galaxies up to z~5

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

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

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

The Molecular Gas in the NGC 6240 Merging Galaxy System at the Highest Spatial Resolution

Star-forming Clumps in Local Luminous Infrared Galaxies

NuSTAR unveils a heavily obscured low-luminosity Active Galactic Nucleus in the Luminous Infrared Galaxy NGC 6286

A massive, distant proto-cluster at z=2.47 caught in a phase of rapid formation?

A Turnover in the Galaxy Main Sequence of Star Formation at $M_{*} \sim 10^{10} M_{\odot}$ for Redshifts $z < 1.3$

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

Kinematic classifications of local interacting galaxies: implications for the merger/disk classifications at high-z

Merger Signatures in the Dynamics of Star-forming Gas

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

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

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

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

A comparison of the morphological properties between local and z~1 infrared luminous galaxies. Are local and high-z (U)LIRGs different?

ALMA Observations of Warm Dense Gas in NGC 1614 --- Breaking of Star Formation Law in the Central kpc

ALMA Observations of Warm Molecular Gas and Cold Dust in NGC 34

Are Dusty Galaxies Blue? Insights on UV Attenuation from Dust-Selected Galaxies

Extended [CII] Emission in Local Luminous Infrared Galaxies

Is There a Maximum Star Formation Rate in High-Redshift Galaxies?

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

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

Warm Molecular Gas in Luminous Infrared Galaxies

A Herschel Survey of the [N II] 205 micron Line in Local Luminous Infrared Galaxies --- The [N II] 205 micron Emission as a Star Formation Rate Indicator

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

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

Explaining the [CII]158um Deficit in Luminous Infrared Galaxies - First Results from a Herschel/PACS Study of the GOALS Sample

Hubble Space Telescope ACS Imaging of the GOALS Sample: Quantitative Structural Properties of Nearby Luminous Infrared Galaxies with L_IR > 10^{11.4} L_sun

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

Mid-Infrared Properties of Nearby Luminous Infrared Galaxies I: Spitzer IRS Spectra for the GOALS Sample

Multi-Wavelength SEDs of Herschel Selected Galaxies in the COSMOS Field

The Build-Up of Nuclear Stellar Cusps in Extreme Starburst Galaxies and Major Mergers

Unveiling The Sigma-Discrepancy II: Revisiting the Evolution of ULIRGs & The Origin of Quasars

A Population of z>2 Far-Infrared Herschel-SPIRE selected Starbursts

A Redshift Survey of Herschel Far-Infrared Selected Starbursts and Implications for Obscured Star Formation

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

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

Identifying Luminous AGN in Deep Surveys: Revised IRAC Selection Criteria

Investigation of Dual Active Nuclei, Outflows, Shock-Heated Gas, and Young Star Clusters in Markarian 266

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

The WIRCam Deep Survey I: Counts, colours and mass-functions derived from near-infrared imaging in the CFHTLS deep fields

C-GOALS: Chandra observations of a complete sample of luminous infrared galaxies from the IRAS Revised Bright Galaxy Survey

Chandra Discovery of a Binary Active Galactic Nucleus in Mrk 739

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

Fe K emission from active galaxies in the COSMOS field

Multi-wavelength GOALS Observations of Star Formation and Active Galactic Nucleus Activity in the Luminous Infrared Galaxy IC 883

Spectroscopy of luminous z>7 galaxy candidates and sources of contamination in z>7 galaxy searches

The location of an active nucleus and the soft X-ray shadowing by a tidal tail in the ULIRG Mrk 273

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

A Far-IR Characterization of 24 micron Selected Galaxies at 0<z<2.5 using Stacking at 70 microns and 160 microns in the COSMOS Field

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

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

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

Bars in early- and late-type disks in COSMOS

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

Herschel observations of water vapour in Markarian 231

No Evolution in the IR-Radio Relation for IR-Luminous Galaxies at z<2 in the COSMOS Field

Radio galaxy feedback in X-ray selected groups from COSMOS: the effect on the ICM

The Buried Starburst in the Interacting Galaxy II Zw 096 as Revealed by the Spitzer Space Telescope

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

The Nuclear Structure in Nearby Luminous Infrared Galaxies: HST NICMOS Imaging of the GOALS Sample

The Spatial Extent of (U)LIRGs in the mid-Infrared I: The Continuum Emission

The VLA-COSMOS Perspective on the IR-Radio Relation. I. New Constraints on Selection Biases and the Non-Evolution of the IR/Radio Properties of Star Forming and AGN Galaxies at Intermediate and High Redshift

A Deep HST H-Band Imaging Survey of Massive Gas-Rich Mergers. II. The QUEST PG QSOs

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

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