Angular power spectra of optical and infrared background anisotropies at wavelengths between 0.5 to 5 $\mu$m are a useful probe of faint sources present during reionization, in addition to faint galaxies and diffuse signals at low redshift. The cross-correlation of these fluctuations with backgrounds at other wavelengths can be used to separate some of these signals. A previous study on the cross-correlation between X-ray and $Spitzer$ fluctuations at 3.6 $\mu$m and 4.5 $\mu$m has been interpreted as evidence for direct collapse blackholes (DCBHs) present at $z > 12$. Here we return to this cross-correlation and study its wavelength dependence from 0.5 to 4.5 $\mu$m using $Hubble$ and $Spitzer$ data in combination with a subset of the 4 Ms $Chandra$ observations in GOODS-S/ECDFS. Our study involves five $Hubble$ bands at 0.6, 0.7, 0.85, 1.25 and 1.6 $\mu$m, and two $Spitzer$-IRAC bands at 3.6 $\mu$m and 4.5 $\mu$m. We confirm the previously seen cross-correlation between 3.6 $\mu$m (4.5 $\mu$m) and X-rays with 3.7$\sigma$ (4.2$\sigma$) and 2.7$\sigma$ (3.7$\sigma$) detections in the soft [0.5-2] keV and hard [2-8] keV X-ray bands, respectively, at angular scales above 20 arcseconds. The cross-correlation of X-rays with $Hubble$ is largely anticorrelated, ranging between the levels of 1.4$-$3.5$\sigma$ for all the $Hubble$ and X-ray bands. This lack of correlation in the shorter optical/NIR bands implies the sources responsible for the cosmic infrared background at 3.6 and 4.5~$\mu$m are at least partly dissimilar to those at 1.6 $\mu$m and shorter.

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.

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 $\mu$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 \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. 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.

We present the cross-correlation between the far-infrared background fluctuations as measured with the Herschel Space Observatory at 250, 350, and 500 {\mu}m and the near-infrared background fluctuations with Spitzer Space Telescope at 3.6 {\mu}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 {\mu}m to 500 {\mu}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 {\mu}m, and (c) intra-halo light, or diffuse stars in dark matter halos, that likely dominates fluctuations at 3.6 {\mu}m. The model is able to reasonably reproduce the auto correlations at each of the far-IR wavelengths and at 3.6 {\mu}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 {\mu}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 {\mu}m to 500 {\mu}m comes from the fact that the galaxies that are primarily contributing to 500 {\mu}m fluctuations peak at a higher redshift than those at 250 {\mu}m.

We present a statistical method based on a maximum likelihood approach to constrain the number counts of extragalactic sources below the nominal flux-density limit of continuum imaging surveys. We extract flux densities from a radio map using positional information from an auxiliary catalogue and show that we can model the number counts of this undetected population down to flux density levels well below the detection threshold of the radio survey. We demonstrate the capabilities that our method will have with future generation wide-area radio surveys by performing simulations over various sky areas. We show that it is possible to accurately constrain the number counts of the simulated distribution down to one tenth of the flux noise rms with just a sky area of 100 deg$^2$. We then test the application of our method using data from the Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST). We extract flux densities from the FIRST map, sensitive to 150 $\mu$Jy/beam (1 $\sigma$), using the positional information from a catalogue in the same field, also acquired at the same frequency, sensitive to 12 $\mu$Jy/beam (1 $\sigma$). Implementing our method, with known source positions, we are able to recover the right differential number counts of the noise-dominated FIRST map fluxes down to a flux density level which is one-tenth the FIRST detection threshold.

We present a measurement of the angular power spectrum of the cosmic far-infrared background (CFIRB) anisotropies in one of the extragalactic fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) at 250, 350 and 500 \mu m bands. Consistent with recent measurements of the CFIRB power spectrum in Herschel-SPIRE maps, we confirm the existence of a clear one-halo term of galaxy clustering on arcminute angular scales with large-scale two-halo term of clustering at 30 arcminutes to angular scales of a few degrees. The power spectrum at the largest angular scales, especially at 250 \mu m, is contaminated by the Galactic cirrus. The angular power spectrum is modeled using a conditional luminosity function approach to describe the spatial distribution of unresolved galaxies that make up the bulk of the CFIRB. Integrating over the dusty galaxy population responsible for the background anisotropies, we find that the cosmic abundance of dust, relative to the critical density, to be between \Omega_dust=10^{-6} and 8 x 10^{-6} in the redshift range z ~ 0-3. This dust abundance is consistent with estimates of the dust content in the Universe using quasar reddening and magnification measurements in the SDSS.

We have developed and characterized an imaging instrument to measure the spatial properties of the diffuse near-infrared extragalactic background light in a search for fluctuations from z > 6 galaxies during the epoch of reionization. The instrument is part of the Cosmic Infrared Background Experiment (CIBER), designed to observe the extragalactic background light above the Earth's atmosphere during a suborbital sounding rocket flight. The imaging instrument incorporates a 2x2 degree field of view, to measure fluctuations over the predicted peak of the spatial power spectrum at 10 arcminutes, and 7"x7" pixels, to remove lower redshift galaxies to a depth sufficient to reduce the low-redshift galaxy clustering foreground below instrumental sensitivity. The imaging instrument employs two cameras with \Delta \lambda / \lambda ~0.5 bandpasses centered at 1.1 and 1.6 microns to spectrally discriminate reionization extragalactic background fluctuations from local foreground fluctuations. CIBER operates at wavelengths where the electromagnetic spectrum of the reionization extragalactic background is thought to peak, and complements fluctuations measurements by AKARI and Spitzer at longer wavelengths. We have characterized the instrument in the laboratory, including measurements of the sensitivity, flat-field response, stray light performance, and noise properties. Several modifications were made to the instrument following a first flight in 2009 February. The instrument performed to specifications in subsequent flights in 2010 July and 2012 March, and the scientific data are now being analyzed.

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $\mu$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Bo\"{o}tes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $\mu$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $\mu$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

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.

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.

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.

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.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.