Search

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

Verifying the fully kinematic nature of the cosmic microwave background (CMB) dipole is of fundamental importance in cosmology. In the standard cosmological model with the Friedman-Lemaitre-Robertson-Walker (FLRW) metric from the inflationary expansion the CMB dipole should be entirely kinematic. Any non-kinematic CMB dipole component would thus reflect the preinflationary structure of spacetime probing the extent of the FLRW applicability. Cosmic backgrounds from galaxies after the matter-radiation decoupling, should have kinematic dipole component identical in velocity with the CMB kinematic dipole. Comparing the two can lead to isolating the CMB non-kinematic dipole. It was recently proposed that such measurement can be done using the near-IR cosmic infrared background (CIB) measured with the currently operating Euclid telescope, and later with Roman. The proposed method reconstructs the resolved CIB, the Integrated Galaxy Light (IGL), from Euclid's Wide Survey and probes its dipole, with a kinematic component amplified over that of the CMB by the Compton-Getting effect. The amplification coupled with the extensive galaxy samples forming the IGL would determine the CIB dipole with an overwhelming signal/noise, isolating its direction to sub-degree accuracy. We develop details of the method for Euclid's Wide Survey in 4 bands spanning 0.6 to 2 mic. We isolate the systematic and other uncertainties and present methodologies to minimize them, after confining the sample to the magnitude range with negligible IGL/CIB dipole from galaxy clustering. These include the required star-galaxy separation, accounting for the extinction correction dipole using the method newly developed here achieving total separation, accounting for the Earth's orbital motion and other systematic effects. (Abridged), Comment: Euclid Key Project paper, A&A, in press

We measured the dipole of the diffuse $\gamma$-ray background (DGB) identifying a highly significant time-independent signal coincidental with that of the Pierre Auger UHECR. The DGB dipole is determined from flux maps in narrow energy bands constructed from 13 years of observations by the Large Area Telescope (LAT) of the {\it Fermi} satellite. The $\gamma$-ray maps were clipped iteratively of sources and foregrounds similar to that done for the cosmic infrared background. The clipped narrow energy band maps were then assembled into one broad energy map out to the given energy starting at $E=2.74$ Gev, where the LAT beam falls below the sky's pixel resolution. Next we consider cuts in Galactic latitude and longitude to probe residual foreground contaminations from the Galactic Plane and Center. In the broad energy range $2.74 < E\leq115.5$ GeV the measured dipoles are stable with respect to the various Galactic cuts, consistent with an extragalactic origin. The $\gamma$-ray sky's dipole/monopole ratio is much greater than that expected from the DGB clustering component and the Compton-Getting effect origin with reasonable velocities. At $\simeq (6.5-7)\%$ it is similar to the Pierre Auger UHECRs with $E_{\rm UHECR}\ge 8$ EeV pointing to a common origin of the two dipoles. However, the DGB flux associated with the found DGB dipole reaches parity with that of the UHECR around $E_{\rm UHECR}\le 1$ EeV, perhaps arguing for a non-cascading mechanism if the DGB dipole were to come from the higher energy UHECRs. The signal/noise of the DGB dipole is largest in the $5-30$ GeV range, possibly suggesting the $\gamma$-photons at these energies are the ones related to cosmic rays., Comment: ApJ Letters, to appear

The Euclid Space Telescope will provide deep imaging at optical and near-infrared wavelengths, along with slitless near-infrared spectroscopy, across ~15,000 sq deg of the sky. Euclid is expected to detect ~12 billion astronomical sources, facilitating new insights into cosmology, galaxy evolution, and various other topics. To optimally exploit the expected very large data set, there is the need to develop appropriate methods and software. Here we present a novel machine-learning based methodology for selection of quiescent galaxies using broad-band Euclid I_E, Y_E, J_E, H_E photometry, in combination with multiwavelength photometry from other surveys. The ARIADNE pipeline uses meta-learning to fuse decision-tree ensembles, nearest-neighbours, and deep-learning methods into a single classifier that yields significantly higher accuracy than any of the individual learning methods separately. The pipeline has `sparsity-awareness', so that missing photometry values are still informative for the classification. Our pipeline derives photometric redshifts for galaxies selected as quiescent, aided by the `pseudo-labelling' semi-supervised method. After application of the outlier filter, our pipeline achieves a normalized mean absolute deviation of ~< 0.03 and a fraction of catastrophic outliers of ~< 0.02 when measured against the COSMOS2015 photometric redshifts. We apply our classification pipeline to mock galaxy photometry catalogues corresponding to three main scenarios: (i) Euclid Deep Survey with ancillary ugriz, WISE, and radio data; (ii) Euclid Wide Survey with ancillary ugriz, WISE, and radio data; (iii) Euclid Wide Survey only. Our classification pipeline outperforms UVJ selection, in addition to the Euclid I_E-Y_E, J_E-H_E and u-I_E,I_E-J_E colour-colour methods, with improvements in completeness and the F1-score of up to a factor of 2. (Abridged), Comment: 37 pages (including appendices), 26 figures; accepted for publication in Astronomy & Astrophysics

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration., Comment: A&A language edited version

While the cosmic microwave background (CMB) dipole is largely assumed entirely kinematic, there appears evidence that a part of it is primordial. Such possibility arises in models implying a tilt, interpreted as a dark flow, across the observable Universe. The kinematic nature of the entire CMB dipole can be probed using the dipole of cosmic backgrounds from galaxies after the last scattering. The near-IR cosmic infrared background (CIB) spectral energy distribution leads to an amplified dipole compared to the CMB. The CIB dipole is affected by galaxy clustering, decreasing with fainter, more distant galaxies, and by Solar System emissions and Galactic dust, which dominate the net CIB cosmological dipole in the optical/near-IR. We propose a technique that enables an accurate measurement of the kinematic near-IR CIB dipole. The CIB, effectively the integrated galaxy light (IGL), would be reconstructed from resolved galaxies in the forthcoming space-borne wide surveys covering four bands 0.9 to 2.5 micron. The galaxies will be sub-selected from the identified magnitude range where the dipole component from galaxy clustering is below the expected kinematic dipole. Using this technique the dipole can be measured in each of the bands at the statistical signal-to-noise S/N>50--100 with the forthcoming Euclid and Roman surveys, isolating CMB dipole's kinematic nature., Comment: MNRAS Letters, in press

