Source author record

Peter Capak

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

52 published item(s)

preprint2020arXiv

ALMA Characterises the Dust Temperature of z ~ 5.5 Star-Forming Galaxies

The infrared spectral energy distributions (SEDs) of main-sequence galaxies in the early universe (z > 4) is currently unconstrained as infrared continuum observations are time consuming and not feasible for large samples. We present Atacama Large Millimetre Array (ALMA) Band 8 observations of four main-sequence galaxies at z ~ 5.5 to study their infrared SED shape in detail. Our continuum data (rest-frame 110$\rm μm$, close to the peak of infrared emission) allows us to constrain luminosity weighted dust temperatures and total infrared luminosities. With data at longer wavelengths, we measure for the first time the emissivity index at these redshifts to provide more robust estimates of molecular gas masses based on dust continuum. The Band 8 observations of three out of four galaxies can only be reconciled with optically thin emission redward of rest-frame 100$\rm μm$. The derived dust peak temperatures at z ~ 5.5 (38$\pm$8K) are elevated compared to average local galaxies, however, 5-10K below what would be predicted from an extrapolation of the trend at $z<4$. This behaviour can be explained by decreasing dust abundance (or density) towards high redshifts, which would cause the infrared SED at the peak to be more optically thin, making hot dust more visible to the external observer. From the 850$\rm μm$ dust continuum, we derive molecular gas masses between $10^{10}$ and $10^{11}\,{\rm M_{\odot}}$ and gas fractions (gas over total mass) of 30-80% (gas depletion times of 100-220Myrs). All in all, our results provide a first measured benchmark SED to interpret future millimetre observations of normal, main-sequence galaxies in the early Universe.

preprint2020arXiv

Bridging between the integrated and resolved main sequence of star formation

The position of galaxies on the stellar mass, star formation rate plane with respect to the star-forming main sequence at each redshift is a convenient way to infer where the galaxy is in its evolution compared to the rest of the population. We use Hubble Space Telescope high resolution images in the GOODS-S field from the CANDELS survey and fit multi wavelength lights in resolution elements of galaxies with stellar population synthesis models. We then construct resolved kpc-scale stellar mass, star formation rate surface density curves for galaxies at z~1. Fitting these resolved main sequence curves with Schechter functions, we parameterize and explain the multi-wavelength structure of galaxies with three variables: phi*, alpha, and M*. For quenched galaxies below the main sequence, we find an average high mass slope (alpha) of the resolved main sequence curves to be ~ -0.4. The scatter of this slope is higher among the lower mass star forming galaxies and those above the main sequence compared to quenched galaxies, due to lack of an evolved bulge. Our findings agree well with an inside-out quenching of star-formation. We find that the knee of the Schechter fits (M*) for galaxies below the main sequence occurs at lower stellar mass surface densities compared to star forming galaxies, which hints at how far quenching has proceeded outwards.

preprint2020arXiv

Spectroscopic Confirmation of a Coma Cluster Progenitor at z ~ 2.2

We report the spectroscopic confirmation of a new protocluster in the COSMOS field at $z$ $\sim$ 2.2, COSMOS Cluster 2.2 (CC2.2), originally identified as an overdensity of narrowband selected H$α$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $μ$m) and $K$ ($\sim$ 1.92-2.40 $μ$m) bands ($\sim$ 1.5 hour each), we confirm 35 unique protocluster members with at least two emission lines detected with S/N $>$ 3. Combined with 12 extra members from the zCOSMOS-deep spectroscopic survey (47 in total), we estimate a mean redshift and a line-of-sight velocity dispersion of $z_{mean}$=2.23224 $\pm$ 0.00101 and $σ_{los}$=645 $\pm$ 69 km s$^{-1}$ for this protocluster, respectively. Assuming virialization and spherical symmetry for the system, we estimate a total mass of $M_{vir}$ $\sim$ $(1-2) \times$10$^{14}$ $M_{\odot}$ for the structure. We evaluate a number density enhancement of $δ_{g}$ $\sim$ 7 for this system and we argue that the structure is likely not fully virialized at $z$ $\sim$ 2.2. However, in a spherical collapse model, $δ_{g}$ is expected to grow to a linear matter enhancement of $\sim$ 1.9 by $z$=0, exceeding the collapse threshold of 1.69, and leading to a fully collapsed and virialized Coma-type structure with a total mass of $M_{dyn}$($z$=0) $\sim$ 9.2$\times$10$^{14}$ $M_{\odot}$ by now. This observationally efficient confirmation suggests that large narrowband emission-line galaxy surveys, when combined with ancillary photometric data, can be used to effectively trace the large-scale structure and protoclusters at a time when they are mostly dominated by star-forming galaxies.

preprint2020arXiv

The ALPINE-ALMA [C II] survey: the luminosity function of serendipitous [C II] line emitters at $z\sim 5$

We present the first [CII] 158 $μ$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2019arXiv

Balmer Break Galaxy Candidates at $z \sim 6$: a Potential View on the Star-Formation Activity at $z \gtrsim 14$

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

preprint2016arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies IV. Physical properties derived from spectral energy distributions

