We derive new CMB anisotropy power spectrum and BBN constraints for evaporating primordial black holes by explicitly solving the electromagnetic particle cascades of emitted particles and the deposition of this emitted energy to the background baryon-photon plasma. We show that the CMB anisotropies can provide stronger constraints compared to BBN and CMB spectral distortions on black holes with masses as small as $M_{\rm BH}=1.1\times 10^{13}$g, a slightly smaller mass than what has been considered in literature until now. We also show that, with more up-to-date data on abundances of deuterium and helium-3, BBN constraints are strengthened significantly.

We calculate the spectral distortions from Hawking evaporation of primordial black holes before the epoch of recombination, taking into account emission of all standard model particles, including quark and gluons, and evolving the resulting particle cascades in the expanding Universe. We show that the constraints on the abundance of primordial black holes are stronger by more than an order of magnitude compared to the previous calculations which take only the primary photon emission into account. We also show that the shapes of the spectral distortions is different from the $y$ or $i$-type distortions and are sensitive to the mass of the primordial black holes. We also extend previous constraints on the decay of long lived unstable particles before recombination to additional decay channels. We show that for dark matter mass $\lesssim$ 1 GeV, the spectral distortion shape is a function of the dark matter mass as well as the decay channel to standard model particles. We also provide new spectral distortion constraints on superconducting cosmic string decay. We explicitly show that consideration of emitted photon spectrum from string decay is not only important for the future experiments but also for already available COBE-FIRAS data.

We propose using the upcoming Cosmic Microwave Background (CMB) ground based experiments to detect the signal of ALPs (Axion like particles) interacting with magnetic fields in galaxy clusters. The conversion between CMB photons and ALPs in the presence of the cluster magnetic field can cause a polarized spectral distortion in the CMB around a galaxy cluster. The strength of the signal depends upon the redshift of the galaxy cluster and will exhibit a distinctive spatial profile around it depending upon the structure of electron density and magnetic field. This distortion produces a different shape from the other known spectral distortions like $y$-type and $\mu$-type and hence are separable from the multi-frequency CMB observation. The spectrum is close to kinematic Sunyaev-Zeldovich (kSZ) signal but can be separated from it using the polarization information. For the future ground-based CMB experiments such as Simons Observatory and CMB-S4, we estimate the measurability of this signal in the presence of foreground contamination, instrument noise and CMB anisotropies. This new avenue can probe the photon-ALP coupling over the ALP mass range from $10^{-13}$ eV to $10^{-12}$ eV with two orders of magnitude better accuracy from CMB-S4 than the current existing bounds.

We detect the diffuse thermal Sunyaev-Zeldovich (tSZ) effect from the gas filaments between the Luminous Red Galaxy (LRG) pairs using a new approach relying on stacking the individual frequency maps. We apply and demonstrate our method on ~88000 LRG pairs in the SDSS DR12 catalogue selected with an improved selection criterion that ensures minimal contamination by the Galactic CO emission as well as the tSZ signal from the clusters of galaxies. We first stack the Planck channel maps and then perform the Internal Linear Combination method to extract the diffuse $y_{\rm sz}$ signal. Our $Stack$ $First$ approach makes the component separation a lot easier as the stacking greatly suppresses the noise and CMB contributions while the dust foreground becomes homogeneous in spectral-domain across the stacked patch. Thus one component, the CMB, is removed while the rest of the foregrounds are made simpler even before component separation algorithm is applied. We obtain the WHIM signal of $y_{\rm whim}=(3.76\pm 0.44)\times 10^{-8}$ in the gas filaments, accounting for the electron overdensity of $\sim 13$. We estimate the detection significance to be $\approx 8.1\sigma$. This excess $y_{\rm sz}$ signal is tracing the warm-hot intergalactic medium and it could account for most of the missing baryons of the Universe. We show that the $Stack$ $First$ approach is more robust to systematics and produces a cleaner signal compared to the methods relying on stacking the $y$-maps to detect weak tSZ signal currently being used by the cosmology community.

Injection of high energy electromagnetic particles around the recombination epoch can modify the standard recombination history and therefore the CMB anisotropy power spectrum. Previous studies have put strong constraints on the amount of electromagnetic energy injection around the recombination era (redshifts $z\lesssim 4500$). However, energy injected in the form of energetic ($>$ keV) visible standard model particles is not deposited instantaneously. The considerable delay between the time of energy injection and the time when all energy is deposited to background baryonic gas and CMB photons, together with the extraordinary precision with which the CMB anisotropies have been measured, means that CMB anisotropies are sensitive to energy that was injected much before the epoch of recombination. We show that the CMB anisotropy power spectrum is sensitive to energy injection even at $z = 10000$, giving stronger constraints compared to big bang nucleosynthesis and CMB spectral distortions. We derive, using Planck CMB data, the constraints on long-lived unstable particles decaying at redshifts $z\lesssim 10000$ (lifetime $\tau_X\gtrsim 10^{11}$s) by explicitly evolving the electromagnetic cascades in the expanding Universe, thus extending previous constraints to lower particle lifetimes. We also revisit the BBN constraints and show that the delayed injection of energy is important for BBN constraints. We find that the constraints can be weaker by a factor of few to almost an order of magnitude, depending on the energy, when we relax the quasi-static or on-the-spot assumptions.