We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al. 2022, SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here we utilize SKYSURF's first sky-surface brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL). Our sky-surface brightness measurements have been verified to an accuracy of better than 1%, which when combined with systematic errors associated with HST, results in sky brightness uncertainties of $\sim$2-4% $\simeq$ 0.005 MJy/sr in each image. We put limits on the amount of diffuse EBL in three near-IR filters (F125W, F140W, and F160W) by comparing our preliminary sky measurements of $> 30,000$ images to Zodiacal light models, carefully selecting the darkest images to avoid contamination from stray light. In addition, we investigate the impact that instrumental thermal emission has on our measurements, finding that it has a limited impact on F125W and F140W measurements, whereas uncertainties in the exact thermal state of HST results in significant uncertainties in the level of astrophysical diffuse light in F160W images. When compared to the Kelsall et al. (1998) Zodiacal model, an isotropic diffuse background of $30$ nW m$^{-2}$ sr$^{-1}$ remains, whereas using the Wright (1998) Zodiacal model results in no discernible diffuse background. Based primarily on uncertainties in the foreground model subtraction, we present limits on the amount of diffuse EBL of 29 nW m$^{-2}$ sr$^{-1}$, 40 nW m$^{-2}$ sr$^{-1}$, and 29 nW m$^{-2}$ sr$^{-1}$ for F125W, F140W, and F160W respectively. While this light is generally isotropic, our modeling at this point does not distinguish between a cosmological origin or a Solar System origin (such as a dim, diffuse, spherical cloud of cometary dust)., Comment: Accepted with Windhorst et al. 2022 to AJ. Main figures are Fig. 10 and 11. Comments welcome!

We give an overview and describe the rationale, methods, and testing of the Hubble Space Telescope (HST) Archival Legacy project "SKYSURF." SKYSURF uses HST's unique capability as an absolute photometer to measure the ~0.2-1.7 $\mu$m sky surface brightness (SB) from 249,861 WFPC2, ACS, and WFC3 exposures in ~1400 independent HST fields. SKYSURF's panchromatic dataset is designed to constrain the discrete and diffuse UV to near-IR sky components: Zodiacal Light (ZL; inner Solar System), Kuiper Belt Objects (KBOs; outer Solar System), Diffuse Galactic Light (DGL), and the discrete plus diffuse Extragalactic Background Light (EBL). We outline SKYSURF's methods to: (1) measure sky-SB levels between its detected objects; (2) measure the integrated discrete EBL, most of which comes from AB$\simeq$17-22 mag galaxies; and (3) estimate how much diffuse light may exist in addition to the extrapolated discrete galaxy counts. Simulations of HST WFC3/IR images with known sky-values and gradients, realistic cosmic ray (CR) distributions, and star plus galaxy counts were processed with nine different algorithms to measure the "Lowest Estimated Sky-SB" (LES) in each image between the discrete objects. The best algorithms recover the inserted LES values within 0.2% when there are no image gradients, and within 0.2-0.4% when there are 5-10% gradients. SKYSURF requires non-standard re-processing of these HST images that includes restoring the lowest sky-level from each visit into each drizzled image. We provide a proof of concept of our methods from the WFC3/IR F125W images, where any residual diffuse light that HST sees in excess of the Kelsall et al. (1998) Zodiacal model prediction does not depend on the total object flux that each image contains. This enables us to present our first SKYSURF results on diffuse light in Carleton et al. (2022)., Comment: Accepted to AJ; see accompanying paper Carleton et al. 2022: arXiv:2205.06347. Comments welcome!

(Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 galaxies within the Euclid Deep Survey. This study is based on ~176,000 real galaxies at z=1-8 in a ~0.7 deg2 area selected from the UltraVISTA ultra-deep survey, and ~96,000 mock galaxies with 25.3$\leq$H<27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from the fiducial, 28-band photometry, and fit spectral energy distributions (SEDs) to various combinations of these simulated data. Our study demonstrates that identifying z>6 with Euclid data alone will be very effective, with a z>6 recovery of 91(88)% for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z=1-5.8 contaminants amongst apparent z>6 galaxies as observed with Euclid alone is 18%, which is reduced to 4(13)% by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimized to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (I-Y)>2.8 and (Y-J)<1.4 colour criteria can separate contaminants from true z>6 galaxies, although these are applicable to only 54% of the contaminants, as many have unconstrained (I-Y) colours. In the most optimistic scenario, these cuts reduce the contamination fraction to 1% whilst preserving 81% of the fiducial z>6 sample. For the faint mock sample, colour cuts are infeasible., Comment: 27 pages, 14 figures, 9 tables; version accepted for publication by A&A

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $\mu$m range, to a 5$\sigma$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime., Comment: 33 pages, 25 figures, accepted for publication in A&A

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic S\'ersic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies., Comment: 22 pages, 21 figures (6 in appendinx). Accepted for publication to A&A. Change in the title for Euclid submission chronology coherence. Few changes for clarification

In the inflation-based cosmology the dark matter (DM) density component starts moving with respect to the universal expansion at $z_{\rm eq}\sim 3,200$ while baryons remain frozen until $z_{\rm rec}\sim 1,100$. It has been suggested that in this case post-linear corrections to the evolution of small fluctuations would result, for the standard $\Lambda$-dominated cold DM (CDM) model, in delayed formation of early objects as supersonic advection flows develop after recombination, so baryons are not immediately captured by the DM gravity on small scales. We develop the hydrodynamical description of such two-component advection and show that, in the supersonic regime, the advection within irrotational fluids is governed by the gradient of the difference of the kinetic energies of the two (DM and baryonic here) components. We then apply this formalism to the case where DM is made up of LIGO-type black holes (BHs) and show that there the advection process on scales relevant for early structure collapse will differ significantly from the earlier discussed (CDM) case because of the additional granulation component to the density field produced during inflation. The advection here will lead efficiently to the common motion of the DM and baryon components on scales relevant for collapse and formation of first luminous sources.This leads to early collapse making easier to explain the existence of supermassive BHs observed in quasars at high $z>7$. The resultant net advection rate reaches minimum around $< 10^9M_\odot$ and subsequently rises to a secondary maximum near the typical mass of $\sim 10^{12}M_\odot$, which may be an important consideration for formation of galaxies at z<(a few)., Comment: Physical Review Letters, to be published

We study the coherence of the near-infrared and X-ray background fluctuations and the X-ray spectral properties of the sources producing it. We use data from multiple Spitzer and Chandra surveys, including the UDS/SXDF surveys, the Hubble Deep Field North, the EGS/AEGIS field, the Chandra Deep Field South and the COSMOS surveys, comprising $\sim$2275 Spitzer/IRAC hours and $\sim$~16 Ms of Chandra data collected over a total area of $\sim$~1~deg$^2$. We report an overall $\sim$5$\sigma$ detection of a cross-power signal on large angular scales $>$ 20$''$ between the 3.6 and 4.5\mum\ and the X-ray bands, with the IR vs [1-2] keV signal detected at 5.2$\sigma$. The [0.5-1] and [2-4] keV bands are correlated with the infrared wavelengths at a $\sim$1$-$3$\sigma$ significance level. The hardest X-ray band ([4-7] keV) alone is not significantly correlated with any infrared wavelengths due to poor photon and sampling statistics. We study the X-ray SED of the cross-power signal. We find that its shape is consistent with a variety of source populations of accreting compact objects, such as local unabsorbed AGNs or high-z absorbed sources. We cannot exclude that the excess fluctuations are produced by more than one population. Because of poor statistics, the current relatively broad photometric bands employed here do not allow distinguishing the exact nature of these compact objects or if a fraction of the fluctuations have instead a local origin., Comment: 14 Pages, 4 figures, Accepted by ApJ, for inquires please contact author N. Cappelluti

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role., Comment: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584

A major scientific goal of JWST is to probe the epoch of re-ionization of the Universe at z above 6, and up to 20 and beyond. At these redshifts, galaxies are just beginning to form and the observable objects are early black holes, supernovae, and cosmic infrared background. The JWST has the necessary sensitivity to observe these targets individually, but a public deep and wide science enabling survey in the wavelength range from 2-5 $\mu$m is needed to discover these black holes and supernovae and to cover the area large enough for cosmic infrared background to be reliably studied. This enabling survey will also discover a large number of other transients and enable sciences such as supernova cosmology up to z $\sim$ 5, star formation history at high redshift through supernova explosions, faint stellar objects in the Milky Way, and galaxy evolution up to z approaching 10. The results of this survey will also serve as an invaluable target feeder for the upcoming era of ELT and SKA., Comment: arXiv admin note: text overlap with arXiv:1710.07005

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing., Comment: Science whitepaper submitted to the Astro2020 Decadal Survey

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes)., Comment: Science whitepaper submitted to the Astro2020 Decadal Survey

