Full-sky maps for gravitational lensing of the cosmic microwave background

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the angle, aiming at gravitational lensing of the cosmic microwave background (CMB), with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to z ∼ 11, producing continuous all-sky maps with arcmin angular resolution. A randomization scheme avoids the repetition of structures along the line of sight, and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcmin, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analysing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees–Sciama effect.     

Observational tests of x-matter models

We study gravitational lensing statistics, matter power spectra and the angular power spectra of the cosmic microwave background (CMB) radiation in x-matter models. We adopt an equation of state of x-matter which can express a wide range of matter from pressureless dust to the cosmological constant. A new ingredient in this model is the sound speed of the x-component, in addition to the equation of state w ₀ =  p _x0/ρ_x0. Except for the cosmological constant case, the perturbations of x-matter itself are considered. Our primary interest is in the effect of non-zero sound speed on the structure formation and the CMB spectra. It is found that there exist parameter ranges where x-matter models are consistent with all current observations. The x-matter generally leaves imprints in the CMB anisotropy and the matter power spectrum, which should be detectable in future observations.     

The Australia Telescope search for cosmic microwave background anisotropy

In an attempt to detect cosmic microwave background (CMB) anisotropy on arcmin scales, we have made an 8.7-GHz image of a sky region with a resolution of 2 arcmin and high surface brightness sensitivity using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration. The foreground discrete-source confusion was estimated from observations with higher resolution at the same frequency and in a scaled array at a lower frequency. Following the subtraction of the foreground confusion, the field shows no features in excess of the instrument noise. This limits the CMB anisotropy flat-band power to Q _flat < 23.6 μ K with 95 per cent confidence; the ATCA filter function (which is available at the website www.atnf.csiro.au/Research/cmbr/cmbr_atca.html) F _l in multipole l -space peaks at l _eff = 4700 and has half-maximum values at l  = 3350 and 6050.     

Testing tidal-torque theory – I. Spin amplitude and direction

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N -body simulation with thousands of well-resolved haloes. The protohaloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the protohaloes match, on average, the spin amplitudes of the virialized haloes of today, if linear growth is assumed until ∼ t ₀/3, or  55–70  per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ∼50° in the TTT prediction at t ₀, such that the linear spatial correlations of spins on scales ≥1  h ⁻¹ Mpc are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second-order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.     

Polarization of synchrotron emission from relativistic reconfinement shocks

We study the polarization properties of relativistic reconfinement shocks with chaotic magnetic fields. Using our hydrodynamical model of their structure, we calculate synthetic polarization maps, longitudinal polarization profiles and discuss the spatially averaged polarization degree as a function of jet half-opening angle  Θ _j   , jet Lorentz factor  Γ _j   and observer inclination angle to the jet axis  θ_obs  . We find that for  θ_obs≲Θ _j   the wave electric vectors are parallel in the vicinity of the structure ends and perpendicular in between, while for  θ_obs > Θ _j   the polarization can only be perpendicular. The spatially averaged polarization degree does not exceed 30 per cent. Parallel average polarization, with polarization degrees lower than 10 per cent, has been found for  θ_obs < Θ _j   under the condition  Γ _j Θ _j > 1  . As earlier works predicted the parallel polarization from relativistic conical shocks, we explain our results by discussing conical shocks with divergent upstream flow.     

Weak lensing of the CMB: extraction of lensing information from the trispectrum

We discuss the four-point correlation function, or the trispectrum in Fourier space, of CMB temperature and polarization anisotropies due to the weak gravitational lensing effect by intervening large scale structure. We discuss the squared temperature power spectrum as a probe of this trispectrum and, more importantly, as an observational approach to extracting the power spectrum of the deflection angle associated with the weak gravitational lensing effect on the CMB. We extend previous discussions on the trispectrum and associated weak lensing reconstruction from CMB data by calculating non-Gaussian noise contributions, beyond the previously discussed dominant Gaussian noise. Non-Gaussian noise contributions are generated by lensing itself and by the correlation between the lensing effect and other foreground secondary anisotropies in the CMB such as the Sunyaev–Zel’dovich (SZ) effect. When the SZ effect is removed from temperature maps using its spectral dependence, we find these additional non-Gaussian noise contributions to be an order of magnitude lower than the dominant Gaussian noise. If the noise-bias due to the dominant Gaussian part of the temperature squared power spectrum is removed, then these additional non-Gaussian contributions provide the limiting noise level for the lensing reconstruction. The temperature squared power spectrum allows a high signal-to-noise extraction of the lensing deflections and a confusion-free separation of the curl (or B-mode) polarization due to inflationary gravitational waves from that due to lensed gradient (or E-mode) polarization. The small angular scale temperature and polarization anisotropy measurements provide a novel approach to weak lensing studies, complementing the approach based on galaxy ellipticities.     