Precision measurement of the cosmological recombination spectrum can provide an entire new window to look at the early universe. We aim to quantify the information hidden in the cosmological recombination spectrum and for this purpose we have developed a new code following the algorithm proposed in Ali-Haimoud and Hirata (2010), Ali-Haimoud (2013). Our code is closely based on the COSMOSPEC code [Chluba and Ali-Haimoud (2016)]. We find, using Fisher information matrix and assuming that the foregrounds can be subtracted by using higher or lower frequency channels and spatial information, that going beyond the detection will need an experiment with sensitivity $25$ times better compared to the proposed experiment PIXIE. Such an experiment will be able to measure the cosmological parameters with a precision that is competitive with the CMB anisotropy experiments. The best constrainted parameter is baryon energy density, $\Omega_{\rm b}$, which can be nailed down with incredible precision in principle. We also show that the shape of the hydrogen lines is connected to the speed of the hydrogen recombination, with the peaks of the recombination lines coinciding with the peak of the recombination rate. In general, the shape of the lines encodes information about the rate of recombination as a function of redshift.

New interactions of neutrinos can stop them from free streaming even after the weak interaction freezeout. This results in a phase shift in the cosmic microwave background (CMB) acoustic peaks which can alleviate the Hubble tension. We demonstrate with Planck CMB and WiggleZ galaxy survey data that this acoustic phase shift, and thus solution to the Hubble tension, can be achieved for neutrinos interacting with dark matter without significantly affecting other observables and without changing the number of relativistic degrees of freedom. We predict potentially observable modification of the CMB B-modes.

We present a new method to determine the probability distribution of the 3-D shapes of galaxy clusters from the 2-D images using stereology. In contrast to the conventional approach of combining different data sets (such as X-rays, Sunyaev-Zeldovich effect and lensing) to fit a 3-D model of a galaxy cluster for each cluster, our method requires only a single data set, such as X-ray observations or Sunyaev-Zeldovich effect observations, consisting of sufficiently large number of clusters. Instead of reconstructing the 3-D shape of an individual object, we recover the probability distribution function (PDF) of the 3-D shapes of the observed galaxy clusters. The shape PDF is the relevant statistical quantity which can be compared with the theory and used to test the cosmological models. We apply this method to publicly available \emph{Chandra} X-ray data of 89 well resolved galaxy clusters. Assuming ellipsoidal shapes, we find that our sample of galaxy clusters is a mixture of prolate and oblate shapes, with a preference for oblateness with the most probable ratio of principle axes 1.4 : 1.3 : 1. The ellipsoidal assumption is not essential to our approach and our method is directly applicable to non-ellipsoidal shapes. Our method is insensitive to the radial density and temperature profiles of the cluster. Our method is sensitive to the changes in shape of the X-ray emitting gas from inner to outer regions and we find evidence for variation in the 3-D shape of the X-ray emitting gas with distance from the centre.

Current constraints on energy injection in the form of energetic particles before the epoch of recombination using CMB spectral distortions assume that all energy goes into $y$ and $\mu$-type distortions. We revisit these constraints with exact calculations of the spectral distortions by evolving the electromagnetic cascades. The actual spectral distortion differs in shape and amplitude from the $y$-type distortion and depends on the energy and nature of injected particles. The constraints on the energy injection processes such as dark matter decay can be relaxed by as much as a factor of 5.

The non-resonant conversion of Cosmic Microwave Background (CMB) photons into scalar as well as light pseudoscalar particles such as axion-like particles (ALPs) in the presence of turbulent magnetic fields can cause a unique, spatially fluctuating spectral distortion in the CMB. We use the publicly available Planck temperature maps for the frequency channels (70-545 GHz) to obtain the first ALP distortion map using $45\%$ clean part of the sky. The $95^{th}$ percentile upper limit on the RMS fluctuation of ALP distortions from the cleanest part of the CMB sky at $15$ arcmin angular resolution is $18.5 \times 10^{-6}$. The RMS fluctuation in the distortion map is also consistent with different combinations of frequency channels and sky-fractions.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Dark matter may have been relativistic and collisional until relatively late times and become cold and collisionless after a phase transition before the matter-radiation equality of the standard $\Lambda$CDM cosmology. We show that such a dark matter has large peculiar velocities due to acoustic oscillations before the phase transition, and evolves ballistically after the phase transition in the collisionless phase until the initial acoustic velocities are redshifted away. We show that this Ballistic Dark Matter (BDM) results in new non-trivial interesting features in the cosmological observables. In particular, the linear matter power spectrum exhibits acoustic oscillations on scales smaller than the Hubble scale at the time of phase transition, and for fast transitions the power at the acoustic peaks in the matter power spectrum exceeds that in a $\Lambda$CDM cosmology. If BDM only forms a part of the total dark matter, an odd vs. even acoustic peak asymmetry becomes prominent. We give an approximate analytical treatment of the linear perturbations in BDM, explaining these features. We also discuss the possibility to constrain BDM using cosmic microwave background and large scale structure data.

It is generally assumed that for energy injection before recombination, all of the injected energy is dissipated as heat in the baryon-photon plasma, giving rise to the $y$-type, $i$-type, and $\mu$-type distortions in the CMB spectrum. We show that this assumption is incorrect when the energy is injected in the form of energetic (i.e. energy much greater than the background CMB temperature) particles. We evolve the electromagnetic cascades, from the injection of high energy particles, in the expanding Universe and follow the non-thermal component of CMB spectral distortions resulting from the interaction of the electromagnetic shower with the background photons, electrons, and ions. The electromagnetic shower loses a substantial fraction of its energy to the CMB spectral distortions before the energy of the particles in the shower has degraded to low enough energies that they can thermalize with the background plasma. This spectral distortion is the result of the interaction of non-thermal energetic electrons in the shower with the CMB and thus has a shape that is substantially different from the $y$-type or $i$-type distortions. The shape of the final \emph{non-thermal relativistic} ($ntr$-type) CMB spectral distortion depends upon the initial energy spectrum of the injected electrons, positrons, and photons and thus has information about the energy injection mechanism e.g. the decay or annihilation channel of the decaying or annihilating dark matter particles. The shape of the spectral distortion is also sensitive to the redshift of energy injection. Our calculations open up a new window into the energy injection at $z\lesssim 2\times 10^5$ which is not degenerate with, and can be distinguished from the low redshift thermal $y$-type distortions.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