The source-subtracted cosmic infrared background (CIB) fluctuations uncovered in deep Spitzer data cannot be explained by known galaxy populations and appear strongly coherent with unresolved cosmic X-ray background (CXB). This suggests that the source-subtracted CIB contains emissions from significantly abundant accreting black holes (BHs). We show that theoretically such populations would have the angular power spectrum which is largely independent of the epochs occupied by these sources, provided they are at z>~ 4, offering an important test of the origin of the new populations. Using the current measurements we reconstruct the underlying soft X-ray CXB from the new sources and show that its fluctuations, while consistent with a high-z origin, have an amplitude that cannot be reached in direct measurements with the foreseeable X-ray space missions. This necessitates application of the methods developed by the authors to future IR and X-ray datasets, which must cover large areas of the sky in order to measure the signal with high precision. The LIBRAE project within ESA's Euclid mission will probe source-subtracted CIB over ~1/2 the sky at three near-IR bands, and its cross-power with unresolved CXB can be measured then from the concurrent eROSITA mission covering the same areas of the sky. We discuss the required methodology for this measurement and evaluate its projected S/N to show the unique potential of this experimental configuration to accurately probe the CXB from the new BH sources and help identify their epochs., Comment: Astrophysical Journal Letters, to be published

We study the source-subtracted near-infrared and X-ray background fluctuations of the COSMOS field using data from the Spitzer SPLASH program ($\sim$1272 hours) and Chandra COSMOS Legacy Survey (4.6 Ms). The new auto power spectra of the cosmic infrared and X-ray background fluctuations reach maximum angular scales of $\sim$ 3000$''$ and $\sim$ 5000$''$, respectively. We measure the cross power spectra between each infrared and X-ray band and calculate the mean power above 20$''$. We find that the soft X-ray band is correlated with 3.6 and 4.5$\mu$m at $\sim$ 4 $\sigma$ significance level. The significance between hard X-ray and the 3.6$\mu$m (4.5$\mu$m) band is $\sim$2.2 $\sigma$ ($\sim$3.8$\sigma$). The combined infrared (3.6 + 4.5$\mu$m) data are correlated with the X-ray data in soft ([0.5-2] keV), hard ([2-7] keV) and broad ([0.5-7] keV) bands at $\sim$5.6$\sigma$, $\sim$4.4$\sigma$ and $\sim$6.6$\sigma$ level, respectively. We compare the new measurements with existing models for the contributions from known populations at $z$$<$7, which are not subtracted. The model predictions are consistent with the measurements but we cannot rule out contributions from other components, such as Direct Collapse Black Holes (DCBH). However, the stacked cross-power spectra, combining other available data, show excess fluctuations about an order of magnitude on average at $\sim$4$\sigma$ confidence at scales within $\sim$300$''$. By studying the X-ray SED of the cross-power signal, assuming no significant variation from the infrared, we find that its shape is consistent with DCBHs., Comment: 20 pages, 20 figures, 4 tables

The cosmic infrared background (CIB) contains emissions accumulated over the entire history of the Universe, including from objects inaccessible to individual telescopic studies. The near-IR (~1-10 mic) part of the CIB, and its fluctuations, reflects emissions from nucleosynthetic sources and gravitationally accreting black holes (BHs). If known galaxies are removed to sufficient depths the source-subtracted CIB fluctuations at near-IR can reveal sources present in the first-stars-era and possibly new stellar populations at more recent times. This review discusses the recent progress in this newly emerging field which identified, with new data and methodology, significant source-subtracted CIB fluctuations substantially in excess of what can be produced by remaining known galaxies. The CIB fluctuations further appear coherent with unresolved cosmic X-ray background (CXB) indicating a very high fraction of BHs among the new sources producing the CIB fluctuations. These observations have led to intensive theoretical efforts to explain the measurements and their properties. While current experimental configurations have limitations in decisively probing these theories, their potentially remarkable implications will be tested in the upcoming CIB measurements with the ESA's Euclid dark energy mission. We describe the goals and methodologies of LIBRAE (Looking at Infrared Background Radiation with Euclid), a NASA-selected project for CIB science with Euclid, which has the potential for transforming the field into a new area of precision cosmology., Comment: Reviews of Modern Physics, to appear

