Cosmological Signal Research Articles

Sparse component separation for accurate cosmic microwave background estimation

The Cosmological Microwave Background (CMB) is of premier importance for the cosmologists to study the birth of our universe. Unfortunately, most CMB experiments such as COBE, WMAP or Planck do not provide a direct measure of the cosmological signal; CMB is mixed up with galactic foregrounds and point sources. For the sake of scientific exploitation, measuring the CMB requires extracting several different astrophysical components (CMB, Sunyaev-Zel'dovich clusters, galactic dust) form multi-wavelength observations. Mathematically speaking, the problem of disentangling the CMB map from the galactic foregrounds amounts to a component or source separation problem. In the field of CMB studies, a very large range of source separation methods have been applied which all differ from each other in the way they model the data and the criteria they rely on to separate components. Two main difficulties are i) the instrument's beam varies across frequencies and ii) the emission laws of most astrophysical components vary across pixels. This paper aims at introducing a very accurate modeling of CMB data, based on sparsity, accounting for beams variability across frequencies as well as spatial variations of the components' spectral characteristics. Based on this new sparse modeling of the data, a sparsity-based component separation method coined Local-Generalized Morphological Component Analysis (L-GMCA) is described. Extensive numerical experiments have been carried out with simulated Planck data. These experiments show the high efficiency of the proposed component separation methods to estimate a clean CMB map with a very low foreground contamination, which makes L-GMCA of prime interest for CMB studies.

On the time arrows, and randomness in cosmological signals

Arrows of time - thermodynamical, cosmological, electromagnetic, quantum mechanical, psychological - are basic properties of Nature. For a quantum system-bath closed system the de-correlated initial conditions and no-memory (Markovian) dynamics are outlined as necessary conditions for the appearance of the thermodynamical arrow. The emergence of the arrow for the system evolving according to non-unitary dynamics due to the presence of the bath, then, is a result of limited observability, and we conjecture the arrow in the observable Universe as determined by the dark sector acting as a bath. The voids in the large scale matter distribution induce hyperbolicity of the null geodesics, with possible observational consequences.

Background reionization history from omniscopes

The measurements of the 21-cm brightness temperature fluctuations from the neutral hydrogen at the Epoch of Reionization should inaugurate the next generation of cosmological observables. In this respect, many works have concentrated on the disambiguation of the cosmological signals from the dominant reionization foregrounds. However, even after perfect foregrounds removal, our ignorance on the background reionization history can significantly affect the cosmological parameter estimation. In particular, the interdependence between the hydrogen ionized fraction, the baryon density and the optical depth to the redshift of observation induce nontrivial degeneracies between the cosmological parameters that have not been considered so far. Using a simple but consistent reionization model, we revisit their expected constraints for a futuristic giant 21-cm omniscope by using for the first time Markov Chain Monte Carlo methods on multiredshift full sky simulated data. Our results agree well with the usual Fisher matrix analysis on the three-dimensional flat sky power spectrum but only when the above-mentioned degeneracies are kept under control. In the opposite situation, Fisher results can be inaccurate. We show that these conditions can be fulfilled by combining cosmic microwave background measurements with multiple observation redshifts probing the beginning of the Epoch of Reionization. This allows a precise reconstruction of the total optical depth, reionization duration and maximal spin temperature. Finally, we discuss the robustness of these results in presence of unresolved ionizing sources. Although most of the standard cosmological parameters remain weakly affected, we find a significant degradation of the background reionization parameter estimation in presence of nuisance ionizing sources.

Linear perturbations in Galileon gravity models

We study the cosmology of Galileon modified gravity models in the linear perturbation regime. We derive the fully covariant and gauge invariant perturbed field equations using two different methods, which give consistent results, and solve them using a modified version of the {\tt CAMB} code. We find that, in addition to modifying the background expansion history and therefore shifting the positions of the acoustic peaks in the cosmic microwave background (CMB) power spectrum, the Galileon field can cluster strongly from early times, and causes the Weyl gravitational potential to grow, rather than decay, at late times. This leaves clear signatures in the low-$l$ CMB power spectrum through the modified integrated Sachs-Wolfe effect, strongly enhances the linear growth of matter density perturbations and makes distinctive predictions for other cosmological signals such as weak lensing and the power spectrum of density fluctuations. The quasi-static approximation is shown to work quite well from small to the near-horizon scales. We demonstrate that Galileon models display a rich phenomenology due to the large parameter space and the sensitive dependence of the model predictions on the Galileon parameters. Our results show that some Galileon models are already ruled out by present data and that future higher significance galaxy clustering, ISW and lensing measurements will place strong constraints on Galileon gravity.