We present a new approach to component separation in multifrequency CMB experiments by formulating the problem as that of partitioning the sky into pixel clusters such that within each pixel cluster the foregrounds have similar spectrum, using only the information available in the data. Only spectral information is used for partitioning, allowing spatially far away pixels to belong to the same cluster if their foreground properties are close. We then apply a modified internal linear combination method to each pixel cluster. Since the foregrounds have similar spectrum within each cluster, the number of components required to describe the foregrounds is smaller compared to all data taken together and simple pixel based ILC algorithm works extremely well. We test our algorithm in the full focal plane simulations provided by the Planck collaboration. We apply our algorithm to the Planck full mission data and compare our CMB maps with the CMB maps released by the Planck collaboration. We show that our CMB maps are clean and unbiased on a larger fraction of the sky, especially at the low Galactic latitudes, compared to publicly available maps released by the Planck collaboration. This is important for constraining beyond the simplest $\Lambda$CDM cosmological models and study of anomalies. Our cleaned CMB maps are made publicly available for use by the cosmology community.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We analyze two volume limited galaxy samples from the SDSS photometric and spectroscopic data to test the isotropy in the local Universe. We use information entropy to quantify the global anisotropy in the galaxy distribution at different length scales and find that the galaxy distribution is highly anisotropic on small scales. The observed anisotropy diminishes with increasing length scales and nearly plateaus out beyond a length scale of 200 Mpc/h in both the datasets. We compare these anisotropies with those predicted by the mock catalogues from the N-body simulations of the Lambda CDM model and find a fairly good agreement with the observations. We find a small residual anisotropy on large scales which decays in a way that is consistent with the linear perturbation theory. The slopes of the observed anisotropy converge to the slopes predicted by the linear theory beyond a length scale of ~ 200 Mpc/h indicating a transition to isotropy. We separately compare the anisotropies observed across the different parts of the sky and find no evidence for a preferred direction in the galaxy distribution.

We revisit the cosmological constraints on resonant and non-resonant conversion of photons to axions in the cosmological magnetic fields. We find that the constraints on photon-axion coupling and primordial magnetic fields are much weaker than previously claimed for low mass axion like particles with masses $m_{\rm a} \lesssim 5\times 10^{-13}\,\text{eV}$. In particular we find that the axion mass range $10^{-14}\,\text{eV} \le m_{\rm a} \le 5\times 10^{-13}\,\text{eV}$ is not excluded by the CMB data contrary to the previous claims. We also examine the photon-axion conversion in the Galactic magnetic fields. Resonant conversion in the large scale coherent Galactic magnetic field results in $100\%$ polarized anisotropic spectral distortions of the CMB for the mass range $10^{-13}\,\text{eV} \lesssim m_{\rm a} \lesssim 10^{-11}\,\text{eV}$. The polarization pattern traces the transverse to line of sight component of the Galactic magnetic field while both the anisotropy in the Galactic magnetic field and electron distribution imprint a characteristic anisotropy pattern in the spectral distortion. Our results apply to scalar as well as pseudoscalar particles. For conversion to scalar particles, the polarization is rotated by $90^{\circ}$ allowing us to distinguish them from the pseudoscalars. For $m_{\rm a} \lesssim 10^{-14} \,\text{eV}$ we have non-resonant conversion in the small scale turbulent magnetic field of the Galaxy resulting in anisotropic but unpolarized spectral distortion in the CMB. These unique signatures are potential discriminants against the isotropic and non-polarized signals such as primary CMB, and $\mu$ and $y$ distortions with the anisotropic nature making it accessible to experiments with only relative calibration like Planck, LiteBIRD, and CORE. We forecast for PIXIE as well as for these experiments using Fisher matrix formalism.

New interactions of neutrinos can stop them from free streaming in the early Universe even after the weak decoupling epoch. This results in the enhancement of the primordial gravitational wave amplitude on small scales compared to the standard $\Lambda$CDM prediction. In this paper we calculate the effect of dark matter neutrino interactions in CMB tensor $B$-modes spectrum. We show that the effect of new neutrino interactions generates a scale or $\ell$ dependent imprint in the CMB $B$-modes power spectrum at $\ell \gtrsim 100$. In the event that primordial $B$-modes are detected by future experiments, a departure from scale invariance, with a blue spectrum, may not necessarily mean failure of simple inflationary models but instead may be a sign of non-standard interactions of relativistic particles. New interactions of neutrinos also induce a phase shift in the CMB B-mode power spectrum which cannot be mimicked by simple modifications of the primordial tensor power spectrum. There is rich information hidden in the CMB $B$-modes spectrum beyond just the tensor to scalar ratio.