We present new measurements of the large scale clustering component of the cross-power spectra of the source-subtracted Spitzer-IRAC Cosmic Infrared Background (CIB) and Chandra-ACIS Cosmic X-ray Background (CXB) surface brightness fluctuations. Our investigation uses data from the Chandra Deep Field South (CDFS), Hubble Deep Field North (HDFN), EGS/AEGIS field and UDS/SXDF surveys, comprising 1160 Spitzer hours and $\sim$ 12 Ms of Chandra data collected over a total area of 0.3 deg$^2$. We report the first ($>$5$\sigma$) detection of a cross-power signal on large angular scales >20 arcsec between [0.5-2] keV and the 3.6 and 4.5$\mu$m bands, at $\sim$5$\sigma$ and 6.3$\sigma$ significance, respectively. The correlation with harder X-ray bands is marginally significant. Comparing the new observations with existing models for the contribution of the known unmasked source population at $z<$7, we find an excess of about an order of magnitude at 5$\sigma$ confidence. We discuss possible interpretations for the origin of this excess in terms of the contribution from accreting early black holes, including both direct collapse black holes and primordial black holes, as well as from scattering in the interstellar medium and intra-halo light., Comment: 7 pages, 3 figures, accepted by ApJ Letters, v2 slight corrections

LIGO's discovery of a gravitational wave from two merging black holes (BHs) of similar masses rekindled suggestions that primordial BHs (PBHs) make up the dark matter (DM). If so, PBHs would add a Poissonian isocurvature density fluctuation component to the inflation-produced adiabatic density fluctuations. For LIGO's BH parameters, this extra component would dominate the small-scale power responsible for collapse of early DM halos at z>10, where first luminous sources formed. We quantify the resultant increase in high-z abundances of collapsed halos that are suitable for producing the first generation of stars and luminous sources. The significantly increased abundance of the early halos would naturally explain the observed source-subtracted near-IR cosmic infrared background (CIB) fluctuations, which cannot be accounted for by known galaxy populations. For LIGO's BH parameters this increase is such that the observed CIB fluctuation levels at 2 to 5 micron can be produced if only a tiny fraction of baryons in the collapsed DM halos forms luminous sources. Gas accretion onto these PBHs in collapsed halos, where first stars should also form, would straightforwardly account for the observed high coherence between the CIB and unresolved cosmic X-ray background in soft X-rays. We discuss modifications possibly required in the processes of first star formation if LIGO-type BHs indeed make up the bulk or all of DM. The arguments are valid only if the PBHs make up all, or at least most, of DM, but at the same time the mechanism appears inevitable if DM is made of PBHs., Comment: Ap.J. Letters, in press

We have performed a specific observational test to measure the effect that the zodiacal light can have on measurements of the spatial fluctuations of the near-IR background. Previous estimates of possible fluctuations caused by zodiacal light have often been extrapolated from observations of the thermal emission at longer wavelengths and low angular resolution, or from IRAC observations of high latitude fields where zodiacal light is faint and not strongly varying with time. The new observations analyzed here target the COSMOS field, at low ecliptic latitude where the zodiacal light intensity varies by factors of $\sim2$ over the range of solar elongations at which the field can be observed. We find that the white noise component of the spatial power spectrum of the background is correlated with the modeled zodiacal light intensity. Roughly half of the measured white noise is correlated with the zodiacal light, but a more detailed interpretation of the white noise is hampered by systematic uncertainties that are evident in the zodiacal light model. At large angular scales ($\gtrsim100"$) where excess power above the white noise is observed, we find no correlation of the power with the modeled intensity of the zodiacal light. This test clearly indicates that the large scale power in the infrared background is not being caused by the zodiacal light., Comment: 17 pp. Accepted for publication in the ApJ

The Euclid space mission, designed to probe evolution of the Dark Energy, will map a large area of the sky at three adjacent near-IR filters, Y, J and H. This coverage will also enable mapping source-subtracted cosmic infrared background (CIB) fluctuations with unprecedented accuracy on sub-degree angular scales. Here we propose methodology, using the Lyman-break tomography applied to the Euclid-based CIB maps, to accurately isolate the history of CIB emissions as a function of redshift from 10 < z < 20, and to identify the baryonic acoustic oscillations (BAOs) at those epochs. To identify the BAO signature, we would assemble individual CIB maps over conservatively large contiguous areas of >~ 400 sq deg. The method can isolate the CIB spatial spectrum by z to sub-percent statistical accuracy. We illustrate this with a specific model of CIB production at high z normalized to reproduce the measured Spitzer-based CIB fluctuation. We show that even if the latter contain only a small component from high-z sources, the amplitude of that component can be accurately isolated with the methodology proposed here and the BAO signatures at z>~ 10 are recovered well from the CIB fluctuation spatial spectrum. Probing the BAO at those redshifts will be an important test of the underlying cosmological paradigm, and would narrow the overall uncertainties on the evolution of cosmological parameters, including the Dark Energy. Similar methodology is applicable to the planned WFIRST mission, where we show that a possible fourth near-IR channel at > 2 micron would be beneficial., Comment: comments welcome

The study of the Cosmic Near-Infrared Background (CIB) light after subtraction of resolved sources can push the limits of current observations and infer the level of galaxy and black hole activity in the early universe. However, disentangling the relative contribution from low- and high-redshift sources is not trivial. Spatial fluctuations of the CIB exhibit a clustering excess at angular scales $\sim 1^\circ$ whose origin has not been conclusively identified. We explore the likelihood that this signal is dominated by emission from galaxies and accreting black holes in the early Universe. We find that, if the first small mass galaxies have a normal IMF, the light of their ageing stars (fossils) integrated over cosmic time contributes a comparable amount to the CIB as their pre-reionization progenitors. However, the measured fluctuation signal is too large to be produced by galaxies at redshifts $z>8$ unless their star formation efficiencies are much larger than those inferred from the observed Lyman-dropout population. In order to produce the observed level of CIB fluctuation without violating constraints from galaxy counts and the electron optical depth of the IGM, minihalos at $z>12$ must form stars with efficiency $f_\star \gtrsim 0.1$ and, although a top-heavy IMF is preferred, have a very low escape fraction of ionizing radiation, $f_{\rm esc}<0.01$. If instead the CIB fluctuations are produced by high-$z$ black holes, one requires vigorous accretion in the early universe reaching $\rho_{\rm acc} \gtrsim 10^5M_\odot{\rm Mpc^{-3}}$ by $z\simeq 10$. This growth must stop by $z \sim 6$ and be significantly obscured not to overproduce the soft cosmic X-ray background (CXB) and its observed coherence with the CIB. We therefore find the range of suitable possibilities at high-$z$ to be narrow, but could possibly be widened by including additional physics and evolution at those epochs., Comment: 13 pages, 11 figures. Accepted for publication in MNRAS

We constrain the deviation of adiabatic evolution of the Universe using the data on the Cosmic Microwave Background (CMB) temperature anisotropies measured by the {\it Planck} satellite and a sample of 481 X-ray selected clusters with spectroscopically measured redshifts. To avoid antenna beam effects, we bring all the maps to the same resolution. We use a CMB template to subtract the cosmological signal while preserving the Thermal Sunyaev-Zeldovich (TSZ) anisotropies; next, we remove galactic foreground emissions around each cluster and we mask out all known point sources. If the CMB black-body temperature scales with redshift as $T(z)=T_0(1+z)^{1-\alpha}$, we constrain deviations of adiabatic evolution to be $\alpha=-0.007\pm 0.013$, consistent with the temperature-redshift relation of the standard cosmological model. This result could suffer from a potential bias $\delta\alpha$ associated with the CMB template, that we quantify it to be $|\delta\alpha|\le 0.02$ and with the same sign than the measured value of $\alpha$, but is free from those biases associated with using TSZ selected clusters; it represents the best constraint to date of the temperature-redshift relation of the Big-Bang model using only CMB data, confirming previous results., Comment: ApJ, in press. Manuscript matches the accepted version: 10 pages, 7 figures, 3 tables

We present new methodology to use cosmic infrared background (CIB) fluctuations to probe sources at 1030., Comment: ApJ, in press. Minor revisions/additions to match the version in proofs

The "dark flow" dipole is a statistically significant dipole found at the position of galaxy clusters in filtered maps of Cosmic Microwave Background (CMB) temperature anisotropies. The dipole measured in WMAP 3, 5 and 7 yr data releases was 1) mutually consistent, 2) roughly aligned with the all-sky CMB dipole and 3) correlated with clusters' X-ray luminosity. We analyzed WMAP 9 yr and the 1st yr Planck data releases using a catalog of 980 clusters outside the Kp0 mask to test our earlier findings. The dipoles measured on these new data sets are fully compatible with our earlier estimates, being similar in amplitude and direction to our previous results and in disagreement with the results of an earlier study by the Planck Collaboration. Further, in Planck datasets dipoles are found independent of frequency, ruling out the Thermal Sunyaev-Zeldovich as the source of the effect. Both, in WMAP and Planck, we find a clear correlation between the dipole measured at the cluster location in filtered maps with the average anisotropy on the original maps, further proving that the dipole is associated with clusters. The dipole signal is dominated by the most massive clusters, with a statistical significance better than 99\%, slightly larger than in WMAP. Since both data sets differ in foreground contributions, instrumental noise and other systematics, the agreement between WMAP and Planck dipoles argues against them being due to systematic effects in either of the experiments., Comment: ApJ, in press