lemomaf: Lensed Mock Map Facility

We present the Lensed Mock Map Facility ( lemomaf ), a tool designed to perform mock weak-lensing measurements on numerically simulated chunks of the Universe. Coupling N -body simulations to a semi-analytical model of galaxy formation, lemomaf can create realistic lensed images and mock catalogues of galaxies, at wavelengths ranging from the ultraviolet to the submillimetre. To demonstrate the power of such a tool, we compute predictions of the source–lens clustering (SLC) effect on the convergence statistics, and quantify the impact of weak lensing on galaxy counts in two different filters. We find that the SLC effect skews the probability density function of the convergence towards low values, with an intensity which strongly depends on the redshift distribution of galaxies. On the other hand, the degree of enhancement or depletion in galaxy counts due to weak lensing is independent of the SLC effect. We discuss the impact on the two-point shear statistics to be measured by future missions like SNAP and LSST . The SLC effect would bias the estimation of σ₈ from two-point statistics up to 5 per cent for a narrow redshift distribution of mean   z ∼ 0.5  , and up to 2 per cent in small angular scales for a redshift distribution of mean   z ∼ 1.5  . We conclude that accurate photometric redshifts for individual galaxies are necessary in order to quantify and isolate the SLC effect.     

Gravitational lensing of anisotropic sources

In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle δ between the light rays emerging from the source, for different lensing configurations. If the lens has a power-law profile   M ∝ R ^γ  , the angle δ initially increases with lens redshift and then either diverges (for a steep profile  γ < 1  ), remains constant (for an isothermal profile  γ= 1  ), or vanishes (for a shallow profile  γ > 1  ) as   z _l→ z _s  . The scaling with lens mass is roughly  δ∝ M ^1/(2−γ)  . The results for an Navarro–Frenk–White (NFW) profile are qualitatively similar to those for a shallow power law, with δ peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle  θ_jet  , there is a probability as high as   P ∼ 0.02–0.07(θ_jet/0.5°)⁻¹  that one of the lensed images may be missed (for  2 ≲ z _s≲ 6  ). Differential absorption within active galactic nuclei (AGNs) could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing.     

Constraints on cosmological parameters from recent measurements of cosmic microwave background anisotropy

A key prediction of cosmological theories for the origin and evolution of structure in the Universe is the existence of a 'Doppler peak' in the angular power spectrum of cosmic microwave background (CMB) fluctuations. We present new results from a study of recent CMB observations which provide the first strong evidence for the existence of a 'Doppler peak' localized in both angular scale and amplitude. This first estimate of the angular position of the peak is used to place a new direct limit on the curvature of the Universe, corresponding to a density of Ω = 0.7^+0.8_−0.5, consistent with a flat universe. Very low-density 'open' universe models are inconsistent with this limit unless there is a significant contribution from a cosmological constant. For a flat standard cold dark matter dominated universe we use our results in conjunction with big bang nucleosynthesis constraints to determine the value of the Hubble constant as H ₀ = 30 − 70 km s⁻¹ Mpc⁻¹ for baryon fractions Ω_b = 0.05 to 0.2. For H ₀ = 50 km s⁻¹ Mpc⁻¹ we find the primordial spectral index of the fluctuations to be n  = 1.1 ± 0.1, in close agreement with the inflationary prediction of n  ≃ 1.0.     

Investigating cosmological weak lensing with the ray-bundle method