We report the detection of thermal Sunyaev-Zeldovich effect (SZ) fluctuations in the intracluster medium (ICM) of Coma cluster observed with ${\it Planck}$. The SZ data links the maximum observable X-ray scale to the large Mpc scale, extending our knowledge of the power spectrum of ICM fluctuations. Deprojecting the 2-d SZ perturbations into 3-d pressure fluctuations, we find an amplitude spectrum which peaks at $\delta P/P = 33\pm 12\%$ and $74\pm19\%$ in the $15'$ and $40'$ radius region, respectively. We perform tests to ensure fluctuations are intrinsic to the cluster and not due to noise contamination. By using high-resolution hydrodynamical models, we improve the ICM turbulence constraints in Coma, finding 3-d Mach number ${\rm Ma_{3d}}= 0.8\pm0.3$ (15' region), increasing to supersonic values at larger radii (40'), and an injection scale $L_{\rm inj}\approx 500$ kpc. Such properties are consistent with driving due to mergers, in particular tied to internal galaxy groups. The large pressure fluctuations show that Coma is in adiabatic mode (mediated by sound waves), rather than isobaric mode (mediated by buoyancy waves). As predicted by turbulence models, the distribution of SZ fluctuations is log-normal with mild non-Gaussianities (heavy tails). The substantial non-thermal pressure support implies hydrostatic mass bias $b_M=-15\%$ to $-45\%$ from the core to the outskirt region, respectively. While total SZ power probes the thermal energy content, the SZ fluctuations constrain the non-thermal deviations important for precision cosmology. The proposed, novel approach can be exploited by multifrequency observations using ground based interferometers and future space CMB missions.

We present an all sky map of the $y$-type distortion calculated from the full mission Planck HFI (High Frequency Instrument) data using the recently proposed approach to component separation based on parametric model fitting and model selection. This simple model selection approach allows us to distinguish between carbon monoxide (CO) line emission and $y$-type distortion, something that is not possible using the internal linear combination based methods. We create a mask to cover the regions of significant CO emission relying on the information in the $\chi^2$ map obtained when fitting for the $y$-distortion and CO emission to the lowest four HFI channels. We revisit the second Planck cluster catalog and try to quantify the quality of the cluster candidates in an approach that is similar in spirit to Aghanim et al. (2014). We find that at least $93\%$ of the clusters in the cosmology sample are free of CO contamination. We also find that $59\%$ of unconfirmed candidates may have significant contamination from molecular clouds. We agree with Planck collaboration (2015) for the worst offenders. We suggest an alternative validation strategy of measuring and subtracting the CO emission from the Planck cluster candidates using radio telescopes thus improving the reliability of the catalog. Our CO mask and annotations to the Planck cluster catalog identifying cluster candidates with possible CO contamination are made publicly available.

We suggest a method for probing global properties of clump populations in Giant Molecular Clouds (GMCs) in the case where these act as X-ray reflection nebulae (XRNe), based on the study of the clumping's overall effect on the reflected X-ray signal, in particular on the Fe K-alpha line's shoulder. We consider the particular case of Sgr B2, one of the brightest and most massive XRN in our Galaxy. We parametrise the gas distribution inside the cloud using a simple clumping model, with the slope of the clump mass function (alpha), the minimum clump mass (m_{min}), the fraction of the cloud's mass contained in clumps (f_{DGMF}), and the mass-size relation of individual clumps as free parameters, and investigate how these affect the reflected X-ray spectrum. In the case of very dense clumps, similar to those presently observed in Sgr B2, these occupy a small volume of the cloud and present a small projected area to the incoming X-ray radiation. We find that these contribute negligibly to the scattered X-rays. Clump populations with volume filling factors of > 10^{-3}, do leave observational signatures, that are sensitive to the clump model parameters, in the reflected spectrum and polarisation. Future high-resolution X-ray observations could therefore complement the traditional optical and radio observations of these GMCs, and prove to be a powerful probe in the study of their internal structure. Finally, clumps in GMCs should be visible both as bright spots and regions of heavy absorption in high resolution X-ray observations. We therefore further study the time-evolution of the X-ray morphology, under illumination by a transient source, as a probe of the 3d distribution and column density of individual clumps by future X-ray observatories.

We use the Planck HFI channel maps to make an all sky map of $\mu$-distortion fluctuations. Our $\mu$-type distortion map is dominated by the $y$-type distortion contamination from the hot gas in the low redshift Universe and we can thus only place upper limits on the $\mu$-type distortion fluctuations. For the amplitude of $\mu$-type distortions on $10'$ scales we get the limit on root mean square (rms) value $\mu_{rms}^{10'}< 6.4\times 10^{-6}$, a limit 14 times stronger than the COBE-FIRAS ($95\%$ confidence) limit on the mean of $< \mu > <90\times 10^{-6}$. Using our maps we also place strong upper limits on the auto angular power spectrum of $\mu$, $C_{\ell}^{\mu\mu}$ and the cross angular power spectrum of $\mu$ with the CMB temperature anisotropies, $C_{\ell}^{\mu T}$. The strongest observational limits are on the largest scales, $\ell(\ell+1)/(2\pi)C_{\ell}^{\mu\mu}|_{\ell=2-26}<(2.3\pm 1.0)\times 10^{-12}$ and $\ell(\ell+1)/(2\pi)C_{\ell}^{\mu T}|_{\ell=2-26}<(2.6\pm 2.6)\times 10^{-12}~{K}$. Our observational limits can be used to constrain new physics which can create spatially varying energy release in the early Universe between redshifts $5\times 10^4\lesssim z\lesssim 2\times 10^6$. We specifically apply our observational results to constrain the primordial non-Gaussianity of the local type, when the source of $\mu$-distortion is Silk damping, for very squeezed configurations with the wavenumber for the short wavelength mode $46 \lesssim k_{S} \lesssim 10^4 ~{Mpc^{-1}}$ and for the long wavelength mode $k_{L}\approx 10^{-3} ~{Mpc^{-1}}$. Our limits on the primordial non-Gaussianity parameters are $f_{NL}<10^5, \tau_{NL}<1.4\times 10^{11}$ for $k_{S}/k_{L}\approx 5\times 10^4- 10^7$. We give a new derivation of the evolution of the $\mu$-distortion fluctuations. We also introduce mixing of Bose-Einstein spectra and $y^{BE}$-type distortions.