The epoch of first star formation and the state of the intergalactic medium (IGM) at that time are not directly observable with current telescopes. The radiation from those early sources is now part of the Cosmic Infrared Background (CIB) and, as these sources ionize the gas around them, the IGM plasma would produce faint temperature anisotropies in the Cosmic Microwave Background (CMB) via the thermal Sunyaev-Zeldovich (TSZ) effect. While these TSZ anisotropies are too faint to be detected, we show that the cross-correlation of maps of source-subtracted CIB fluctuations from {\it Euclid}, with suitably constructed microwave maps at different frequencies can probe the physical state of the gas during reionization and test/constrain models of the early CIB sources. We identify the frequency-combined CMB-subtracted microwave maps from space and ground-based instruments to show that they can be cross-correlated with the forthcoming all-sky {\it Euclid} CIB maps to detect the cross-power at scales $\sim 5'-60'$ with the signal/noise of up to $S/N\sim 4-8$ depending on the contribution to the Thomson optical depth during those pre-reionization epochs ($\Delta \tau\simeq 0.05$) and the temperature of IGM (up to $\sim10^4$K). Such a measurement would offer a new window to explore emergence and physical properties of these first light sources., Comment: ApJ Lett, in press

A spatial clustering signal has been established in Spitzer/IRAC measurements of the unresolved Cosmic near-Infrared Background (CIB) out to large angular scales, ~1 deg. This CIB signal, while significantly exceeding the contribution from the remaining known galaxies, was further found to be coherent at a highly statistically significant level with the unresolved soft Cosmic X-ray Background (CXB). This measurement probes the unresolved CXB to very faint source levels using deep near-IR source subtraction. We study contributions from extragalactic populations at low to intermediate redshifts to the measured positive cross-power signal of the CIB fluctuations with the CXB. We model the X-ray emission from AGN, normal galaxies and hot gas residing in virialized structures, calculating their CXB contribution including their spatial coherence with all infrared emitting counterparts. We use a halo model framework to calculate the auto and cross-power spectra of the unresolved fluctuations based on the latest constraints of the halo occupation distribution and the biasing of AGN, galaxies and diffuse emission. At small angular scales (<1'), the 4.5mic vs 0.5-2 keV coherence can be explained by shot noise from galaxies and AGN. However, at large angular scales (~10 arcmin) we find that the net contribution from the modeled populations is only able to account for ~3% of the measured CIBxCXB cross-power. The discrepancy suggests that the CIBxCXB signal originates from the same unknown source population producing the CIB clustering signal out to ~1 deg., Comment: 18 pages, 12 figures. Accepted for publication in ApJ

Recent measurements of large-scale peculiar velocities from the cumulative kinematic Sunyaev-Zel'dovich effect using WMAP data and X-ray selected clusters from ROSAT have identified a bulk flow of galaxy clusters at $\sim 600-1,000$ km s$^{-1}$ on scales of $\sim0.5-1$ Gpc, roughly aligned with the all-sky Cosmic Microwave Background (CMB) dipole. The flow is inferred from the detection of a residual dipole generated by CMB fluctuations exclusively in the direction of galaxy clusters, and measured within apertures containing zero monopole. Consistent with this interpretation, the amplitude of the dipole correlates with the X-ray luminosity of the clusters. To enable this measurement, the CMB data need to be filtered to remove the primary CMB, thereby increasing the data's signal-to-noise ratio. Filtering cannot imprint a signal with the mentioned properties at cluster positions; however, an inadequately designed filter can greatly suppress such a signal. We show here that recent studies that failed to detect a large-scale flow with various filters indeed adopted flawed implementations; when correctly implemented, these alternative filters lead to results that are in fact consistent with the Dark Flow signal. The discrepancies can be traced to the likely presence of residual dipoles caused by the thermal Sunyaev-Zel'dovich effect, to assumptions about cluster profiles incompatible with the data as well as failure to compute dipoles at the zero monopole aperture. PLANCK maps, with their large frequency coverage and a 217 GHz channel, will be instrumental to probe bulk flows, to remove spurious dipole signals and to help identify filtering schemes appropriate for this measurement., Comment: Revised arguments in section 3, results unchanged

In order to understand the nature of the sources producing the recently uncovered cosmic infrared background (CIB) fluctuations, we study cross-correlations between the fluctuations in the source-subtracted CIB from Spitzer/IRAC data and the unresolved cosmic X-ray background from deep Chandra observations. Our study uses data from the EGS/AEGIS field, where both data sets cover an 8'x45' region of the sky. Our measurement is the cross-power spectrum between the IR and X-ray data. The cross-power signal between the IRAC maps at 3.6 um and 4.5 um and the Chandra [0.5-2] keV data has been detected, at angular scales 20'', with an overall significance of 3.8 sigma and 5.6 sigma, respectively. At the same time we find no evidence of significant cross-correlations at the harder Chandra bands. The cross-correlation signal is produced by individual IR sources with 3.6 um and 4.5 um magnitudes m_AB 25-26 and [0.5-2] keV X-ray fluxes 7e-17 cgs.We determine that at least 15%-25% of the large scale power of the CIB fluctuations is correlated with the spatial power spectrum of the X-ray fluctuations. If this correlation is attributed to emission from accretion processes at both IR and X-ray wavelengths, this implies a much higher fraction of accreting black holes than among the known populations. We discuss the various possible origins for the cross-power signal and show that neither local foregrounds nor the known remaining normal galaxies and active galactic nuclei can reproduce the measurements. These observational results are an important new constraint on theoretical modeling of the near-IR CIB fluctuations., Comment: 28 pages, 9 Figures, Accepted by ApJ

Spatial Fluctuations in the Cosmic Infrared Background have now been measured out to sub-degree scales showing a strong clustering signal from unresolved sources. We attempt to explain these measurement by considering faint galaxy populations at z<6 as the underlying sources for this signal using 233 measured UV, optical and NIR luminosity functions (LF) from a variety of surveys covering a wide range of redshifts. We populate the lightcone and calculate the total emission redshifted into the near-IR bands in the observer frame and recover the observed optical and near-IR galaxy counts to a good accuracy. Using a halo model for the clustering of galaxies with an underlying LCDM density field, we find that fluctuations from known galaxy populations are unable to account for the large scale CIB clustering signal seen by HST/NICMOS, Spitzer/IRAC and AKARI/IRC and continue to diverge out to larger angular scales. Our purely empirical reconstruction shows that known galaxy populations are not responsible for the bulk of the fluctuation signal seen in the measurements and suggests an unknown population of very faint and highly clustered sources dominating the signal., Comment: 3 pages, 2 figures. To appear in proceedings of First Stars IV meeting (Kyoto, Japan; 2012)

We reconstruct \gamma-ray opacity of the Universe out to z<3-4 using an extensive library of 342 observed galaxy luminosity function surveys extending to high redshifts. We cover the whole range from UV to mid-IR (0.15-25mic) providing for the first time a robust empirical calculation of the \gamma\gamma-optical depth out to several TeV. Here, we use the same database as Helgason et al. 2012 where the EBL was reconstructed from luminosity functions out to 4.5mic and was shown to recover observed galaxy counts to high accuracy. We extend our earlier library of LFs to 25mic such that it covers the energy range of pair production with \gamma-rays 1) in the entire Fermi/LAT energy range, and 2) at higher TeV energies probed by ground-based Cherenkov telescopes. In the absence of significant contributions to the cosmic diffuse background from unknown populations, such as the putative Population III era sources, the Universe appears to be largely transparent to \gamma-rays at all Fermi/LAT energies out to z~2 whereas becoming opaque to TeV photons already at z<0.2 and reaching \tau~10 at z=1. Comparing with the currently available Fermi/LAT GRB and blazar data shows that there is room for significant emissions originating in the first stars era., Comment: 5 pages, 5 figures. Published in ApJ Letters

