Source author record

Asantha Cooray

Asantha Cooray 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.IM gr-qc hep-ph astro-ph.EP astro-ph.HE hep-th physics.soc-ph

Catalog footprint

What is connected

128works

10topics

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.

preprint2022arXiv

Anisotropies of Cosmic Optical and Near-IR Background from China Space Station Telescope (CSST)

Anisotropies of the cosmic optical background (COB) and cosmic near-IR background (CNIRB) are capable of addressing some of the key questions in cosmology and astrophysics. In this work, we measure and analyze the angular power spectra of the simulated COB and CNIRB in the ultra-deep field of the China Space Station Telescope (CSST-UDF). The CSST-UDF covers about 9 square degrees, with magnitude limits ~28.3, 28.2, 27.6, 26.7 AB mag for point sources with 5-sigma detection in the r (0.620 um), i (0.760 um), z (0.915 um), and y (0.965 um) bands, respectively. According to the design parameters and scanning pattern of the CSST, we generate mock data, merge images and mask the bright sources in the four bands. We obtain four angular power spectra from l=200 to 2,000,000 (from arcsecond to degree), and fit them with a multi-component model including intrahalo light (IHL) using the Markov chain Monte Carlo (MCMC) method. We find that the signal-to-noise ratio (SNR) of the IHL is larger than 8 over the range of angular scales that are useful for astrophysical studies (l~10,000-400,000). Comparing to previous works, the constraints on the model parameters are improved by factors of 3~4 in this study, which indicates that the CSST-UDF survey can be a powerful probe on the cosmic optical and near-IR backgrounds.

preprint2022arXiv

Gas Phase Metallicities of Local Ultra-Luminous Infrared Galaxies Follow Normal Star-Forming Galaxies

Despite advances in observational data, theoretical models, and computational techniques to simulate key physical processes in the formation and evolution of galaxies, the stellar mass assembly of galaxies still remains an unsolved problem today. Optical spectroscopic measurements appear to show that the gas-phase metallicities of local ultra-luminous infrared galaxies (ULIRGs) are significantly lower than those of normal star-forming galaxies. This difference has resulted in the claim that ULIRGs are fueled by metal-poor gas accretion from the outskirts\cite{Mannucci10}. Here we report on a new set of gas-phase metallicity measurements making use of the far-infrared spectral lines of [O{\sc iii}]52 $μ$m, [O{\sc iii}]88 $μ$m, and [N{\sc iii}]57 $μ$m instead of the usual optical lines. Photoionization models have resulted in a metallicity diagnostic based on these three lines that break the electron density degeneracy and reduce the scatter of the correlation significantly. Using new data from SOFIA and archival data from Herschel Space Observatory, we find that local ULIRGs lie on the mass-metallicity relation of star-forming galaxies and have metallicities comparable to other galaxies with similar stellar masses and star formation rates. The lack of a departure suggests that ULIRGs follow the same mass assembly mechanism as luminous star-forming galaxies and $\sim 0.3$ dex under-abundance in metallicities derived from optical lines is a result of heavily obscured metal-rich gas which has a negligible effect when using the FIR line diagnostics.

preprint2022arXiv

Snowmass 2021 CMB-S4 White Paper

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.

preprint2021arXiv

A Long Stream of Metal-Poor Cool Gas around a Massive Starburst Galaxy at z = 2.67

We present the first detailed dissection of the circumgalactic medium (CGM) of massive starburst galaxies at z > 2. Our target is a submillimeter galaxy (SMG) at z = 2.674 that has a star formation rate of 1200 $M_\odot$/yr and a molecular gas reservoir of $1.3\times10^{11} M_\odot$. We characterize its CGM with two background QSOs at impact parameters of 93 kpc and 176 kpc. We detect strong HI and metal-line absorption near the redshift of the SMG towards both QSOs, each consisting of three main subsystems spanning over 1500 km/s. The absorbers show remarkable kinematic and metallicity coherence across a separation of 86 kpc. In particular, the cool gas in the CGM of the SMG exhibits high HI column densities ($\log N_{\rm HI}/{\rm cm}^{-2} = 20.2, 18.6$), low metallicities ([M/H] $\approx$ -2.0), and similar radial velocities ($\approx$ -300 km/s). While the HI column densities match previous results on the CGM around QSOs at z > 2, the metallicities are lower by more than an order of magnitude, making it an outlier in the line width$-$metallicity relation of damped Ly$α$ absorbers. The large physical extent, the velocity coherence, the high surface density, and the low metallicity are all consistent with the cool, inflowing, and near-pristine gas streams predicted to penetrate hot massive halos at z > 1.5. We estimate a total gas accretion rate of ~100 $M_\odot$/yr from three such streams, which falls short of the star formation rate but is consistent with simulations. At this rate, it takes about a gigayear to acquire the molecular gas reservoir of the central starburst.

preprint2021arXiv

Measurements of the Zodiacal Light Absolute Intensity through Fraunhofer Absorption Line Spectroscopy with CIBER

Scattered sunlight from the interplanetary dust (IPD) cloud in our Solar system presents a serious foreground challenge for spectro-photometric measurements of the Extragalactic Background Light (EBL). In this work, we report on measurements of the absolute intensity of the Zodiacal Light (ZL) using the novel technique of Fraunhofer line spectroscopy on the deepest 8542 Angstrom line of the near-infrared CaII absorption triplet. The measurements are performed with the Narrow Band Spectrometer (NBS) aboard the Cosmic Infrared Background Experiment (CIBER) sounding rocket instrument. We use the NBS data to test the accuracy of two ZL models widely cited in the literature; the Kelsall and Wright models, which have been used in foreground removal analyses that produce high and low EBL results respectively. We find a mean reduced chi squared of 3.5 for the Kelsall model and a chi squared of 2.0 for the Wright model. The best description of our data is provided by a simple modification to the Kelsall model which includes a free ZL offset parameter. This adjusted model describes the data with a reduced chi squared of 1.5 and yields an inferred offset amplitude of 46 +- 19 nW m^-2 sr^-1 extrapolated to 12500 Angstroms. These measurements elude to the potential existence of a dust cloud component in the inner Solar system whose intensity does not strongly modulate with the Earth's motion around the Sun.

preprint2021arXiv

Modelling high-resolution ALMA observations of strongly lensed dustystar forming galaxies detected by Herschel

We present modelling of ~0.1arcsec resolution Atacama Large Millimetre/sub-millimeter Array imaging of seven strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Four of these systems are galaxy-galaxy scale strong lenses, with the remaining three being group-scale lenses. Through careful modelling of visibilities, we infer the mass profiles of the lensing galaxies and by determining the magnification factors, we investigate the intrinsic properties and morphologies of the lensed sub-millimetre sources. We find that these sub-millimetre sources all have ratios of star formation rate to dust mass that is consistent with or in excess of the mean ratio for high-redshift sub-millimetre galaxies and low redshift ultra-luminous infrared galaxies. The contribution to the infrared luminosity from possible AGN is not quantified and so could be biasing our star formation rates to higher values. The majority of our lens models have mass density slopes close to isothermal, but some systems show significant differences.

preprint2021arXiv

Polarization spectrum of near infrared zodiacal light observed with CIBER

We report the first measurement of the zodiacal light (ZL) polarization spectrum in the near-infrared between 0.8 and 1.8 $μ$m. Using the low-resolution spectrometer (LRS) on board the Cosmic Infrared Background Experiment (CIBER), calibrated for absolute spectrophotometry and spectropolarimetry, we acquire long-slit polarization spectral images of the total diffuse sky brightness towards five fields. To extract the ZL spectrum, we subtract contribution of other diffuse radiation, such as the diffuse galactic light (DGL), the integrated star light (ISL), and the extragalactic background light (EBL). The measured ZL polarization spectrum shows little wavelength dependence in the near-infrared and the degree of polarization clearly varies as a function of the ecliptic coordinates and solar elongation. Among the observed fields, the North Ecliptic Pole shows the maximum degree of polarization of $\sim$ 20$\%$, which is consistent with an earlier observation from the Diffuse Infrared Background Experiment (DIRBE) aboard on the Cosmic Background Explorer (COBE). The measured degree of polarization and its solar elongation dependence are reproduced by the empirical scattering model in the visible band and also by the Mie scattering model for large absorptive particles, while the Rayleigh scattering model is ruled out. All of our results suggest that the interplanetary dust is dominated by large particles.

preprint2020arXiv

Cosmological constraints from line intensity mapping with interlopers

Understanding the formation and evolution of the Universe is crucial for cosmological studies, and the line intensity mapping provides a powerful tool for this kind of study. We propose to make use of multipole moments of redshift-space line intensity power spectrum to constrain the cosmological and astrophysical parameters, such as the equation of state of dark energy, massive neutrinos, primordial non-Gaussianity, and star formation rate density. As an example, we generate mock data of multipole power spectra for H-alpha 6563AA, [OIII] 5007AA and [OII] 3727AA measured by SPHEREx experiment at z=1 considering contaminations from interloper lines, and use Markov Chain Monte Carlo (MCMC) method to constrain the parameters in the model. We find a good fitting result of the parameters compared to their fiducial values, which means that the multipole power spectrum can effectively distinguish signal and interloper lines, and break the degeneracies between parameters, such as line mean intensity and bias. We also explore the cross power spectrum with CSST (Chinese Space Station Telescope) spectroscopic galaxy survey in the constraints. Since more accurate fitting results can be obtained by including measurements of the emission lines at higher redshifts out to z=3 at least and cross-correlations between emission lines can be involved, the line intensity mapping is expected to offer excellent results in future cosmological and astrophysical studies.

preprint2020arXiv

Emergence of an Ultra-Red Ultra-Massive Galaxy Cluster Core at $z=4$

Recent simulations and observations of massive galaxy cluster evolution predict that the majority of stellar mass build up happens within cluster members by $z=2$, before cluster virialization. Protoclusters rich with dusty, star-forming galaxies (DSFGs) at $z>3$ are the favored candidate progenitors for these massive galaxy clusters at $z\sim0$. We present here the first study analyzing stellar emission along with cold dust and gas continuum emission in a spectroscopically confirmed $z=4.002$ protocluster core rich with DSFGs, the Distant Red Core (DRC). We combine new HST and Spitzer data with existing Gemini, Herschel, and ALMA observations to derive individual galaxy-level properties, and compare them to coeval field and other protocluster galaxies. All of the protocluster members are massive ($>10^{10}$ M$_\odot$), but not significantly more so than their coeval field counterparts. Within uncertainty, all are nearly indistinguishable from galaxies on the star-forming vs. stellar mass main-sequence relationship, and on the star formation efficiency plane. Assuming no future major influx of fresh gas, we estimate that these gaseous DSFGs will deplete their gas reservoirs in $\sim300$ Myr, becoming the massive quiescent ellipticals dominating cluster cores by $z\sim3$. Using various methodologies, we derive a total $z=4$ halo mass of $\sim10^{14}$ M$_\odot$, and estimate that the DRC will evolve to become an ultra-massive cluster core of mass $\gtrsim10^{15}$ M$_\odot$ by $z=0$.

preprint2020arXiv

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.

preprint2020arXiv

The Case for Probe-class NASA Astrophysics Missions

Astrophysics spans an enormous range of questions on scales from individual planets to the entire cosmos. To address the richness of 21st century astrophysics requires a corresponding richness of telescopes spanning all bands and all messengers. Much scientific benefit comes from having the multi-wavelength capability available at the same time. Most of these bands,or measurement sensitivities, require space-based missions. Historically, NASA has addressed this need for breadth with a small number of flagship-class missions and a larger number of Explorer missions. While the Explorer program continues to flourish, there is a large gap between Explorers and strategic missions. A fortunate combination of new astrophysics technologies with new, high capacity, low dollar-per-kg to orbit launchers, and new satellite buses allow for cheaper missions with capabilities approaching strategic mission levels. NASA has recognized these developments by calling for Probe-class mission ideas for mission studies, spanning most of the electromagnetic spectrum from GeV gamma-rays to the far infrared, and the new messengers of neutrinos and ultra-high energy cosmic rays. The key insight from the Probes exercise is that order-of-magnitude advances in science performance metrics are possible across the board for initial total cost estimates in the range 500M-1B dollars.

preprint2020arXiv

Three Dusty Star Forming Galaxies at $z\sim1.5$: Mergers and Disks on the Main Sequence

The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.

preprint2019arXiv

Photometric Redshift Estimation with Galaxy Morphology using Self-Organizing Maps

We use multi-band optical and near-infrared photometric observations of galaxies in the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) to predict photometric redshifts using artificial neural networks. The multi-band observations span over 0.39 microns to 8.0 microns for a sample of about 1000 galaxies in the GOODS-S field for which robust size measurements are available from Hubble Space Telescope Wide Field Camera 3 observations. We use Self Organizing Maps (SOMs) to map the multi dimensional photometric and galaxy size observations while taking advantage of existing spectroscopic redshifts at 0 < z < 2 for independent training and testing sets. We show that use of photometric and morphological data led to redshift estimates comparable to redshift measurements from SED modeling and from self-organizing maps without morphological measurements.

preprint2017arXiv

Tomographic Intensity Mapping versus Galaxy Surveys: Observing the Universe in H-alpha emission with new generation instruments

The H-alpha line emission is an important probe for a number of fundamental quantities in galaxies, including their number density, star formation rate (SFR) and overall gas content. A new generation of low-resolution intensity mapping probes, e.g. SPHEREx and CDIM, will observe galaxies in H-alpha emission over a large fraction of the sky from the local Universe till a redshift of z ~ 6 to 10, respectively. This will also be the target line for observations by the high-resolution Euclid and WFIRST instruments in the z ~ 0.7 - 2 redshift range. In this paper, we estimate the intensity and power spectra of the H-alpha line in the z ~ 0 - 5 redshift range using observed line luminosity functions (LFs), when possible, and simulations, otherwise. We estimate the significance of our predictions by accounting for the modelling uncertainties (e.g. SFR, extinction, etc.) and observational contamination. We find that Intensity Mapping (IM) surveys can make a statistical detection of the full H-alpha emission between z ~ 0.8 - 5. Moreover, we find that the high-frequency resolution and the sensitivity of the planned CDIM surveys allow for the separation of H-alpha emission from several interloping lines. We explore ways to use the combination of these line intensities to probe galaxy properties. As expected, our study indicates that galaxy surveys will only detect bright galaxies that contribute up to a few percent of the overall H-alpha intensity. However, these surveys will provide important constraints on the high end of the H-alpha LF and put strong constraints on the AGN LF.

preprint2016arXiv

Axion decay and anisotropy of near-IR extragalactic background light