We use the published Planck and SPT cluster catalogs and recently published $y$-distortion maps to put strong observational limits on the contribution of the fluctuating part of the $y$-type distortions to the $y$-distortion monopole. Our bounds are $5.4\times 10^{-8} < \langle y\rangle < 2.2\times 10^{-6}$. Our upper bound is a factor of 6.8 stronger than the currently best upper $95\%$ confidence limit from COBE-FIRAS of $\langle y\rangle <15\times 10^{-6}$. In the standard cosmology, large scale structure is the only source of such distortions and our limits therefore constrain the baryonic physics involved in the formation of the large scale structure. Our lower limit, from the detected clusters in the Planck and SPT catalogs, also implies that a Pixie-like experiment should detect the $y$-distortion monopole at $>27$-$\sigma$. The biggest sources of uncertainty in our upper limit are the monopole offsets between different HFI channel maps that we estimate to be $<10^{-6}$.

We present an efficient algorithm for the least squares parameter fitting optimized for component separation in multi-frequency CMB experiments. We sidestep some of the problems associated with non-linear optimization by taking advantage of the quasi-linear nature of the foreground model. We demonstrate our algorithm, linearized iterative least-squares (LIL), on the publicly available Planck sky model FFP6 simulations and compare our result with the other algorithms. We work at full Planck resolution and show that degrading the resolution of all channels to that of the lowest frequency channel is not necessary. Finally we present results for the publicly available Planck data. Our algorithm is extremely fast, fitting 6 parameters to 7 lowest Planck channels at full resolution (50 million pixels) in less than 160 CPU-minutes (or few minutes running in parallel on few tens of cores). LIL is therefore easily scalable to future experiments which may have even higher resolution and more frequency channels. We also naturally propagate the uncertainties in different parameters due to noise in the maps as well as degeneracies between the parameters to the final errors on the parameters using Fisher matrix. One indirect application of LIL could be a front-end for Bayesian parameter fitting to find the maximum of the likelihood to be used as the starting point for the Gibbs sampling. We show for rare components, such as the carbon-monoxide emission, present in small fraction of sky, the optimal approach should combine parameter fitting with model selection. LIL may also be useful in other astrophysical applications which satisfy the quasi-linearity criteria.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

We predict a thin (scale height $\sim$ 80 pc) diffuse component of the Galactic Ridge X-ray emission (GRXE) arising from the scattering of the radiation of bright X-ray binaries (XBs) by the interstellar medium (ISM). The morphology of this scattered component is expected to trace the clumpy molecular and HI clouds. We calculate this contribution to the GRXE from known Galactic XBs assuming that they are all persistent. The known XBs sample is incomplete, however, because it is flux limited and spans the lifetime of X-ray astronomy ($\sim 50$ years). We therefore also use a simulated sample of sources, to estimate the diffuse emission we should expect in an optimistic case assuming that the X-ray luminosity of our Galaxy is on average similar to that of other galaxies. In the calculations we also take into account the enhancement of the total scattering cross-section due to coherence effects in the elastic scattering from multi-electron atoms and molecules. This scattered emission can be distinguished from the contribution of low X-ray luminosity stars by the presence of narrow fluorescent K-$\alpha$ lines of Fe, Si, and other abundant elements present in the ISM and by directly resolving the contribution of low X-ray luminosity stars. We find that within $1^\circ$ latitude of the Galactic plane the scattered emission contributes on average $10-30\%$ of the GRXE flux in the case of known sources and over $50\%$ in the case of simulated sources. In the latter case, the scattered component is found to even dominate the stellar emission in certain parts of the Galactic plane. X-rays with energies $\gtrsim 1$ keV from XBs should also penetrate deep inside the HI and molecular clouds, where they are absorbed and heat the ISM. We find that this heating rate dominates the heating by cosmic rays (assuming a solar neighborhood energy density) in a considerable part of the Galaxy. [abridged]

We present an approximate calculation of the full Bayesian posterior probability distribution for the local non-Gaussianity parameter $f_{\text{nl}}$ from observations of cosmic microwave background anisotropies within the framework of information field theory. The approximation that we introduce allows us to dispense with numerically expensive sampling techniques. We use a novel posterior validation method (DIP test) in cosmology to test the precision of our method. It transfers inaccuracies of the calculated posterior into deviations from a uniform distribution for a specially constructed test quantity. For this procedure we study toy cases that use one- and two-dimensional flat skies, as well as the full spherical sky. We find that we are able to calculate the posterior precisely under a flat-sky approximation, albeit not in the spherical case. We argue that this is most likely due to an insufficient precision of the used numerical implementation of the spherical harmonic transform, which might affect other non-Gaussianity estimators as well. Furthermore, we present how a nonlinear reconstruction of the primordial gravitational potential on the full spherical sky can be obtained in principle. Using the flat-sky approximation, we find deviations for the posterior of $f_{\text{nl}}$ from a Gaussian shape that become more significant for larger values of the underlying true $f_{\text{nl}}$. We also perform a comparison to the well-known estimator of Komatsu et al. [Astrophys. J. 634, 14 (2005)] and finally derive the posterior for the local non-Gaussianity parameter $g_{\text{nl}}$ as an example of how to extend the introduced formalism to higher orders of non-Gaussianity.