In this paper we investigate the power spectrum of the unresolved 0.5-2 keV CXB with deep Chandra 4 Ms observations in the CDFS. We measured a signal which, on scales >30", is significantly higher than the Shot-Noise and is increasing with the angular scale. We interpreted this signal as the joint contribution of clustered undetected sources like AGN, Galaxies and Inter-Galactic-Medium (IGM). The power of unresolved cosmic sources fluctuations accounts for \sim 12% of the 0.5-2 keV extragalactic CXB. Overall, our modeling predicts that \sim 20% of the unresolved CXB flux is made by low luminosity AGN, \sim 25% by galaxies and \sim 55% by the IGM (Inter Galactic Medium). We do not find any direct evidence of the so called Warm Hot Intergalactic Medium (i.e. matter with 10^5K7.5 and compared it with SMBH evolution models., Comment: 15 pages, 9 figures, accepted by MNRAS

We study the power of PLANCK data to constrain deviations of the Cosmic Microwave Background black body temperature from adiabatic evolution using the thermal Sunyaev-Zeldovich anisotropy induced by clusters of galaxies. We consider two types of data sets: the cosmological signal is removed in the Time Ordered Information or is removed from the final maps; and two different statistical estimators, based on the ratio of temperature anisotropies at two different frequencies and on a fit to the spectral variation of the cluster signal with frequency. To test for systematics, we construct a template from clusters drawn from a hydro-simulation included in the pre-launch Planck Sky Model. We demonstrate that, using a proprietary catalog of X-ray selected clusters with measured redshifts, electron densities and X-ray temperatures, we can constrain deviations of adiabatic evolution, measured by the parameter $\alpha$ in the redshift scaling $T(z)=T_0(1+z)^{1-\alpha}$, with an accuracy of $\sigma_\alpha=0.011$ in the most optimal case and with $\sigma_\alpha=0.016$ for a less optimal case. These results represent a factor 2-3 improvement over similar measurements carried out using quasar spectral lines and a factor 6-20 with respect to earlier results using smaller cluster samples., Comment: 12 pages in ApJ style

In this review we present a comprehensive discussion of peculiar velocity field measured recently on very large scales with a novel method using X-ray galaxy clusters as tracers. The measurement is based on the kinematic component of the Sunyaev-Zeldovich (KSZ) effect produced by Compton scattering of cosmic microwave background (CMB) photons off the hot intracluster gas, and uses a large catalog of X-ray selected clusters and all-sky CMB maps obtained with the WMAP satellite. The method probes the dipole of the CMB temperature field evaluated at the cluster positions and within the apertures in which the CMB monopole contribution vanishes, thereby isolating the signal remaining from the KSZ effect produced by coherently moving clusters. The detection of a highly significant dipole out to the depth of at least ~ 800 Mpc casts doubt on the notion that gravitational instability from the observed mass distribution is the sole -- or even dominant -- cause of the detected motions. Rather it appears that the flow may extend across the entire observable Universe. Possible implications include the possibility to constrain the primeval preinflationary structure of space-time and its landscape, and/or the need for modifications of presently known physics (e.g. arising from a higher-dimensional structure of gravity). We review these possibilities in light of the measurements described here and specifically discuss the prospects of future measurements and the issues they should resolve. We address the consistency of these large-scale velocity measurements with those obtained on smaller scales by studies using galaxies as tracers, and resolve the discrepancies with two recent claims based on modified CMB analysis schemes., Comment: Physics Reports, draft version - comments welcome! Manuscript with all figures at full resolution is available from http://www.kashlinsky.info/bulkflows/physicsreport

