The Hubble constant $H_{0}$ and matter density $\Omega_{m}$ of the Universe are measured using the latest $\gamma$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $\gamma$-ray attenuation along the line of sight depends on the expansion rate and matter content of the Universe. This novel strategy results in a value of $H_{0}=67.4_{-6.2}^{+6.0}$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.14_{-0.07}^{+0.06}$. These estimates are independent and complementary to those based on the distance ladder, cosmic microwave background (CMB), clustering with weak lensing, and strong lensing data. We also produce a joint likelihood analysis of our results from $\gamma$ rays and these from more mature methodologies, excluding the CMB, yielding a combined value of $H_{0}=66.6\pm 1.6$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.29\pm 0.02$.

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 Extragalactic Background Light (EBL) can be probed via the absorption imprint it leaves in the spectra of gamma-ray sources ($\gamma\gamma \rightarrow e^-e^+$). We recently developed a dedicated technique to reconstruct the EBL, and its evolution with redshift, from $\gamma$ ray optical depth data using a large sample of blazars detected by the $Fermi$ Large Area Telescope. Here, we extend this dataset to the TeV regime using ground-based Cherenkov observations of 38 blazars and report the first homogeneous measurement of the EBL spectral intensity covering the ultraviolet to infrared wavelengths ($\sim$0.1-100$\mathrm{\mu m}$). A minimal EBL throughout the wavelength range with respect to integrated galaxy light is found, allowing little additional unresolved emission from faint or truly diffuse populations setting an upper limit of $\lesssim 4~{\rm nW\cdot m^{-2}sr^{-1}}$ at 1.4 ${\rm \mu m}$. In particular, the cosmic optical background (COB) at $z=0$ is found to be $27.8_{-2.0}^{+2.1}~{\rm nW\cdot m^{-2}sr^{-1}}$. This work lays the foundation for accurate gamma-ray measurements of the EBL across its whole spectral range using a combination of GeV and TeV data.

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.

The light emitted by all galaxies over the history of the Universe produces the extragalactic background light (EBL) at ultraviolet, optical, and infrared wavelengths. The EBL is a source of opacity for $\gamma$ rays via photon-photon interactions, leaving an imprint in the spectra of distant $\gamma$-ray sources. We measure this attenuation using {739} active galaxies and one gamma-ray burst detected by the {\it Fermi} Large Area Telescope. This allows us to reconstruct the evolution of the EBL and determine the star-formation history of the Universe over 90\% of cosmic time. Our star-formation history is consistent with independent measurements from galaxy surveys, peaking at redshift $z\sim2$. Upper limits of the EBL at the epoch of re-ionization suggest a turnover in the abundance of faint galaxies at $z\sim 6$.

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.

We use high resolution hydrodynamical simulations to study the contribution to the X-ray background from high-$z$ energetic sources, such as X-ray binaries, accreting nuclear black holes and shock heated interstellar medium. Adopting the model discussed in Eide et al. (2018), we find that these X-ray sources during the Epoch of Reionization (EoR) contribute less than a few percent of the unresolved X-ray background. The same sources contribute to less than $\sim$2\% of the measured angular power spectrum of the fluctuations of the X-ray background. The outputs of radiative transfer simulations modeling the EoR are used to evaluate the cross-correlations of X-ray background with the 21~cm signal from neutral hydrogen. Such correlation could be used to confirm the origin of the 21 cm signal, as well as give information on the properties of the X-ray sources during the EoR. We find that the correlations are positive during the early stages of reionization when most of the hydrogen is neutral, while they become negative when the intergalactic medium gets highly ionized, with the transition from positive to negative depending on both the X-ray model and the scale under consideration. With {\tt SKA} as the reference instrument for the 21~cm experiment, the predicted S/N for such correlations is $<1$ if the corresponding X-ray survey is only able to resolve and remove X-ray sources with observed flux $>10^{-15}\,\rm erg\, cm^{-2} \, s^{-1}$, while the cumulative S/N from $l=1000$ to $10^{4}$ at $x_{\rm HI}=0.5$ is $\sim 5$ if sources with observed flux $>10^{-17}\,\rm erg\, cm^{-2} \, s^{-1}$ are detected.