Using the ray-bundle method for calculating gravitational lens magnifications, we outline a method by which the magnification probability may be determined specifically in the weak lensing limit for cosmological models obtained from N -body simulations.
16 different models are investigated, which are variations on three broad classes of cold dark matter model: the standard model with  (Ω₀, λ ₀)=(1.0,0.0)  , the open model with  (Ω₀, λ ₀)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω₀, λ ₀)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H ₀, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ ₈, are studied. It is shown that there is no signature of these parameters in the weak lensing magnification distributions. The magnification probability distributions are also shown to be independent of the numerical parameters such as the lens mass and simulation box size in the N -body simulations.     

Sunyaev–Zel'dovich cluster survey simulations for Planck

We examine the ability of the future Planck mission to provide a catalogue of galaxy clusters observed via their Sunyaev–Zel'dovich (SZ) distortion in the cosmic microwave background (CMB). For this purpose we produce full-sky SZ maps based on N -body simulations and scaling relations between cluster properties for several cosmological models. We extrapolate the N -body simulations by a mass function to high redshifts in order to obtain a realistic SZ background. The simulated Planck observations include, besides the thermal and kinematic SZ effects, contributions from the primordial CMB, extragalactic point sources as well as Galactic dust, free–free and synchrotron emission. A harmonic-space maximum-entropy method is used to separate the SZ signal from contaminating components in combination with a cluster detection algorithm based on thresholding and flux integration to identify clusters and to obtain their fluxes. We estimate a survey sensitivity limit (depending on the quality of the recovered cluster flux) and provide cluster survey completeness and purity estimates. We find that, given our modelling and detection algorithm, Planck will reliably detect at least several thousands of clusters over the full sky. The exact number depends on the particular cosmological model (up to 10 000 cluster detections in a concordance ΛCDM model with  σ₈= 0.9  ). We show that the Galaxy does not significantly affect the cluster detection. Furthermore, the dependence of the thermal SZ power spectrum on the matter variance on scales of  8 h ⁻¹  Mpc and the quality of its reconstruction by the employed method are investigated. Our simulations suggest that the Planck cluster sample will not only be useful as a basis for follow-up observations, but also will have the ability to provide constraints on cosmological parameters.     

Modelling the first probable two-plane lens system B2114+022: reproducing two compact radio cores A and D

We test possible lensing scenarios of the JVAS system B2114+022, in which two galaxies at different redshifts ('G1' at z ₁=0.3157 and 'G2' at z ₂=0.5883) are found within 2 arcsec of quadruple radio sources. For our investigation, we use possible lensing constraints derived from a wealth of data on the radio sources obtained with VLA, MERLIN, VLBA and EVN as well as HST imaging data on the two galaxies, which were presented recently in Augusto et al. In the present study, we focus on reproducing the widest separated, observationally similar radio components A and D as lensed images. We first treat G2 (which is the more distant one from the geometric centre) as a shear term, and then consider two-plane lensing explicitly including G2's potential at the z ₂ plane as the first case of two-plane lens modelling. Our modelling results not only support the hypothesis that the system includes gravitationally lensed images of a higher-redshift extragalactic object, but they also show that the explicit inclusion of G2's potential at the second lens plane is necessary in order to fit the data with astrophysically plausible galaxy parameters. Finally, we illustrate a natural consequence of a two-plane lens system, namely the prediction of distortion as well as shift and stretching of G2's isophotes by G1's potential, which can in principle be measured by subtracting out G1's light distribution in an image of high signal-to-noise ratio and good angular resolution, especially a multicolour one.     

Statistics of the Sunyaev–Zel'dovich effect power spectrum

Using large numbers of simulations of the microwave sky, incorporating the cosmic microwave background (CMB) and the Sunyaev–Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10 000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 per cent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager (CBI) and Berkeley–Illinois–Maryland Association (BIMA) experiments. Our results indicate that the observed excess could be explained using a lower value of σ₈ than previously suggested, however the effect is not enough to match  σ₈= 0.825  . The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.     

Evolution of the dark matter distribution with three-dimensional weak lensing

