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.

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.

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.

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.

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).

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We present new methodology to use cosmic infrared background (CIB) fluctuations to probe sources at 10<z<30 from a JWST/NIRCam configuration that will isolate known galaxies to 28 AB mag at 0.5--5 micron. At present significant mutually consistent source-subtracted CIB fluctuations have been identified in the Spitzer and Akari data at 2--5 micron, but we demonstrate internal inconsistencies at shorter wavelengths in the recent CIBER data. We evaluate CIB contributions from remaining galaxies and show that the bulk of the high-z sources will be in the confusion noise of the NIRCam beam, requiring CIB studies. The accurate measurement of the angular spectrum of the fluctuations and probing the dependence of its clustering component on the remaining shot noise power would discriminate between the various currently proposed models for their origin and probe the flux distribution of its sources. We show that the contribution to CIB fluctuations from remaining galaxies is large at visible wavelengths for the current instruments precluding probing the putative Lyman-break of the CIB fluctuations. We demonstrate that with the proposed JWST configuration such measurements will enable probing the Lyman break. We develop a Lyman-break tomography method to use the NIRCam wavelength coverage to identify or constrain, via the adjacent two-band subtraction, the history of emissions over 10<z<30 as the Universe comes out of the 'Dark Ages'. We apply the proposed tomography to the current Spitzer/IRAC measurements at 3.6 and 4.5 micron, to find that it already leads to interestingly low upper limit on emissions at z>30.

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.

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.

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.

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.

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^5K<T<10^7K and density contrast {\delta} <1000), but we estimated that it could produce about 1/7 of the unresolved CXB. We placed an upper limit to the space density of postulated X-ray-emitting early black hole at z>7.5 and compared it with SMBH evolution models.

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.

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.

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.

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.

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.

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.

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.

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.

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)

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.

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.

(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.

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.

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.

Using WMAP 3-year data at the locations of close to $\sim 700$ X-ray selected clusters we have detected the amplitude of the thermal Sunyaev-Zeldovich (TSZ) effect at the 15$\sigma$ level, the highest statistical significance reported so far. Owing to the large size of our cluster sample, we are able to detect the corresponding CMB distortions out to large cluster-centric radii. The region over which the TSZ signal is detected is, on average, four times larger in radius than the X-ray emitting region, extending to $\sim 3h_{70}^{-1}$Mpc. We show that an isothermal $\beta$ model does not fit the electron pressure at large radii; instead, the baryon profile is consistent with the Navarro-Frenk-White profile, expected for dark matter in the concordance $\Lambda$CDM model. The X-ray temperature at the virial radius of the clusters falls by a factor $\sim 3-4$ from the central value, depending on the cluster concentration parameter. Our results suggest that cluster dynamics at large radii is dominated by dark matter and is well described by Newtonian gravity.

We present the latest results on CIB fluctuations from early epochs from deep Spitzer data. The results show the existence of significant CIB fluctuations at the IRAC wavelengths (3.6 to 8 mic) which remain after removing galaxies down to very faint levels. These fluctuations must arise from populations with a significant clustering component, but only low levels of the shot noise. There are no correlations between the source-subtracted IRAC maps and the corresponding fields observed with the HST ACS at optical wavelengths. Taken together, these data imply that 1) the sources producing the CIB fluctuations are individually faint with S_\nu< a few nJy at 3.6 and 4.5 mic; 2) have different clustering pattern than the more recent galaxy populations; 3) are located within the first 0.7 Gyr (unless these fluctuations can somehow be produced by - so far unobserved - local galaxies of extremely low luminosity and with the unusual for local populations clustering pattern), 4) produce contribution to the net CIB flux of at least 1-2 nW/m^2/sr at 3.6 and 4.5 mic and must have mass-to-light ratio significantly below the present-day populations, and 5) they have angular density of ~ a few per arcsec^2 and are in the confusion of the present day instruments, but can be individually observable with JWST.