The extragalactic background light (EBL) is comprised of the cumulative radiation from all galaxies and active galactic nuclei over the cosmic history. In addition to point sources, EBL also contains information from diffuse sources of radiation. The angular power spectra of the near-infrared intensities could contain additional signals and a complete understanding of the nature of the IR background is still lacking in the literature. Here we explore the constraints that can be placed on particle decays, especially candidate dark matter models involving axions that trace dark matter halos of galaxies. Axions with a mass around a few eV will decay via two photons with wavelengths in the near-IR band, and will leave a signature in the IR background intensity power spectrum. Using recent power spectra measurements from the Hubble Space Telescope (HST) and Cosmic Infrared Background Experiment (CIBER), we find that the 0.6 to 1.6 micron power spectra can be explained by axions with masses around 4 eV. The total axion abundance Omega_a~0.05, and it is comparable to the baryon density of the Universe. The suggested mean axion mass and abundance are not ruled out by existing cosmological observations. Interestingly, the axion model with a mass distribution is preferred by the data, which cannot be explained by the standard quantum chromodynamics (QCD) theory and needs further discussion.

preprint2016arXiv

Cosmic Dawn Intensity Mapper

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

preprint2016arXiv

Extragalactic Background Light: Measurements and Applications

This review covers the measurements related to the extragalactic background light (EBL) intensity from gamma-rays to radio in the electromagnetic spectrum over 20 decades in the wavelength. The cosmic microwave background (CMB) remains the best measured spectrum with an accuracy better than 1%. The measurements related to the cosmic optical background (COB), centered at 1 microns, are impacted by the large zodiacal light associated with interplanetary dust in the inner Solar system. The best measurements of COB come from an indirect technique involving Gamma-ray spectra of bright blazars with an absorption feature resulting from pair-production off of COB photons. The cosmic infrared background (CIB) peaking at around 100 microns established an energetically important background with an intensity comparable to the optical background. This discovery paved the path for large aperture far-infrared and sub-millimeter observations resulting in the discovery of dusty, starbursting galaxies. Their role in galaxy formation and evolution remains an active area of research in modern-day astrophysics. The extreme UV background remains mostly unexplored and will be a challenge to measure due to the high Galactic background and absorption of extragalactic photons by the intergalactic medium at these EUV/soft X-ray energies. We also summarize our understanding of the spatial anisotropies and angular power spectra of intensity fluctuations. We motivate a precise direct measurement of the COB between 0.1 to 5 microns using a small aperture telescope observing either from the outer Solar system, at distances of 5 AU or more, or out of the ecliptic plane. Other future applications include improving our understanding of the background at TeV energies and spectral distortions of CMB and CIB.

preprint2016arXiv

Multi-wavelength lens reconstruction of a $\textit{Planck}$ $\&$ $\textit{Herschel}$-detected star-bursting galaxy

We present a source-plane reconstruction of a ${\it Herschel}$ and ${\it Planck}$-detected gravitationally-lensed dusty star-forming galaxy (DSFG) at $z=1.68$ using {\it Hubble}, Sub-millimeter Array (SMA), and Keck observations. The background sub-millimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at $z=0.997$ and appears as an arc of length $\sim 15^{\prime \prime}$ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of a multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model SEDs leading to a magnification-corrected star formation rate of $390 \pm 60$ M$_{\odot}$ yr$^{-1}$ and a stellar mass of $1.1 \pm 0.4\times 10^{11}$ M$_{\odot}$. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of $0.43 \pm 0.13$, $350 \pm 200$ M$_{\odot}$ pc$^{-2}$, and $\sim 12 \pm 7~$M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, respectively. The ratio of star formation rate surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGS.

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

Spectral Line De-confusion in an Intensity Mapping Survey

Spectral line intensity mapping has been proposed as a promising tool to efficiently probe the cosmic reionization and the large-scale structure. Without detecting individual sources, line intensity mapping makes use of all available photons and measures the integrated light in the source confusion limit, to efficiently map the three-dimensional matter distribution on large scales as traced by a given emission line. One particular challenge is the separation of desired signals from astrophysical continuum foregrounds and line interlopers. Here we present a technique to extract large-scale structure information traced by emission lines from different redshifts, embedded in a three-dimensional intensity mapping data cube. The line redshifts are distinguished by the anisotropic shape of the power spectra when projected onto a common coordinate frame. We consider the case where high-redshift [CII] lines are confused with multiple low-redshift CO rotational lines. We present a semi-analytic model for [CII] and CO line estimates based on the cosmic infrared background measurements, and show that with a modest instrumental noise level and survey geometry, the large-scale [CII] and CO power spectrum amplitudes can be successfully extracted from a confusion-limited data set, without external information. We discuss the implications and limits of this technique for possible line intensity mapping experiments.

preprint2016arXiv

Star Formation in Luminous Quasars at 2<z<3

We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.

preprint2016arXiv

Toward a tomographic analysis of the cross-correlation between Planck CMB lensing and H-ATLAS galaxies

We present an improved and extended analysis of the cross-correlation between the map of the Cosmic Microwave Background (CMB) lensing potential derived from the \emph{Planck} mission data and the high-redshift galaxies detected by the \emph{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) in the photometric redshift range $z_{\rm ph} \ge 1.5$. We compare the results based on the 2013 and 2015 \textit{Planck} datasets, and investigate the impact of different selections of the H-ATLAS galaxy samples. Significant improvements over our previous analysis have been achieved thanks to the higher signal-to-noise ratio of the new CMB lensing map recently released by the \textit{Planck} collaboration. The effective galaxy bias parameter, $b$, for the full galaxy sample, derived from a joint analysis of the cross-power spectrum and of the galaxy auto-power spectrum is found to be $b = 3.54^{+0.15}_{-0.14}$. Furthermore, a first tomographic analysis of the cross-correlation signal is implemented, by splitting the galaxy sample into two redshift intervals: $1.5 \le z_{\rm ph} < 2.1$ and $z_{\rm ph}\ge 2.1$. A statistically significant signal was found for both bins, indicating a substantial increase with redshift of the bias parameter: $b=2.89\pm0.23$ for the lower and $b=4.75^{+0.24}_{-0.25}$ for the higher redshift bin. Consistently with our previous analysis we find that the amplitude of the cross correlation signal is a factor of $1.45^{+0.14}_{-0.13}$ higher than expected from the standard $Λ$CDM model for the assumed redshift distribution. The robustness of our results against possible systematic effects has been extensively discussed although the tension is mitigated by passing from 4 to 3$σ$.

preprint2015arXiv

ALMA maps the Star-Forming Regions in a Dense Gas Disk at z~3

We exploit long-baseline ALMA sub-mm observations of the lensed star-forming galaxy SDP 81 at z=3.042 to investigate the properties of inter-stellar medium on scales of 50-100pc. The kinematics of the CO gas within this system are well described by a rotationally-supported disk with an inclination-corrected rotation speed, v=320+/-20km/s and a dynamical mass of M=(3.5+/-1.0)x10^10Mo within a radius of 1.5 kpc. The disk is gas rich and unstable, with a Toomre parameter, Q=0.30+/-0.10 and so should collapse in to star-forming regions with Jeans length L_J~130pc. We identify five star-forming regions within the ISM on these scales and show that their scaling relations between luminosity, line-widths and sizes are significantly offset from those typical of molecular clouds in local Galaxies (Larson's relations). These offsets are likely to be caused by the high external hydrostatic pressure for the interstellar medium (ISM), P/kB=(40+/-20)x10^7K/cm3, which is ~10,000x higher than the typical ISM pressure in the Milky Way. The physical conditions of the star-forming ISM and giant molecular clouds appears to be similar to the those found in the densest environments in the local Universe, such as those in the Galactic center.

preprint2015arXiv

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

preprint2015arXiv

Astronomy Job Crisis

Astronomers in CANDELS outline changes for the academic system to promote a smooth transition for junior scientists from academia to industry.

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

Dust energy balance study of two edge-on spiral galaxies in the Herschel-ATLAS survey

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

preprint2015arXiv

Extinction and nebular line properties of a Herschel-selected lensed dusty starburst AT z=1.027

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $μ$m, resulted in detections of H$α$+[NII], H$β$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$α$/H$β$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $Hα$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation.

preprint2015arXiv

Far-infrared observations of an unbiased sample of gamma-ray burst host galaxies

Gamma-ray bursts (GRBs) are the most energetic phenomena in the Universe; believed to result from the collapse and subsequent explosion of massive stars. Even though it has profound consequences for our understanding of their nature and selection biases, little is known about the dust properties of the galaxies hosting GRBs. We present analysis of the far-infrared properties of an unbiased sample of 20 \textit{BeppoSAX} and \textit{Swift} GRB host galaxies (at an average redshift of $z\,=\,3.1$) located in the {\it Herschel} Astrophysical Terahertz Large Area Survey, the {\it Herschel} Virgo Cluster Survey, the {\it Herschel} Fornax Cluster Survey, the {\it Herschel} Stripe 82 Survey and the {\it Herschel} Multi-tiered Extragalactic Survey, totalling $880$ deg$^2$, or $\sim 3$\% of the sky in total. Our sample selection is serendipitous, based only on whether the X-ray position of a GRB lies within a large-scale {\it Herschel} survey -- therefore our sample can be considered completely unbiased. Using deep data at wavelengths of 100\,--\,500$\,μ$m, we tentatively detected 1 out of 20 GRB hosts located in these fields. We constrain their dust masses and star formation rates (SFRs), and discuss these in the context of recent measurements of submillimetre galaxies and ultraluminous infrared galaxies. The average far-infrared flux of our sample gives an upper limit on SFR of $<114\,{\rm M}\odot\,\mbox{yr}^{-1}$. The detection rate of GRB hosts is consistent with that predicted assuming that GRBs trace the cosmic SFR density in an unbiased way, i.e. that the fraction of GRB hosts with $\mbox{SFR}>500\,{\rm M}\odot\,\mbox{yr}^{-1}$ is consistent with the contribution of such luminous galaxies to the cosmic star formation density.

preprint2015arXiv

H-ATLAS/GAMA: Quantifying the Morphological Evolution of the Galaxy Population Using Cosmic Calorimetry

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars.

preprint2015arXiv

Planck Trispectrum Constraints on Primordial Non-Gaussianity at Cubic Order

Non-Gaussianity of the primordial density perturbations provides an important measure to constrain models of inflation. At cubic order the non-Gaussianity is captured by two parameters $τ_{\rm NL}$ and $g_{\rm NL}$ that determine the amplitude of the density perturbation trispectrum. Here we report measurements of the kurtosis power spectra of the cosmic microwave background (CMB) temperature as mapped by Planck by making use of correlations between square temperature-square temperature and cubic temperature-temperature anisotropies. In combination with noise simulations, we find the best joint estimates to be $τ_{\rm{NL}}=0.3 \pm 0.9 \times 10^4$ and $g_{\rm{NL}}=-1.2 \pm 2.8 \times 10^5$. If $τ_{\rm NL}=0$, we find $g_{\rm NL}= -1.3\pm 1.8 \times 10^5$.

preprint2015arXiv

Prospects for detecting CII emission during the Epoch of Reionization

We produce simulations of emission of the atomic CII line in large sky fields in order to determine the current prospects for mapping this line during the high redshift Epoch of Reionization. We estimate the CII line intensity, redshift evolution and spatial fluctuations using observational relations between CII emission and the SFR in a galaxy for the frequency range of 200 GHz to 300 GHz. We obtained a frequency averaged intensity of CII emission of ${\rm I_{\rm CII}=(4 \pm 2)\times10^{2}\, Jy\, \rm sr^{-1}}$ in the redshift range $z\, \sim\, 5.3\, -\, 8.5$. Observations of CII emission in this frequency range will suffer contamination from emission lines at lower redshifts, in particular from the CO rotation lines. For the relevant frequency range we estimated the CO contamination (originated in emission from galaxies at $z\, <\, 2.5$), using simulations, to be ${\rm I_{\rm CO} \approx 10^{3}\, Jy \, sr^{-1}}$ and independently confirmed the result based in observational relations. We generated maps as a function of angle and frequency using detailed simulations of the CII and CO emission across several redshifts in order to properly take into account the observational pipeline and light cone effects. In order to reduce the foreground contamination we found that we should mask galaxies below redshifts $\sim 2.5$ with a CO flux in one of the CO(J:2-1) to CO(J:6-5) lines higher than ${\rm 5\times 10^{-22}\, W\ m^{-2}}$ or a AB magnitude lower than ${\rm m_{\rm K}\, =\, 22}$. We estimate that the additional continuum contamination is of the order of ${\rm 10^{5}\, Jy\, sr^{-1}}$. It is also considered the possibility of cross correlating foreground lines with galaxies in order to probe the intensity of the foregrounds.

preprint2015arXiv

Spitzer Imaging of Strongly-Lensed Herschel-Selected Dusty Star Forming Galaxies

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

preprint2015arXiv

Ultraviolet Luminosity Density of the Universe During the Epoch of Reionization

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $μ$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log ρ_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1σ)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

preprint2014arXiv

Cross-Correlation of Near and Far-Infrared Background Anisotropies as Traced by Spitzer and Herschel

We present the cross-correlation between the far-infrared background fluctuations as measured with the Herschel Space Observatory at 250, 350, and 500 μm and the near-infrared background fluctuations with Spitzer Space Telescope at 3.6 μm. The cross-correlation between far and near-IR background anisotropies are detected such that the correlation coefficient at a few to ten arcminute angular scales decreases from 0.3 to 0.1 when the far-IR wavelength increases from 250 μm to 500 μm. We model the cross-correlation using a halo model with three components: (a) far-IR bright or dusty star-forming galaxies below the masking depth in Herschel maps, (b) near-IR faint galaxies below the masking depth at 3.6 μm, and (c) intra-halo light, or diffuse stars in dark matter halos, that likely dominates fluctuations at 3.6 μm. The model is able to reasonably reproduce the auto correlations at each of the far-IR wavelengths and at 3.6 μm and their corresponding cross-correlations. While the far and near-IR auto-correlations are dominated by faint dusty, star-forming galaxies and intra-halo light, respectively, we find that roughly half of the cross-correlation between near and far-IR backgrounds is due to the same galaxies that remain unmasked at 3.6 μm. The remaining signal in the cross-correlation is due to intra-halo light present in the same dark matter halos as those hosting the same faint and unmasked galaxies. In this model, the decrease in the cross-correlation signal from 250 μm to 500 μm comes from the fact that the galaxies that are primarily contributing to 500 μm fluctuations peak at a higher redshift than those at 250 μm.

preprint2014arXiv

Dusty Star-Forming Galaxies at High Redshift