We characterise the physical nature of a 1.1 mm-selected, flux-limited, and interferometrically followed up sample of SMGs in COSMOS. We used the MAGPHYS code to fit the multiwavelength (UV-radio) SEDs of 16 of the target SMGs. We also constructed the pure radio SEDs of our SMGs using three different radio bands (325 MHz, 1.4 GHz, and 3 GHz). Moreover, since two SMGs in our sample, AzTEC1 and AzTEC3, benefit from previous CO line observations, we studied their properties in more detail. We found that 63% of our target SMGs lie above the galaxy main-sequence by more than a factor of 3, and hence are starbursts. The 3 GHz radio sizes we have previously measured for the target SMGs were compared with the present stellar mass estimates, and we found that the z>3 SMGs are fairly consistent with the mass-size relationship of z~2 compact, quiescent galaxies (cQGs). The median IR-radio correlation parameter is found to be q=2.27, which is lower than measured locally (median q=2.64). AzTEC1 is found to have a sub-Eddington SFR surface density (by a factor of 2.6), while AzTEC3 appears to be an Eddington-limited starburster. The gas reservoir in these two high-z SMGs would be exhausted in only ~86 and 19 Myr at the current SFR, respectively. A comparison of the MAGPHYS-based properties of our SMGs with those of equally bright ALESS SMGs suggests that the two populations share fairly similar physical characteristics, including the q parameter. A hint of negative correlation is found between the 3 GHz size and the level of starburstiness, and hence cosmic-ray electrons in more compact starbursts might be more susceptible to free-free absorption. Some of the derived low and high q values (compared to the local median) could be the result of a specific merger/post-starburst phase of galaxy evolution. Overall, our results support the scenario where z>3 SMGs evolve into today's giant ellipticals.

preprint2016arXiv

A Tight Relation Between N/O Ratio and Galaxy Stellar Mass Can Explain the Evolution of Strong Emission Line Ratios with Redshift

The offset of high redshift star-forming galaxies in the [OIII]/H$β$ versus [NII]/H$α$ (O3N2) diagram in comparison with the local star-forming galaxy sequence is now well established. The physical origin of the shift is the subject of some debate, and has important implications for metallicity measurements based on strong lines at all redshifts. To investigate the origin of the O3N2 offset, we use a sample of ~100,000 star-forming galaxies from SDSS DR12 to probe the empirical correlations between emission line diagnostics and measurable galaxy physical properties. In particular, we examine how surface density of star formation, ionization parameter, nitrogen-to-oxygen (N/O) ratio, and stellar mass drive position in two key diagnostic diagrams: O3N2 and [OIII]/H$β$ versus [SII]/H$α$ (O3S2). We show that, at fixed [OIII]/H$β$, galaxies falling closer to the high-redshift locus in O3N2 have higher star formation density, stellar mass and N/O ratios. We also find a tight correspondence in the distributions of stellar mass and N/O in the diagnostic diagrams. This relation, spanning a range of galaxy evolutionary states, suggests that the N/O-$M_{*}$ relation is more fundamental than the N/O-metallicity relation. We argue that a tight N/O-$M_{*}$ relation is well-motivated physically, and that the observed correlation of N/O with O/H in the local universe is primarily a reflection of the existence of the mass-metallicity relation. Because the mass-metallicity relation evolves much more rapidly with redshift than N/O-$M_{*}$, the N/O ratios of high redshift galaxies are significantly elevated in comparison with local galaxies with the same gas-phase O/H. The O3N2 shift and elevated N/O ratios observed in high redshift galaxies therefore come about as a natural consequence of the N/O-$M_{*}$ relation combined with the evolution of the mass-metallicity relation.

preprint2016arXiv

An Empirical Approach to Cosmological Galaxy Survey Simulation: Application to SPHEREx Low-Resolution Spectroscopy