Silk damping at redshifts 1.5 x 10^4 < z < 2 x 10^6 erases CMB anisotropies on scales corresponding to the comoving wavenumbers 8 < k < 10^4 Mpc^-1 (10^5 < \ell < 10^8). This dissipated energy is gained by the CMB monopole, creating distortions from a blackbody in the CMB spectrum of the \mu-type and the i-type. We study, using Fisher matrices, the constraints we can get from measurements of these spectral distortions on the primordial power spectrum from future experiments such as Pixie, and how these constraints change as we change the frequency resolution and the sensitivity of the experiment. We show that the additional information in the shape of the $i$-type distortions, in combination with the \mu-type distortions, allows us to break the degeneracy between the amplitude and the spectral index of the power spectrum on these scales and leads to much tighter constraints. We quantify the information contained in both the \mu-type distortions and the i-type distortions taking into account the partial degeneracy with the y-type distortions and the temperature of the blackbody part of the CMB. We also calculate the constraints possible on the primordial power spectrum when the spectral distortion information is combined with the CMB anisotropies measured by the WMAP, SPT, ACT and Planck experiments.

Spectral features in the CMB energy spectrum contain a wealth of information about the physical processes in the early Universe, z < 2 x 10^6. The CMB spectral distortions are complementary to all other probes of cosmology. In fact, most of the information contained in the CMB spectrum is inaccessible by any other means. This review outlines the main physics behind the spectral features in the CMB throughout the history of the Universe, concentrating on the distortions which are inevitable and must be present at a level observable by the next generation of proposed CMB experiments. The spectral distortions considered here include spectral features from cosmological recombination, resonant scattering of CMB by metals during reionization which allows us to measure their abundances, y-type distortions during and after reionization and \mu-type and i-type (intermediate between \mu and y) distortions created at redshifts z > 1.5 x 10^4.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

y-type spectral distortions of the cosmic microwave background allow us to detect clusters and groups of galaxies, filaments of hot gas and the non-uniformities in the warm hot intergalactic medium. Several CMB experiments (on small areas of sky) and theoretical groups (for full sky) have recently published y-type distortion maps. We propose to search for two artificial hot spots in such y-type maps resulting from the incomplete subtraction of the effect of the motion induced dipole on the cosmic microwave background sky. This dipole introduces, at second order, additional temperature and y-distortion anisotropy on the sky of amplitude few \mu K which could potentially be measured by Planck HFI and Pixie experiments and can be used as a source of cross channel calibration by CMB experiments. This y-type distortion is present in every pixel and is not the result of averaging the whole sky. This distortion, calculated exactly from the known linear dipole, can be subtracted from the final y-type maps, if desired.

Silk damping in the early Universe, before and during recombination, erases anisotropies in the cosmic microwave background (CMB) on small scales. This power, which disappears from anisotropies, appears in the monopole as y-type, i-type and \mu-type distortions. The observation of the CMB spectral distortions will thus make available to us the information about the primordial power spectrum on scales corresponding to the comoving wavenumbers $8< k < 10^4 Mpc^{-1}$ increasing our total view of inflation, when combined with CMB anisotropies, to span 17 e-folds. These distortions can be understood simply as mixing of blackbodies of different temperatures and the subsequent comptonization of the resulting distortions.

In this article we focus on the astrophysical results and the related cosmological implications derived from recent microwave surveys, with emphasis to those coming from the Planck mission. We critically discuss the impact of systematics effects and the role of methods to separate the cosmic microwave background signal from the astrophysical emissions and each different astrophysical component from the others. We then review of the state of the art in diffuse emissions, extragalactic sources, cosmic infrared back- ground, and galaxy clusters, addressing the information they provide to our global view of the cosmic structure evolution and for some crucial physical parameters, as the neutrino mass. Finally, we present three different kinds of scientific perspectives for fundamental physics and cosmology offered by the analysis of on-going and future cosmic microwave background projects at different angular scales dedicated to anisotropies in total intensity and polarization and to absolute temperature.

We calculate numerical solutions and analytic approximations for the intermediate-type spectral distortions. Detection of a \mu-type distortion (saturated comptonization) in the CMB will constrain the time of energy injection to be at a redshift 2x10^6> z > 2x10^5, while a detection of a y-type distortion (minimal comptonization) will mean that there was heating of CMB at redshift z< 1.5x10^4. We point out that the partially comptonized spectral distortions, generated in the redshift range 1.5x10^4 < z x 2x10^5, are much richer in information than the pure y and \mu-type distortions. The spectrum created during this period is intermediate between y and \mu-type distortions and depends sensitively on the redshift of energy injection. These intermediate-type distortions cannot be mimicked by a mixture of y and \mu-type distortions at all frequencies and vice versa. The measurement of these intermediate-type CMB spectral distortions has the possibility to constrain precisely not only the amount of energy release in the early Universe but also the mechanism, for example, particle annihilation and Silk damping can be distinguished from particle decay. The intermediate-type distortion templates and software code using these templates to calculate the CMB spectral distortions for user-defined energy injection rate are made publicly available.

