BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Michael L. Brown appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan.
This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements.
We present approaches to quickly simulate systematics affecting CMB observations, including the effects of the scanning strategy. Using summary properties of the scan we capture features of full time ordered data (TOD) simulations, allowing maps and power spectra to be generated at much improved speed for a number of systematics - the cases we present experienced speed ups of 3-4 orders of magnitude when implementing the map-based approaches. We demonstrate the effectiveness of the approaches at capturing the salient features of the scan by directly comparing to full TOD simulations - seeing agreement at sub-percent levels of accuracy. We simulate the effects of differential gain, pointing, and ellipticity to show the effectiveness of the approaches, but note that one could extend these techniques to other systematics. We finally show how to apply these fast map-based simulations of systematic effects to a full focal plane showing their ability to incorporate thousands of detectors as seen in modern CMB experiments.
Contamination by polarized foregrounds is one of the biggest challenges for future polarized cosmic microwave background (CMB) surveys and the potential detection of primordial $B$-modes. Future experiments, such as Simons Observatory (SO) and CMB-S4, will aim at very deep observations in relatively small ($f_{\rm sky} \sim 0.1$) areas of the sky. In this work, we investigate the forecasted performance, as a function of the survey field location on the sky, for regions over the full sky, balancing between polarized foreground avoidance and foreground component separation modeling needs. To do this, we simulate observations by an SO-like experiment and measure the error bar on the detection of the tensor-to-scalar ratio, $σ(r)$, with a pipeline that includes a parametric component separation method, the Correlated Component Analysis, and the use of the Fisher information matrix. We forecast the performance over 192 survey areas covering the full sky and also for optimized low-foreground regions. We find that modeling the spectral energy distribution of foregrounds is the most important factor, and any mismatch will result in residuals and bias in the primordial $B$-modes. At these noise levels, $σ(r)$ is not especially sensitive to the level of foreground contamination, provided the survey targets the least-contaminated regions of the sky close to the Galactic poles.
The CMB $B$-mode polarisation signal -- both the primordial gravitational wave signature and the signal sourced by lensing -- is subject to many contaminants from systematic effects. Of particular concern are systematics that result in mixing of signals of different ``spin'', particularly leakage from the much larger spin-0 intensity signal to the spin-2 polarisation signal. We present a general formalism, which can be applied to arbitrary focal plane setups, that characterises signals in terms of their spin. We provide general expressions to describe how spin-coupled signals observed by the detectors manifest at map-level, in the harmonic domain, and in the power spectra, focusing on the polarisation spectra -- the signals of interest for upcoming CMB surveys. We demonstrate the presence of a previously unidentified cross-term between the systematic and the intrinsic sky signal in the power spectrum, which in some cases can be the dominant source of contamination. The formalism is not restricted to intensity to polarisation leakage but provides a complete elucidation of all leakage including polarisation mixing, and applies to both full and partial (masked) sky surveys, thus covering space-based, balloon-borne, and ground-based experiments. Using a pair-differenced setup, we demonstrate the formalism by using it to completely characterise the effects of differential gain and pointing systematics, incorporating both intensity leakage and polarisation mixing. We validate our results with full time ordered data simulations. Finally, we show in an Appendix that an extension of simple binning map-making to include additional spin information is capable of removing spin-coupled systematics during the map-making process.
We investigate whether a Gaussian likelihood is sufficient to obtain accurate parameter constraints from a Euclid-like combined tomographic power spectrum analysis of weak lensing, galaxy clustering and their cross-correlation. Testing its performance on the full sky against the Wishart distribution, which is the exact likelihood under the assumption of Gaussian fields, we find that the Gaussian likelihood returns accurate parameter constraints. This accuracy is robust to the choices made in the likelihood analysis, including the choice of fiducial cosmology, the range of scales included, and the random noise level. We extend our results to the cut sky by evaluating the additional non-Gaussianity of the joint cut-sky likelihood in both its marginal distributions and dependence structure. We find that the cut-sky likelihood is more non-Gaussian than the full-sky likelihood, but at a level insufficient to introduce significant inaccuracy into parameter constraints obtained using the Gaussian likelihood. Our results should not be affected by the assumption of Gaussian fields, as this approximation only becomes inaccurate on small scales, which in turn corresponds to the limit in which any non-Gaussianity of the likelihood becomes negligible. We nevertheless compare against N-body weak lensing simulations and find no evidence of significant additional non-Gaussianity in the likelihood. Our results indicate that a Gaussian likelihood will be sufficient for robust parameter constraints with power spectra from Stage IV weak lensing surveys.
Future CMB experiments will require exquisite control of systematics in order to constrain the $B$-mode polarisation power spectrum. One class of systematics that requires careful study is instrumental systematics. The potential impact of such systematics is most readily understood by considering analysis pipelines based on pair differencing. In this case, any differential gain, pointing or beam ellipticity between the two detectors in a pair can result in intensity leakage into the $B$-mode spectrum, which needs to be controlled to a high precision due to the much greater magnitude of the total intensity signal as compared to the $B$-mode signal. One well known way to suppress such systematics is through careful design of the scan-strategy, in particular making use of any capability to rotate the instrument about its pointing (boresight) direction. Here, we show that the combination of specific choices of such partial boresight rotation angles with redundancies present in the scan strategy is a powerful approach for suppressing systematic effects. This mitigation can be performed in analysis in advance of map-making and, in contrast to other approaches (e.g. deprojection or filtering), results in no signal loss. We demonstrate our approach explicitly with time ordered data simulations relevant to next-generation ground-based CMB experiments, using deep and wide scan strategies appropriate for experiments based in Chile. These simulations show a reduction of multiple orders of magnitude in the spurious $B$-mode signal arising from differential gain and differential pointing systematics.
The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.
We present optical and near-infrared imaging covering a $\sim$1.53 deg$^2$ region in the Super-Cluster Assisted Shear Survey (SuperCLASS) field, which aims to make the first robust weak lensing measurement at radio wavelengths. We derive photometric redshifts for $\approx$176,000 sources down to $i^\prime_{\rm AB}\sim24$ and present photometric redshifts for 1.4 GHz $e$-MERLIN and VLA detected radio sources found in the central 0.26 deg$^{2}$. We compile an initial catalog of 149 radio sources brighter than S$_{1.4}>75$ $μ$Jy and find their photometric redshifts span $0<z_{\rm phot}<4$ with radio luminosities between $10^{21}-10^{25}$ W Hz$^{-1}$, with medians of $\langle z \rangle =0.55$ and $\langle L_{1.4}\rangle =1.9\times10^{23}$ W Hz$^{-1}$ respectively. We find 95% of the \uJy\ radio source sample (141/149) have SEDs best fit by star-forming templates while 5% (8/149) are better fit by AGN. Spectral indices are calculated for sources with radio observations from VLA and GMRT at 325 MHz, with an average spectral slope of $α=0.59\pm0.04$. Using the full photometric redshift catalog we construct a density map at the redshift of the known galaxy clusters, $z=0.20\pm0.08$. Four of the five clusters are prominently detected at $>7 σ$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers.
We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.
We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.
We present a new model of the microwave sky in polarization that can be used to simulate data from CMB polarization experiments. We exploit the most recent results from the Planck satellite to provide an accurate description of the diffuse polarized foreground synchrotron and thermal dust emission. Our model can include the two mentioned foregrounds, and also a constructed template of Anomalous Microwave Emission (AME). Several options for the frequency dependence of the foregrounds can be easily selected, to reflect our uncertainties and to test the impact of different assumptions. Small angular scale features can be added to the foreground templates to simulate high-resolution observations. We present tests of the model outputs to show the excellent agreement with Planck and WMAP data. We determine the range within which the foreground spectral indices can be varied to be consistent with the current data. We also show forecasts for a high-sensitivity, high-resolution full-sky experiment such as the Cosmic ORigin Explorer (COrE). Our model is released as a python script that is quick and easy to use, available at \url{http://www.jb.man.ac.uk/~chervias}.
Galaxy shapes are not randomly oriented, rather they are statistically aligned in a way that can depend on formation environment, history and galaxy type. Studying the alignment of galaxies can therefore deliver important information about the physics of galaxy formation and evolution as well as the growth of structure in the Universe. In this review paper we summarise key measurements of galaxy alignments, divided by galaxy type, scale and environment. We also cover the statistics and formalism necessary to understand the observations in the literature. With the emergence of weak gravitational lensing as a precision probe of cosmology, galaxy alignments have taken on an added importance because they can mimic cosmic shear, the effect of gravitational lensing by large-scale structure on observed galaxy shapes. This makes galaxy alignments, commonly referred to as intrinsic alignments, an important systematic effect in weak lensing studies. We quantify the impact of intrinsic alignments on cosmic shear surveys and finish by reviewing practical mitigation techniques which attempt to remove contamination by intrinsic alignments.
The shapes of galaxies are not randomly oriented on the sky. During the galaxy formation and evolution process, environment has a strong influence, as tidal gravitational fields in the large-scale structure tend to align nearby galaxies. Additionally, events such as galaxy mergers affect the relative alignments of both the shapes and angular momenta of galaxies throughout their history. These "intrinsic galaxy alignments" are known to exist, but are still poorly understood. This review will offer a pedagogical introduction to the current theories that describe intrinsic galaxy alignments, including the apparent difference in intrinsic alignment between early- and late-type galaxies and the latest efforts to model them analytically. It will then describe the ongoing efforts to simulate intrinsic alignments using both N-body and hydrodynamic simulations. Due to the relative youth of this field, there is still much to be done to understand intrinsic galaxy alignments and this review summarises the current state of the field, providing a solid basis for future work.
We investigate the correlations in galaxy shapes between optical and radio wavelengths using archival observations of the COSMOS field. Cross-correlation studies between different wavebands will become increasingly important for precision cosmology as future large surveys may be dominated by systematic rather than statistical errors. In the case of weak lensing, galaxy shapes must be measured to extraordinary accuracy (shear systematics of $< 0.01\%$) in order to achieve good constraints on dark energy parameters. By using shape information from overlapping surveys in optical and radio bands, robustness to systematics may be significantly improved without loss of constraining power. Here we use HST-ACS optical data, VLA radio data, and extensive simulations to investigate both our ability to make precision measurements of source shapes from realistic radio data, and to constrain the intrinsic astrophysical scatter between the shapes of galaxies as measured in the optical and radio wavebands. By producing a new image from the VLA-COSMOS L-band radio visibility data that is well suited to galaxy shape measurements, we are able to extract precise measurements of galaxy position angles. Comparing to corresponding measurements from the HST optical image, we set a lower limit on the intrinsic astrophysical scatter in position angles, between the optical and radio bands, of $σ_α> 0.212π$ radians (or $38.2^{\circ}$) at a $95\%$ confidence level.
We construct forecasts for cosmological parameter constraints from weak gravitational lensing surveys involving the Square Kilometre Array (SKA). Considering matter content, dark energy and modified gravity parameters, we show that the first phase of the SKA (SKA1) can be competitive with other Stage III experiments such as the Dark Energy Survey (DES) and that the full SKA (SKA2) can potentially form tighter constraints than Stage IV optical weak lensing experiments, such as those that will be conducted with LSST, WFIRST-AFTA or Euclid-like facilities. Using weak lensing alone, going from SKA1 to SKA2 represents improvements by factors of $\sim10$ in matter, $\sim10$ in dark energy and $\sim5$ in modified gravity parameters. We also show, for the first time, the powerful result that comparably tight constraints (within $\sim5\%$) for both Stage III and Stage IV experiments, can be gained from cross-correlating shear maps between the optical and radio wavebands, a process which can also eliminate a number of potential sources of systematic errors which can otherwise limit the utility of weak lensing cosmology.
We construct a pipeline for simulating weak lensing cosmology surveys with the Square Kilometre Array (SKA), taking as inputs telescope sensitivity curves; correlated source flux, size and redshift distributions; a simple ionospheric model; source redshift and ellipticity measurement errors. We then use this simulation pipeline to optimise a 2-year weak lensing survey performed with the first deployment of the SKA (SKA1). Our assessments are based on the total signal-to-noise of the recovered shear power spectra, a metric that we find to correlate very well with a standard dark energy figure of merit. We first consider the choice of frequency band, trading off increases in number counts at lower frequencies against poorer resolution; our analysis strongly prefers the higher frequency Band 2 (950-1760 MHz) channel of the SKA-MID telescope to the lower frequency Band 1 (350-1050 MHz). Best results would be obtained by allowing the centre of Band 2 to shift towards lower frequency, around 1.1 GHz. We then move on to consider survey size, finding that an area of 5,000 square degrees is optimal for most SKA1 instrumental configurations. Finally, we forecast the performance of a weak lensing survey with the second deployment of the SKA. The increased survey size (3$π$\,steradian) and sensitivity improves both the signal-to-noise and the dark energy metrics by two orders of magnitude.
We develop and apply the position angle-only shear estimator of Whittaker et al. (2014) to realistic galaxy images. This is done by demonstrating the method on the simulations of the GREAT3 challenge (Mandelbaum et al. 2014b), which include contributions from anisotropic PSFs. We measure the position angles of the galaxies using three distinct methods - the integrated light method, quadrupole moments of surface brightness, and using model-based ellipticity measurements provided by IM3SHAPE. A weighting scheme is adopted to address biases in the position angle measurements which arise in the presence of an anisotropic PSF. Biases on the shear estimates, due to measurement errors on the position angles and correlations between the measurement errors and the true position angles, are corrected for using simulated galaxy images and an iterative procedure. The properties of the simulations are estimated using the deep field images provided as part of the challenge. A method is developed to match the distributions of galaxy fluxes and half-light radii from the deep fields to the corresponding distributions in the field of interest. We recover angle-only shear estimates with a performance close to current well-established model and moments-based methods for all three angle measurement techniques. The Q-values for all three methods are found to be Q~400. The code is freely available online at http://www.jb.man.ac.uk/~mbrown/angle_only_shear/.
Radio continuum surveys have, in the past, been of restricted use in cosmology. Most studies have concentrated on cross-correlations with the cosmic microwave background to detect the integrated Sachs-Wolfe effect, due to the large sky areas that can be surveyed. As we move into the SKA era, radio continuum surveys will have sufficient source density and sky area to play a major role in cosmology on the largest scales. In this chapter we summarise the experiments that can be carried out with the SKA as it is built up through the coming decade. We show that the SKA can play a unique role in constraining the non-Gaussianity parameter to σ(f_NL) ~ 1, and provide a unique handle on the systematics that inhibit weak lensing surveys. The SKA will also provide the necessary data to test the isotropy of the Universe at redshifts of order unity and thus evaluate the robustness of the cosmological principle.Thus, SKA continuum surveys will turn radio observations into a central probe of cosmological research in the coming decades.
We measure the cosmic shear power spectrum on large angular scales by cross-correlating the shapes of ~9 million galaxies measured in the optical SDSS survey with the shapes of ~2.7x10^5 radio galaxies measured by the overlapping VLA-FIRST survey. Our measurements span the multipole range 10 < l < 130, corresponding to angular scales 2deg < θ < 20deg. On these scales, the shear maps from both surveys suffer from significant systematic effects that prohibit a measurement of the shear power spectrum from either survey alone. Conversely we demonstrate that a power spectrum measured by cross-correlating the two surveys is unbiased. We measure an E-mode power spectrum from the data that is inconsistent with zero signal at the 99\% confidence (~2.7σ) level. The odd-parity B-mode signal and the EB cross- correlation are both found to be consistent with zero (within 1σ). These constraints are obtained after a careful error analysis that accounts for uncertainties due to cosmic variance, random galaxy shape noise and shape measurement errors, as well as additional errors associated with the observed large-scale systematic effects in the two surveys. Our constraints are consistent with the expected signal in the concordance cosmological model assuming recent estimates of the cosmological parameters from the Planck satellite, and literature values for the median redshifts of the SDSS and FIRST galaxy populations. The cross-power spectrum approach described in this paper represents a powerful technique for mitigating shear systematics and will be ideal for extracting robust results, with the exquisite control of systematics required, from future cosmic shear surveys with the SKA, LSST, Euclid and WFIRST-AFTA.
Over the past few years two of the largest and highest fidelity experiments conceived have been approved for construction: Euclid is an ESA M-Class mission that will map three-quarters of the extra galactic sky with Hubble Space Telescope resolution optical and NIR imaging, and NIR spectroscopy, its scientific aims (amongst others) are to create a map of the dark Universe and to determine the nature of dark energy. The Square Kilometre Array (SKA) has similar scientific aims (and others) using radio wavelength observations. The two experiments are synergistic in several respects, both through the scientific objectives and through the control of systematic effects. SKA Phase-1 and Euclid will be commissioned on similar timescales offering an exciting opportunity to exploit synergies between these facilities.
We give an overview of complementarity and synergy in cosmology between the Square Kilometre Array and future survey projects in other wavelengths. In the SKA era, precision cosmology will be limited by systematic errors and cosmic variance, rather than statistical errors. However, combining and/or cross-correlating multi-wavelength data, from the SKA to the cosmic microwave background, optical/infrared and X-ray, substantially reduce these limiting factors. In this chapter, we summarize future survey projects and show highlights of complementarity and synergy, which can be very powerful to probe major cosmological problems such as dark energy, modified gravity and primordial non-Gaussianity.
We discuss methods for performing weak lensing using radio observations to recover information about the intrinsic structural properties of the source galaxies. Radio surveys provide unique information that can benefit weak lensing studies, such as HI emission, which may be used to construct galaxy velocity maps, and polarized synchrotron radiation; both of which provide information about the unlensed galaxy and can be used to reduce galaxy shape noise and the contribution of intrinsic alignments. Using a proxy for the intrinsic position angle of an observed galaxy, we develop techniques for cleanly separating weak gravitational lensing signals from intrinsic alignment contamination in forthcoming radio surveys. Random errors on the intrinsic orientation estimates introduce biases into the shear and intrinsic alignment estimates. However, we show that these biases can be corrected for if the error distribution is accurately known. We demonstrate our methods using simulations, where we reconstruct the shear and intrinsic alignment auto and cross-power spectra in three overlapping redshift bins. We find that the intrinsic position angle information can be used to successfully reconstruct both the lensing and intrinsic alignment power spectra with negligible residual bias.
This document was submitted as supporting material to an Engineering Change Proposal (ECP) for the Square Kilometre Array (SKA). This ECP requests gridded visibilities as an extra imaging data product from the SKA, in order to enable bespoke analysis techniques to measure source morphologies to the accuracy necessary for precision cosmology with radio weak lensing. We also discuss the properties of an SKA weak lensing data set and potential overlaps with other cosmology science goals.
Weak gravitational lensing measurements are traditionally made at optical wavelengths where many highly resolved galaxy images are readily available. However, the Square Kilometre Array (SKA) holds great promise for this type of measurement at radio wavelengths owing to its greatly increased sensitivity and resolution over typical radio surveys. The key to successful weak lensing experiments is in measuring the shapes of detected sources to high accuracy. In this document we describe a simulation pipeline designed to simulate radio images of the quality required for weak lensing, and will be typical of SKA observations. We provide as input, images with realistic galaxy shapes which are then simulated to produce images as they would have been observed with a given radio interferometer. We exploit this pipeline to investigate various stages of a weak lensing experiment in order to better understand the effects that may impact shape measurement. We first show how the proposed SKA1-Mid array configurations perform when we compare the (known) input and output ellipticities. We then investigate how making small changes to these array configurations impact on this input-outut ellipticity comparison. We also demonstrate how alternative configurations for SKA1-Mid that are smaller in extent, and with a faster survey speeds produce similar performance to those originally proposed. We then show how a notional SKA configuration performs in the same shape measurement challenge. Finally, we describe ongoing efforts to utilise our simulation pipeline to address questions relating to how applicable current (mostly originating from optical data) shape measurement techniques are to future radio surveys. As an alternative to such image plane techniques, we lastly discuss a shape measurement technique based on the shapelets formalism that reconstructs the source shapes directly from the visibility data.
We apply the Correlated Component Analysis (CCA) method on simulated data of the Square Kilometre Array, with the aim of accurately cleaning the 21 cm reionization signal from diffuse foreground contamination. The CCA has been developed for the Cosmic Microwave Background, but the application of the Fourier-domain implementation of this method to the reionization signal is straightforward. The CCA is a parametric method to estimate the frequency behaviour of the foregrounds from the data by using second-order statistics. We test its performance on foreground simulations of increasing complexity, designed to challenge the parametric models adopted. We also drop the assumption of spectral smoothness that most of the methods rely upon. We are able to clean effectively the simulated data across the explored frequency range (100-200 MHz) for all the foreground simulations. This shows that the CCA method is very promising for EoR component separation.
We develop a method for performing a weak lensing analysis using only measurements of galaxy position angles. By analysing the statistical properties of the galaxy orientations given a known intrinsic ellipticity distribution, we show that it is possible to obtain estimates of the shear by minimizing a $χ^2$ statistic. The method is demonstrated using simulations where the components of the intrinsic ellipticity are taken to be Gaussian distributed. Uncertainties in the position angle measurements introduce a bias into the shear estimates which can be reduced to negligible levels by introducing a correction term into the formalism. We generalize our approach by developing an algorithm to obtain direct shear estimators given any azimuthally symmetric intrinsic ellipticity distribution. We introduce a method of measuring the position angles of the galaxies from noisy pixelized images, and propose a method to correct for biases which arise due to pixelization and correlations between measurement errors and galaxy ellipticities. We also develop a method to constrain the sample of galaxies used to obtain an estimate of the intrinsic ellipticity distribution such that fractional biases in the resulting shear estimates are below a given threshold value. We demonstrate the angle only method by applying it to simulations where the ellipticities are taken to follow a log-normal distribution. We compare the performance of the position angle only method with the standard method based on full ellipticity measurements by reconstructing lensing convergence maps from both numerical simulations and from the CFHTLenS data. We find that the difference between the convergence maps reconstructed using the two methods is consistent with noise.
We present a finely-binned tomographic weak lensing analysis of the Canada-France-Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, with accurate shear and photometric redshifts for a galaxy sample with a median redshift of zm =0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six redshift bins, each spanning the angular range of 1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological probes: the cosmic microwave background with data from WMAP7, baryon acoustic oscillations with data from BOSS, and a prior on the Hubble constant from the HST distance ladder. This leads to constraints on the normalisation of the matter power spectrum sigma_8 = 0.799 +/- 0.015 and the matter density parameter Omega_m = 0.271 +/- 0.010 for a flat Lambda CDM cosmology. For a flat wCDM cosmology we constrain the dark energy equation of state parameter w = -1.02 +/- 0.09. We also provide constraints for curved Lambda CDM and wCDM cosmologies. We find the intrinsic alignment contamination to be galaxy-type dependent with a significant intrinsic alignment signal found for early-type galaxies, in contrast to the late-type galaxy sample for which the intrinsic alignment signal is found to be consistent with zero.
We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.
In a recent paper (Brown & Battye 2011), we proposed the use of integrated polarization measurements of background galaxies in radio weak gravitational lensing surveys and investigated the potential impact on the statistical measurement of cosmic shear. Here we extend this idea to reconstruct maps of the projected dark matter distribution, or lensing convergence field. The addition of polarization can, in principle, greatly reduce shape noise due to the intrinsic dispersion in galaxy ellipticities. We show that maps reconstructed using this technique in the radio band can be competitive with those derived using standard lensing techniques which make use of many more galaxies. In addition, since the reconstruction noise is uncorrelated between these standard techniques and the polarization technique, their comparison can serve as a powerful check for systematics and their combination can reduce noise further. We examine the convergence reconstruction which could be achieved with two forthcoming facilities: (i) a deep survey, covering 1.75 square degrees using the e-MERLIN instrument currently being commissioned in the UK and (ii) the high resolution, deep wide field surveys which will eventually be conducted with the Square Kilometre Array.
We investigate the observational signatures of three models of the early Universe in the $B$-mode polarization of the Cosmic Microwave Background (CMB) radiation. In addition to the standard single field inflationary model, we also consider the constraints obtainable on the loop quantum cosmology model (from Loop Quantum Gravity) and on cosmic strings, expected to be copiously produced during the latter stages of Brane inflation. We first examine the observational features of the three models, and then use current $B$-mode polarization data from the BICEP and QUaD experiments to constrain their parameters. We also examine the detectability of the primordial $B$-mode signal predicted by these models and forecast the parameter constraints achievable with future CMB polarization experiments. We find that: (a) these three models of the early Universe predict different features in the CMB B-mode polarization power spectrum, which are potentially distinguishable from the CMB experiments; (b) since $B$-mode polarization measurements are mostly unaffected by parameter degeneracies, they provide the cleanest probe of these early Universe models; (c) using the BICEP and QUaD data we obtain the following parameter constraints: $r=0.02^{+0.31}_{-0.26}$ ($1σ$ for the tensor-to-scalar ratio in the single field inflationary model); $m < 1.36\times 10^{-8} \text{M}_{\text{pl}}$ and $k_{*} < 2.43 \times 10^{-4} \text{Mpc}^{-1}$ ($1σ$ for the mass and scale parameters in the loop quantum cosmology model); and $Gμ< 5.77 \times 10^{-7}$ ($1σ$ for the cosmic string tension); (d) future CMB observations (both satellite missions and forthcoming sub-orbital experiments) will provide much more rigorous tests of these early Universe models.
We propose a new technique for weak gravitational lensing in the radio band making use of polarization information. Since the orientation of a galaxy's polarized emission is both unaffected by lensing and is related to the galaxy's intrinsic orientation, it effectively provides information on the unlensed galaxy position angle. We derive a new weak lensing estimator which exploits this effect and makes full use of both the observed galaxy shapes and the estimates of the intrinsic position angles as provided by polarization. Our method has the potential to both reduce the effects of shot noise, and to reduce to negligible levels, in a model-independent way, all effects of intrinsic galaxy alignments. We test our technique on simulated weak lensing skies, including an intrinsic alignment contaminant consistent with recent observations, in three overlapping redshift bins. Adopting a standard weak lensing analysis and ignoring intrinsic alignments results in biases of 5-10% in the recovered power spectra and cosmological parameters. Applying our new estimator to one tenth the number of galaxies used for the standard case, we recover both power spectra and the input cosmology with similar precision as compared to the standard case and with negligible residual bias, even in the presence of a substantial (astrophysical) scatter in the relationship between the observed orientation of the polarized emission and the intrinsic orientation. Assuming a reasonable polarization fraction for star-forming galaxies, and no cosmological conspiracy in the relationship between polarization direction and intrinsic morphology, our estimator should prove a valuable tool for weak lensing analyses of forthcoming radio surveys, in particular, deep wide field surveys with e-MERLIN, MeerKAT and ASKAP and ultimately, definitive radio lensing surveys with the SKA.
We investigate the impact of both slow and fast polarization modulation strategies on the science return of upcoming ground-based experiments aimed at measuring the B-mode polarization of the CMB. Using simulations of the Clover experiment, we compare the ability of modulated and un-modulated observations to recover the signature of gravitational waves in the polarized CMB sky in the presence of a number of anticipated systematic effects. The general expectations that fast modulation is helpful in mitigating low-frequency detector noise, and that the additional redundancy in the projection of the instrument's polarization sensitivity directions onto the sky when modulating reduces the impact of instrumental polarization, are borne out by our simulations. Neither low-frequency polarized atmospheric fluctuations nor systematic errors in the polarization sensitivity directions are mitigated by modulation. Additionally, we find no significant reduction in the effect of pointing errors by modulation. For a Clover-like experiment, pointing jitter should be negligible but any systematic mis-calibration of the polarization coordinate reference system results in significant E-B mixing on all angular scales and will require careful control. We also stress the importance of combining data from multiple detectors in order to remove the effects of common-mode systematics (such as 1/f atmospheric noise) on the measured polarization signal. Finally we compare the performance of our simulated experiment with the predicted performance from a Fisher analysis. We find good agreement between the Fisher predictions and the simulations except for the very largest scales where the power spectrum estimator we have used introduces additional variance to the B-mode signal recovered from our simulations.