Highly accurate models of the galaxy population over cosmological volumes are necessary in order to predict the performance of upcoming cosmological missions. We present a data-driven model of the galaxy population constrained by deep 0.1-8 $\rm μm$ imaging and spectroscopic data in the COSMOS survey, with the immediate goal of simulating the spectroscopic redshift performance of the proposed SPHEREx mission. SPHEREx will obtain over the full-sky $R\sim41$ spectrophotometry at moderate spatial resolution ($\sim6"$) over the wavelength range 0.75-4.18 $\rm μm$ and $R\sim135$ over the wavelength range 4.18-5 $\rm μm$. We show that our simulation accurately reproduces a range of known galaxy properties, encapsulating the full complexity of the galaxy population and enables realistic, full end-to-end simulations to predict mission performance. Finally, we discuss potential applications of the simulation framework to future cosmology missions and give a description of released data products.

preprint2016arXiv

Effects of Local Environment and Stellar Mass on Galaxy Quenching out to z~3

We study the effects of local environment and stellar mass on galaxy properties using a mass complete sample of quiescent and star-forming systems in the COSMOS field at $z\lesssim$ 3. We show that at $z\lesssim$ 1, the median star-formation rate (SFR) and specific SFR (sSFR) of all galaxies depend on environment, but they become independent of environment at $z\gtrsim$ 1. However, we find that only for \textit{star-forming} galaxies, the median SFR and sSFR are similar in different environments, regardless of redshift and stellar mass. We find that the quiescent fraction depends on environment at $z\lesssim$ 1, and on stellar mass out to $z\sim$ 3. We show that at $z\lesssim$ 1, galaxies become quiescent faster in denser environments and that the overall environmental quenching efficiency increases with cosmic time. Environmental and mass quenching processes depend on each other. At $z\lesssim$ 1, denser environments more efficiently quench galaxies with higher masses (log($M/M_{\odot}$)$\gtrsim$ 10.7), possibly due to a higher merger rate of massive galaxies in denser environments, and that mass quenching is more efficient in denser regions. We show that the overall mass quenching efficiency ($ε_{mass}$) for more massive galaxies (log($M/M_{\odot}$)$\gtrsim$ 10.2) rises with cosmic time until $z\sim$ 1 and flattens out since then. However, for less massive galaxies, the rise in $ε_{mass}$ continues to the present time. Our results suggest that environmental quenching is only relevant at $z\lesssim$ 1, likely a fast process, whereas mass quenching is the dominant mechanism at $z\gtrsim$ 1, with a possible stellar feedback physics.

preprint2016arXiv

Large scale structure around a z=2.1 cluster

The most prodigious starburst galaxies are absent in massive galaxy clusters today, but their connection with large scale environments is less clear at $z\gtrsim2$. We present a search of large scale structure around a galaxy cluster core at $z=2.095$ using a set of spectroscopically confirmed galaxies. We find that both color-selected star-forming galaxies (SFGs) and dusty star-forming galaxies (DSFGs) show significant overdensities around the $z=2.095$ cluster. A total of 8 DSFGs (including 3 X-ray luminous active galactic nuclei, AGNs) and 34 SFGs are found within a 10 arcmin radius (corresponds to $\sim$15 cMpc at $z\sim2.1$) from the cluster center and within a redshift range of $Δz=0.02$, which leads to galaxy overdensities of $δ_{\rm DSFG}\sim12.3$ and $δ_{\rm SFG}\sim2.8$. The cluster core and the extended DSFG- and SFG-rich structure together demonstrate an active cluster formation phase, in which the cluster is accreting a significant amount of material from large scale structure while the more mature core may begin to virialize. Our finding of this DSFG-rich structure, along with a number of other protoclusters with excess DSFGs and AGNs found to date, suggest that the overdensities of these rare sources indeed trace significant mass overdensities. However, it remains puzzling how these intense star formers are triggered concurrently. Although an increased probability of galaxy interactions and/or enhanced gas supply can trigger the excess of DSFGs, our stacking analysis based on 850 $μ$m images and morphological analysis based on rest-frame optical imaging do not show such enhancements of merger fraction and gas content in this structure.

preprint2016arXiv

Metallicity and Age of the Stellar Stream Around the Disk Galaxy NGC 5907

Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron data with visible-light images allows us to use a long wavelength baseline to estimate the metallicity and age of the stellar population along a ~60 kpc long segment of the stream. We have fitted the stellar spectral energy distribution (SED) with a single-burst stellar population synthesis model and we use it to distinguish between the proposed satellite accretion and minor/major merger formation models of the stellar stream around this galaxy. We conclude that a massive minor merger (stellar mass ratio of at least 1:8) can best account for the metallicity of -0.3 inferred along the brightest parts of the stream.

preprint2016arXiv

Repeatability and Accuracy of Exoplanet Eclipse Depths Measured with Post-Cryogenic Spitzer

We examine the repeatability, reliability, and accuracy of differential exoplanet eclipse depth measurements made using the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope during the post-cryogenic mission. We have re-analyzed an existing 4.5 μm data set, consisting of 10 observations of the XO-3b system during secondary eclipse, using seven different techniques for removing correlated noise. We find that, on average, for a given technique, the eclipse depth estimate is repeatable from epoch to epoch to within 156 parts per million (ppm). Most techniques derive eclipse depths that do not vary by more than a factor 3 of the photon noise limit. All methods but one accurately assess their own errors: for these methods, the individual measurement uncertainties are comparable to the scatter in eclipse depths over the 10 epoch sample. To assess the accuracy of the techniques as well as to clarify the difference between instrumental and other sources of measurement error, we have also analyzed a simulated data set of 10 visits to XO-3b, for which the eclipse depth is known. We find that three of the methods (BLISS mapping, Pixel Level Decorrelation, and Independent Component Analysis) obtain results that are within three times the photon limit of the true eclipse depth. When averaged over the 10 epoch ensemble, 5 out of 7 techniques come within 60 ppm of the true value. Spitzer exoplanet data, if obtained following current best practices and reduced using methods such as those described here, can measure repeatable and accurate single eclipse depths, with close to photon-limited results.

preprint2016arXiv

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey: Report of a Community Workshop Examining Extragalactic, Galactic, Stellar and Planetary Science

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $μ$m [with R$\sim$41.4] and 4.18 and 5.00 $μ$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

preprint2016arXiv

The Impossibly Early Galaxy Problem

The current hierarchical merging paradigm and $Λ$CDM predict that the $z \sim 4-8$ universe should be a time in which the most massive galaxies are transitioning from their initial halo assembly to the later baryonic evolution seen in star-forming galaxies and quasars. However, no evidence of this transition has been found in many high redshift galaxy surveys including CFHTLS, CANDELS and SPLASH, the first studies to probe the high-mass end at these redshifts. Indeed, if halo mass to stellar mass ratios estimated at lower-redshift continue to $z \sim 6-8$, CANDELS and SPLASH report several orders of magnitude more $M \sim 10^{12-13} M_\odot$ halos than are possible to have formed by those redshifts, implying these massive galaxies formed impossibly early. We consider various systematics in the stellar synthesis models used to estimate physical parameters and possible galaxy formation scenarios in an effort to reconcile observation with theory. Although known uncertainties can greatly reduce the disparity between recent observations and cold dark matter merger simulations, even taking the most conservative view of the observations, there remains considerable tension with current theory.

preprint2016arXiv

The Local [CII] 158 um Emission Line Luminosity Function

We present, for the first time, the local [CII] 158 um emission line luminosity function measured using a sample of more than 500 galaxies from the Revised Bright Galaxy Sample (RBGS). [CII] luminosities are measured from the Herschel PACS observations of the Luminous Infrared Galaxies in the Great Observatories All-sky LIRG Survey (GOALS) and estimated for the rest of the sample based on the far-IR luminosity and color. The sample covers 91.3% of the sky and is complete at S_60 um > 5.24 Jy. We calculated the completeness as a function of [CII] line luminosity and distance, based on the far-IR color and flux densities. The [CII] luminosity function is constrained in the range ~10^(7-9) (Lo) from both the 1/V_max and a maximum likelihood methods. The shape of our derived [CII] emission line luminosity function agrees well with the IR luminosity function. For the CO(1-0) and [CII] luminosity functions to agree, we propose a varying ratio of [CII]/CO(1-0) as a function of CO luminosity, with larger ratios for fainter CO luminosities. Limited [CII] high redshift observations as well as estimates based on the IR and UV luminosity functions are suggestive of an evolution in the [CII] luminosity function similar to the evolution trend of the cosmic star formation rate density. Deep surveys using ALMA with full capability will be able to confirm this prediction.

preprint2015arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies - I. Multiwavelength identifications and redshift distribution

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4<S/N<=4.5 if a multiwavelength counterpart within 1.5" is also present, the total number of detections in our survey is 22. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or multiple components. Using optical to NIR photometric redshifts, available spectroscopic redshifts, and redshifts estimated from the radio-to-submm spectral index we infer a median redshift of $\tilde{z}=3.20\pm0.25$ for our sample. To study the overall multiplicity and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest AzTEC SMGs detected at 1.1 mm, AzTEC1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. We find that the median redshift for the 15 brightest AzTEC SMGs ($\tilde{z}=3.05\pm0.44$) is consistent with that for AzTEC16-30. This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux density. For the combined AzTEC1-30 sample we derive a median redshift of $\tilde{z}=3.17\pm0.27$, consistent with previous results based on mm-selected samples. We further infer that within the combined AzTEC1-30 sample $\sim25\pm9\%$ of sources separate into multiple components.

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

Cosmology with the SPHEREX All-Sky Spectral Survey

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<λ<$4.1$μ$m and R=150 spectra for 4.1$<λ<$4.8$μ$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

preprint2015arXiv

Mapping the Galaxy Color-Redshift Relation: Optimal Photometric Redshift Calibration Strategies for Cosmology Surveys

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

preprint2014arXiv

Dust Attenuation in High Redshift Galaxies -- 'Diamonds in the Sky'

We use observed optical to near infrared spectral energy distributions (SEDs) of 266 galaxies in the COSMOS survey to derive the wavelength dependence of the dust attenuation at high redshift. All of the galaxies have spectroscopic redshifts in the range z = 2 to 6.5. The presence of the CIV absorption feature, indicating that the rest-frame UV-optical SED is dominated by OB stars, is used to select objects for which the intrinsic, unattenuated spectrum has a well-established shape. Comparison of this intrinsic spectrum with the observed broadband photometric SED then permits derivation of the wavelength dependence of the dust attenuation. The derived dust attenuation curve is similar in overall shape to the Calzetti curve for local starburst galaxies. We also see the 2175 Å~bump feature which is present in the Milky Way and LMC extinction curves but not seen in the Calzetti curve. The bump feature is commonly attributed to graphite or PAHs. No significant dependence is seen with redshift between sub-samples at z = 2 - 4 and z = 4 - 6.5. The 'extinction' curve obtained here provides a firm basis for color and extinction corrections of high redshift galaxy photometry.

preprint2014arXiv

Inter-comparison of Radio-Loudness Criteria for Type 1 AGNs in the XMM-COSMOS Survey

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24μm}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500Å})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(νL_ν(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24μm}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.

preprint2014arXiv

Star Formation at $4 < z < 6$ From the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH)

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