We present a direct detection of the growth of large-scale structure, using weak gravitational lensing and photometric redshift data from the COMBO-17 survey. We use deep R -band imaging of two  0.5 × 0.5 deg²  fields, affording shear estimates for over 52 000 galaxies; we combine these with photometric redshift estimates from our 17-band survey, in order to obtain a 3D shear field. We find theoretical models for evolving matter power spectra and correlation functions, and fit the corresponding shear correlation functions to the data as a function of redshift. We detect the evolution of the power at the 4.7σ level given reasonable priors, and measure the rate of evolution for  0 < z < 1  . We also fit correlation functions to our 3D data as a function of cosmological parameters σ₈ and  Ω_Λ  . We find joint constraints on  Ω_Λ  and σ₈, demonstrating an improvement in accuracy by ≃40 per cent over that available from 2D weak lensing for the same area.     

Fast power spectrum estimation

We present a new algorithm to rapidly and optimally compute power spectra. This new algorithm is based on a generalization of iterative multigrid, and has computational cost     , compared to the standard brute force approach which costs     . The procedure retains this speed on the full sky and for ill-conditioned matrices. It is applicable to galaxy power spectra, cosmic microwave background (CMB), polarization and weak lensing data. We present a mathematical convergence analysis, and performance results.     

Weak gravitational flexion

Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy–galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy–galaxy lensing, examining isothermal sphere and Navarro–Frenk–White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in two-dimensional projection on the sky, and in three dimensions in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence–flexion cross-correlations, and we find that the signal is an effective probe of structure on non-linear scales.     

Statistical distribution of gravitational-lensing excursion angles: winding ways to us from the deep Universe

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, quasi-stellar objects and also the cosmic microwave background temperature map), is known to consist of two components: a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a ΛCDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating these two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive haloes with   M > 10¹⁴  h ⁻¹ M_⊙  at   z < 1  and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arcmin, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at subarcminute scales.     

Weak lensing of the CMB by large-scale structure

Several recent papers have studied lensing of the CMB by large-scale structures, which probes the projected matter distribution from z=10³ to z0. This interest is motivated in part by upcoming high resolution, high sensitivity CMB experiments, such as APEX/SZ, ACT, SPT or Planck, which should be sensitive to lensing. In this paper, we examine the reconstruction of the large-scale dark matter distribution from lensed CMB temperature anisotropies. We go beyond previous work in using numerical simulations to include higher order, non-Gaussian effects and find that the convergence and its power spectrum are biased, with the bias increasing with the angular resolution. We also study the contamination by the kinetic Sunyaev–Zel'dovich signal, which is spectrally indistinguishable from lensed CMB anisotropies, and find that it leads to an overestimate of the convergence. We finish by estimating the sensitivity of the previously cited experiments and find that all of them could detect the lensing effect, but would be biased at around the 10% level.     

Data analysis of continuous gravitational wave: all-sky search and study of templates

We have studied the problem of an all-sky search in reference to a continuous gravitational wave particularly for such sources whose wave-form are known in advance. We have made an analysis of the number of templates required for matched-filter analysis as applicable to these sources. We have employed the concept of fitting factor (FF); treating the source location as the parameters of the signal manifold and have studied the matching of the signal with templates corresponding to different source locations. We have investigated the variation of FF with source location and have noticed a symmetry in template parameters,  θ_T  and  φ_T  . It has been found that the two different template values in source location, each in  θ_T  and  φ_T  , have the same FF. We have also computed the number of templates required assuming the noise power spectral density   S _n ( f )  to be flat. It is observed that higher FF requires an exponentially increasing, large number of templates.     

Variance and skewness in the FIRST survey

We investigate the large-scale clustering of radio sources in the FIRST 1.4-GHz survey by analysing the distribution function ( counts in cells ). We select a reliable sample from the the FIRST catalogue, paying particular attention to the problem of how to define single radio sources from the multiple components listed. We also consider the incompleteness of the catalogue. We estimate the angular two-point correlation function w (θ), the variance Ψ₂ and skewness Ψ₃ of the distribution for the various subsamples chosen on different criteria. Both w (θ) and Ψ₂ show power-law behaviour with an amplitude corresponding to a spatial correlation length of r ₀ ∼ 10  h ⁻¹Mpc. We detect significant skewness in the distribution, the first such detection in radio surveys. This skewness is found to be related to the variance through Ψ₃ =  S ₃(Ψ₂)^α, with α = 1.9 ± 0.1, consistent with the non-linear gravitational growth of perturbations from primordial Gaussian initial conditions. We show that the amplitude of variance and the skewness are consistent with realistic models of galaxy clustering.