MEASURING THE REDSHIFT DEPENDENCE OF THE COSMIC MICROWAVE BACKGROUND MONOPOLE TEMPERATURE WITH PLANCK DATA

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.

Imaging neutral hydrogen on large scales during the Epoch of Reionization with LOFAR

The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.

MEASURING BARYON ACOUSTIC OSCILLATIONS ON 21 cm INTENSITY FLUCTUATIONS AT MODERATE REDSHIFTS

After reionization, emission in the 21 cm hyperfine transition provides a direct probe of neutral hydrogen distributed in galaxies. Different from galaxy redshift surveys, observation of baryon acoustic oscillations in the cumulative 21 cm emission may offer an attractive method for constraining dark energy properties at moderate redshifts. Keys to this program are techniques to extract the faint cosmological signal from various contaminants, such as detector noise and continuum foregrounds. In this paper, we investigate the possible systematic and statistical errors in the acoustic scale estimates using ground-based radio interferometers. Based on the simulated 21 cm interferometric measurements, we analyze the performance of a Fourier-space, light-of-sight algorithm in subtracting foregrounds, and further study the observing strategy as a function of instrumental configurations. Measurement uncertainties are presented from a suite of simulations with a variety of parameters, in order to have an estimate of what behaviors will be accessible in the future generation of hydrogen surveys. We find that 10 separate interferometers, each of which contains $\sim 300$ dishes, observes an independent patch of the sky and produces an instantaneous field-of-view of $\sim 100$ $\rm deg^2$, can be used to make a significant detection of acoustic features over a period of a few years. Compared to optical surveys, the broad bandwidth, wide field-of-view and multi-beam observation are all unprecedented capabilities of low-frequency radio experiments.

The potential for very high-frequency gravitational wave detection

The science case for observing gravitational waves at frequencies in the millihertz–kilohertz range using LIGO, VIRGO, GEO600 or LISA is very strong and the first results are expected at these frequencies. However, as gravitational wave astronomy progresses beyond the first detections, other frequency bands may be worth exploring. Early predictions of gravitational wave emission from discrete sources at very much higher frequencies (megahertz and above) have been published and more recent studies of cosmological signals from inflation, Kaluza–Klein modes from gravitational interactions in brane worlds and plasma instabilities surrounding violent astrophysical events, are all possible sources. This communication examines current observational possibilities and the detector technology required to make meaningful observations at these frequencies.

STATISTICAL DESCRIPTION OF SYNCHROTRON INTENSITY FLUCTUATIONS: STUDIES OF ASTROPHYSICAL MAGNETIC TURBULENCE

We provide a theoretical description of synchrotron fluctuations arising from magnetic turbulence. We derive an expression that relates the correlation of synchrotron fluctuations for an arbitrary index of relativistic electrons with the correlations of squared fluctuations of magnetic field component perpendicular to the line of sight. The latter correlations we study assuming that the turbulence is axisymmetric. We obtain general relations valid for an arbitrary model of magnetic turbulence and analyse the relations for particular example of magnetic turbulence that is supported by numerical simulations. We predict that the synchrotron intensity fluctuations are anisotropic with larger correlation present along the direction of magnetic field. This anisotropy is dominated by the quadrupole component with the ratio between quadrupole and monopole parts being sensitive to the compressibility of underlying turbulence. Our work opens avenues for quantitative studies of magnetic turbulence in our galaxy and beyond using synchrotron emission. It also outlines the directions of how synchrotron foreground emission can be separated from cosmological signal, i.e. from CMB or highly redshifted HI emission. For the sake of completeness we also provide the expressions for the synchrotron polarization (Stocks parameters and their combinations) for the model of axisymmetric magnetic turbulence.

The impact of high spatial frequency atmospheric distortions on weak-lensing measurements