Cosmic infrared background (CIB) contains emission from epochs inaccessible to current telescopic studies, such as the era of the first stars. We discuss theoretical expectations for the CIB contributions from the early population of massive stars. We then present the latest results from the ongoing project by our team (Kashlinsky, Arendt, Mather & Moseley 2005,2007a,b,c,) to measure CIB fluctuations from early epochs using deep Spitzer data. The results show the existence of significant CIB fluctuations at the IRAC wavelengths (3.6 to 8 mic) which remain after removing galaxies down to very faint levels. These fluctuations must arise from populations that have a significant clustering component, but only low levels of the shot noise. Furthermore, there are no correlations between the source-subtracted IRAC maps and the corresponding fields observed with the HST ACS at optical wavelengths. Taken together, these data imply that 1) the sources producing the CIB fluctuations are individually faint with flux < a few nJy at 3.6 and 4.5 mic; 2) are located within the first 0.7 Gyr (unless these fluctuations can somehow be produced by - so far unobserved - local galaxies of extremely low luminosity and with the unusual for local populations clustering pattern), 3) they produce contribution to the net CIB flux of at least 1-2 nW/m^2/sr at 3.6 and 4.5 mic and must have mass-to-light ratio significantly below the present-day populations, and 4) they have angular density of ~ a few per arcsec^2 and are in the confusion of the present day instruments, but can be individually observable with JWST.

We study the possible contribution of optical galaxies detected with the {\it Hubble} ACS instrument to the near-IR cosmic infrared (CIB) fluctuations in deep {\it Spitzer} images. The {\it Spitzer} data used in this analysis are obtained in the course of the GOODS project from which we select four independent $10^\prime\times10^\prime$ regions observed at both 3.6 and 4.5 \um. ACS source catalogs for all of these areas are used to construct maps containing only their emissions in the ACS $B, V, i, z$-bands. We find that deep Spitzer data exhibit CIB fluctuations remaining after removal of foreground galaxies of a very different clustering pattern at both 3.6 and 4.5 \um than the ACS galaxies could contribute. We also find that there are very good correlations between the ACS galaxies and the {\it removed} galaxies in the Spitzer maps, but practically no correlations remain with the residual Spitzer maps used to identify the CIB fluctuations. These contributions become negligible on larger scales used to probe the CIB fluctuations arising from clustering. This means that the ACS galaxies cannot contribute to the large-scale CIB fluctuations found in the residual Spitzer data. The absence of their contributions also means that the CIB fluctuations arise at $z\gsim 7.5$ as the Lyman break of their sources must be redshifted past the longest ACS band, or the fluctuations have to originate in the more local but extremely low luminosity galaxies.

Cosmic infrared background (CIB) includes emissions from objects inaccessible to current telescopic studies, such as the putative Population III, the first stars. Recently, strong direct evidence for significant CIB levels produced by the first stars came from CIB fluctuations discovered in deep Spitzer images. Such CIB levels should have left a unique absorption feature in the spectra of high-z GRBs and blazars as suggested in [4]. This is observable with GLAST sources at z>2 and measuring this absorption will give important information on energetics and constituents of the first stars era.

We point out that in their analysis of the deep Spitzer images, Cooray et al (2006) perform Fourier transform on maps which have very few pixels left (only 20 to 30 percent). For such deeply cut maps one cannot reliably compute large-scale map properties using Fourier transforms. Instead the maps must be analyzed via the correlation function, C(\theta), which is immune to mask effects. We find, when computing C(\theta) for their maps, that removing ACS/HST galaxies does not lead to appreciable change in the correlation properties of the remaining diffuse emission. We then demonstrate with simulations that the power spectrum of CIB fluctuations prior to removal of the ACS galaxies reproduces C(\theta) in the maps from which the ACS galaxies have been removed. This implies that these galaxies cannot be responsible for the CIB fluctuations detected in Kashlinsky et al (2005, 2007), contrary to the claims of Cooray et al (2006).

Cosmic infrared background fluctuations may contain measurable contribution from objects inaccessible to current telescopic studies, such as the first stars and other luminous objects in the first Gyr of the Universe's evolution. In an attempt to uncover this contribution we have analyzed the GOODS data obtained with the Spitzer IRAC instrument, which are deeper and cover larger scales than the Spitzer data we have previously analyzed. Here we report these new measurements of the cosmic infrared background (CIB) fluctuations remaining after removing cosmic sources to fainter levels than before. The remaining anisotropies on scales > 0.5 arcmin have a significant clustering component with a low shot-noise contribution. We show that these fluctuations cannot be accounted for by instrumental effects, nor by the Solar system and Galactic foreground emissions and must arise from extragalactic sources.