Mixing of blackbodies with different temperatures creates a spectral distortion which, at lowest order, is a y-type distortion, indistinguishable from the thermal y-type distortion produced by the scattering of CMB photons by hot electrons residing in clusters of galaxies. This process occurs in the radiation-pressure dominated early Universe, when the primordial perturbations excite standing sound waves on entering the sound horizon. Photons from different phases of the sound waves, having different temperatures, diffuse through the electron-baryon plasma and mix together. This diffusion, with the length defined by Thomson scattering, dissipates sound waves and creates spectral distortions in the CMB. Of the total dissipated energy, 2/3 raises the average temperature of the blackbody part of spectrum, while 1/3 creates a distortion of y-type. It is well known that at redshifts 10^5< z< 2x10^6, comptonization rapidly transforms y-distortions into a Bose-Einstein spectrum. The chemical potential of the Bose-Einstein spectrum is again 1/3 the value we would get if all the dissipated energy was injected into a blackbody spectrum but no extra photons were added. We study the mixing of blackbody spectra, emphasizing the thermodynamic point of view, and identifying spectral distortions with entropy creation. This allows us to obtain the main results connected with the dissipation of sound waves in the early Universe in a very simple way. We also show that mixing of blackbodies in general, and dissipation of sound waves in particular, leads to creation of entropy.

Silk damping of primordial small-scale perturbations in the photon-baryon fluid due to diffusion of photons inevitably creates spectral distortions in the CMB. With the proposed CMB experiment PIXIE it might become possible to measure these distortions and thereby constrain the primordial power spectrum at comoving wavenumbers 50 Mpc^{-1} < k < 10^4 Mpc^{-1}. Since primordial fluctuations in the CMB on these scales are completely erased by Silk damping, these distortions may provide the only way to shed light on otherwise unobservable aspects of inflationary physics. A consistent treatment of the primordial dissipation problem requires going to second order in perturbation theory, while thermalization of these distortions necessitates consideration of second order in Compton scattering energy transfer. Here we give a full 2x2 treatment for the creation and evolution of spectral distortions due to the acoustic dissipation process, consistently including the effect of polarization and photon mixing in the free streaming regime. We show that 1/3 of the total energy (9/4 larger than previous estimates) stored in small-scale temperature perturbations imprints observable spectral distortions, while the remaining 2/3 only raises the average CMB temperature, an effect that is unobservable. At high redshift dissipation is mainly mediated through the quadrupole anisotropies, while after recombination peculiar motions are most important. During recombination the damping of the higher multipoles is also significant. We compute the average distortion for several examples using CosmoTherm, analyzing their dependence on parameters of the primordial power spectrum. For one of the best fit WMAP7 cosmologies, with n_S=1.027 and n_run=-0.034, the cooling of baryonic matter practically compensates the heating from acoustic dissipation in the mu-era. (abridged)

The blackbody spectrum of CMB was created in the blackbody photosphere at redshifts z>2x10^6. At these early times, the Universe was dense and hot enough that complete thermal equilibrium between baryonic matter (electrons and ions) and photons could be established. Any perturbation away from the blackbody spectrum was suppressed exponentially. New physics, for example annihilation and decay of dark matter, can add energy and photons to CMB at redshifts z>10^5 and result in a Bose-Einstein spectrum with a non-zero chemical potential ($\mu$). Precise evolution of the CMB spectrum around the critical redshift of z~2x10^6 is required in order to calculate the $\mu$-type spectral distortion and constrain the underlying new physics. Although numerical calculation of important processes involved (double Compton process, comptonization and bremsstrahlung) is not difficult, analytic solutions are much faster and easier to calculate and provide valuable physical insights. We provide precise (better than 1%) analytic solutions for the decay of $\mu$, created at an earlier epoch, including all three processes, double Compton, Compton scattering on thermal electrons and bremsstrahlung in the limit of small distortions. This is a significant improvement over the existing solutions with accuracy ~10% or worse. We also give a census of important sources of energy injection into CMB in standard cosmology. In particular, calculations of distortions from electron-positron annihilation and primordial nucleosynthesis illustrate in a dramatic way the strength of the equilibrium restoring processes in the early Universe. Finally, we point out the triple degeneracy in standard cosmology, i.e., the $\mu$ and $y$ distortions from adiabatic cooling of baryons and electrons, Silk damping and annihilation of thermally produced WIMP dark matter are of similar order of magnitude (~ 10^{-8}-10^{-10}).

The difference in the adiabatic indices of photons and non-relativistic baryonic matter in the early Universe causes the electron temperature to be slightly lower than the radiation temperature. Thermalization of photons with a colder plasma results in the accumulation of photons in the Rayleigh-Jeans tail, aided by stimulated recoil, while the higher frequency spectrum tries to approach Planck spectrum at the electron temperature $T_{\gamma}^{final}=\Te<T_{\gamma}^{initial}$; i.e., Bose-Einstein condensation of photons occurs. We find new solutions of the Kompaneets equation describing this effect. No actual condensate is, in reality, possible since the process is very slow and photons drifting to low frequencies are efficiently absorbed by bremsstrahlung and double Compton processes. The spectral distortions created by Bose-Einstein condensation of photons are within an order of magnitude (for the present range of allowed cosmological parameters), with exactly the same spectrum but opposite in sign, of those created by diffusion damping of the acoustic waves on small scales corresponding to comoving wavenumbers $45< k< 10^4\, Mpc^{-1}$. The initial perturbations on these scales are completely unobservable today due to their being erased completely by Silk damping. There is partial cancellation of these two distortions, leading to suppression of $\mu$ distortions expected in the standard model of cosmology. The net distortion depends on the scalar power index $n_S$ and its running $d n_S/d\ln k$, and may vanish for special values of parameters, for example, for a running spectrum with, $n_S=1,d n_S/d\ln k=-0.038$. We arrive at an intriguing conclusion: even a null result, non-detection of $\mu$-type distortion at a sensitivity of $10^{-9}$, gives a quantitative measure of the primordial small-scale power spectrum.