High precision cosmology with weak gravitational lensing requires a precise measure of the Point Spread Function across the imaging data where the accuracy to which high spatial frequency variation can be modelled is limited by the stellar number density across the field. We analyse dense stellar fields imaged at the Canada-France-Hawaii Telescope to quantify the degree of high spatial frequency variation in ground-based imaging Point Spread Functions and compare our results to models of atmospheric turbulence. The data shows an anisotropic turbulence pattern with an orientation independent of the wind direction and wind speed. We find the amplitude of the high spatial frequencies to decrease with increasing exposure time as $t^{-1/2}$, and find a negligibly small atmospheric contribution to the Point Spread Function ellipticity variation for exposure times $t>180$ seconds. For future surveys analysing shorter exposure data, this anisotropic turbulence will need to be taken into account as the amplitude of the correlated atmospheric distortions becomes comparable to a cosmological lensing signal on scales less than $\sim 10$ arcminutes. This effect could be mitigated, however, by correlating galaxy shear measured on exposures imaged with a time separation greater than 50 seconds, for which we find the spatial turbulence patterns to be uncorrelated.

SUBTRACTING FOREGROUNDS FROM INTERFEROMETRIC MEASUREMENTS OF THE REDSHIFTED 21 cm EMISSION

The ability to subtract foreground contamination from low-frequency observations is crucial to reveal the underlying 21 cm signal. The traditional line-of-sight methods can deal with the removal of diffuse emission and unresolved point sources, but not bright point sources. In this paper, we introduce a foreground cleaning technique in Fourier space, which allows us to handle all such foregrounds simultaneously and thus sidestep any special treatments to bright point sources. Its feasibility is tested with a simulated data cube for the 21 CentiMeter Array experiment. This data cube includes more realistic models for the 21 cm signal, continuum foregrounds, detector noise and frequency-dependent instrumental response. We find that a combination of two weighting schemes can be used to protect the frequency coherence of foregrounds: the uniform weighting in the uv plane and the inverse-variance weighting in the spectral fitting. The visibility spectrum is therefore well approximated by a quartic polynomial along the line of sight. With this method, we demonstrate that the huge foreground contamination can be cleaned out effectively with residuals on the order of \sim10 mK, while the spectrally smooth component of the cosmological signal is also removed, bringing about systematic underestimate in the extracted power spectrum primarily on large scales.

How well can we measure and understand foregrounds with 21-cm experiments?