We extend previous measurements of cosmic infrared background (CIB) fluctuations to ~ 1 deg using new data from the Spitzer Extended Deep Survey. Two fields, with depths of ~12 hr/pixel over 3 epochs, are analyzed at 3.6 and 4.5 mic. Maps of the fields were assembled using a self-calibration method uniquely suitable for probing faint diffuse backgrounds. Resolved sources were removed from the maps to a magnitude limit of AB mag ~ 25, as indicated by the level of the remaining shot noise. The maps were then Fourier-transformed and their power spectra were evaluated. Instrumental noise was estimated from the time-differenced data, and subtracting this isolates the spatial fluctuations of the actual sky. The power spectra of the source-subtracted fields remain identical (within the observational uncertainties) for the three epochs indicating that zodiacal light contributes negligibly to the fluctuations. Comparing to 8 mic power spectra shows that Galactic cirrus cannot account for the fluctuations. The signal appears isotropically distributed on the sky as required for an extragalactic origin. The CIB fluctuations continue to diverge to > 10 times those of known galaxy populations on angular scales out to < 1 deg. The low shot noise levels remaining in the diffuse maps indicate that the large scale fluctuations arise from the spatial clustering of faint sources well below the confusion noise. The spatial spectrum of these fluctuations is in reasonable agreement with an origin in populations clustered according to the standard cosmological model (LCDM) at epochs coinciding with the first stars era., Comment: ApJ, to be published

We model fluctuations in the Cosmic Infrared Background (CIB) arising from known galaxy populations using 233 measured UV, optical and NIR luminosity functions (LF) from a variety of surveys spanning a wide range of redshifts. We compare best-fit Schechter parameters across the literature and find clear indication of evolution with redshift. Providing fitting formulae for the multi-band evolution of the LFs out to z~5, we calculate the total emission redshifted into the near-IR bands in the observer frame and recover the observed optical and near-IR galaxy counts to a good accuracy. Our empirical approach, in conjunction with a halo model describing the clustering of galaxies, allows us to compute the fluctuations of the unresolved CIB and compare the models to current measurements. We find that fluctuations from known galaxy populations are unable to account for the large scale CIB clustering signal seen by Spitzer/IRAC and AKARI/IRC and continue to diverge out to larger angular scales. This holds true even if the LFs are extrapolated out to faint magnitudes with a steep faint-end slope all the way to z=8. We also show that removing resolved sources to progressively fainter magnitude limits, isolates CIB fluctuations to increasingly higher redshifts. Our empirical approach suggests that known galaxy populations are not responsible for the bulk of the fluctuation signal seen in the measurements and favors a very faint population of highly clustered sources., Comment: 20 pages, 12 figures. Typos in V2 corrected. Accepted for publication in ApJ

Conventional interpretation of the observed cosmic microwave background (CMB) dipole is that all of it is produced by local peculiar motions. Alternative explanations requiring part of the dipole to be primordial have received support from measurements of large-scale bulk flows. A test of the two hypothesis is whether other cosmic dipoles produced by collapsed structures later than last scattering coincide with the CMB dipole. One background is the cosmic infrared background (CIB) whose absolute spectrum was measured to ~30% by the COBE satellite. Over the 100 to 500 {\mu}m wavelength range its spectral energy distribution can provide a probe of its alignment with CMB. This is tested with the COBE FIRAS dataset which is available for such a measurement because of its low noise and frequency resolution important for Galaxy subtraction. Although the FIRAS instrument noise is in principle low enough to determine the CIB dipole, the Galactic foreground is sufficiently close spectrally to keep the CIB dipole hidden. A similar analysis is performed with DIRBE, which - because of the limited frequency coverage - provides a poorer a dataset. We discuss strategies for measuring the CIB dipole with future instruments to probe the tilt and apply it to the Planck, Herschel and the proposed Pixie missions. We demonstrate that a future FIRAS-like instrument with instrument noise a factor of ~10 lower than FIRAS would make a statistically significant measurement of the CIB dipole. We find that the Planck and Herschel data sets will not allow a robust CIB dipole measurement. The Pixie instrument promises a determination of the CIB dipole and its alignment with either the CMB dipole or the dipole galaxy acceleration vector., Comment: 9 pages 9 figures

We present new results on the "dark flow" from a measurement of the dipole in the distribution of peculiar velocities of galaxy clusters, applying the methodology proposed and developed by us earlier. Our latest measurement is conducted using new, low-noise 7-yr WMAP data as well as an all-sky sample of X-ray selected galaxy clusters compiled exclusively from published catalogs. Our analysis of the CMB signature of the kinematic Sunyaev-Zeldovich (SZ) effect finds a statistically significant dipole at the location of galaxy clusters. The residual dipole outside the cluster regions is small, rendering our overall measurement 3-4 sigma significant. The amplitude of the dipole correlates with cluster properties, being larger for the most X-ray luminous clusters, as required if the signal is produced by the SZ effect. Since it is measured at zero monopole, the dipole can not be due to the thermal SZ effect. Our results are consistent with those obtained earlier by us from 5-yr WMAP data and using a proprietary cluster catalog. In addition, they are robust to quadrupole removal, demonstrating that quadrupole leakage contributes negligibly to the signal. The lower noise of the 7-yr WMAP also allows us, for the first time, to obtain tentative empirical confirmation of our earlier conjecture that the adopted filtering flips the sign of the KSZ effect for realistic clusters and thus of the deduced direction of the flow. The latter is consistent with our earlier measurement in both the amplitude and direction. Assuming the filtering indeed flips the sign of the KSZ effect from the clusters, the direction agrees well also with the results of independent work using galaxies as tracers at lower distances. We make all maps and cluster templates derived by us from public data available to the scientific community to allow independent tests of our method and findings., Comment: ApJ, in press. Replaced with accepted version. The data needed for these results are at http://www.kashlinsky.info/bulkflows/data_public