Far-infrared and submillimeter wavelength surveys have now established the important role of dusty, star-forming galaxies (DSFGs) in the assembly of stellar mass and the evolution of massive galaxies in the Universe. The brightest of these galaxies have infrared luminosities in excess of 10$^{13}$ L$_{\odot}$ with implied star-formation rates of thousands of solar masses per year. They represent the most intense starbursts in the Universe, yet many are completely optically obscured. Their easy detection at submm wavelengths is due to dust heated by ultraviolet radiation of newly forming stars. When summed up, all of the dusty, star-forming galaxies in the Universe produce an infrared radiation field that has an equal energy density as the direct starlight emission from all galaxies visible at ultraviolet and optical wavelengths. The bulk of this infrared extragalactic background light emanates from galaxies as diverse as gas-rich disks to mergers of intense starbursting galaxies. Major advances in far-infrared instrumentation in recent years, both space-based and ground-based, has led to the detection of nearly a million DSFGs, yet our understanding of the underlying astrophysics that govern the start and end of the dusty starburst phase is still in nascent stage. This review is aimed at summarizing the current status of DSFG studies, focusing especially on the detailed characterization of the best-understood subset (submillimeter galaxies, who were summarized in the last review of this field over a decade ago, Blain et al., 2002), but also the selection and characterization of more recently discovered DSFG populations. We review DSFG population statistics, their physical properties including dust, gas and stellar contents, their environments, and current theoretical models related to the formation and evolution of these galaxies.

preprint2014arXiv

Foreground contamination in Ly-alpha intensity mapping during the epoch of reionization

The intensity mapping of Ly-alpha emission during the epoch of reionization (EoR) will be contaminated by foreground emission lines from lower redshifts. We calculate the mean intensity and power spectrum of Ly-alpha emission at z~7, and estimate the uncertainties according to the relevant astrophysical processes. We find that the low-redshift emission lines from 6563 A H-alpha, 5007 A [OIII] and 3727 A [OII] will be the strong contaminants on the observed Ly-alpha power spectrum. We make use of both the star formation rate (SFR) and luminosity functions (LF) to estimate the mean intensity and power spectra of the three foreground lines at z~0.5 for H-alpha, z~0.9 for [OIII] and z~1.6 for [OII], as they will contaminate the Ly-alpha emission at z~7. The [OII] line is found to be the strongest. We analyze the masking of the bright survey pixels with a foreground line above some line intensity threshold as a way to reduce the contamination in the intensity mapping survey. We find that the foreground contamination can be neglected if we remove the pixels with fluxes above 1.4x10^-20 W/m^2.

preprint2014arXiv

HerMES: The Rest-Frame UV Emission and A Lensing Model for the z=6.34 Luminous Dusty Starburst Galaxy HFLS3

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3.

preprint2014arXiv

Herschel-ATLAS and ALMA: HATLAS J142935.3-002836, a lensed major merger at redshift 1.027

[Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]).

preprint2014arXiv

Herschel-ATLAS/GAMA:How does the far-IR luminosity function depend on galaxy group properties?

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

preprint2014arXiv

High-resolution, 3D radiative transfer modeling : I. The grand-design spiral galaxy M51

Context: Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered driven merely by young stars and, consequently, often applied to trace the star formation rate in galaxies. Recent studies have argued that the old stellar population might anticipate a non-negligible fraction of the radiative dust heating. Aims: In this work, we aim to analyze the contribution of young (< 100 Myr) and old (~ 10 Gyr) stellar populations to radiative dust heating processes in the nearby grand-design spiral galaxy M51 using radiative transfer modeling. High-resolution 3D radiative transfer (RT) models are required to describe the complex morphologies of asymmetric spiral arms and clumpy star-forming regions and model the propagation of light through a dusty medium. Methods: In this paper, we present a new technique developed to model the radiative transfer effects in nearby face-on galaxies. We construct a high-resolution 3D radiative transfer model with the Monte-Carlo code SKIRT accounting for the absorption, scattering and non-local thermal equilibrium (NLTE) emission of dust in M51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 μm, GALEX FUV, Hα and MIPS 24 μm wavebands, assuming an exponential vertical distribution with an appropriate scale height. The dust geometry is constrained through the far-ultraviolet (FUV) attenuation, which is derived from the observed total-infrared-to-far-ultraviolet luminosity ratio. The stellar luminosity, star formation rate and dust mass have been scaled to reproduce the observed stellar spectral energy distribution (SED), FUV attenuation and infrared SED. (abridged)

preprint2014arXiv

On the Origin of Near-Infrared Extragalactic Background Light Anisotropy

Extragalactic background light (EBL) anisotropy traces variations in the total production of photons over cosmic history, and may contain faint, extended components missed in galaxy point source surveys. Infrared EBL fluctuations have been attributed to primordial galaxies and black holes at the epoch of reionization (EOR), or alternately, intra-halo light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report new EBL anisotropy measurements from a specialized sounding rocket experiment at 1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely explained by IHL emission. The measured fluctuations are associated with an EBL intensity that is comparable to the background from known galaxies measured through number counts, and therefore a substantial contribution to the energy contained in photons in the cosmos.

preprint2014arXiv

Optical-faint, Far-infrared-bright Herschel Sources in the CANDELS Fields: Ultra-Luminous Infrared Galaxies at z>1 and the Effect of Source Blending

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

preprint2014arXiv

The applicability of FIR fine-structure lines as Star Formation Rate tracers over wide ranges of metallicities and galaxy types

We analyze the applicability of far-infrared fine-structure lines [CII] 158 micron, [OI] 63 micron and [OIII] 88 micron to reliably trace the star formation rate (SFR) in a sample of low-metallicity dwarf galaxies from the Herschel Dwarf Galaxy Survey and compare with a broad sample of galaxies of various types and metallicities in the literature. We study the trends and scatter in the relation between the SFR (as traced by GALEX FUV and MIPS 24 micron) and far-infrared line emission, on spatially resolved and global galaxy scales, in dwarf galaxies. We assemble far-infrared line measurements from the literature and infer whether the far-infrared lines can probe the SFR (as traced by the total-infrared luminosity) in a variety of galaxy populations. In metal-poor dwarfs, the [OI] and [OIII] lines show the strongest correlation with the SFR with an uncertainty on the SFR estimates better than a factor of 2, while the link between [CII] emission and the SFR is more dispersed (uncertainty factor of 2.6). The increased scatter in the SFR-L([CII]) relation towards low metal abundances, warm dust temperatures, large filling factors of diffuse, highly ionized gas suggests that other cooling lines start to dominate depending on the density and ionization state of the gas. For the literature sample, we evaluate the correlations for a number of different galaxy populations. The [CII] and [OI] lines are considered to be reliable SFR tracers in starburst galaxies, recovering the star formation activity within an uncertainty of factor 2. [Abridged]

preprint2014arXiv

The co-evolution of black hole growth and star formation from a cross-correlation analysis between quasars and the cosmic infrared background

We present the first cross-correlation measurement between Sloan Digital Sky Survey (SDSS) Type 1 quasars and the cosmic infrared background (CIB) measured by Herschel. The distribution of the quasars at 0.15<z<3.5 covers the redshift range where we expect most of the CIB to originate. We detect the sub-mm emission of the quasars, which dominates on small scales, as well as correlated emission from dusty star-forming galaxies (DSFGs) dominant on larger scales. The mean sub-mm flux densities of the DR7 quasars (median redshift <z>=1.4) is $11.1^{+1.6}_{-1.4}$, $7.1^{+1.6}_{-1.3}$ and $3.6^{+1.4}_{-1.0}$ mJy at 250, 350 and 500 microns, respectively, while the mean sub-mm flux densities of the DR9 quasars (<z>=2.5) is $5.7^{+0.7}_{-0.6}$, $5.0^{+0.8}_{-0.7}$ and $1.8^{+0.5}_{-0.4}$ mJy. We find that the correlated sub-mm emission includes both the emission from satellite DSFGs in the same halo as the central quasar and the emission from DSFGs in separate halos (correlated with the quasar-hosting halo). The amplitude of the one-halo term is ~10 times smaller than the sub-mm emission of the quasars, implying the the satellites have a lower star-formation rate than the quasars. The satellite fraction for the DR7 quasars is $0.008^{+0.008}_{-0.005}$ and the host halo mass scale for the central and satellite quasars is $10^{12.36\pm0.87}$ M$_{\odot}$ and $10^{13.60\pm0.38}$ M$_{\odot}$, respectively. The satellite fraction of the DR9 quasars is $0.065^{+0.021}_{-0.031}$ and the host halo mass scale for the central and satellite quasars is $10^{12.29\pm0.62}$ M$_{\odot}$ and $10^{12.82\pm0.39}$ M$_{\odot}$, respectively. Thus, the typical halo environment of the SDSS Type 1 quasars is found to be similar to that of DSFGs, which supports the generally accepted view that dusty starburst and quasar activity are evolutionarily linked.

preprint2014arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA). VI. The distribution and properties of molecular cloud associations in M31

In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.

preprint2013arXiv

Constraints on Spatial Variations in the Fine-Structure constant from Planck

We use the Cosmic Microwave Background (CMB) anisotropy data from Planck to constrain the spatial fluctuations of the fine-structure constant α. Through Thompson scattering of CMB photons, spatial anisotropies of αlead to higher-order correlations in the CMB anisotropies. We use a quadratic estimator based on the four-point correlation function of the CMB temperature anisotropy to extract the angular power spectrum of the spatial variation of the fine-structure constant projected along the line of sight at the last scattering surface. At tens of degree angular scales and above, we constrain the rms fluctuations of the fine structure constant to be δα/α_0= (1.34 +/- 5.82) x 10^-2 at the 95% confidence level with respect to the standard value α_0. We find no evidence for a spatially varying αat a redshift of 10^3.

preprint2013arXiv

Disentangling a group of lensed submm galaxies at z~2.9

MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.

preprint2013arXiv

Far-infrared Fine-Structure Line Diagnostics of Ultraluminous Infrared Galaxies

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

preprint2013arXiv

Herschel SPIRE-FTS Observations of Excited CO and [CI] in the Antennae (NGC 4038/39): Warm and Cold Molecular Gas

We present Herschel SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205μm$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63μm$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $α_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$ in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63μm$ to be $\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength $G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.

preprint2013arXiv

Herschel-ATLAS/GAMA: What determines the far-infrared properties of radio-galaxies?

We perform a stacking analysis of H-ATLAS data in order to obtain isothermal dust temperatures and rest-frame luminosities at 250um (L250), for 1599 radio sources over the H-ATLAS P1 GAMA area. The radio sample is generated using a combination of NVSS data and K-band UKIDSS-LAS data, over 0.01<z<0.8. The FIR properties of the sample are investigated as a function of L_1.4 GHz, z, projected radio-source size and radio spectral index. In order to search for stellar-mass dependent relations, we split the parent sample into those sources which are below and above 1.5 L_{K}^{*}. Correcting for stellar mass and redshift, we find no relation between the L250 and L_1.4GHz of radio AGN. This implies that a galaxy's nominal radio luminosity has little or no bearing on the SFR and/or dust mass content of the host system. The L250 of both the radio detected and non radio-detected galaxies (defined as those sources not detected at 1.4GHz but detected in the SDSS with r'<22) rises with increasing z. For mass and colour matched samples, sub-1.5 L_{K}^{*} and super-1.5 L_{K}^{*} radio-detected galaxies have 0.89\pm0.18 and 0.49\pm0.12 times the L250 of their non radio-detected counterparts. We explain these results in terms of the hotter, denser and richer halo environments massive radio-galaxies maintain and are embedded in. These results indicate that all massive radio galaxies (>1.5 L_{K}^{*}) may have systematically lower FIR luminosities (~25%) than their colour-matched non radio-detected counterparts. Compact radio sources (<30kpc) are associated with higher L250s and dust temperatures than their more extended (>30kpc) counterparts. The higher dust temperature suggests that this may be attributed to enhanced SFRs, but whether this is directly or indirectly due to radio activity (e.g. jet induced or merger-driven SF) is as yet unknown.

preprint2013arXiv

Herschel/SPIRE Sub-millimeter Spectra of Local Active Galaxies

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. We detect CO transitions from J_up = 4 to 12, as well as the two [CI] fine structure lines at 492 and 809 GHz and the [NII] 461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n_H2 ~ 10^3.2 - 10^3.9 cm^-3 and T_kin ~ 300 - 800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, however it is possibly related to the star-formation activity in these galaxies. Our modeling of the [CI] emission suggests that it is produced in cold (T_kin < 30 K) and dense (n_H2 > 10^3 cm^-3) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J=1-0 transition at 1232 GHz is detected in absorption in UGC05101 and in emission in NGC7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC7130.

preprint2013arXiv

New Approaches to Probing Minkowski Functionals

We generalize the concept of the ordinary skew-spectrum to probe the effect of non-Gaussianity on the morphology of Cosmic Microwave Background (CMB) maps in several domains: in real-space (where they are commonly known as cumulant-correlators), and in harmonic and needlet bases. The essential aim is to retain more information than normally contained in these statistics, in order to assist in determining the source of any measured non-Gaussianity, in the same spirit as Munshi & Heavens (2010) skew-spectra were used to identify foreground contaminants to the CMB bispectrum in Planck data. Using a perturbative series to construct the Minkowski Functionals (MFs), we provide a pseudo-Cl based approach in both harmonic and needlet representations to estimate these spectra in the presence of a mask and inhomogeneous noise. Assuming homogeneous noise we present approx- imate expressions for error covariance for the purpose of joint estimation of these spectra. We present specific results for four different models of primordial non-Gaussianity local, equilateral, orthogonal and enfolded models, as well as non-Gaussianity caused by unsubtracted point sources. Closed form results of next-order corrections to MFs too are obtained in terms of a quadruplet of kurt-spectra. We also use the method of modal decomposition of the bispectrum and trispectrum to reconstruct the MFs as an alternative method of reconstruction of morphological properties of CMB maps. Finally, we introduce the odd-parity skew-spectra to probe the odd-parity bispectrum and its impact on the morphology of the CMB sky. Although developed for the CMB, the generic results obtained here can be useful in other areas of cosmology.

preprint2013arXiv

The extragalactic background light from the measurements of the attenuation of high-energy gamma-ray spectrum

The attenuation of high-energy gamma-ray spectrum due to the electron-positron pair production against the extragalactic background light (EBL) provides an indirect method to measure the EBL of the universe. We use the measurements of the absorption features of the gamma-rays from blazars as seen by Fermi Gamma-ray Space Telescope to explore the EBL flux density and constrain the EBL spectrum, star formation rate density (SFRD) and photon escape fraction from galaxies out to z=6. Our results are basically consistent with the existing determinations of the quantities. We find a larger photon escape fraction at high redshifts, especially at z=3, compared to the result from the recent Ly-alpha measurements. Our SFRD result is consistent with the data from both gamma-ray burst and UV observations in 1-sigma level. However, the average SFRD we obtain at z>~3 matches the gamma-ray data better than the UV data. Thus our SFRD result at z>~6 favors that it is sufficiently high enough to reionize the universe.