Before it becomes a sensitive probe of the epoch of reionization, the dark ages and fundamental physics, 21-cm tomography must successfully contend with the issue of foreground contamination. Broad-band foreground sources are expected to be roughly 4 orders of magnitude larger than any cosmological signals, so precise foreground models will be necessary. Such foreground models often contain a large number of parameters, reflecting the complicated physics that governs foreground sources. In this paper, we concentrate on spectral modelling (neglecting, for instance, bright point source removal from spatial maps) and show that 21-cm tomography experiments will likely not be able to measure these parameters without large degeneracies, simply because the foreground spectra are so featureless and generic. However, we show that this is also an advantage, because it means that the foregrounds can be characterized to a high degree of accuracy once a small number of parameters (likely three or four, depending on one's instrumental specifications) are measured. This provides a simple understanding for why 21-cm foreground subtraction schemes are able to remove most of the contaminants by suppressing just a small handful of simple spectral forms. In addition, this suggests that the foreground modelling process should be relatively simple and will likely not be an impediment to the foreground subtraction schemes that are necessary for a successful 21-cm tomography experiment.

Polarized foreground removal at low radio frequencies using rotation measure synthesis: uncovering the signature of hydrogen reionization

Measurement of redshifted 21-cm emission from neutral hydrogen promises to be the most effective method for studying the reionisation history of hydrogen and, indirectly, the first galaxies. These studies will be limited not by raw sensitivity to the signal, but rather, by bright foreground radiation from Galactic and extragalactic radio sources and the Galactic continuum. In addition, leakage due to gain errors and non-ideal feeds conspire to further contaminate low-frequency radio obsevations. This leakage leads to a portion of the complex linear polarisation signal finding its way into Stokes I, and inhibits the detection of the non-polarised cosmological signal from the epoch of reionisation. In this work, we show that rotation measure synthesis can be used to recover the signature of cosmic hydrogen reionisation in the presence of contamination by polarised foregrounds. To achieve this, we apply the rotation measure synthesis technique to the Stokes I component of a synthetic data cube containing Galactic foreground emission, the effect of instrumental polarisation leakage, and redshifted 21-cm emission by neutral hydrogen from the epoch of reionisation. This produces an effective Stokes I Faraday dispersion function for each line of sight, from which instrumental polarisation leakage can be fitted and subtracted. Our results show that it is possible to recover the signature of reionisation in its late stages (z ~ 7) by way of the 21-cm power spectrum, as well as through tomographic imaging of ionised cavities in the intergalactic medium.

Symmetry of the CMB sky as a new test of its statistical isotropy. Non cosmological octupole?

In this article we propose a novel test for statistical anisotropy of the CMB ΔT(n̂ = (θ,ϕ)). The test is based on the fact, that the Galactic foregrounds have a remarkably strong symmetry with respect to their antipodal points S1: = n̂↔−n̂,n̂ = (θ,ϕ) and with respect to the Galactic plane, while the cosmological signal should not be symmetric or asymmetric under these transitions.We have applied the test for the octupole component of the WMAP ILC 7 map, by looking at a3,1 and a3,3, and their ratio to a3,2 both for real and imaginary values. We find abnormal symmetry of the octupole component at the level of 0.58%, compared to Monte Carlo simulations. By using the analysis of the phases of the octupole we found remarkably strong cross-correlations between the phases of the kinematic dipole and the ILC 7 octupole, in full agreement with previous results. We further test the multipole range 2 < l < 100, by investigating the ratio between the l+m = even and l+m = odd parts of power spectra. We compare the results to simulations of a Gaussian random sky, and find significant departure from the statistically isotropic and homogeneous case, for a very broad range of multipoles. We found that for the most prominent peaks of our estimator, the phases of the corresponding harmonics are coherent with phases of the octupole. We believe, our test would be very useful for detections of various types of residuals of the foreground and systematic effects at a very broad range of multipoles 2 ⩽ l ⩽ 1500−3000 for the forthcoming PLANCK CMB map, before any conclusions about primordial non-Gaussianity and statistical anisotropy of the CMB.

SIMPLE FOREGROUND CLEANING ALGORITHM FOR DETECTING PRIMORDIALB-MODE POLARIZATION OF THE COSMIC MICROWAVE BACKGROUND

We reconsider the pixel-based, "template" polarized foreground removal method within the context of a next-generation, low-noise, low-resolution (0.5 degree FWHM) space-borne experiment measuring the cosmological B-mode polarization signal in the cosmic microwave background (CMB). This method was put forward by the Wilkinson Microwave Anisotropy Probe (WMAP) team and further studied by Efstathiou et al. We need at least 3 frequency channels: one is used for extracting the CMB signal, whereas the other two are used to estimate the spatial distribution of the polarized dust and synchrotron emission. No external template maps are used. We extract the tensor-to-scalar ratio (r) from simulated sky maps consisting of CMB, noise (2 micro K arcmin), and a foreground model, and find that, even for the simplest 3-frequency configuration with 60, 100, and 240 GHz, the residual bias in r is as small as Delta r~0.002. This bias is dominated by the residual synchrotron emission due to spatial variations of the synchrotron spectral index. With an extended mask with fsky=0.5, the bias is reduced further down to <0.001.

Extracting cosmological signals from foregrounds in deep mm maps of the sky

The high Galactic latitude sky at millimeter and submm wavelengths contains significant cosmological information about the early Universe (in terms of the cosmic microwave background) but also the process of structure formation in the Universe from the far infrared background produced by early galaxies and the Sunyaev-Zeldovich effect in clusters of galaxies. As the Planck mission will produce full sky maps in this frequency range, deeper maps of selected low-foregrounds patches of the sky can produce complementary and important information. Here we analyze the performance of a balloon-borne survey covering a 10^\circ x 10^\circ patch of the sky with a few arcminute resolution and very high pixel sensitivity. We simulate the different components of the mm/submm sky (i.e., CMB anisotropies, SZ effect, radio and infrared sources, far infrared background, and interstellar dust) using current knowledge about each of them. We then combine them, adding detector noise, to produce detailed simulated observations in four observational bands ranging from 130 to 500 GHz. Finally, we analyze the simulated maps and estimate the performance of the instrument in extracting the relevant information about each of the components. We find that the CMB angular power spectrum is accurately recovered up to l ~ 3000. Using the Sunyaev-Zel'dovich effect, most of the galaxy clusters present in our input map are detected (60% efficiency overall). Our results also show that much stronger constrains can be placed on far infrared background models.

21cmfast: a fast, seminumerical simulation of the high-redshift 21-cm signal

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

Omniscopes: Large area telescope arrays with onlyNlogNcomputational cost

We show that the class of antenna layouts for telescope arrays allowing cheap analysis hardware (with correlator cost scaling as N log N rather than N^2 with the number of antennas N) is encouragingly large, including not only previously discussed rectangular grids but also arbitrary hierarchies of such grids, with arbitrary rotations and shears at each level. We show that all correlations for such a 2D array with an n-level hierarchy can be efficiently computed via a Fast Fourier Transform in not 2 but 2n dimensions. This can allow major correlator cost reductions for science applications requiring exquisite sensitivity at widely separated angular scales, for example 21cm tomography (where short baselines are needed to probe the cosmological signal and long baselines are needed for point source removal), helping enable future 21cm experiments with thousands or millions of cheap dipole-like antennas. Such hierarchical grids combine the angular resolution advantage of traditional array layouts with the cost advantage of a rectangular Fast Fourier Transform Telescope. We also describe an algorithm for how a subclass of hierarchical arrays can efficiently use rotation synthesis to produce global sky maps with minimal noise and a well-characterized synthesized beam.

The cosmological free-free signal from galaxy groups and clusters

Using analytical models and cosmological N-body simulations, we study the free-free radio emission from ionized gas in clusters and groups of galaxies. The results obtained with the simulations are compared with analytical predictions based on the mass function and scaling relations. Earlier works based on analytical models have shown that the average free-free signal from small haloes (galaxies) during and after the reionization time could be detected with future experiments as a distortion of the CMB spectrum at low frequencies ($\nu <$ 5 GHz). We focus on the period after the reionization time (from redshift $z=0$ up to $z=7$) and on haloes that are more massive than in previous works (groups and clusters). We show how the average signal from haloes with $M > 10^{13} h^{-1} M_{\odot}$ is less than 10% the signal from the more abundant and colder smaller mass haloes. However, the individual signal from the massive haloes could be detected with future experiments opening the door for a new window to study the intracluster medium.

HOW TO IDENTIFY AND SEPARATE BRIGHT GALAXY CLUSTERS FROM THE LOW-FREQUENCY RADIO SKY

In this work we simulate the $50-200$ MHz radio sky that is constrained in the field of view ($5^{\circ}$ radius) of the 21 Centimeter Array (21CMA), by carrying out Monte-Carlo simulations to model redshifted cosmological reionization signals and strong contaminating foregrounds, including emissions from our Galaxy, galaxy clusters, and extragalactic point sources. As an improvement of previous works, we consider in detail not only random variations of morphological and spectroscopic parameters within the ranges allowed by multi-band observations, but also evolution of radio halos in galaxy clusters, assuming that relativistic electrons are re-accelerated in the ICM in merger events and lose energy via both synchrotron emission and inverse Compton scattering with CMB photons. By introducing a new approach designed on the basis of independent component analysis (ICA) and wavelet detection algorithm, we prove that, with a cumulative observation of one month with the 21CMA array, about $80\%$ of galaxy clusters with central brightness temperatures of $> 10~{\rm K}$ at 65 MHz can be safely identified and separated from the overwhelmingly bright foreground. We find that the morphological and spectroscopic distortions are extremely small as compared to the input simulated clusters, and the reduced $\chi^2$ of brightness temperature profiles and spectra are controlled to be $\lesssim 0.5$ and $\lesssim 1.3$, respectively. These results robustly indicate that in the near future a sample of dozens of bright galaxy clusters will be disentangled from the foreground in 21CMA observations, the study of which will greatly improve our knowledge about cluster merger rates, electron acceleration mechanisms in cluster radio halos, and magnetic field in the ICM.