We discuss interpretation of the cosmic infrared background (CIB) anisotropies detected by us recently in the Spitzer IRAC based measurements. The fluctuations are approximately isotropic on the sky consistent with their cosmological origin. They remain after removal of fairly faint intervening sources and must arise from a population which has a strong CIB clustering component with only a small shot-noise level. We discuss the constraints the data place on the luminosities, epochs and mass-to-light ratios of the indvidual sources producing them. Assuming the concordance LambdaCDM cosmology the measurements imply that the luminous sources producing them lie at cosmic times < 1 Gyr and were individually much brighter per unit mass than the present stellar populations.

We review the recent measurements on the cosmic infrared background (CIB) and their implications for the physics of the first stars era, including Population III. The recently obtained CIB results range from the direct measurements of CIB fluctuations from distant sources using deep Spitzer data to strong upper limits on the near-IR CIB from blazar spectra. This allows to compare the Population III models with the CIB data to gain direct insight into the era of the first stars and the formation and evolution of Population III and the microphysics of the feedback processes in the first halos of collapsing material. We also discuss the cosmological confusion resulting from these CIB sources and the prospects for resolving them individually with NASA's upcoming space instruments such as the JWST.

We discuss the contribution of Population III stars to the near-IR (NIR) cosmic infrared background (CIB) and its effect on spectra of high-$z$ high-energy gamma-ray bursts (GRBs). It is shown that if Population III were massive stars, the claimed NIR CIB excess will be reproduced if only ~ 4+/-2% of all baryons went through these stars. Regardless of the precise value of the NIR CIB produced by them, they would leave enough photons to provide a large optical depth for high-energy photons from high-z GRBs. Observations of such GRBs are expected following the planned launch of NASA's GLAST mission. The presence or absence of such damping in the spectra of high-$z$ GRBs will then provide important information on the emissions from the Population III. The location of this cutoff may also serve as an indicator of the GRB's redshift.

Cosmic infrared background (CIB) contains information about galaxy luminosities over the entire history of the Universe and can be a powerful diagnostic of the early populations otherwise inaccessible to telescopic studies. Its measurements are very difficult because of the strong IR foregrounds from the Solar system and the Galaxy. Nevertheless, substantial recent progress in measuring the CIB and its structure has been made. The measurements now allow to set significant constraints on early galaxy evolution and, perhaps, even detect the elusive Population III era. We discuss briefly the theory behind the CIB, review the latest measurements of the CIB and its structure, and discuss their implications for detecting and/or constraining the first stars and their epochs.

The deepest space and ground-based observations find metal-enriched galaxies at cosmic times when the Universe was <1 Gyr old. These stellar populations had to be preceded by the metal-free first stars, Population III. Recent cosmic microwave background polarization measurements indicate that stars started forming early when the Universe was $\lsim200$ Myr old. Theoretically it is now thought that Population III stars were significantly more massive than the present metal-rich stellar populations. While such sources will not be individually detectable by existing or planned telescopes, they would have produced significant cosmic infrared background radiation in the near- infrared, whose fluctuations reflect the conditions in the primordial density field. Here we report a measurement of diffuse flux fluctuations after removing foreground stars and galaxies. The anisotropies exceed the instrument noise and the more local foregrounds and can be attributed to emission from massive Population III stars, providing observational evidence of an era dominated by these objects.

The Cosmic Infrared Background (CIB) reflects the sum total of galactic luminosities integrated over the entire age of the universe. From its measurements one can deduce information about objects at epochs currently inaccessible to telescopic studies. This review discusses the state of the current CIB measurements and the (mostly space-based) instruments with which these measurements have been made, the obstacles (the various foreground emissions)and the physics behind the CIB and its structure. Theoretical discussion of the CIB levels can now be normalized to the standard cosmological model narrowing down theoretical uncertainties. We review the information behind and theoretical modeling of both the mean (isotropic) levels of the CIB and their fluctuations. The CIB is divided into three broad bands: near-IR, mid-IR and far-IR. For each of the bands we review the main contributors to the CIB flux and the epochs at which the bulk of the flux originates. We also discuss the data on the various quantities relevant for correct interpretation of the CIB levels: the star-formation history, the present-day luminosity function measurements, resolving the various galaxy contributors to the CIB, etc. The integrated light of all galaxies in the deepest near-IR galaxy counts to date fails to match the observed mean level of the CIB, probably indicating a significant high-redshift contribution to the CIB. Additionally, Population III stars should have left a strong and measurable signature via their contribution to the cosmic infrared background (CIB) anisotropies for a wide range of their formation scenarios, and measuring the excess CIB anisotropies coming from high z would provide direct information on that epoch.