preprint2014arXiv

The VLA-COSMOS Survey: V. 324 MHz continuum observations

We present 90 cm VLA imaging of the COSMOS field, comprising a circular area of 3.14 square degrees at 8.0"x6.0" angular resolution with an average rms of 0.5 mJy/beam. The extracted catalog contains 182 sources (down to 5.5sigma), 30 of which are multi-component sources. Using Monte Carlo artificial source simulations we derive the completeness of the catalog, and we show that our 90 cm source counts agree very well with those from previous studies. Using X-ray, NUV-NIR and radio COSMOS data to investigate the population mix of our 90 cm radio sample, we find that our sample is dominated by active galactic nuclei (AGN). The average 90-20 cm spectral index (S_nu~nu**alpha, where S_nu is the flux density at frequency nu, and alpha the spectral index) of our 90 cm selected sources is -0.70, with an interquartile range of -0.90 to -0.53. Only a few ultra-steep-spectrum sources are present in our sample, consistent with results in the literature for similar fields. Our data do not show clear steepening of the spectral index with redshift. Nevertheless, our sample suggests that sources with spectral indices steeper than -1 all lie at z>1, in agreement with the idea that ultra-steep-spectrum radio sources may trace intermediate-redshift galaxies (z>1).

preprint2014arXiv