One of the important light elements created during the big bang nucleosynthesis is Be-7 which then decays to Li-7 by electron capture when recombination becomes effective but well before the Saha equilibrium recombination is reached. This means that Be-7 should wait until its recombination epoch even though the half-life of the hydrogenic beryllium atom is only 106.4 days. We calculate when the conversion from primordial Be-7 to Li-7 occurs taking into account the population of the hyperfine structure sublevels and solving the kinetic equations for recombination, photoionization and conversion rate. We also calculate the energies and the spectrum of narrow neutrino doublet lines resulting from Be-7 decay.

Reply to Flambaum and Porsev comment arXiv:1004.2540 on "21 cm radiation - a new probe of variation in the fine structure constant" arXiv:astro-ph/0701752

Inhomogeneous recombination can give rise to perturbations in the electron number density which can be a factor of five larger than the perturbations in baryon density. We do a thorough analysis of the second order anisotropies generated in the cosmic microwave background (CMB) due to perturbations in the electron number density. We show that solving the second order Boltzmann equation for photons is equivalent to solving the first + second order Boltzmann equations and then taking the second order part of the solution. We find the approximate solution to the photon Boltzmann hierarchy in l modes and show that the contributions from inhomogeneous recombination to the second order monopole, dipole and quadrupole are numerically small. We also point out that perturbing the electron number density in the first order tight coupling and damping solutions for the monopole, dipole and quadrupole is not equivalent to solving the second order Boltzmann equations for inhomogeneous recombination. Finally we confirm our result in a previous paper that inhomogeneous recombination gives rise to a local type non-Gaussianity parameter f_{NL}~ -1. The signal to noise for the detection of the temperature bispectrum generated by inhomogeneous recombination is ~ 1 for an ideal full sky experiment measuring modes up to l_{max}=2500.

We calculate the cosmic microwave background (CMB) bispectrum due to inhomogeneous reionization. We calculate all the terms that can contribute to the bispectrum that are products of first order terms on all scales in conformal Newtonian gauge. We also correctly account for the de-correlation between the matter density and initial conditions using perturbation theory up to third order. We find that the bispectrum is of local type as expected. For a reasonable model of reionization, in which the Universe is completely ionized by redshift z_{ri} ~ 8 with optical depth to the last scattering surface \tau_0=0.087 the signal to noise for detection of the CMB temperature bispectrum is S/N ~ 0.1 and confusion in the estimation of primordial non-Gaussianity is f_{NL} ~ -0.1. For an extreme model with z_{ri} ~ 12.5, \tau_0=0.14 we get S/N ~ 0.5 and f_{NL} ~ -0.2.

New low frequency radio telescopes currently being built open up the possibility of observing the 21-cm radiation before the Epoch of Reionization in the future, in particular at redshifts 200 > z > 30, also known as the dark ages. At these high redshifts, Cosmic Microwave Background (CMB) radiation is absorbed by neutral hydrogen at its 21-cm hyperfine transition. This redshifted 21-cm signal thus carries information about the state of the early Universe and can be used to test fundamental physics. We study the constraints these observations can put on the variation of fundamental constants. We show that the 21-cm radiation is very sensitive to the variations in the fine structure constant and can in principle place constraints comparable to or better than the other astrophysical experiments (fractional change < 10^ {-5}). Making such observations will require radio telescopes of collecting area 10 - 10^6 sq. km compared to 1 sq. km of current telescopes. These observations will thus provide independent constraints on the fine structure constant at high redshifts, observations of quasars being the only alternative. More importantly the 21-cm absorption of CMB is the only way to probe the redshift range between recombination and reionization.

The perturbations in the electron number density during recombination contributes to the Cosmic Microwave Background bispectrum through second order terms. Perturbations in the electron density can be a factor of ~5 larger than the baryon density fluctuations on large scales as shown in the calculations by Novosyadlyj. This raises the possibility that the contribution to bispectrum arising from perturbations in the optical depth may be non-negligible. We calculate this bispectrum and find it to peak for squeezed triangles and of peak amplitude of the order of primordial non-Gaussianity of local type with fNL of 0.05 ~ -1 depending on the l-modes being considered. This is because the shape of the bispectrum is different from the primordial one although it peaks for squeezed configurations, similar to the local type primordial non-Gaussianity.

We calculate the contribution of cosmic strings arising from a phase transition in the early universe, or cosmic superstrings arising from brane inflation, to the cosmic 21 cm power spectrum at redshifts z > 30. Future experiments can exploit this effect to constrain the cosmic string tension Gu and probe virtually the entire brane inflation model space allowed by current observations. Although current experiments with a collecting area of ~ 1 km^2 will not provide any useful constraints, future experiments with a collecting area of 10^4-10^6 km^2 covering the cleanest 10% of the sky can in principle constrain cosmic strings with tension Gu > 10^(-10) to 10^(-12) (superstring/phase transition mass scale >10^13 GeV).

We investigate the effect of variation in the value of the fine structure constant at high redshifts (recombination > z > 30) on the absorption of the cosmic microwave background (CMB) at 21 cm hyperfine transition of the neutral atomic hydrogen. We find that the 21 cm signal is very sensitive to the variations in the fine structure constant and it is so far the only probe of the fine structure constant in this redshift range. A change in the value of the fine structure constant by 1% changes the mean brightness temperature decrement of the CMB due to 21 cm absorption by > 5% over the redshift range z < 50. There is an effect of similar magnitude on the amplitude of the fluctuations in the brightness temperature. The redshift of maximum absorption also changes by ~ 5%.