Following the successful mapping of the last scattering surface by WMAP and balloon experiments, the epoch of the first stars, when Population III stars formed, is emerging as the next cosmological frontier. It is not clear what these stars' properties were, when they formed or how long their era lasted before leading to the stars and galaxies we see today. We show that these questions can be answered with the current and future measurements of the near-IR cosmic infrared background (CIB). Theoretical arguments suggest that Population III stars were very massive and short-lived stars that formed at $z\sim 10-20$ at rare peaks of the density field in the cold-dark-matter Universe. Because Population III stars probably formed individually in small mini-halos, they are not directly accessible to current telescopic studies. We show that these stars left a strong and measurable signature via their contribution to the CIB anisotropies for a wide range of their formation scenarios. The excess in the recently measured near-IR CIB anisotropies over that from normal galaxies can be explained by contribution from early Population III stars. These results imply that Population III were indeed very massive stars and their epoch started at $z\sim 20$ and lasted past $z\lsim 13$. We show the importance of accurately measuring the CIB anisotropies produced by Population III with future space-based missions.

We introduce a new statistic to measure more accurately the cosmic sound speed of clusters of galaxies at different redshifts. This statistic is evaluated by cross-correlating cosmic microwave background (CMB) fluctuations caused by the Sunyaev-Zel'dovich effect from observed clusters of galaxies with their redshifts. When clusters are distributed in redshift bins of narrow width, one could measure the mean squared cluster peculiar velocity with an error \sigma_{C_S^2}\lsim (300{\rm km/s})^2. This can be done around z>0.3 with clusters of flux above 200 mJy which will be detected by PLANCK, coupled with high resolution microwave images to eliminate the cosmological part of the CMB fluctuations. The latter can be achieved with observations by the planned ALMA array or the NSF South Pole telescope and other surveys. By measuring the cosmic sound speed and the bulk flow in, e.g., 4 spheres of ~ 100h^{-1}Mpc at z=0.3, we could have a direct measurement of the matter density 0.21<\Omega_m<0.47 at 95 % confidence level.