preprint2013arXiv

The rapid assembly of an elliptical galaxy of 400 billion solar masses at a redshift of 2.3

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

preprint2013arXiv

The roles of star formation and AGN activity of IRS sources in the HerMES fields

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

preprint2012arXiv

A Conditional Luminosity Function Model of the Cosmic Far-Infrared Background Anisotropy Power Spectrum

The cosmic far-infrared background (CFIRB) is expected to be generated by faint, dusty star-forming galaxies during the peak epoch of galaxy formation. The anisotropy power spectrum of the CFIRB captures the spatial distribution of these galaxies in dark matter halos and the spatial distribution of dark matter halos in the large-scale structure. Existing halo models of CFIRB anisotropy power spectrum are either incomplete or lead to halo model parameters that are inconsistent with the galaxy distribution selected at other wavelengths. Here we present a conditional luminosity function approach to describe the far-IR bright galaxies. We model the 250 um luminosity function and its evolution with redshift and model-fit the CFIRB power spectrum at 250 um measured by the Herschel Space Observatory. We introduce a redshift dependent duty-cycle parameter so that we are able to estimate the typical duration of the dusty star formation process in the dark matter halos as a function of redshifts. We find the duty cycle of galaxies contributing to the far-IR background is 0.3 to 0.5 with a dusty star-formation phase lasting for \sim0.3-1.6 Gyrs. This result confirms the general expectation that the far-IR background is dominated by star-forming galaxies in an extended phases, not bright starbursts that are driven by galaxy mergers and last \sim10-100 Myrs. The halo occupation number for satellite galaxies has a power-law slope that is close to unity over 0<z<4. We find that the minimum halo mass for dusty, star-forming galaxies with L_250>10^{10} L_Sun is 2\times10^{11}M_Sun and 3\times 10^{10}M_Sun at z=1 and 2, respectively. Integrating over the galaxy population with L_250>10^{9} L_Sun, we find that the cosmic density of dust residing in the dusty, star-forming galaxies responsible for the background anisotropies Ω_{dust}\sim3\times10^{-6} to 2\times10^{-5}.

preprint2012arXiv

A measurement of the intrahalo light fraction with near-infrared background anisotropies

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

preprint2012arXiv

A Population of Dust-Rich Quasars at z ~ 1.5

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

preprint2012arXiv

Amplitudes of thermal and kinetic Sunyaev-Zel'dovich signals from small-scale CMB anisotropies

While the arcminute-scale Cosmic Microwave Background (CMB) anisotropies are due to secondary effects, point sources dominate the total anisotropy power spectrum. At high frequencies the point sources are primarily in the form of dusty, star-forming galaxies. Both Herschel and Planck have recently measured the anisotropy power spectrum of cosmic infrared background (CIB) generated by dusty, star-forming galaxies from degree to sub-arcminute angular scales, including the non-linear clustering of these galaxies at multipoles of 3000 to 6000 relevant to CMB secondary anisotropy studies. We scale the CIB angular power spectra to CMB frequencies and interpret the combined WMAP-7 year and arcminute-scale Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) CMB power spectra measurements to constrain the Sunyaev-Zel'dovich (SZ) effects. Allowing the CIB clustering amplitude to vary, we constrain the amplitudes of thermal and kinetic SZ power spectra at 150 GHz.

preprint2012arXiv

Can dust emission be used to map the interstellar medium in high-redshift galaxies? Results from the Herschel Reference Survey