We analyze the uncertainties and possible systematics associated with the "Dark Flow" measurements using the cumulative Sunyaev-Zeldovich effect combined with all-sky catalogs of clusters of galaxies. Filtering of all-sky cosmic microwave background (CMB) maps is required to remove the intrinsic cosmological signal down to the limit imposed by cosmic variance. Contributions to the errors come from the remaining cosmological signal, that integrates down with the number of clusters, and the instrumental noise, that scales with the number of pixels; the latter decreases with integration time and is subdominant for the Wilkinson Microwave Anisotropy Probe 5-year data. It is proven both analytically and numerically that the errors for 5-year WMAP data are $\simeq 15\sqrt{3/N_{\rm clusters}} \mu$K per dipole component. The measured bulk flow velocity is significant at up to the $>99.95%$ confidence level, depending on the cluster sample. We discuss different methods to compute error bars and demonstrate that they have biases that would overpredict the errors, as is the case in a recent reanalysis of our earlier results. If the signal is caused by systematic effects present in the data, such systematics must have a dipole pattern, which correlates with cluster X-ray luminosity and is present only at cluster positions. Only contributions from the Sunyaev-Zeldovich effect could provide such contaminants via several potential effects. We discuss such candidates apart from the bulk-motion of the cluster samples and demonstrate that their contributions to our measurements are negligible. (abridged), Comment: ApJ, in press

We present new measurements of the large-scale bulk flows of galaxy clusters based on 5-year WMAP data and a significantly expanded X-ray cluster catalogue. Our method probes the flow via measurements of the kinematic Sunyaev-Zeldovich (SZ) effect produced by the hot gas in moving clusters. It computes the dipole in the cosmic microwave background (CMB) data at cluster pixels, which preserves the SZ component while integrating down other contributions. Our improved catalog of over 1,000 clusters enables us to further investigate possible systematic effects and, thanks to a higher median cluster redshift, allows us to measure the bulk flow to larger scales. We present a corrected error treatment and demonstrate that the more X-ray luminous clusters, while fewer in number, have much larger optical depth, resulting in a higher dipole and thus a more accurate flow measurement. This results in the observed correlation of the dipole derived at the aperture of zero monopole with the monopole measured over the cluster central regions. This correlation is expected if the dipole is produced by the SZ effect and cannot be caused by unidentified systematics (or primary cosmic microwave background anisotropies). We measure that the flow is consistent with approximately constant velocity out to at least 800 Mpc. The significance of the measured signal peaks around 500 Mpc, most likely because the contribution from more distant clusters becomes progressively more diluted by the WMAP beam. We can, however, at present not rule out that these more distant clusters simply contribute less to the overall motion., Comment: Ap.J. (Letters), Mar 20/2010. Replaced with published version

This paper provides a detailed description of the data reduction and analysis procedures that have been employed in our previous studies of spatial fluctuation of the cosmic infrared background (CIB) using deep Spitzer IRAC observations. The self-calibration we apply removes a strong instrumental signal from the fluctuations which would otherwise corrupt our results. The procedures and results for masking bright sources, and modeling faint sources down to levels set by the instrumental noise are presented. Various tests are performed to demonstrate that the resulting power spectra of these fields are not dominated by instrumental or procedural effects. These tests indicate that the large scale (>~30') fluctuations that remain in the deepest fields are not directly related to the galaxies that are bright enough to be individually detected. We provide the parameterization of these power spectra in terms of separate instrument noise, shot noise, and power law components. Our measurements of spatial fluctuations of the CIB intensity indicate the mean emission from the objects producing the fluctuations is quite low (>~1 nW m-2 sr-1 at 3-5 micron), and thus consistent with current gamma-ray absorption constraints. The source of the fluctuations may be high-z Population III objects, or a more local component of very low luminosity objects with clustering properties that differ from the resolved galaxies. Finally, we discuss the prospects of the upcoming space-based surveys to directly measure the epochs inhabited by the populations producing these source-subtracted CIB fluctuations, and to isolate the individual fluxes of these populations., Comment: 76 pages, 39 Postscript figures. Submitted to ApJS. (Abstract abridged.); 2010 ApJS, in press (Jan)

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude., Comment: 7 pages; Science White Paper for the US Astro 2010-2020 Decadal Survey. If interested in further community-wide efforts on this topic please contact the first author

Peculiar velocities of clusters of galaxies can be measured by studying the fluctuations in the cosmic microwave background (CMB) generated by the scattering of the microwave photons by the hot X-ray emitting gas inside clusters. While for individual clusters such measurements result in large errors, a large statistical sample of clusters allows one to study cumulative quantities dominated by the overall bulk flow of the sample with the statistical errors integrating down. We present results from such a measurement using the largest all-sky X-ray cluster catalog combined to date and the 3-year WMAP CMB data. We find a strong and coherent bulk flow on scales out to at least > 300 h^{-1} Mpc, the limit of our catalog. This flow is difficult to explain by gravitational evolution within the framework of the concordance LCDM model and may be indicative of the tilt exerted across the entire current horizon by far-away pre-inflationary inhomogeneities., Comment: Ap.J. (Letters), in press. 20 Oct issue (Vol. 686)

This paper presents detailed analysis of large-scale peculiar motions derived from a sample of ~ 700 X-ray clusters and cosmic microwave background (CMB) data obtained with WMAP. We use the kinematic Sunyaev-Zeldovich (KSZ) effect combining it into a cumulative statistic which preserves the bulk motion component with the noise integrated down. Such statistic is the dipole of CMB temperature fluctuations evaluated over the pixels of the cluster catalog (Kashlinsky & Atrio-Barandela 2000). To remove the cosmological CMB fluctuations the maps are Wiener-filtered in each of the eight WMAP channels (Q, V, W) which have negligible foreground component. Our findings are as follows: The thermal SZ (TSZ) component of the clusters is described well by the Navarro-Frenk-White profile expected if the hot gas traces the dark matter in the cluster potential wells. Such gas has X-ray temperature decreasing rapidly towards the cluster outskirts, which we demonstrate results in the decrease of the TSZ component as the aperture is increased to encompass the cluster outskirts. We then detect a statistically significant dipole in the CMB pixels at cluster positions. Arising exclusively at the cluster pixels this dipole cannot originate from the foreground or instrument noise emissions and must be produced by the CMB photons which interacted with the hot intracluster gas via the SZ effect. The dipole remains as the monopole component, due to the TSZ effect, vanishes within the small statistical noise out to the maximal aperture where we still detect the TSZ component. We demonstrate with simulations that the mask and cross-talk effects are small for our catalog and contribute negligibly to the measurements. The measured dipole thus arises from the KSZ effect produced by the coherent large scale bulk flow motion., Comment: Minor changes to match the published version - Ap.J., 1 Feb 2009 issue