We implement and further refine the recently proposed method (Kashlinsky, Hern\'andez-Monteagudo & Atrio-Barandela, 2001 - KHA) for a time efficient extraction of the power spectrum from future cosmic microwave background (CMB) maps. The method is based on the clustering properties of peaks and troughs of the Gaussian CMB sky. The procedure takes only ${1/2}[f(\nu)]^2N^2$ steps where $f(\nu)$ is the fraction of pixels with $|\delta T|\geq\nu$ standard deviations in the map of $N$ pixels. We use the new statistic introduced in KHA, $\xi_\nu$, which characterizes spatial clustering of the CMB sky peaks of progressively increasing thresholds. The tiny fraction of the remaining pixels (peaks and troughs) contains the required information on the CMB power spectrum of the entire map. The threshold $\nu$ is the only parameter that determines the accuracy of the final spectrum. We performed detailed numerical simulations for parameters of the two-year WMAP and Planck CMB sky data including cosmological signal, inhomogeneous noise and foreground residuals. In all cases we find that the method can recover the power spectrum out to the Nyquist scale of the experiment channel. We discuss how the error bars scale with $\nu$ allowing to decide between accuracy and speed. The method can determine with significant accuracy the CMB power spectrum from the upcoming CMB maps in only $\sim(10^{-5}-10^{-3})\times N^2$ operations.

We report first results for the cosmic infrared background (CIB) fluctuations at 1.25, 1.65 and 2.17 micron obtained from long exposures constructed from 2MASS standard star fields. We have co-added and analyzed scans from one such field with a total exposure time > 1 hour, and removed sources and other artifacts. The stars and galaxies were clipped out to K_s~19^m leaving only high-z galaxies (or possibly local low-surface-brightness systems). The residual component of the diffuse emission on scales from a few arc-sec to a few arc-min has a power-law slope consistent with emission produced by clustered galaxies. The noise (and residual artifacts) contribution to the signal is small and the colors of the signal are very different from Galactic stars or air-glow. We therefore identify the signal as CIB fluctuations from the faint unresolved galaxies. We show that the present-day galaxies with no evolution would produce a significant deficit in the observed CIB fluctuations. Thus the dominant contribution to the observed signal must come from high z and may indicate high rates of star formation at those epochs.

This is one of two papers in which we report the detection of structure in the cosmic infrared background (CIB) between 1.25 - 2.2 micron through the use of data from the Two Micron Sky Survey (2MASS). This paper concentrates on data assembly, analysis and the estimate of the various foreground contributions; the companion paper (Kashlinsky, Odenwald, Mather, Skrutskie, Cutri 2002, hereafter KOMSC) presents the cosmological results for the CIB fluctuations and their implications. By using repeated observations of a specific calibration star field, we were able to achieve integration times in excess of 3900 seconds compared to the 7.8 seconds in the standard 2MASS data product. This yielded a point source detection limit (3 \sigma) of +18.5^m in K_s band. The resulting co-added images were processed to remove point sources to a limiting surface brightness of +20^m/arcsec$^2 or 40 nW/m^2/sr. The remaining maps contained over 90% of the pixels and were Fourier transformed to study the spatial structure of the diffuse background light. After removing resolved sources and other artifacts, we find that the power spectrum of the final images has a power-law distribution consistent with clustering by distant galaxies. We estimate here the contributions to this signal from Galactic foregrounds, atmospheric OH-glow, zodiacal light and instrument noise, all of which are small and of different slopes. Hence, this supports the KOMSC identification of the signal as coming from the CIB fluctuations produced by distant clustered galaxies.

To understand the cosmic history of element synthesis it will be important to obtain extinction-free measures of the heavy element contents of high-redshift objects and to chart two monumental events: the collapse of the first metal-free clouds to form stars, and the initial seeding of the universe with dust. The information needed to achieve these objectives is uniquely available in the far-infrared/submillimeter (FIR/SMM) spectral region. Following the Decadal Report and anticipating the development of the Single Aperture Far-IR (SAFIR) telescope and FIR/SMM interferometry, we estimate the measurement capabilities of a large-aperture, background-limited FIR/SMM observatory and an interferometer on a boom, and discuss how such instruments could be used to measure the element synthesis history of the universe.

The question "How did we get here and what will the future bring?" captures the human imagination and the attention of the National Academy of Science's Astronomy and Astrophysics Survey Commitee (AASC). Fulfillment of this "fundamental goal" requires astronomers to have sensitive, high angular and spectral resolution observations in the far-infrared/submillimeter (far-IR/sub-mm) spectral region. With half the luminosity of the universe and vital information about galaxy, star and planet formation, observations in this spectral region require capabilities similar to those currently available or planned at shorter wavelengths. In this paper we summarize the scientific motivation, some mission concepts and technology requirements for far-IR/sub-mm space interferometers that can be developed in the 2010-2020 timeframe.

We present a new method for time-efficient and accurate extraction of the power spectrum from future cosmic microwave background (CMB) maps based on properties of peaks and troughs of the Gaussian CMB sky. We construct a statistic describing their angular clustering - analogously to galaxies, the 2-point angular correlation function, $\xi_\nu(\theta)$. We show that for increasing peak threshold, $\nu$, the $\xi_\nu(\theta)$ is strongly amplified and becomes measurable for $\nu\geq$1 on angular scales $\leq 10^\circ$. Its amplitude at every scale depends uniquely on the CMB temperature correlation function, $C(\theta)$, and thus the measured $\xi_\nu$ can be uniquely inverted to obtain $C(\theta)$ and its Legendre transform, the power spectrum of the CMB field. Because in this method the CMB power spectrum is deduced from high peaks/troughs of the CMB field, the procedure takes only $[f(\nu)]^2N^2$ operations where $f(\nu)$ is the fraction of pixels with $|\delta T|\geq\nu$ standard deviations in the map of $N$ pixels and is e.g. 0.045 and 0.01 for $\nu$=2 and 2.5 respectively. We develop theoretical formalism for the method and show with detailed simulations, using MAP mission parameters, that this method allows to determine very accurately the CMB power spectrum from the upcoming CMB maps in only $\sim(10^{-4}-10^{-3})\times N^2$ operations.

We propose a new method to measure the possible large-scale bulk flows in the Universe from the cosmic microwave background (CMB) maps from the upcoming missions, MAP and Planck. This can be done by studying the statistical properties of the CMB temperature field at many X-ray cluster positions. At each cluster position, the CMB temperature fluctuation will be a combination of the Sunyaev-Zeldovich (SZ) kinematic and thermal components, the cosmological fluctuations and the instrument noise term. When averaged over many such clusters the last three will integrate down, whereas the first one will be dominated by a possible bulk flow component. In particular, we propose to use all-sky X-ray cluster catalogs that should (or could) be available soon from X-ray satellites, and then to evaluate the dipole component of the CMB field at the cluster positions. We show that for the MAP and Planck mission parameters the dominant contributions to the dipole will be from the terms due to the SZ kinematic effect produced by the bulk flow (the signal we seek) and the instrument noise (the noise in our signal). Computing then the expected signal-to-noise ratio for such measurement, we get that at the 95 % confidence level the bulk flows on scales >100h^{-1} Mpc can be probed down to the amplitude of $< 200$ km/sec with the MAP data and down to only 30 km/sec with the Planck mission.

The cosmic infrared background (CIB) radiation is the cosmic repository for energy release throughout the history of the universe. Using the all-sky data from the COBE DIRBE instrument at wavelengths 1.25 - 100 mic we attempt to measure the CIB fluctuations. In the near-IR, foreground emission is dominated by small scale structure due to stars in the Galaxy. There we find a strong correlation between the amplitude of the fluctuations and Galactic latitude after removing bright foreground stars. Using data outside the Galactic plane ($|b| > 20\deg$) and away from the center ($90\deg< l <270\deg$) we extrapolate the amplitude of the fluctuations to cosec$|b|=0$. We find a positive intercept of $\delta F_{\rm rms} = 15.5^{+3.7}_{-7.0},5.9^{+1.6}_{-3.7}, 2.4^{+0.5}_{-0.9}, 2.0^{+0.25}_{-0.5}$ nW/m2/sr at 1.25, 2.2,3.5 and 4.9 mic respectively, where the errors are the range of 92% confidence limits. For color subtracted maps between band 1 and 2 we find the isotropic part of the fluctuations at $7.6^{+1.2}_{-2.4}$ nW/m2/sr. Based on detailed numerical and analytic models, this residual is not likely to originate from the Galaxy, our clipping algorithm, or instrumental noise. We demonstrate that the residuals from the fit used in the extrapolation are distributed isotropically and suggest that this extra variance may result from structure in the CIB. For $2\deg< \theta < 15^\deg$, a power-spectrum analysis yields firm upper limits of $(\theta/5^\deg) \times\delta F_{\rm rms} (\theta) < $ 6, 2.5, 0.8, 0.5 nW/m2/sr at 1.25, 2.2, 3.5 and 4.9 mic respectively. From 10-100 mic, the upper limits <1 nW/m2/sr.

Accurate synthetic models of stellar populations are constructed and used in evolutionary models of stellar populations in forming galaxies. Following their formation, the late type galaxies are assumed to follow the Schmidt law for star formation, while early type galaxies are normalized to the present-day fundamental plane relations assumed to mimic the metallicity variations along their luminosity sequence. We then compute predictions of these models for the observational data at early epochs for various cosmological parameters $\Omega, \Omega_\Lambda$ and $H_0$. We find good match to the metallicity data from the damped $L_\alpha$ systems and the evolution of the luminosity density out to $z\simeq 1$. Likewise, our models provide good fits for low values of $\Omega$ to the deep number counts of galaxies in all bands where data is available; this is done without assuming existence of extra populations of galaxies at high $z$. Our models also match the data on the redshift distribution of galaxy counts in $B$ and $K$ bands. We compute the predicted mean levels and angular distribution of the cosmic infrared background produced from the early evolution of galaxies. The predicted fluxes and fluctuations are still below the current observational limits, but not by a large factor. Finally, we find that the recent detection of the diffuse extragalactic light in the visible bands requires for our models high redshift of galaxy formation, $z_f \geq$(3-4); otherwise the produced flux of the extragalactic light at optical bands exceeds the current observational limits.

It is shown how data on the cluster correlation function can be used in order to reconstruct the density of the pregalactic density field on the cluster mass scale. The method is applied to the data on the cluster correlation amplitude -- richness dependence. The spectrum of the recovered density field has the same shape as the density field derived from data on the galaxy correlation function which is measured as function of linear scales. Matching the two amplitudes relates the mass to the comoving scale it contains and thereby leads to a direct determination of $\Omega$. The resultant density parameter turns out to be $\Omega$=0.25.

In this paper we evaluate the spectrum of the pregalactic density field on scales $1h^{-1}Mpc < r < 100h^{-1}$Mpc from a variety of astronomical data. APM data on $w(\theta)$ in six narrow magnitude is used, after correcting to possible evolutionary effects, to constrain the spectrum of galaxy clustering on scales $10h^{-1}Mpc < r < 50-100h^{-1}Mpc$. Fitting power spectra of CDM models to the data at all depths requires $\Omega h=0.2$ if the primordial index $n=1$ and $\Omega h=0.3$ if the spectrum is tilted with $n=0.7$. Then we compare the peculiar velocity field predicted by the APM spectrum of galaxy (light) distribution with the actual velocity data. The two fields are consistent and the comparison suggests that the bias factor is scale independent with $\Omega^{0.6}/b\simeq$(0.2-0.3). The next dataset used comes from the cluster correlation data. We calculate in detail the amplification of the cluster correlation function due to gravitational clustering and use the data on both the slope of the cluster correlation function and its amplitude-richness dependence. Cluster masses are normalized using the Coma cluster. We find that CDM models are hard to reconcile with all the three datasets: APM data on $w(\theta)$, the data on cluster correlation function, and the data on the latter's amplitude-richness dependence. We show that the data on the amplitude-richness dependence can be used directly to obtain the spectrum of the pregalactic density field. Applying the method to the data, we recover the density field on scales between 5 and 25$h^{-1}$Mpc whose slope is in good agreement with the APM data on the same scales. Requiring the two amplitudes to be the same, fixes the value of $\Omega$ to be 0.3 in agreement with observations of the dynamics of the Coma cluster. Finally we use the data

We discuss a new method of estimating the cosmic infrared background (CIB) from the spatial properties of infrared maps and give the limits on the CIB from applying it to the COBE/DIRBE maps. The strongest limits are obtained at mid- to far-IR where foregrounds are bright, but smooth. If the CIB comes from matter clustered like galaxies, the smoothness of the maps implies CIB levels less than $\sim$(10-15) nW/m$^2$/sr over this wavelength range.

The recently discovered high redshift galaxy, 53W091, with accurate age measurements (Dunlop et al 1996) provides a measure of the small-scale power of the primordial density field and, as we show, a crucial test of the inflation-inspired models. It allows $\Omega$=1 cosmologies only for low values of $H_0$, but then pushes formation of that galaxy to redshift much greater than allowed for by the cold-dark-matter density field. Cold-dark-matter (CDM) models with cosmological constant ($\Lambda$) and low $\Omega$ would decrease the redshift at which this galaxy has collapsed. However, in CDM models decreasing $\Omega$ suppresses the small scale power in the density field and this effect turns out to be dominant. We estimate the mass of the galaxy and show that it represents a very rare and unlikely event in the density field of such models. Similar problems would occur in other modifications of the CDM cosmogonies.

The epoch of galaxy formation provides an important additional test of cosmological theories. Cold-dark-matter (CDM) models with cosmological constant ($\Lambda$) are designed to account for the observed excess power in galaxy distribution, but at the same time suppress the power on small scales pushing galaxy formation to recent epochs. We point out that the recently discovered high redshift galaxy, 53W091, with accurate age measurements (Dunlop et al 1996) provides a particularly important test of these models. In the flat $\Lambda$-dominated Universe, the redshift of formation of 53W091 decreases with decreasing $\Omega$. However, in the modified CDM models decreasing $\Omega$ suppresses the small scale power in the density field and this effect turns out to be dominant. We estimate the mass of the galaxy and show that it represents a very rare and unlikely event in the density field of such models. Similar problems would occur in other modifications of the CDM cosmogonies.

We measure the smoothness of the infrared sky using the COBE DIRBE maps, and obtain interesting limits on the production of the diffuse cosmic infrared background (CIB) light by matter clustered like galaxies. The predicted fluctuations of the CIB with the DIRBE beam size of 0.7\deg\ are of the order of 10\%, and the maps are smooth at the level of $\delta\nu I_\nu \sim $ a few \nwm2sr rms from 2.2 to 100 \um. The lowest numbers are achieved at mid- to far-IR where the foreground is bright but smooth; they are $\sqrt{C(0)}$$\leq$$(1-1.5)$ \nwm2sr at $\lambda=$ 10-100 \um. If the CIB comes from clustered matter evolving according to typical scenarios, then the smoothness of the maps implies CIB levels less than $\sim (10-15) $ \nwm2sr over this wavelength range.

The cosmic infrared background from galaxies should have a minimum fluctuation of the order of 10\% on angular scales probed by the DIRBE (Diffuse Infrared Background Experiment) on board COBE (Cosmic Background Explorer). We present the theoretical basis for the analysis and discuss the all-sky all-wavelengths survey that we conducted. We further show that a linear combination of maps at different wavelengths can greatly reduce the fluctuations produced by foreground stars in the near-IR, while not eliminating the fluctuations of the background from high redshift galaxies. The method is potentially very powerful, especially at wavelengths in the mid- to far-IR where the foreground is bright but smooth.

This paper is devoted to studying the CIB through its correlation properties. We studied the limits on CIB anisotropy in the near IR (1.25, 2.2, and 3.5 \um, or $J,\;K,\;L$) bands at a scale of 0.7\deg\ using the COBE\footnote{ The National Aeronautics and Space Administration/Goddard Space Flight Center (NASA/GSFC) is responsible for the design, development, and operation of the {\it COBE}. Scientific guidance is provided by the {\it COBE} Science Working Group. GSFC is also responsible for the development of the analysis software and for the production of the mission data sets.} Diffuse Infrared Background Experiment (DIRBE) data. In single bands we obtain the upper limits on the zero-lag correlation signal $C(0)= \langle(\nu \delta I_\nu)^2\rangle < 3.6 \times 10^{-16},\; 5.1 \times 10^{-17},\; 5.7 \times 10^{-18}$ \w2m4sr2 for the $J,K,L$ bands respectively. The DIRBE data exhibit a clear color between the various bands with a small dispersion. On the other hand most of the CIB is expected to come from redshifted galaxies and thus should have different color properties. We use this observation to develop a `color subtraction' method of linear combinations of maps at two different bands. This method is expected to suppress the dominant fluctuations from foreground stars and nearby galaxies, while not reducing (or perhaps even amplifying) the extragalactic contribution to $C(0)$. Applying this technique gives significantly lower and more isotropic limits.

We outline a new method for estimating the cosmic infrared background using the spatial and spectral correlation properties of infrared maps. The cosmic infrared background from galaxies should have a minimum fluctuation of the order of 10\% on angular scales of the order of 1\deg. We show that a linear combination of maps at different wavelengths can greatly reduce the fluctuations produced by foreground stars, while not eliminating the fluctuations of the background from high redshift galaxies. The method is potentially very powerful, especially at wavelengths where the foreground is bright but smooth.

We study the microwave background anisotropy due to superhorizon-size perturbations (the Grischuk-Zel'dovich effect) in open universes with negative spatial curvature. Using COBE results on the low-order temperature multipole moments, we find that if the homogeneity of the observable Universe arises from an early epoch of inflation, the present density parameter cannot differ from unity by more than the observed quadrupole anisotropy, $|1-\Omega_0| \alt Q \simeq 5\times 10^{-6}$. Thus, inflation models with low $\Omega_0$ either do not fit the microwave background observations or they do not solve the horizon problem.

We review how the various large-scale data constrain cosmological parameters and, consequently, theories for the origin of large-scale structure in the Universe. We discuss the form of the power spectrum implied by the correlation data of galaxies and argue by comparing the velocity field implied by the distribution of light with the observed velocity flows that the bias parameter, $b$, is likely to be constant in the linear regime. This then allows one to estimate the density parameter, $\Omega$, and $b$ directly from the \underline{data} on $\xi(r)$ and the velocity fields. We show that it is consistent with low values of $\Omega^{0.6}/b$. We discuss the ways to normalise the optical data at $z\sim0$ directly to the COBE (or other microwave background) data. The data on high-$z$ \underline{galaxies} allows one to further constrain the shape of the \underline{primordial} power spectrum at scales which are non-linear today ($< 8h^{-1}$Mpc) and we discuss the consistency of the data with inflationary models normalised to the large-scale structure observations.