It has often been suggested that an alternative to the standard CO/21-cm method for estimating the mass of the interstellar medium (ISM) in a galaxy might be to estimate the mass of the ISM from the continuum dust emission. In this paper, we investigate the potential of this technique using Herschel observations of ten galaxies in the Herschel Reference Survey and in the Herschel Virgo Cluster Survey. We show that the emission detected by Herschel is mostly from dust that has a temperature and emissivity index similar to that of dust in the local ISM in our galaxy, with the temperature generally increasing towards the centre of each galaxy. We calibrate the dust method using the CO and 21-cm observations to provide an independent estimate of the mass of hydrogen in each galaxy, solving the problem of the uncertain `X factor' for the molecular gas by minimizing the dispersion in the ratio of the masses estimated using the two methods. With the calibration for the dust method and the estimate of the X-factor produced in this way, the dispersion in the ratio of the two gas masses is 30%, which gives an upper limit on the fundamental accuracy of the dust method. The calibration we obtain for the dust method is very similar to an independent Herschel measurement for M31 and to the calibration for the Milky Way from Planck measurements.

preprint2012arXiv

CANDELS: The Contribution of the Observed Galaxy Population to Cosmic Reionization

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 6<z<8. These galaxies were selected from new deep near-infrared HST imaging from the CANDELS, HUDF09 and ERS programs. In contrast to the majority of previous analyses, which assume that the distribution of galaxy ultraviolet (UV) luminosities follows a Schechter distribution, and that the distribution continues to luminosities far below our observable limit, we investigate the contribution to reionization from galaxies which we can observe, free from these assumptions. We find that the observable population of galaxies can sustain a fully reionized IGM at z=6, if the average ionizing photon escape fraction (f_esc) is ~30%. A number of previous studies have measured UV luminosity densities at these redshifts that vary by 5X, with many concluding that galaxies could not complete reionization by z=6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high values of f_esc were present. The observed UV luminosity density from our observed galaxy samples at z=7-8 is not sufficient to maintain a fully reionized IGM unless f_esc>50%. Combining our observations with constraints on the emission rate of ionizing photons from Ly-alpha forest observations at z=6, we can constrain f_esc<34% (2-sigma) if the observed galaxies are the only contributors to reionization, or <13% (2-sigma) if the luminosity function extends to M_UV = -13. These escape fractions are sufficient to complete reionization by z=6. These constraints imply that the volume ionized fraction of the IGM becomes less than unity at z>7, consistent with a number of complementary reionization probes. If faint galaxies dominate reionization, future JWST observations will probe deep enough to see them, providing an indirect constraint on the ionizing photon escape fraction [abridged].

preprint2012arXiv

GOODS-Herschel and CANDELS: The Morphologies of Ultraluminous Infrared Galaxies at z~2

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

preprint2012arXiv

HerMES: Candidate Gravitationally Lensed Galaxies and Lensing Statistics at Submillimeter Wavelengths

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

preprint2012arXiv

Intensity Mapping of Lyman-alpha Emission During the Epoch of Reionization

We calculate the absolute intensity and anisotropies of the Lyman-alpha radiation field present during the epoch of reionization. We consider emission from both galaxies and the intergalactic medium (IGM) and take into account the main contributions to the production of Lyman-alpha photons: recombinations, collisions, continuum emission from galaxies and scattering of Lyman-n photons in the IGM. We find that the emission from individual galaxies dominates over the IGM with a total Lyman-alpha intensity (times frequency) of about (1.43-3.57)x10^{-8} erg s^{-1} cm^{-2} sr^{-1} at a redshift of 7. This intensity level is low so it is unlikely that the Lyman-αbackground during reionization can be established by an experiment aiming at an absolute background light measurement. Instead we consider Lyman-αintensity mapping with the aim of measuring the anisotropy power spectrum which has rms fluctuations at the level of 1 x 10^{-16} [erg s^[-1} cm^{-2} sr^{-1}]^2 at a few Mpc scales. These anisotropies could be measured with a spectrometer at near-IR wavelengths from 0.9 to 1.4 μm with fields in the order of 0.5 to 1 sq. degrees. We recommend that existing ground-based programs using narrow band filters also pursue intensity fluctuations to study statistics on the spatial distribution of faint Lyman-αemitters. We also discuss the cross-correlation signal with 21 cm experiments that probe HI in the IGM during reionization. A dedicated sub-orbital or space-based Lyman-αintensity mapping experiment could provide a viable complimentary approach to probe reionization, when compared to 21 cm experiments, and is likely within experimental reach.

preprint2012arXiv

Spatially resolved stellar, dust and gas properties of the post-interacting Whirlpool Galaxy system

Using infrared imaging from the Herschel Space Observatory, observed as part of the VNGS, we investigate the spatially resolved dust properties of the interacting Whirlpool galaxy system (NGC 5194 and NGC 5195), on physical scales of 1 kpc. Spectral energy distribution modelling of the new infrared images in combination with archival optical, near- through mid-infrared images confirms that both galaxies underwent a burst of star formation ~370-480 Myr ago and provides spatially resolved maps of the stellar and dust mass surface densities. The resulting average dust-to-stellar mass ratios are comparable to other spiral and spheroidal galaxies studied with Herschel, with NGC 5194 at log M(dust)/M(star)= -2.5+/-0.2 and NGC 5195 at log M(dust)/M(star)= -3.5+/-0.3. The dust-to-stellar mass ratio is constant across NGC 5194 suggesting the stellar and dust components are coupled. In contrast, the mass ratio increases with radius in NGC 5195 with decreasing stellar mass density. Archival mass surface density maps of the neutral and molecular hydrogen gas are also folded into our analysis. The gas-to-dust mass ratio, 94+/-17, is relatively constant across NGC 5194. Somewhat surprisingly, we find the dust in NGC 5195 is heated by a strong interstellar radiation field, over 20 times that of the ISRF in the Milky Way, resulting in relatively high characteristic dust temperatures (~30 K). This post-starburst galaxy contains a substantial amount of low-density molecular gas and displays a gas-to-dust ratio (73+/-35) similar to spiral galaxies. It is unclear why the dust in NGC 5195 is heated to such high temperatures as there is no star formation in the galaxy and its active galactic nucleus is 5-10 times less luminous than the one in NGC 5194, which exhibits only a modest enhancement in the amplitude of its ISRF.

preprint2012arXiv

Spitzer IRAC identification of Herschel-ATLAS SPIRE sources

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

preprint2012arXiv

Submillimetre line spectrum of the Seyfert galaxy NGC1068 from the Herschel-SPIRE Fourier Transform Spectrometer

The first complete submillimetre spectrum (190-670um) of the Seyfert 2 galaxy NGC1068 has been observed with the SPIRE Fourier Transform Spectrometer onboard the {\it Herschel} Space Observatory. The sequence of CO lines (Jup=4-13), lines from water, the fundamental rotational transition of HF, two o-H_2O+ lines and one line each from CH+ and OH+ have been detected, together with the two [CI] lines and the [NII]205um line. The observations in both single pointing mode with sparse image sampling and in mapping mode with full image sampling allow us to disentangle two molecular emission components, one due to the compact circum-nuclear disk (CND) and one from the extended region encompassing the star forming ring (SF-ring). Radiative transfer models show that the two CO components are characterized by density of n(H_2)=10^4.5 and 10^2.9 cm^-3 and temperature of T=100K and 127K, respectively. The comparison of the CO line intensities with photodissociation region (PDR) and X-ray dominated region (XDR) models, together with other observational constraints, such as the observed CO surface brightness and the radiation field, indicate that the best explanation for the CO excitation of the CND is an XDR with density of n(H_2) 10^4 cm^-3 and X-ray flux of 9 erg s^-1 cm^-2, consistent with illumination by the active galactic nucleus, while the CO lines in the SF-ring are better modeled by a PDR. The detected water transitions, together with those observed with the \her \sim PACS Spectrometer, can be modeled by an LVG model with low temperature (T_kin \sim 40K) and high density (n(H_2) in the range 10^6.7-10^7.9 cm^-3).

preprint2012arXiv

The dust energy balance in the edge-on spiral galaxy NGC 4565

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

preprint2012arXiv

The Impact of Assuming Flatness in the Determination of Neutrino Properties from Cosmological Data

Cosmological data have provided new constraints on the number of neutrino species and the neutrino mass. However these constraints depend on assumptions related to the underlying cosmology. Since a correlation is expected between the number of effective neutrinos N_{eff}, the neutrino mass \sum m_ν, and the curvature of the universe Ω_k, it is useful to investigate the current constraints in the framework of a non-flat universe. In this paper we update the constraints on neutrino parameters by making use of the latest cosmic microwave background (CMB) data from the ACT and SPT experiments and consider the possibility of a universe with non-zero curvature. We first place new constraints on N_{eff} and Ω_k, with N_{eff} = 4.03 +/- 0.45 and 10^3 Ω_k = -4.46 +/- 5.24. Thus, even when Ω_k is allowed to vary, N_{eff} = 3 is still disfavored with 95% confidence. We then investigate the correlation between neutrino mass and curvature that shifts the 95% upper limit of \sum m_ν< 0.45 eV to \sum m_ν< 0.95 eV. Thus, the impact of assuming flatness in neutrino cosmology is significant and an essential consideration with future experiments.

preprint2012arXiv

The Infrared Properties of Sources Matched in the WISE all-sky and Herschel-ATLAS Surveys

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

preprint2012arXiv

The Near-IR Background Intensity and Anisotropies During The Epoch of Reionization

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The absolute intensity and the anisotropies of the near-IR background provide an observational probe of the first-light galaxies and their spatial distribution. We estimate the extragalactic background light intensity during reionization by accounting for the stellar and nebular emission from first-light galaxies. We require the UV photon density from these galaxies to generate a reionization history that is consistent with the optical depth to electron scattering from cosmic microwave background measurements. We also require the bright-end luminosity function of galaxies in our models to reproduce the measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 0.1 to 0.3 nW m$^{-2}$ sr$^{-1}$ at the peak of its spectrum at $\sim$ 1.1 $μ$m. We also discuss the anisotropy power spectrum of the near-IR background using a halo model to describe the galaxy distribution. We compare our predictions for the anisotropy power spectrum to existing measurements from deep near-IR imaging data from {\it Spitzer}/IRAC, {\it Hubble}/NICMOS, and {\it AKARI}. The predicted rms fluctuations at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing measurements. While strong arguments have been made that the measured fluctuations do not have an origin involving faint low-redshift galaxies, we find that the existing measurements are also incompatible with an origin during the era of reionization. The measured near-IR background anisotropies remain unexplained and could be associated with an unidentified non-astrophysical origin.

preprint2011arXiv

Constraints on the SZ Power Spectrum on Degree Angular Scales in WMAP Data

The Sunyaev-Zel'dovich (SZ) effect has a distinct spectral signature that allows its separation from fluctuations in the cosmic microwave background (CMB) and foregrounds. Using CMB anisotropies measured in Wilkinson Microwave Anisotropy Probe's five-year maps, we constrain the SZ fluctuations at large, degree angular scales corresponding to multipoles in the range from 10 to 400. We provide upper bounds on SZ fluctuations at multipoles greater than 50, and find evidence for a hemispherically asymmetric signal at ten degrees angular scales. The amplitude of the detected signal cannot be easily explained with the allowed number density and temperature of electrons in the Galactic halo. We have failed to explain the excess signal as a residual from known Galactic foregrounds or instrumental uncertainties such as 1/f-noise.

preprint2011arXiv

Future weak lensing constraints in a dark coupled universe

Coupled cosmologies can predict values for the cosmological parameters at low redshifts which may differ substantially from the parameters values within non-interacting cosmologies. Therefore, low redshift probes, as the growth of structure and the dark matter distribution via galaxy and weak lensing surveys constitute a unique tool to constrain interacting dark sector models. We focus here on weak lensing forecasts from future Euclid and LSST-like surveys combined with the ongoing Planck cosmic microwave background experiment. We find that these future data could constrain the dimensionless coupling to be smaller than a few $\times 10^{-2}$. The coupling parameter $ξ$ is strongly degenerate with the cold dark matter energy density $Ω_{c}h^2$ and the Hubble constant $H_0$.These degeneracies may cause important biases in the cosmological parameter values if in the universe there exists an interaction among the dark matter and dark energy sectors.

preprint2011arXiv

Intensity Mapping of the [CII] Fine Structure Line during the Epoch of Reionization

The atomic CII fine-structure line is one of the brightest lines in a typical star-forming galaxy spectrum with a luminosity ~ 0.1% to 1% of the bolometric luminosity. It is potentially a reliable tracer of the dense gas distribution at high redshifts and could provide an additional probe to the era of reionization. By taking into account of the spontaneous, stimulated and collisional emission of the CII line, we calculate the spin temperature and the mean intensity as a function of the redshift. When averaged over a cosmologically large volume, we find that the CII emission from ionized carbon in individual galaxies is larger than the signal generated by carbon in the intergalactic medium (IGM). Assuming that the CII luminosity is proportional to the carbon mass in dark matter halos, we also compute the power spectrum of the CII line intensity at various redshifts. In order to avoid the contamination from CO rotational lines at low redshift when targeting a CII survey at high redshifts, we propose the cross-correlation of CII and 21-cm line emission from high redshifts. To explore the detectability of the CII signal from reionization, we also evaluate the expected errors on the CII power spectrum and CII-21 cm cross power spectrum based on the design of the future milimeter surveys. We note that the CII-21 cm cross power spectrum contains interesting features that captures physics during reionization, including the ionized bubble sizes and the mean ionization fraction, which are challenging to measure from 21-cm data alone. We propose an instrumental concept for the reionization CII experiment targeting the frequency range of $\sim$ 200 to 300 GHz with 1, 3 and 10 meter apertures and a bolometric spectrometer array with 64 independent spectral pixels with about 20,000 bolometers.

preprint2011arXiv

Lensed galaxies in CANDELS

We present results from a search for gravitationally lensed galaxies present in the Hubble Space Telescope (HST)/Wide Field Camera-3 (WFC3) images of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). We present one bona fide lens system in UDS and two compact lens candidates in the GOODS-S field. The lensing system in UDS involves two background galaxies, one at z=1.847 lensed to an arc and a counterimage, and the second at a photometric redshift of z=2.32^{+0.10}_{-0.06} lensed to a double image. We reconstruct the lensed sources in the source plane and find in each of the two cases the sources can be separated to a pair of galaxies. The sources responsible for the arc are compact with effective radii of 0.3 to 0.4 kpc in WFC3 J_{125}-band and a total stellar mass and a star-formation rate of 2.1_{-0.4}^{+2.4} times 10^7 M_sun and 2.3_{-1.7}^{+ 0.6} M_sun yr^{-1}, respectively.The abnormally high H_{160}-band flux of this source is likely due to OIII emission lines with a rest-frame equivalent width about 700 Angstroms for OIII 5007 Angstroms. The sources responsible for the double image have corresponding values of about 0.4 to 0.5 kpc, 1.4_{-0.8}^{+1.9} times 10^9 M_sun, and 8.7_{-7.0}^{+11.1} M_sun yr^{-1}. Once completed CANDELS is expected to contain about 15 lensing systems and will allow statistical studies on both lensing mass profiles and z ~ 2 lensed galaxies.

preprint2011arXiv

Measuring Light from the Epoch of Reionization with CIBER, the Cosmic Infrared Background Experiment

Ultraviolet emission from the first generation of stars in the Universe ionized the intergalactic medium in a process which was completed by z~6; the wavelength of these photons has been redshifted by (1+z) into the near infrared today and can be measured using instruments situated above the Earth's atmosphere. First flying in February 2009, the Cosmic Infrared Background Experiment (CIBER) comprises four instruments housed in a single reusable sounding rocket borne payload. CIBER will measure spatial anisotropies in the extragalactic IR background caused by cosmological structure from the epoch of reionization using two broadband imaging instruments, make a detailed characterization of the spectral shape of the IR background using a low resolution spectrometer, and measure the absolute brightness of the Zodical light foreground with a high resolution spectrometer in each of our six science fields. This paper presents the scientific motivation for CIBER and details of its first two flights, including a review of the published scientific results from the first flight and an outlook for future reionization science with CIBER data.

preprint2011arXiv

Observations of Arp 220 using Herschel-SPIRE: An Unprecedented View of the Molecular Gas in an Extreme Star Formation Environment

We present Herschel SPIRE-FTS observations of Arp~220, a nearby ULIRG. The FTS continuously covers 190 -- 670 microns, providing a good measurement of the continuum and detection of several molecular and atomic species. We detect luminous CO (J = 4-3 to 13-12) and water ladders with comparable total luminosity; very high-J HCN absorption; OH+, H2O+, and HF in absorption; and CI and NII. Modeling of the continuum yields warm dust, with T = 66 K, and an unusually large optical depth of ~5 at 100 microns. Non-LTE modeling of the CO shows two temperature components: cold molecular gas at T ~ 50 K and warm molecular gas at T ~1350 K. The mass of the warm gas is 10% of the cold gas, but dominates the luminosity of the CO ladder. The temperature of the warm gas is in excellent agreement with H2 rotational lines. At 1350 K, H2 dominates the cooling (~20 L_sun/M_sun) in the ISM compared to CO (~0.4 L_sun/M_sun). We found that only a non-ionizing source such as the mechanical energy from supernovae and stellar winds can excite the warm gas and satisfy the energy budget of ~20 L_sun/M_sun. We detect a massive molecular outflow in Arp 220 from the analysis of strong P-Cygni line profiles observed in OH+, H2O+, and H2O. The outflow has a mass > 10^{7} M_sun and is bound to the nuclei with velocity < 250 km/s. The large column densities observed for these molecular ions strongly favor the existence of an X-ray luminous AGN (10^{44} ergs/s) in Arp 220.

preprint2011arXiv

Probing Reionization with Intensity Mapping of Molecular and Fine Structure Lines

We propose observations of the molecular gas distribution during the era of reionization. At $z \sim $ 6 to 8, $^{12}$CO$(J=1-0)$ line intensity results in a mean brightness temperature of about 0.5 $μ$K with rms fluctuations of 0.1 $μ$K at 1 to 10 Mpc spatial scales, corresponding to 30 arcminute angular scales. This intensity fluctuations can be mapped with an interferometer, similar to existing and planned 21-cm background experiments, but operating at $\sim$ 12 to 17 GHz. We discuss the feasibility of detecting the cross-correlation between HI and CO molecular gas since such a cross-correlation has the advantage that it will be independent of systematics and most foregrounds in each of the 21-cm and CO$(1-0)$ line experiments. Additional instruments tuned to higher-order transitions of the CO molecule or an instrument operating with high spectral resolution at mm-wavelengths targeting 158 $μ$m CII could further improve the reionization studies with molecular gas. The combined 21-cm and CO line observations has the potential to establish the relative distribution of gas in the inter-galactic medium and molecular gas that are clumped in individual first-light galaxies that are closely connected to the formation of massive stars in these galaxies.

preprint2011arXiv

Sub-millimetre galaxies reside in dark matter halos with masses greater than 3x10^11 solar masses

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

preprint2011arXiv

The Average Physical Properties and Star Formation Histories of the UV-Brightest Star-Forming Galaxies at z~3.7

[Abridged] We investigate the average physical properties and star formation histories of the most UV-luminous star-forming galaxies at z~3.7. Our results are derived from analyses of the average spectral energy distributions (SEDs), constructed from stacked optical to infrared photometry, of a sample of the 1,902 most UV-luminous star-forming galaxies found in 5.3 square degrees of the NOAO Deep Wide-Field Survey. We bin the sample according to UV luminosity, and find that the shape of the average SED in the rest-frame optical and infrared is fairly constant with UV luminosity: i.e., more UV luminous galaxies are, on average, also more luminous at longer wavelengths. In the rest-UV, however, the spectral slope (measured at 0.13-0.28 um) rises steeply with the median UV luminosity from -1.8 at L L* to -1.2 in the brightest bin (L~4-5L*). We use population synthesis analyses to derive the average physical properties of these galaxies and find that: (1) L_UV, and thus star formation rates (SFRs), scale closely with stellar mass such that more UV-luminous galaxies are also more massive; (2) The median ages indicate that the stellar populations are relatively young (200-400 Myr) and show little correlation with UV luminosity; and (3) More UV-luminous galaxies are dustier than their less-luminous counterparts, such that L~4-5L* galaxies are extincted up to A(1600)=2 mag while L L* galaxies have A(1600)=0.7-1.5 mag. Based on these observations, we argue that the average star formation histories of UV-luminous galaxies are better described by models in which SFR increases with time in order to simultaneously reproduce the tight correlation between the observed SFR and stellar mass, and the universally young ages of these galaxies. We demonstrate the potential of measurements of the SFR-M* relation at multiple redshifts to discriminate between simple models of star formation histories.

preprint2011arXiv

The OH line contamination of 21 cm intensity fluctuation measurements for z=1~4

The large-scale structure of the Universe can be mapped with unresolved intensity fluctuations of the 21 cm line. The power spectrum of the intensity fluctuations has been proposed as a probe of the baryon acoustic oscillations at low to moderate redshifts with interferometric experiments now under consideration. We discuss the contamination to the low-redshift 21 cm intensity power spectrum generated by the 18 cm OH line since the intensity fluctuations of the OH line generated at a slightly higher redshift contribute to the intensity fluctuations observed in an experiment. We assume the OH megamaser luminosity is correlated with the star formation rate, and use the simulation to estimate the OH signal and the spatial anisotropies. We also use a semi-analytic simulation to predict the 21 cm power spectrum. At z=1 to 3, we find that the OH contamination could reach 0.1% to 1% of the 21 cm rms fluctuations at the scale of the first peak of the baryon acoustic oscillation. When z>3 the OH signal declines quickly, so that the contamination on the 21-cm becomes negligible at high redshifts.

preprint2010arXiv

A Constraint On the Integrated Mass Power Spectrum out to z = 1100 from Lensing of the Cosmic Microwave Background

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

preprint2010arXiv

A Halo Model of Local IRAS Galaxies Selected at 60 Micron Using Conditional Luminosity Functions

Using conditional luminosity functions (CLFs) which encode the luminosity distribution of galaxies as a function of halo mass, we construct a halo model of IRAS galaxies selected at 60 micron. An abundance matching technique is used to link galaxy luminosity to the host halo mass. The shape of the mass - light relation at 60 micron is different from those derived at r-, K- and B-band. This is because the 60 micron LF can not be fitted by a Schechter function with a sharp exponential cutoff. We then seek the parameters in the CLFs that best fit the LF and power spectrum. We find that the predicted galaxy bias as a function of L60 from the best-fit model agrees well with the clustering measurements. At the faint end of the LF where quiescent star-forming galaxies dominate, most IRAS galaxies are central galaxies in halos of M >~ 10^{10} h^{-1} M_sun but a non-negligible fraction are satellites typically hosted in more massive halos. The majority of IRAS galaxies with L60 >~ 10^{10} h^{-2} L_sun are M82 type starbursts which are central galaxies hosted in halos of M >~ 10^{12.5} h^{-1} M_sun. In comparison, optical galaxies generally reside in much more massive halos. The rate of change in L60 (an indicator of recent star formation) as a function of halo mass at M >~ 10^{12.5} h^{-1} M_sun is much larger that d L_optical / dM or d L_NIR / dM indicating the existence of physical mechanisms which are very efficient in converting cold gas into stars, possibly dynamical effects arising from interactions or mergers. We further calculate the space density of major mergers for halos massive enough to host ultraluminous infrared galaxies (ULIRGs) using the mean merger rate derived from the Millennium simulations. Compared to the space density of local ULIRGs, it implies that either the majority of major mergers at z~0 do not lead to ULIRGs or the ULIRG phase is relatively short.

preprint2010arXiv

A Measurement of Cubic-Order Primordial Non-Gaussianity (g_{NL} and τ_{NL}) With WMAP 5-Year Data

We measure two higher-order power spectra involving weighted cubic and squared temperature anisotropy maps from WMAP 5-year data to study the trispectrum generated by primordial non-Gaussianity. Using these measurements combined with Gaussian and noise simulations, we constrain the cubic order non-Gaussianity parameters τ_{NL}, and g_{NL}. With V+W-band data out to l_{max}=600, we find -7.4 < g_{\rm NL}/10^5 < 8.2 and -0.6 < τ_{\rm NL}/10^4 < 3.3 improving the previous COBE-based limit on τ_{\rm NL} < 10^8 nearly four orders of magnitude with WMAP. We find that the ratio of trispectrum to bispectrum amplitude as captured by the ratio of τ_{\rm Nl}/(6f_{\rm NL}/5)^2 ranges from -3 to 21 at the 95% confidence level.

preprint2010arXiv

CMB Constraints on Primordial non-Gaussianity from the Bispectrum (f_{NL}) and Trispectrum (g_{NL} and τ_{NL}) and a New Consistency Test of Single-Field Inflation

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (f_{NL}) and fourth-order (g_{NL} and τ_{NL}) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of τ_{NL} and g_{NL} with WMAP 5-year data, our study provides the first forecasts for future constraints on g_{NL}. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio A_{NL} = τ_{NL}/(6f_{NL}/5)^2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to l_{max}=600 we constrain the cubic order non-Gaussianity parameters τ_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < τ_{NL}/10^4 < 3.3 improving the previous COBE-based limit on τ_{NL} < 10^8 nearly four orders of magnitude with WMAP.

preprint2010arXiv

CMB Neutrino Mass Bounds and Reionization

Current cosmic microwave background (CMB) bounds on the sum of the neutrino masses assume a sudden reionization scenario described by a single parameter that determines the onset of reionization. We investigate the bounds on the neutrino mass in a more general reionization scenario based on a principal component approach. We found the constraint on the sum of the neutrino masses from CMB data can be relaxed by a $\sim$40% in a generalized reionization scenario. Moreover, the amplitude of the r.m.s. mass fluctuations $σ_8$ is also considerably lower providing a better consistency with a low amplitude of the Sunyaev-Zel'dovich signal.

preprint2010arXiv

Constraints on primordial non-Gaussianity from WMAP7 and Luminous Red Galaxies power spectrum and forecast for future surveys

We place new constraints on the primordial local non-Gaussianity parameter f_NL using recent Cosmic Microwave Background anisotropy and galaxy clustering data. We model the galaxy power spectrum according to the halo model, accounting for a scale dependent bias correction proportional to f_NL/k^2. We first constrain f_NL in a full 13 parameters analysis that includes 5 parameters of the halo model and 7 cosmological parameters. Using the WMAP7 CMB data and the SDSS DR4 galaxy power spectrum, we find f_NL=171\pm+140 at 68% C.L. and -69<f_NL<+492 at 95% C.L.. We discuss the degeneracies between f_NL and other cosmological parameters. Including SN-Ia data and priors on H_0 from Hubble Space Telescope observations we find a stronger bound: -35<f_NL<+479 at 95% C.L.. We also fit the more recent SDSS DR7 halo power spectrum data finding, for a Λ-CDM+f_NL model, f_NL=-93\pm128 at 68% C.L. and -327<f_{NL}<+177 at 95% C.L.. We finally forecast the constraints on f_NL from future surveys as EUCLID and from CMB missions as Planck showing that their combined analysis could detect f_NL\sim 5.

preprint2010arXiv

Cosmology with Photometric Surveys of Type Ia Supernovae

We discuss the extent to which photometric measurements alone can be used to identify Type Ia supernovae (SNIa) and to determine redshift and other parameters of interest for cosmological studies. We fit the light curve data of the type expected from a survey such as the one planned with Large Synoptic Survey Telescope (LSST) and also to remove the contamination from the core-collapse supernovae to SNIa samples. We generate 1000 SNIa mock flux data for each of the LSST filters based on existing design parameters, then use a Markov Chain Monte-Carlo (MCMC) analysis to fit for the redshift, apparent magnitude, stretch factor and the phase of the SNIa. We find that the model fitting works adequately well when the true SNe redshift is below 0.5, while at $z < 0.2$ the accuracy of the photometric data is almost comparable with spectroscopic measurements of the same sample. We discuss the contamination of Type Ib/c (SNIb/c) and Type II supernova (SNII) on the SNIa data set. We find it is easy to distinguish the SNII through the large $χ^2$ mismatch when fitting to photometric data with Ia light curves. This is not the case for SNIb/c. We implement a statistical method based on the Bayesian estimation in order to statistically reduce the contamination from SNIb/c for cosmological parameter measurements from the whole SNe sample. The proposed statistical method also evaluate the fraction of the SNIa in the total SNe data set, which provides a valuable guide to establish the degree of contamination.

preprint2010arXiv

HerMES: Halo Occupation Number and Bias Properties of Dusty Galaxies from Angular Clustering Measurements

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

preprint2010arXiv

Herschel-ATLAS: Statistical Properties of Galactic Cirrus in the GAMA-9 Hour Science Demonstration Phase Field

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $β$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = â H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

preprint2010arXiv

Higher-order Statistics of Weak Lensing Shear and Flexion

Owing to their more extensive sky coverage and tighter control on systematic errors, future deep weak lensing surveys should provide a better statistical picture of the dark matter clustering beyond the level of the power spectrum. In this context, the study of non-Gaussianity induced by gravity can help tighten constraints on the background cosmology by breaking parameter degeneracies, as well as throwing light on the nature of dark matter, dark energy or alternative gravity theories. Analysis of the shear or flexion properties of such maps is more complicated than the simpler case of the convergence due to the spinorial nature of the fields involved. Here we develop analytical tools for the study of higher-order statistics such as the bispectrum (or trispectrum) directly using such maps at different source redshift. The statistics we introduce can be constructed from cumulants of the shear or flexions, involving the cross-correlation of squared and cubic maps at different redshifts. Typically, the low signal-to-noise ratio prevents recovery of the bispectrum or trispectrum mode by mode. We define power spectra associated with each multi- spectra which compresses some of the available information of higher order multispectra. We show how these can be recovered from a noisy observational data even in the presence of arbitrary mask, which introduces mixing between Electric (E-type) and Magnetic (B-type) polarization, in an unbiased way. We also introduce higher order cross-correlators which can cross-correlate lensing shear with different tracers of large scale structures.

preprint2010arXiv

Impact of general reionization scenarios on extraction of inflationary parameters

Determination of whether the Harrison--Zel'dovich spectrum for primordial scalar perturbations is consistent with observations is sensitive to assumptions about the reionization scenario. In light of this result, we revisit constraints on inflationary models using more general reionization scenarios. While the bounds on the tensor-to-scalar ratio are largely unmodified, when different reionization schemes are addressed, hybrid models are back into the inflationary game. In the general reionization picture, we reconstruct both the shape and amplitude of the inflaton potential. We find a broader spectrum of potential shapes when relaxing the simple reionization restriction. An upper limit of $10^{16}$ GeV to the amplitude of the potential is found, regardless of the assumptions on the reionization history.

preprint2010arXiv

Primordial Non-Gaussianity from a Joint Analysis of Cosmic Microwave Background Temperature and Polarization

We explore a systematic approach to the analysis of primordial non-Gaussianity using fluctuations in temperature and polarization of the Cosmic Microwave Background (CMB). Following Munshi & Heavens (2009), we define a set of power-spectra as compressed forms of the bispectrum and trispectrum derived from CMB temperature and polarization maps; these spectra compress the information content of the corresponding full multispectra and can be useful in constraining early Universe theories. We generalize the standard pseudo-C_l estimators in such a way that they apply to these spectra involving both spin-0 and spin-2 fields, developing explicit expressions which can be used in the practical implementation of these estimators. While these estimators are suboptimal, they are nevertheless unbiased and robust hence can provide useful diagnostic tests at a relatively small computational cost. We next consider approximate inverse-covariance weighting of the data and construct a set of near-optimal estimators based on that approach. Instead of combining all available information from the entire set of mixed bi- or trispectra, i.e multispectra describing both temperature and polarization information, we provide analytical constructions for individual estimators, associated with particular multispectra. The bias and scatter of these estimators can be computed using Monte-Carlo techniques. Finally, we provide estimators which are completely optimal for arbitrary scan strategies and involve inverse covariance weighting; we present the results of an error analysis performed using a Fisher-matrix formalism at both the one-point and two-point level.

preprint2010arXiv

Testing General Relativity with Current Cosmological Data

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large scale structure and the deflection of light by that structure. We clarify the relations between several different model independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2010arXiv

The Impact of Secondary non-Gaussianities in the CMB on Cosmological Parameter Estimation

We consider corrections to the underlying cosmology due to secondary contributions from weak gravitational lensing, the integrated Sachs-Wolfe effect, and the Sunyaev-Zel'dovich effect contained in the trispectrum. We incorporate these additional contributions to the covariance of a binned angular power spectrum of temperature anisotropies in the analysis of current and prospective data sets. Although recent experiments such as ACBAR and CBI are not particularly sensitive to these additional non-Gaussian effects, the interpretation of Planck and CMBPol anisotropy spectra will require an accounting of non-Gaussian covariance leading to a degradation in cosmological parameter estimates by up to 20% and 30%, respectively.

preprint2010arXiv

The Near-Infrared Sky Surveyor

[NIRSS is one of three concepts that contributed to the Wide-Field Infrared Survey Telescope (WFIRST) mission advocated by the Decadal Survey.] Operating beyond the reaches of the Earth's atmosphere, free of its limiting absorption and thermal background, the Near-Infrared Sky Surveyor (NIRSS) will deeply map the entire sky at near-infrared wavelengths, thereby enabling new and fundamental discoveries ranging from the identification of extrasolar planets to probing the reionization epoch by identifying thousands of quasars at z>10. NIRSS will directly address the NASA scientific objective of studying cosmic origins by using a 1.5-meter telescope to reach full-sky 0.2 uJy (25.6 mag AB) sensitivities in four passbands from 1 to 4 microns in a 4-yr mission. At the three shorter passbands (1 - 2.5 microns), the proposed depth is comparable to the deepest pencil-beam surveys done to date and is 3000 times more sensitive than the only previous all-sky near-infrared survey, 2MASS. At the longest passband (3.5 micron), which is not feasible from the ground, NIRSS will be 500 times more sensitive than WISE. NIRSS fills a pivotal gap in our knowledge of the celestial sphere, is a natural complement to WISE, and is well matched to the next generation of deep (0.1 uJy), wide-area (>2 pi ster), ground-based optical surveys (LSST and Pan-Starrs). With the high thermal backgrounds of ground-based infrared observations, a near-infrared full sky survey at sub-uJy sensitivity is only feasible from space.

preprint2010arXiv

Weak lensing and dark energy: the impact of dark energy on nonlinear dark matter clustering

We examine the influence of percent-level dark energy corrections to the nonlinear matter power spectrum on constraints of the dark energy equation of state from future weak lensing probes. We explicitly show that a poor approximation (off by > 10%) to the nonlinear corrections causes a > 1 sigma bias on the determination of the dark energy equation of state. Future weak lensing surveys must therefore incorporate dark energy modifications to the nonlinear matter power spectrum accurate to the percent-level, to avoid introducing significant bias in their measurements. For the WMAP5 cosmology, the more accurate power spectrum is more sensitive to dark energy properties, resulting in a factor of two improvement in dark energy equation of state constraints. We explore the complementary constraints on dark energy from future weak lensing and supernova surveys. A space-based, JDEM-like survey measures the equation of state in five independent redshift bins to ~10%, while this improves to ~5% for a wide-field ground-based survey like LSST. These constraints are contingent upon our ability to control weak lensing systematic uncertainties to the sub-percent level.

preprint2009arXiv

A Measurement of Primordial Non-Gaussianity Using WMAP 5-Year Temperature Skewness Power Spectrum

We constrain the primordial non-Gaussianity parameter of the local model f_{NL} using the skewness power spectrum associated with the two-to-one cumulant correlator of cosmic microwave background temperature anisotropies. This bispectrum-related power spectrum was constructed after weighting the temperature map with the appropriate window functions to form an estimator that probes the multipolar dependence of the underlying bispectrum associated with the primordial non-Gaussianity. We also estimate a separate skewness power spectrum sensitive more strongly to unresolved point sources. When compared to previous attempts at measuring the primordial non-Gaussianity with WMAP data, our estimators have the main advantage that we do not collapse information to a single number. When model fitting the two-to-one skewness power spectrum we make use of bispectra generated by the primordial non-Gaussianity, radio point sources, and lensing-secondary correlation. We analyze Q, V and W-band WMAP 5-year data using the KQ75 mask out to l_{max}=600. Using V and W-band data and marginalizing over model parameters related to point sources and lensing-secondary bispectrum, our overall and preferred constraint on f_{NL} is 11.0\pm23.7 at the 68% confidence level (-36.4 < f_{NL} < 58.4 at 95% confidence). We find no evidence for a non-zero value of f_{NL} even marginally at the 1σlevel.

preprint2009arXiv

CMB Lensing Constraints on Dark Energy and Modified Gravity Scenarios

Weak gravitational lensing leaves a characteristic imprint on the cosmic microwave background temperature and polarization angular power spectra. Here we investigate the possible constraints on the integrated lensing potential from future CMB angular spectra measurements expected from Planck and EPIC. We find that Planck and EPIC will constrain the amplitude of the integrated projected potential responsible for lensing at 6% and 1% level, respectively with very little sensitivity to the shape of the lensing potential. We discuss the implications of such a measurement in constraining dark energy and modified gravity scalar-tensor theories. We then discuss the impact of a wrong assumption on the weak lensing potential amplitude on cosmological parameter inference.

preprint2009arXiv

Cosmic Infrared Background ExpeRiment (CIBER): A Probe of Extragalactic Background Light from Reionization

The Cosmic Infrared Background ExpeRiment (CIBER) is a rocket-borne absolute photometry imaging and spectroscopy experiment optimized to detect signatures of first-light galaxies present during reionization in the unresolved IR background. CIBER-I consists of a wide-field two-color camera for fluctuation measurements, a low-resolution absolute spectrometer for absolute EBL measurements, and a narrow-band imaging spectrometer to measure and correct scattered emission from the foreground zodiacal cloud. CIBER-I was successfully flown on February 25th, 2009 and has one more planned flight in early 2010. We propose, after several additional flights of CIBER-I, an improved CIBER-II camera consisting of a wide-field 30 cm imager operating in 4 bands between 0.5 and 2.1 microns. It is designed for a high significance detection of unresolved IR background fluctuations at the minimum level necessary for reionization. With a FOV 50 to 2000 times largerthan existing IR instruments on satellites, CIBER-II will carry out the definitive study to establish the surface density of sources responsible for reionization.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

New Optimised Estimators for the Primordial Trispectrum

Cosmic microwave background studies of non-Gaussianity involving higher-order multispectra can distinguish between early universe theories that predict nearly identical power spectra. However, the recovery of higher-order multispectra is difficult from realistic data due to their complex response to inhomogeneous noise and partial sky coverage, which are often difficult to model analytically. A traditional alternative is to use one-point cumulants of various orders, which collapse the information present in a multispectrum to one number. The disadvantage of such a radical compression of the data is a loss of information as to the source of the statistical behaviour. A recent study by Munshi & Heavens (2009) has shown how to define the skew spectrum (the power spectra of a certain cubic field, related to the bispectrum) in an optimal way and how to estimate it from realistic data. The skew spectrum retains some of the information from the full configuration-dependence of the bispectrum, and can contain all the information on non-Gaussianity. In the present study, we extend the results of the skew spectrum to the case of two degenerate power-spectra related to the trispectrum. We also explore the relationship of these power-spectra and cumulant correlators previously used to study non-Gaussianity in projected galaxy surveys or weak lensing surveys. We construct nearly optimal estimators for quick tests and generalise them to estimators which can handle realistic data with all their complexity in a completely optimal manner. We show how these higher-order statistics and the related power spectra are related to the Taylor expansion coefficients of the potential in inflation models, and demonstrate how the trispectrum can constrain both the quadratic and cubic terms.

preprint2009arXiv

No Evidence for Dark Energy Dynamics from a Global Analysis of Cosmological Data

We use a variant of principal component analysis to investigate the possible temporal evolution of the dark energy equation of state, w(z). We constrain w(z) in multiple redshift bins, utilizing the most recent data from Type Ia supernovae, the cosmic microwave background, baryon acoustic oscillations, the integrated Sachs-Wolfe effect, galaxy clustering, and weak lensing data. Unlike other recent analyses, we find no significant evidence for evolving dark energy; the data remains completely consistent with a cosmological constant. We also study the extent to which the time-evolution of the equation of state would be constrained by a combination of current- and future-generation surveys, such as Planck and the Joint Dark Energy Mission.

preprint2009arXiv

Non-Gaussianity in WMAP Data Due to the Correlation of CMB Lensing Potential with Secondary Anisotropies

We measure the skewness power spectrum of the Cosmic Microwave Background (CMB) anisotropies optimized for a detection of the secondary bispectrum generated by the correlation of the CMB lensing potential with integrated Sachs-Wolfe effect and the Sunyaev-Zel'dovich effect. The covariance of our measurements is generated by Monte-Carlo simulations of Gaussian CMB fields with noise properties consistent with Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data. When interpreting multi-frequency measurements we also take into account the confusion resulting from unresolved radio point sources. We analyze Q, V and W-band WMAP 5-year raw and foreground-cleaned maps using the KQ75 mask out to l_max=600. We find no significant evidence for a non-zero non-Gaussian signal from the lensing-secondary correlation in all three bands and we constrain the overall amplitude of the cross power spectrum between CMB lensing potential and the sum of SZ and ISW fluctuations to be 0.42 \pm 0.86 and 1.19 \pm 0.86 in combined V and W-band raw and foreground-cleaned maps provided by the WMAP team, respectively. The point source amplitude at the bispectrum level measured with this skewness power spectrum is higher than previous measurements of point source non-Gaussianity. We also consider an analysis where we also account for the primordial non-Gaussianity in addition to lensing-secondary bispectrum and point sources. The focus of this paper is on secondary anisotropies. Consequently the estimator is not optimised for primordial non-Gaussianity and the limit we find on local non-Gaussianity from the foreground-cleaned V+W maps is f_NL = -13 \pm 62, when marginalized over point sources and lensing-ISW/SZ constributions to the total bispectrum.

preprint2009arXiv

Secondary non-Gaussianity and Cross-Correlation Analysis

We develop optimised estimators of two sorts of power spectra for fields defined on the sky, in the presence of partial sky coverage. The first is the cross-power spectrum of two fields on the sky; the second is the skew spectrum of three fields. These can probe the Integrated Sachs Wolfe Effect (ISW) at large angular scales and the Sunyaev Zéldovich (SZ) effect from hot gas in clusters at small angular scales. The skew spectrum, recently introduced by Munshi & Heavens (2009), is an optimised statistic which can be tuned to study a particular form of non-Gaussianity, such as may arise in the early Universe, but which retains information on the nature of non-Gaussianity. In this paper we develop the mathematical formalism for the skew spectrum of 3 different fields. When applied to the CMB, this allows us to explore the contamination of the skew spectrum by secondary sources of CMB fluctuations. Considering the three-point function, the study of the bispectrum provides valuable information regarding cross-correlation of secondaries with lensing of CMB with much higher significance compared to just the study involving CMB sky alone. We use our analytical models to study specific cases of cosmological interest which include cross-correlating CMB with various large scale tracers to probe ISW and SZ effects for cross spectral analysis and use the formalism to study the signal-to-noise ratio for detection of the weak lensing of the CMB by cross-correlating it with different tracers as well as point sources for CMB experiments such as Planck (abridged).

preprint2008arXiv

CMBPol Mission Concept Study: Gravitational Lensing

Gravitational lensing of the cosmic microwave background by large-scale structure in the late universe is both a source of cosmological information and a potential contaminant of primordial gravity waves. Because lensing imprints growth of structure in the late universe on the CMB, measurements of CMB lensing will constrain parameters to which the CMB would not otherwise be sensitive, such as neutrino mass. If the instrumental noise is sufficiently small (<~ 5 uK-arcmin), the gravitational lensing contribution to the large-scale B-mode will be the limiting source of contamination when constraining a stochastic background of gravity waves in the early universe, one of the most exciting prospects for future CMB polarization experiments. High-sensitivity measurements of small-scale B-modes can reduce this contamination through a lens reconstruction technique that separates the lensing and primordial contributions to the B-mode on large scales. A fundamental design decision for a future CMB polarization experiment such as CMBpol is whether to have coarse angular resolution so that only the large-scale B-mode (and the large-scale E-mode from reionization) is measured, or high resolution to additionally measure CMB lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, and required control of beam systematics.

preprint2008arXiv

Implications of Two Type Ia Supernova Populations for Cosmological Measurements

Recent work suggests that Type Ia supernovae (SNe) are composed of two distinct populations: prompt and delayed. By explicitly incorporating properties of host galaxies, it may be possible to target and eliminate systematic differences between these two putative populations. However, any resulting {\em post}-calibration shift in luminosity between the components will cause a redshift-dependent systematic shift in the Hubble diagram. Utilizing an existing sample of 192 SNe Ia, we find that the average luminosity difference between prompt and delayed SNe is constrained to be $(4.5 \pm 8.9)%$. If the absolute difference between the two populations is 0.025 mag, and this is ignored when fitting for cosmological parameters, then the dark energy equation of state (EOS) determined from a sample of 2300 SNe Ia is biased at $\sim1σ$. By incorporating the possibility of a two-population systematic, this bias can be eliminated. However, assuming no prior on the strength of the two-population effect, the uncertainty in the best-fit EOS is increased by a factor of 2.5, when compared to the equivalent sample with no underlying two-population systematic. To avoid introducing a bias in the EOS parameters, or significantly degrading the measurement accuracy, it is necessary to control the post-calibration luminosity difference between prompt and delayed SN populations to better than 0.025 mag.

preprint2008arXiv

Measuring dark energy spatial inhomogeneity with supernova data

The gravitational lensing distortion of distant sources by the large-scale distribution of matter in the Universe has been extensively studied. In contrast, very little is known about the effects due to the large-scale distribution of dark energy. We discuss the use of Type Ia supernovae as probes of the spatial inhomogeneity and anisotropy of dark energy. We show that a shallow, almost all-sky survey can limit rms dark energy fluctuations at the horizon scale down to a fractional energy density of ~10^-4

preprint2007arXiv

Searching for primordial black hole dark matter with pulsar timing arrays

We discuss the possibility of detecting the presence of primordial black holes (PBHs), such as those that might account for galactic dark matter, using modification of pulsar timing residuals when PBHs pass within ~1000 AU and impart impulse accelerations to the Earth. With this technique, PBHs with masses around 10^{25} g (~0.1 lunar mass) can be detected. Currently, the constraints on the abundance of such dark matter candidates are weak. A 30 year-long monitoring campaign with the proposed Square Kilometer Array (SKA) can rule out a PBH fraction more than ~1/10 in the solar neighborhood in the form of dark matter with mass ~10^{25} g.

preprint2007arXiv

The Warm Spitzer Mission: Prospects for Studies of the Distant Universe

IRAC excels at detecting distant objects. Due to a combination of the shapes of the spectral energy distributions of galaxies and the low background achieved from space, IRAC reaches greater depth in comparable exposure time at 3.6 and 4.5 micron than any ground- or space-based facility currently can at 2.2 micron. Furthermore, the longer wavelengths probed by IRAC enable studies of the rest-frame optical and near-infrared light of galaxies and AGN to much higher redshift than is possible from the ground. This white paper explores the merits of different survey strategies for studying the distant universe during the warm mission. A three-tiered approach serves a wide range of science goals and uses the spacecraft effectively: 1) an ultra-deep survey of ~0.04 square degrees to a depth of ~250 hrs (in conjunction with an HST/WFC3 program), to study the Universe at 7<z<14; 2) a survey of ~2 square degrees to the GOODS depth of 20 hrs, to identify luminous galaxies at z>6 and characterize the relation between the build-up of dark matter halos and their constituent galaxies at 2<z<6, and 3) a 500 square degree survey to the SWIRE depth of 120 s, to systematically study large scale structure at 1<z<2 and characterize high redshift AGN. One or more of these programs could conceivably be implemented by the SSC, following the example of the Hubble Deep Field campaigns. As priorities in this field continuously shift it is also crucial that a fraction of the exposure time remains unassigned, thus enabling science that will reflect the frontiers of 2010 and beyond rather than those of 2007.

preprint2006arXiv

Angular Diameter Distance Measurement with Galaxy Clustering in the Multipole Space

The shape of the angular power spectrum of galaxies in the linear regime is defined by the horizon size at the matter-radiation equality. When calibrated by cosmic microwave background measurements, the shape of the clustering spectrum can be used as a standard ruler to estimate angular diameter distance as a function of redshift at which galaxy clustering is measured. We apply the proposed cosmological test of Cooray et al. (2001) to a recent set of luminous red galaxy angular clustering spectra from the Sloan Digital Sky Survey between redshifts of 0.2 and 0.6. Using the overall shape of the clustering power spectrum in the linear regime, w e measure comoving angular diameter distances to eight redshift bins by marginalizing over the bias factors that determine the overall amplitude of the clustering spectrum in each of the bins. The Hubble constant consistent with these distance estimates is 68.5^+6.7_-6.1 km s^-1 Mpc-1 at the 68% confidence level. We comment on the expected improvements with future surveys and the potential to measure dark energy parameters with this method.

preprint2006arXiv

Deciphering Inflation with Gravitational Waves: Cosmic Microwave Background Polarization vs. Direct Detection with Laser Interferometers

A detection of the primordial gravitational wave background is considered to be the ``smoking-gun '' evidence for inflation. While super-horizon waves are probed with cosmic microwave background (CMB) polarization, the relic background will be studied with laser interferometers. The long lever arm spanned by the two techniques improves constraints on the inflationary potential and validation of consistency relations expected under inflation. If gravitational waves with a tensor-to-scalar amplitude ratio greater than 0.01 are detected by the CMB, then a direct detection experiment with a sensitivity consistent with current concept studies should be pursued vigorously. If no primordial tensors are detected by the CMB, a direct detection experiment to understand the simplest form of inflation must have a sensitivity improved by two to three orders of magnitude over current plans.

preprint2006arXiv

Modifications to the Cosmic 21-cm Background Frequency Spectrum by Scattering via electrons in Galaxy Clusters

The cosmic 21-cm background frequency spectrum related to the spin-flip transition of neutral Hydrogen present during and before the era of reionization is rich in features associated with physical processes that govern transitions between the two spin states. The intervening electrons in foreground galaxy clusters inversely Compton scatter the 21-cm background spectrum and modify it just as the cosmic microwave background (CMB) spectrum is modified by inverse-Compton scattering. Towards typical galaxy clusters at low redshifts, the resulting modification is a few tenths milli-Kelvin correction to the few tens milli-Kelvin temperature of 21-cm signal relative to that of the cosmic microwave background black body spectrum. The modifications are mostly associated with sharp changes in the cosmic 21-cm background spectrum such as due to the onset of a Lyman-$α$ radiation field or heating of neutral gas. Though low frequency radio interferometers that are now planned for 21-cm anisotropy measurements are insensitive to the mean 21-cm spectrum, differential observations of galaxy clusters with these interferometers can be utilized to indirectly establish global features in the 21-cm frequency spectrum. We discuss the feasibility to detect the spectrum modified by clusters and find that for upcoming interferometers, while a detection towards an individual cluster is challenging, one can average signals over a number of clusters, selected based on the strength of the Sunyave-Zel'dovich effect at high radio frequencies involving CMB scattering alone, to establish the mean 21-cm spectrum.

preprint2006arXiv

The Born and Lens-Lens Corrections to Weak Gravitational Lensing Angular Power Spectra

We revisit the estimation of higher order corrections to the angular power spectra of weak gravitational lensing. Extending a previous calculation of Cooray and Hu, we find two additional terms to the fourth order in potential perturbations of large-scale structure corresponding to corrections associated with the Born approximation and the neglect of line-of-sight coupling of two foreground lenses in the standard first order result. These terms alter the convergence ($κκ$), the lensing shear E-mode ($εε$), and their cross-correlation ($κε$) power spectra on large angular scales, but leave the power spectra of the lensing shear B-mode ($ββ$) and rotational ($ωω$) component unchanged as compared to previous estimates. The new terms complete the calculation of corrections to weak lensing angular power spectra associated with both the Born approximation and the lens-lens coupling to an order in which the contributions are most significant. Taking these features together, we find that these corrections are unimportant for any weak lensing survey, including for a full sky survey limited by cosmic variance.

preprint2005arXiv

Cluster Merger Variance and the Luminosity Gap Statistic

The presence of multiple luminous galaxies in clusters can be explained by the finite time over which a galaxy sinks to the center of the cluster and merges with the the central galaxy. The simplest measurable statistic to quantify the dynamical age of a system of galaxies is the luminosity (magnitude) gap, which is the difference in photometric magnitude between the two most luminous galaxies. We present a simple analytical estimate of the luminosity gap distribution in groups and clusters as a function of dark matter halo mass. The luminosity gap is used to define "fossil" groups; we expect the fraction of fossil systems to exhibit a strong and model-independent trend with mass: ~1-3% of massive clusters and ~5-40% of groups should be fossil systems. We also show that, on cluster scales, the observed intrinsic scatter in the central galaxy luminosity-halo mass relation can be ascribed to dispersion in the merger histories of satellites within the cluster. We compare our predictions to the luminosity gap distribution in a sample of 730 clusters in the Sloan Digital Sky Survey C4 Catalog and find good agreement. This suggests that theoretical excursion set merger probabilities and the standard theory of dynamical segregation are valid on cluster scales.

preprint2005arXiv

Dissipationless Merging and the Assembly of Central Galaxies

We reanalyze the galaxy-mass correlation function measured by the Sloan Digital Sky Survey to obtain host dark matter halo masses at galaxy and galaxy group scales. We extend the data to galaxy clusters in the 2MASS catalog and study the relation between central galaxy luminosity and halo mass. While the central galaxy luminosity scales as ~M^{0.7-0.8} at low masses, the relation flattens to ~M^{<0.3} above ~4x10^{13} M_sun. The total luminosity of galaxies in the halo, however, continues to grow as a power-law ~M^{0.8-0.9}. Starting from the hypothesis that the central galaxies grow by merging ("galactic cannibalism"), we develop a simple model for the evolution of their luminosities as a consequence of the accretion of satellite galaxies. The luminosity-mass relation flattens when the time scale on which dynamical friction induces orbital decay in the satellite galaxies exceeds the age of the dark matter halo. Then, the growth of the central galaxy is suppressed as it can cannibalize only the rare, massive satellite galaxies. The model takes the dependence of the total luminosity of galaxies in a halo on its mass and the global galaxy luminosity function as input, and reproduces the observed central galaxy luminosity-mass relation over three decades in halo mass, (10^{12}-10^{15}) M_sun. The success of the model suggests that gas cooling and subsequent star formation did not play an important role in the final assembly of central galaxies from sub-L_star precursors.

preprint2005arXiv

Problems with Pencils: Lensing Covariance of Supernova Distance Measurements

While luminosity distances from Type Ia supernovae (SNe) provide a powerful probe of cosmological parameters, the accuracy with which these distances can be measured is limited by cosmic magnification due to gravitational lensing by the intervening large-scale structure. Spatial clustering of foreground mass fluctuations leads to correlated errors in distance estimates from SNe. By including the full covariance matrix of supernova distance measurements, we show that a future survey covering more than a few square degrees on the sky, and assuming a total of ~2000 SNe, will be largely unaffected by covariance noise. ``Pencil beam'' surveys with small fields of view, however, will be prone to the lensing covariance, leading to potentially significant degradations in cosmological parameter estimates. For a survey with 30 arcmin mean separation between SNe, lensing covariance leads to a ~45% increase in the expected errors in dark energy parameters compared to fully neglecting lensing, and a ~20% increase compared to including just the lensing variance. Given that the lensing covariance is cosmology dependent and cannot be mapped out sufficiently accurately with direct weak lensing observations, surveys with small mean SN separation must incorporate the effects of lensing covariance, including its dependence on the cosmological parameters.

preprint2004arXiv

Polarization Signals of the 21 cm Background from the Era of Reionization

While emission and absorption lines of the 21 cm spin-flip transition of neutral Hydrogen are intrinsically unpolarized, a magnetic field creates left- and right-handed polarized components through the Zeeman effect. Here we consider the resulting polarization of the redshifted 21 cm background from the intergalactic medium before reionization. The polarization is detectable in regions with a strong gradient in the mean brightness temperature. In principle, this can open a new window on the evolution of intergalactic magnetic fields. One possible approach is an extended integration of an individual target during this era, such as the Mpc-scale HII regions inferred to surround quasars at z ~ 6.5. The differential intensity between the two polarization states can be used as a probe of the magnetic field at the edge of the HII region. We estimate that the SKA could (ignoring systematics) detect B \~ 200 (10) \muG coherent over several kiloparsecs with an observational bandwidth of 100 (2) kHz. Beyond individual sources, the statistical properties of wide-field 21 cm polarization maps, such as the angular power spectrum, can be used to constrain the large-scale magnetic field. In this case, the SKA can detect $B \sim 100 μ$G fields coherent over many megaparsecs. The magnetic field can be measured in any epoch over which the 21 cm background changes rapidly (for example because the ionized fraction or spin temperature change). Although the resulting constraints with SKA are relatively weak compared to theoretical expectations, they nevertheless offer a unique direct probe of magnetic fields in the high-redshift universe.

preprint2004arXiv

Uncorrelated Estimates of Dark Energy Evolution

Type Ia supernova data have recently become strong enough to enable, for the first time, constraints on the time variation of the dark energy density and its equation of state. Most analyses, however, are using simple two or three-parameter descriptions of the dark energy evolution, since it is well known that allowing more degrees of freedom introduces serious degeneracies. Here we present a method to produce uncorrelated and nearly model-independent band power estimates of the equation of state of dark energy and its density as a function of redshift. We apply the method to recently compiled supernova data. Our results are consistent with the cosmological constant scenario, in agreement with other analyses that use traditional parameterizations, though we find marginal (2-sigma) evidence for w(z) < -1 at z < 0.2. In addition to easy interpretation, uncorrelated, localized band powers allow intuitive and powerful testing of the constancy of either the energy density or equation of state. While we have used relatively coarse redshift binning suitable for the current set of about 150 supernovae, this approach should reach its full potential in the future, when applied to thousands of supernovae found from ground and space, combined with complementary information from other cosmological probes.

preprint2003arXiv

CMB-induced Cluster Polarization as a Cosmological Probe

Scattering of the temperature anisotropy quadrupole by free electrons in galaxy clusters leads to a secondary polarization signal in the cosmic microwave background (CMB) fluctuations. At low redshifts, the temperature quadrupole contains a significant contribution from the integrated Sachs-Wolfe (ISW) effect associated with the growth of density fluctuations. Using polarization data from a sample of clusters over a wide range in redshift, one can statistically establish the presence of the ISW effect and determine its redshift evolution. Given the strong dependence of the ISW effect on the background cosmology, cluster polarization can eventually be used as a powerful probe of dark energy. As a further application, we also discuss how it might be used to understand the potential lack of power on large scales.

preprint2003arXiv

Lensing Reconstruction with CMB Temperature and Polarization

Weak gravitational lensing by intervening large-scale structure induces a distinct signature in the cosmic microwave background (CMB) that can be used to reconstruct the weak-lensing displacement map. Estimators for individual Fourier modes of this map can be combined to produce an estimator for the lensing-potenial power spectrum. The naive estimator for this quantity will be biased upwards by the uncertainty associated with reconstructing individual modes; we present an iterative scheme for removing this bias. The variance and covariance of the lensing-potenial power spectrum estimator are calculated and evaluated numerically in a $Λ$CDM universe for Planck and future polarization-sensitive CMB experiments.

preprint2002arXiv

Halo Models of Large Scale Structure

We review the formalism and applications of the halo-based description of nonlinear gravitational clustering. In this approach, all mass is associated with virialized dark matter halos; models of the number and spatial distribution of the halos, and the distribution of dark matter within each halo, are used to provide estimates of how the statistical properties of large scale density and velocity fields evolve as a result of nonlinear gravitational clustering. We first describe the model, and demonstrate its accuracy by comparing its predictions with exact results from numerical simulations of nonlinear gravitational clustering. We then present several astrophysical applications of the halo model: these include models of the spatial distribution of galaxies, the nonlinear velocity, momentum and pressure fields, descriptions of weak gravitational lensing, and estimates of secondary contributions to temperature fluctuations in the cosmic microwave background.

preprint2002arXiv

Is the Cosmic Microwave Background Circularly Polarized?

The primordial anisotropies of the cosmic microwave background (CMB) are linearly polarized via Compton-scattering. The Faraday conversion process during the propagation of polarized CMB photons through regions of the large-scale structure containing magnetized relativistic plasma, such as galaxy clusters, will lead to a circularly polarized contribution. Though the resulting Stokes-V parameter is of order 10^-9 at frequencies of 10 GHz, the contribution can potentially reach the total Stokes-U at low frequencies due to the cubic dependence on the wavelength. In future, the detection of circular polarization of CMB can be used as a potential probe of the physical properties associated with relativistic particle populations in large-scale structures.

preprint2002arXiv

Second Moment of Halo Occupation Number

The halo approach to large scale structure provides a physically motivated model to understand clustering properties of galaxies. An important aspect of the halo model involves a description on how galaxies populate dark matter halos or what is now called the halo occupation distribution. We discuss a way in which clustering information, especially in the non-linear regime, can be used to determine moments of this halo occupation number. We invert the non-linear part of the real space power spectrum from the PSCz survey to determine the second moment of the halo occupation distribution in a model independent manner. The precise measurement of higher order correlations can eventually be used to determine successive higher order moments of this distribution.

preprint2002arXiv

Second Order Corrections to Weak Lensing by Large-Scale Structure

We calculate corrections to the power spectrum predictions of weak lensing by large scale structure due to higher order effects in the gravitational potential. Using a perturbative approach to third order in transverse displacements, we calculate a second order correction to the angular power spectra of E and B mode shear and convergence resulting from dropping the so-called Born approximation, where one integrates along the unperturbed photon path. We also consider a correction to the power spectra from the coupling between lenses at different redshifts. Both effects generate B-mode shear and the latter also causes a net rotation of the background galaxy images. We show all these corrections are at least two orders of magnitude below the convergence or E-mode power and hence relevant only to future ultra high precision measurements. These analytical calculations are consistent with previous numerical estimates and validate the use of current large scale structure weak lensing predictions for cosmological studies and future use of B-modes as a monitor of systematic effects.

preprint2002arXiv

Small-Scale Cosmic Microwave Background Polarization from Reionization

We discuss fluctuations in the cosmic microwave background (CMB) polarization due to scattering from reionized gas at low redshifts. Polarization is produced by re-scattering of the primordial temperature anisotropy quadrupole and of the kinematic quadrupole that arises from gas motion transverse to the line of sight. We show that both effects produce equal E- and B-parity polarization, and are, in general, several orders of magnitude below the dominant polarization contributions at the last scattering surface to E-modes or the gravitational-lensing contribution to B-modes at intermediate redshifts. These effects are also several orders of magnitude below the B polarization due to lensing even after subtraction with higher-order correlations, and are thus too small to constitute a background for searches for the polarization signature of inflationary gravitational waves.

Asantha Cooray

What is connected

Connect this record

See the researcher in context

Building this map preview

128 published item(s)

Anisotropies of Cosmic Optical and Near-IR Background from China Space Station Telescope (CSST)

Gas Phase Metallicities of Local Ultra-Luminous Infrared Galaxies Follow Normal Star-Forming Galaxies

Snowmass 2021 CMB-S4 White Paper

A Long Stream of Metal-Poor Cool Gas around a Massive Starburst Galaxy at z = 2.67

Measurements of the Zodiacal Light Absolute Intensity through Fraunhofer Absorption Line Spectroscopy with CIBER

Modelling high-resolution ALMA observations of strongly lensed dustystar forming galaxies detected by Herschel

Polarization spectrum of near infrared zodiacal light observed with CIBER

Cosmological constraints from line intensity mapping with interlopers

Emergence of an Ultra-Red Ultra-Massive Galaxy Cluster Core at $z=4$

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

The Case for Probe-class NASA Astrophysics Missions

Three Dusty Star Forming Galaxies at $z\sim1.5$: Mergers and Disks on the Main Sequence

Photometric Redshift Estimation with Galaxy Morphology using Self-Organizing Maps

Tomographic Intensity Mapping versus Galaxy Surveys: Observing the Universe in H-alpha emission with new generation instruments

Axion decay and anisotropy of near-IR extragalactic background light

Cosmic Dawn Intensity Mapper

Extragalactic Background Light: Measurements and Applications

Multi-wavelength lens reconstruction of a $\textit{Planck}$ $\&$ $\textit{Herschel}$-detected star-bursting galaxy

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

Spectral Line De-confusion in an Intensity Mapping Survey

Star Formation in Luminous Quasars at 2<z<3

Toward a tomographic analysis of the cross-correlation between Planck CMB lensing and H-ATLAS galaxies

ALMA maps the Star-Forming Regions in a Dense Gas Disk at z~3

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Astronomy Job Crisis

Cosmology with the SPHEREX All-Sky Spectral Survey

Dust energy balance study of two edge-on spiral galaxies in the Herschel-ATLAS survey

Extinction and nebular line properties of a Herschel-selected lensed dusty starburst AT z=1.027

Far-infrared observations of an unbiased sample of gamma-ray burst host galaxies

H-ATLAS/GAMA: Quantifying the Morphological Evolution of the Galaxy Population Using Cosmic Calorimetry

Planck Trispectrum Constraints on Primordial Non-Gaussianity at Cubic Order

Prospects for detecting CII emission during the Epoch of Reionization

Spitzer Imaging of Strongly-Lensed Herschel-Selected Dusty Star Forming Galaxies

Ultraviolet Luminosity Density of the Universe During the Epoch of Reionization

Cross-Correlation of Near and Far-Infrared Background Anisotropies as Traced by Spitzer and Herschel

Dusty Star-Forming Galaxies at High Redshift

Foreground contamination in Ly-alpha intensity mapping during the epoch of reionization

HerMES: The Rest-Frame UV Emission and A Lensing Model for the z=6.34 Luminous Dusty Starburst Galaxy HFLS3

Herschel-ATLAS and ALMA: HATLAS J142935.3-002836, a lensed major merger at redshift 1.027

Herschel-ATLAS/GAMA:How does the far-IR luminosity function depend on galaxy group properties?

High-resolution, 3D radiative transfer modeling : I. The grand-design spiral galaxy M51

On the Origin of Near-Infrared Extragalactic Background Light Anisotropy

Optical-faint, Far-infrared-bright Herschel Sources in the CANDELS Fields: Ultra-Luminous Infrared Galaxies at z>1 and the Effect of Source Blending

The applicability of FIR fine-structure lines as Star Formation Rate tracers over wide ranges of metallicities and galaxy types

The co-evolution of black hole growth and star formation from a cross-correlation analysis between quasars and the cosmic infrared background

The Herschel Exploitation of Local Galaxy Andromeda (HELGA). VI. The distribution and properties of molecular cloud associations in M31

Constraints on Spatial Variations in the Fine-Structure constant from Planck

Disentangling a group of lensed submm galaxies at z~2.9

Far-infrared Fine-Structure Line Diagnostics of Ultraluminous Infrared Galaxies

Herschel SPIRE-FTS Observations of Excited CO and [CI] in the Antennae (NGC 4038/39): Warm and Cold Molecular Gas

Herschel-ATLAS/GAMA: What determines the far-infrared properties of radio-galaxies?

Herschel/SPIRE Sub-millimeter Spectra of Local Active Galaxies

New Approaches to Probing Minkowski Functionals

The extragalactic background light from the measurements of the attenuation of high-energy gamma-ray spectrum

The rapid assembly of an elliptical galaxy of 400 billion solar masses at a redshift of 2.3

The roles of star formation and AGN activity of IRS sources in the HerMES fields

A Conditional Luminosity Function Model of the Cosmic Far-Infrared Background Anisotropy Power Spectrum

A measurement of the intrahalo light fraction with near-infrared background anisotropies

A Population of Dust-Rich Quasars at z ~ 1.5

Amplitudes of thermal and kinetic Sunyaev-Zel'dovich signals from small-scale CMB anisotropies

Can dust emission be used to map the interstellar medium in high-redshift galaxies? Results from the Herschel Reference Survey

CANDELS: The Contribution of the Observed Galaxy Population to Cosmic Reionization

GOODS-Herschel and CANDELS: The Morphologies of Ultraluminous Infrared Galaxies at z~2

HerMES: Candidate Gravitationally Lensed Galaxies and Lensing Statistics at Submillimeter Wavelengths

Intensity Mapping of Lyman-alpha Emission During the Epoch of Reionization

Spatially resolved stellar, dust and gas properties of the post-interacting Whirlpool Galaxy system

Spitzer IRAC identification of Herschel-ATLAS SPIRE sources

Submillimetre line spectrum of the Seyfert galaxy NGC1068 from the Herschel-SPIRE Fourier Transform Spectrometer

The dust energy balance in the edge-on spiral galaxy NGC 4565

The Impact of Assuming Flatness in the Determination of Neutrino Properties from Cosmological Data

The Infrared Properties of Sources Matched in the WISE all-sky and Herschel-ATLAS Surveys

The Near-IR Background Intensity and Anisotropies During The Epoch of Reionization

Constraints on the SZ Power Spectrum on Degree Angular Scales in WMAP Data

Future weak lensing constraints in a dark coupled universe