Weighing the Giants II: Improved Calibration of Photometry from Stellar Colors and Accurate Photometric Redshifts

We present improved methods for using stars found in astronomical exposures to calibrate both star and galaxy colors as well as to adjust the instrument flat field. By developing a spectroscopic model for the SDSS stellar locus in color-color space, synthesizing an expected stellar locus, and simultaneously solving for all unknown zeropoints when fitting to the instrumental locus, we increase the calibration accuracy of stellar locus matching. We also use a new combined technique to estimate improved flat-field models for the Subaru SuprimeCam camera, forming `star flats' based on the magnitudes of stars observed in multiple positions or through comparison with available SDSS magnitudes. These techniques yield galaxy magnitudes with reliable color calibration (< 0.01 - 0.02 mag accuracy) that enable us to estimate photometric redshift probability distributions without spectroscopic training samples. We test the accuracy of our photometric redshifts using spectroscopic redshifts z_s for ~5000 galaxies in 27 cluster fields with at least five bands of photometry, as well as galaxies in the COSMOS field, finding sigma((z_p - z_s)/(1 + z_s)) ~ 0.03 for the most probable redshift z_p. We show that the full posterior probability distributions for the redshifts of galaxies with five-band photometry exhibit good agreement with redshifts estimated from thirty-band photometry in the COSMOS field. The growth of shear with increasing distance behind each galaxy cluster shows the expected redshift-distance relation for a flat Lambda-CDM cosmology. Photometric redshifts and calibrated colors are used in subsequent papers to measure the masses of 51 galaxy clusters from their weak gravitational shear. We make our Python code for stellar locus matching available at http://big-macs-calibrate.googlecode.com; the code requires only a catalog and filter functions.

preprint2013arXiv

A comparative analysis of virial black-hole mass estimates of moderate-luminosity active galactic nuclei using Subaru/FMOS

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

preprint2013arXiv

Characterisation of SCUBA-2 450um and 850um-selected Galaxies in the COSMOS Field

We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.

preprint2013arXiv

Encoding the infrared excess (IRX) in the NUVrK color diagram for star-forming galaxies

We present an empirical method of assessing the star formation rate (SFR) of star-forming galaxies based on their locations in the rest-frame color-color diagram (NUV-r) vs (r-K). By using the Spitzer 24 micron sample in the COSMOS field (~16400 galaxies with 0.2 < z < 1.3) and a local GALEX-SDSS-SWIRE sample (~700 galaxies with z < 0.2), we show that the mean infrared excess <IRX>= < L_IR / L_UV > can be described by a single vector, NRK, that combines the two colors. The calibration between <IRX> and NRK allows us to recover the IR luminosity, L_IR, with an accuracy of ~0.21 dex for the COSMOS sample and ~0.27 dex for the local one. The SFRs derived with this method agree with the ones based on the observed (UV+IR) luminosities and on the spectral energy distribution fitting for the vast majority (~85 %) of the star-forming population. Thanks to a library of model galaxy SEDs with realistic prescriptions for the star formation history, we show that we need to include a two-component dust model (i.e., birth clouds and diffuse ISM) and a full distribution of galaxy inclinations in order to reproduce the behavior of the <IRX> stripes in the NUVrK diagram. In conclusion, the NRK method, based only on rest-frame UV and optical colors available in most of the extragalactic fields, offers a simple alternative of assessing the SFR of star-forming galaxies in the absence of far-IR or spectral diagnostic observations.

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

Spectral Energy Distributions of Type 1 AGN in XMM-COSMOS Survey II - Shape Evolution

The mid-infrared to ultraviolet (0.1 -- 10 $μm$) spectral energy distribution (SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field ``Cosmic Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub-sample of 200 radio-quiet quasars with black hole mass estimates and host galaxy corrections, we studied their mean SEDs as a function of a broad range of redshift, bolometric luminosity, black hole mass and Eddington ratio, and compared them with the Elvis et al. (1994, E94) type 1 AGN mean SED. We found that the mean SEDs in each bin are closely similar to each other, showing no statistical significant evidence of dependence on any of the analyzed parameters. We also measured the SED dispersion as a function of these four parameters, and found no significant dependencies. The dispersion of the XMM-COSMOS SEDs is generally larger than E94 SED dispersion in the ultraviolet, which might be due to the broader ``window function'' for COSMOS quasars, and their X-ray based selection.

preprint2013arXiv

The role of galaxy interaction in the SFR-M relation: characterizing morphological properties of Herschel-selected galaxies at 0.2<z<1.5

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

preprint2013arXiv

Weak Lensing Calibrated M-T Scaling Relation of Galaxy Groups in the COSMOS Field

The scaling between X-ray observables and mass for galaxy clusters and groups is instrumental for cluster based cosmology and an important probe for the thermodynamics of the intracluster gas. We calibrate a scaling relation between the weak lensing mass and X-ray spectroscopic temperature for 10 galaxy groups in the COSMOS field, combined with 55 higher mass clusters from the literature. The COSMOS data includes HST imaging and redshift measurements of 46 source galaxies per square arcmin, enabling us to perform unique weak lensing measurements of low mass systems. Our sample extends the mass range of the lensing calibrated M-T relation an order of magnitude lower than any previous study, resulting in a power-law slope of 1.48$^{+0.13}_{-0.09}$. The slope is consistent with the self-similar model, predictions from simulations, and observations of clusters. However, X-ray observations relying on mass measurements derived under the assumption of hydrostatic equilibrium have indicated that masses at group scales are lower than expected. Both simulations and observations suggest that hydrostatic mass measurements can be biased low. Our external weak lensing masses provides the first observational support for hydrostatic mass bias at group level, showing an increasing bias with decreasing temperature and reaching a level of 30-50% at 1 keV.

preprint2012arXiv

Evolution of the Quasar Luminosity Function Over 3 < z < 5 in the COSMOS Survey Field

We investigate the high-redshift quasar luminosity function (QLF) down to an apparent magnitude of I(AB) = 25 in the Cosmic Evolution Survey (COSMOS). Careful analysis of the extensive COSMOS photometry and imaging data allows us to identify and remove stellar and low-redshift contaminants, enabling a selection that is nearly complete for type-1 quasars at the redshifts of interest. We find 155 likely quasars at z > 3.1, 39 of which have prior spectroscopic confirmation. We present our sample in detail and use these confirmed and likely quasars to compute the rest-frame UV QLF in the redshift bins 3.1 < z < 3.5 and 3.5 < z < 5. The space density of faint quasars decreases by roughly a factor of four from z \sim 3.2 to z \sim 4, with faint-end slopes of β \sim -1.7 at both redshifts. The decline in space density of faint optical quasars at z > 3 is similar to what has been found for more luminous optical and X-ray quasars. We compare the rest-frame UV luminosity functions found here with the X-ray luminosity function at z > 3, and find that they evolve similarly between z \sim 3.2 and z \sim 4; however, the different normalizations imply that roughly 75% of X-ray bright active galactic nuclei (AGN) at z \sim 3 - 4 are optically obscured. This fraction is higher than found at lower redshift and may imply that the obscured, type-2 fraction continues to increase with redshift at least to z \sim 4. Finally, the implications of the results derived here for the contribution of quasars to cosmic reionization are discussed.

preprint2012arXiv

Lyα Emission from High Redshift Sources in COSMOS

We investigate spectroscopically measured Lyα equivalent widths and escape fractions of 244 sources of which 95 are Lyman Break Galaxies (LBGs) and 106 Lyman Alpha Emitters (LAEs) at z~4.2, z~4.8, and z~5.6 selected from intermediate and narrow band observations. The sources were selected from the Cosmic Evolution Survey (COSMOS), and observed with the DEIMOS spectrograph. We find that the distribution of equivalent widths shows no evolution with redshift for both the LBG selected sources and the intermediate/narrowband LAEs. We also find that the Lyα escape fraction of intermediate and narrow band LAEs is on average higher and has a larger variation than the escape fraction of LBG selected sources. The escape fraction does not show a dependence with redshift. Similar to what has been found for LAEs at low redshifts, the sources with the highest extinctions show the lowest escape fractions. The range of escape fractions increases with decreasing extinction. This is evidence that the dust extinction is the most important factor affecting the escape of Lyα photons, but at low extinctions other factors such as HI covering fraction and gas kinematics can be just as effective at inhibiting the escape of Lyα photons.

preprint2012arXiv

Magnification by Galaxy Group Dark Matter Halos

We report on the detection of gravitational lensing magnification by a population of galaxy groups, at a significance level of 4.9 sigma. Using X-ray selected groups in the COSMOS 1.64 deg^2 field, and high-redshift Lyman break galaxies as sources, we measure a lensing-induced angular cross-correlation between the samples. After satisfying consistency checks that demonstrate we have indeed detected a magnification signal, and are not suffering from contamination by physical overlap of samples, we proceed to implement an optimally weighted cross-correlation function to further boost the signal to noise of the measurement. Interpreting this optimally weighted measurement allows us to study properties of the lensing groups. We model the full distribution of group masses using a composite-halo approach, considering both the singular isothermal sphere and Navarro-Frenk-White profiles, and find our best fit values to be consistent with those recovered using the weak-lensing shear technique. We argue that future weak-lensing studies will need to incorporate magnification along with shear, both to reduce residual systematics and to make full use of all available source information, in an effort to maximize scientific yield of the observations.

preprint2012arXiv

Near-Infrared Survey of the GOODS-North Field: Search for Luminous Galaxy Candidates at z=>6.5

We present near-infrared (NIR; J & Ks) survey of the Great Observatories Origins Deep Survey-North (GOODS-N) field. The publicly available imaging data were obtained using the MOIRCS instrument on the 8.2m Subaru and the WIRCam instrument on the 3.6m Canada-France-Hawaii Telescope (CFHT). These observations fulfill a serious wavelength gap in the GOODS-N data - i.e., lack of deep NIR observations. We combine the Subaru/MOIRCS and CFHT/WIRCam archival data to generate deep J and Ks band images, covering the full GOODS-N field (~169 sq. arcmin) to an AB magnitude limit of ~25 mag (3sigma). We applied z'-band dropout color selection criteria, using the NIR data generated here. We have identified two possible Lyman Break Galaxy (LBG) candidates at z\gtrsim6.5 with J\lesssim24.5. The first candidate is a likely LBG at z\sim6.5 based on a weak spectral feature tentatively identified as Lyalpha line in the deep Keck/DEIMOS spectrum, while the second candidate is a possible LBG at z\sim7 based on its photometric redshift. These z'-dropout objects, if confirmed, are among the brightest such candidates found so far. At z\gtrsim6.5, their star formation rate is estimated as 100-200 solar mass per year. If they continue to form stars at this rate, they assemble a stellar mass of ~5x10^10 solar mass after about 400 million years, becoming the progenitors of massive galaxies observed at z\sim5. We study the implication of the z'-band dropout candidates discovered here, in constraining the bright-end of the luminosity function and understanding the nature of high redshift galaxies.

preprint2012arXiv

SDSS 0956+5128: A Broad-line Quasar with Extreme Velocity Offsets

We report on the discovery of a Type 1 quasar, SDSS 0956+5128, with a surprising combination of extreme velocity offsets. SDSS 0956+5128 is a broad-lined quasar exhibiting emission lines at three substantially different redshifts: a systemic redshift of z ~ 0.714 based on narrow emission lines, a broad MgII emission line centered 1200 km/s bluer than the systemic velocity, at z ~ 0.707, and broad Hαand Hβemission lines centered at z ~ 0.690. The Balmer line peaks are 4100 km/s bluer than the systemic redshift. There are no previously known objects with such an extreme difference between broad MgII and broad Balmer emission. The two most promising explanations are either an extreme disk emitter or a high-velocity black hole recoil. However, neither explanation appears able to explain all of the observed features of SDSS 0956+5128, so the object may provide a challenge to our general understanding of quasar physics.

preprint2011arXiv

Galaxies in X-ray Groups I: Robust Membership Assignment and the Impact of Group Environments on Quenching

Understanding the mechanisms that lead dense environments to host galaxies with redder colors, more spheroidal morphologies, and lower star formation rates than field populations remains an important problem. As most candidate processes ultimately depend on host halo mass, accurate characterizations of the local environment, ideally tied to halo mass estimates and spanning a range in halo mass and redshift are needed. In this work, we present and test a rigorous, probabalistic method for assigning galaxies to groups based on precise photometric redshifts and X-ray selected groups drawn from the COSMOS field. The groups have masses in the range 10^13 < M_200c/M_sun < 10^14 and span redshifts 0<z<1. We characterize our selection algorithm via tests on spectroscopic subsamples, including new data obtained at the VLT, and by applying our method to detailed mock catalogs. We find that our group member galaxy sample has a purity of 84% and completeness of 92% within 0.5 R200c. We measure the impact of uncertainties in redshifts and group centering on the quality of the member selection with simulations based on current data as well as future imaging and spectroscopic surveys. As a first application of our new group member catalog which will be made publicly available, we show that member galaxies exhibit a higher quenched fraction compared to the field at fixed stellar mass out to z~1, indicating a significant relationship between star formation and environment at group scales. We also address the suggestion that dusty star forming galaxies in such groups may impact the high-l power spectrum of the cosmic microwave background and find that such a population cannot explain the low power seen in recent SZ measurements.

preprint2011arXiv

Measuring the Geometry of the Universe from Weak Gravitational Lensing behind Galaxy Groups in the HST COSMOS survey

Gravitational lensing can provide pure geometric tests of the structure of space-time, for instance by determining empirically the angular diameter distance-redshift relation. This geometric test has been demonstrated several times using massive clusters which produce a large lensing signal. In this case, matter at a single redshift dominates the lensing signal, so the analysis is straightforward. It is less clear how weaker signals from multiple sources at different redshifts can be stacked to demonstrate the geometric dependence. We introduce a simple measure of relative shear which for flat cosmologies separates the effect of lens and source positions into multiplicative terms, allowing signals from many different source-lens pairs to be combined. Applying this technique to a sample of groups and low-mass clusters in the COSMOS survey, we detect a clear variation of shear with distance behind the lens. This represents the first detection of the geometric effect using weak lensing by multiple, low-mass systems. The variation of distance with redshift is measured with sufficient precision to constrain the equation of state of the universe under the assumption of flatness, equivalent to a detection of a dark energy component Omega_X at greater than 99% confidence for an equation-of-state parameter -2.5 < w < -0.1. For the case w = -1, we find a value for the cosmological constant density parameter Omega_Lambda = 0.85+0.044-0.19 (68% C.L.), and detect cosmic acceleration (q_0 < 0) at the 98% C.L.. We consider the systematic uncertainties associated with this technique and discuss the prospects for applying it in forthcoming weak-lensing surveys.

preprint2011arXiv

New constraints on the evolution of the stellar-to-dark matter connection: a combined analysis of galaxy-galaxy lensing, clustering, and stellar mass functions from z=0.2 to z=1

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

preprint2011arXiv

SpecPro: An Interactive IDL Program for Viewing and Analyzing Astronomical Spectra

We present an interactive IDL program for viewing and analyzing astronomical spectra in the context of modern imaging surveys. SpecPro's interactive design lets the user simultaneously view spectroscopic, photometric, and imaging data, allowing for rapid object classification and redshift determination. The spectroscopic redshift can be determined with automated cross-correlation against a variety of spectral templates or by overlaying common emission and absorption features on the 1-D and 2-D spectra. Stamp images as well as the spectral energy distribution (SED) of a source can be displayed with the interface, with the positions of prominent photometric features indicated on the SED plot. Results can be saved to file from within the interface. In this paper we discuss key program features and provide an overview of the required data formats.

preprint2010arXiv

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

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

preprint2010arXiv

The 2008 extreme outburst of the young eruptive variable star Ex Lupi

In early 2008, the young low-mass star EX Lupi, the prototype of the EXor class of eruptive variables, optically brightened by over five magnitudes for a period of 7 months. The previous time a change of such amplitude had been observed in EX Lup was over 50 years ago. In this Letter we present new optical and near-IR high resolution spectroscopy of EX~Lup during the 2008 outburst. We investigate the physical characteristics of the outburst both soon after it began and some four months later, and consider the energetics and kinematics observed. Emission line strengths, widths, and profiles changed significantly between the two observations. Also, modeling of the 2.2935 um CO overtone bandhead emission suggests that an inner gap in the circumstellar gas disk around the star may be present and it is from the inner edge of the gas disk that the CO overtone emission probably arises. We derive a mass accretion luminosity and rate during the extreme outburst of ~2+-0.5~Lsun and ~2+-0.5x10^-7 Msun yr^-1, respectively, which suggests that this outburst was indeed one of the strongest witnessed in EX Lup, yet not as strong as those observed in FU Orionis stars.

preprint2010arXiv

The bulk of the black hole growth since z~1 occurs in a secular universe: No major merger-AGN connection

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

preprint2010arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S^4G)

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.

preprint2009arXiv

Environment of MAMBO galaxies in the COSMOS field

Submillimeter galaxies (SMG) represent a dust-obscured high-redshift population undergoing massive star formation activity. Their properties and space density have suggested that they may evolve into spheroidal galaxies residing in galaxy clusters. In this paper, we report the discovery of compact (~10"-20") galaxy overdensities centered at the position of three SMGs detected with the Max-Planck Millimeter Bolometer camera (MAMBO) in the COSMOS field. These associations are statistically significant. The photometric redshifts of galaxies in these structures are consistent with their associated SMGs; all of them are between z=1.4-2.5, implying projected physical sizes of ~170 kpc for the overdensities. Our results suggest that about 30% of the radio-identified bright SMGs in that redshift range form in galaxy density peaks in the crucial epoch when most stars formed.

preprint2009arXiv

Massive Galaxies in COSMOS: Evolution of Black hole versus bulge mass but not versus total stellar mass over the last 9 Gyrs?

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

preprint2009arXiv

The Nature of Optically Dull Active Galactic Nuclei in COSMOS

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

preprint2008arXiv

The COSMOS AGN Spectroscopic Survey I: XMM Counterparts

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

preprint2007arXiv

Dark matter maps reveal cosmic scaffolding

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

preprint2007arXiv

Evolution of the Bar Fraction in COSMOS: Quantifying the Assembly of the Hubble Sequence

We have analyzed the redshift-dependent fraction of galactic bars over 0.2<z<0.84 in 2,157 luminous face-on spiral galaxies from the COSMOS 2-square degree field. Our sample is an order of magnitude larger than that used in any previous investigation, and is based on substantially deeper imaging data than that available from earlier wide-area studies of high-redshift galaxy morphology. We find that the fraction of barred spirals declines rapidly with redshift. Whereas in the local Universe about 65% of luminous spiral galaxies contain bars (SB+SAB), at z ~0.84 this fraction drops to about 20%. Over this redshift range the fraction of strong (SB) bars drops from about 30% to under 10%. It is clear that when the Universe was half its present age, the census of galaxies on the Hubble sequence was fundamentally different from that of the present day. A major clue to understanding this phenomenon has also emerged from our analysis, which shows that the bar fraction in spiral galaxies is a strong function of stellar mass, integrated color and bulge prominence. The bar fraction in very massive, luminous spirals is about constant out to z ~ 0.84 whereas for the low mass, blue spirals it declines significantly with redshift beyond z=0.3. There is also a slight preference for bars in bulge dominated systems at high redshifts which may be an important clue towards the co-evolution of bars, bulges and black holes. Our results thus have important ramifications for the processes responsible for galactic downsizing, suggesting that massive galaxies matured early in a dynamical sense, and not just as a result of the regulation of their star formation rate.

preprint2002arXiv

A Redshift z = 6.56 Galaxy Behind the Cluster Abell 370

We report the discovery of a redshift z=6.56 galaxy lying behind the cluster Abell 370. The object HCM 6A was found in a narrowband imaging survey using a 118 Angstrom bandpass filter centered at 9152 Ang in the LRIS camera on the 10 m Keck II Telescope. Candidate Lyman alpha emitters were identified by the equivalent width of the emission and the absence of lower wavelength flux in ultradeep broadband images. HCM6A is the first galaxy to be confirmed at redshift z>6, and has W_lambda(observed)=190 Angstroms, flux = 2.7 x 10^-17 erg cm^-2 s^-1. Spectra obtained with LRIS confirm the emission line and the continuum break across the line, and show an asymmetric line profile with steep fall-off on the blue side. Deep Subaru near-infrared CISCO images in J, H and K' which extend the sampled continuum to longer wavelengths give a consistent estimate of the continuum flux density in these line-free regions of 2.6+/-0.7 x 10^-30 erg cm^-2 s^-1 Hz^-1. The line width and strength, asymmetric profile, and very deep spectral break are only consistent with the interpretation of the line as a redshifted Lyman alpha feature. From the detailed lensing model of this cluster, we estimate a lensing amplification of 4.5 for this galaxy, which is located slightly over an arcminute from the center of the cluster, for an unlensed flux of 6.5 x 10^-18 erg cm^-2 s^-1. The presence of such a galaxy suggests that the reionizing epoch is beyond z=6.6.