We study cross-correlations of the kinetic Sunyaev-Zel'dovich effect (kSZ) and 21 cm signals during the epoch of reionisation (EoR) to measure the effects of patchy reionisation. Since the kSZ effect is proportional to the line-of-sight velocity, the kSZ-21 cm cross correlation suffers from cancellation at small angular scales. We thus focus on the correlation between the kSZ-squared field (kSZ$^2$) and 21 cm signals. When the global ionisation fraction is low ($x_e\lesssim 0.7$), the kSZ$^2$ fluctuation is dominated by rare ionised bubbles which leads to an anti-correlation with the 21 cm signal. When $0.8\lesssim x_e<1$, the correlation is dominated by small pockets of neutral regions, leading to a positive correlation. However, at very high redshifts when $x_e<0.15$, the spin temperature fluctuations change the sign of the correlation from negative to positive, as weakly ionised regions can have strong 21 cm signals in this case. To extract this correlation, we find that Wiener filtering is effective in removing large signals from the primary CMB anisotropy. The expected signal-to-noise ratios for a $\sim$10-hour integration of upcoming Square Kilometer Array data cross-correlated with maps from the current generation of CMB observatories with 3.4~$\mu$K arcmin noise and 1.7~arcmin beam over 100~deg$^2$ are 51, 60, and 37 for $x_e=0.2$, 0.5, and 0.9, respectively.

The extragalactic background light (EBL), from ultraviolet to infrared wavelengths, is predominantly due to emission from stars, accreting black holes and reprocessed light due to Galactic dust. The EBL can be studied through the imprint it leaves, via $\gamma$-$\gamma$ absorption of high-energy photons, in the spectra of distant $\gamma$-ray sources. The EBL has been probed through the search for the attenuation it produces in the spectra of BL Lacertae (BL Lac) objects and individual $\gamma$-ray bursts (GRBs). GRBs have significant advantages over blazars for the study of the EBL especially at high redshifts. Here we analyze a combined sample of twenty-two GRBs, detected by the Fermi Large Area Telescope between 65 MeV and 500 GeV. We report a marginal detection (at the ~2.8 $\sigma$ level) of the EBL attenuation in the stacked spectra of the source sample. This measurement represents a first constraint of the EBL at an effective redshift of ~1.8. We combine our results with prior EBL constraints and conclude that Fermi-LAT is instrumental to constrain the UV component of the EBL. We discuss the implications on existing empirical models of EBL evolution.

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.

We show that measurements of the fluctuations in the near-infrared background (NIRB) from the AKARI satellite can be explained by faint galaxy populations at low redshifts. We demonstrate this using reconstructed images from deep galaxy catalogs (HUGS/S-CANDELS) and two independent galaxy population models. In all cases, we find that the NIRB fluctuations measured by AKARI are consistent with faint galaxies and there is no need for a contribution from unknown populations. We find no evidence for a steep Rayleigh-Jeans spectrum for the underlying sources as previously reported. The apparent Rayleigh-Jeans spectrum at large angular scales is likely a consequence of galaxies being removed systematically to deeper levels in the longer wavelength channels.

Emission from high-$z$ galaxies must unquestionably contribute to the near-infrared background (NIRB). However, this contribution has so far proven difficult to isolate even after subtracting the resolved galaxies to deep levels. Remaining NIRB fluctuations are dominated by unresolved low-$z$ galaxies on small angular scales, and by an unidentified component with unclear origin on large scales ($\approx 1000''$). In this paper, by analyzing mock maps generated from semi-numerical simulations and empirically determined $L_{\rm UV} - M_{\rm h}$ relations, we find that fluctuations associated with galaxies at $5 < z < 10$ amount to several percent of the unresolved NIRB flux fluctuations. We investigate the properties of this component for different survey areas and limiting magnitudes. In all cases, we show that this signal can be efficiently, and most easily at small angular scales, isolated by cross-correlating the source-subtracted NIRB with Lyman Break Galaxies (LBGs) detected in the same field by {\tt HST} surveys. This result provides a fresh insight into the properties of reionization sources.

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.

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.

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.

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.