GalSim is a collaborative, open-source project aimed at providing an image simulation tool of enduring benefit to the astronomical community. It provides a software library for generating images of ...astronomical objects such as stars and galaxies in a variety of ways, efficiently handling image transformations and operations such as convolution and rendering at high precision. We describe the GalSim software and its capabilities, including necessary theoretical background. We demonstrate that the performance of GalSim meets the stringent requirements of high precision image analysis applications such as weak gravitational lensing, for current datasets and for the Stage IV dark energy surveys of the Large Synoptic Survey Telescope, ESA’s Euclid mission, and NASA’s WFIRST-AFTA mission. The GalSim project repository is public and includes the full code history, all open and closed issues, installation instructions, documentation, and wiki pages (including a Frequently Asked Questions section). The GalSim repository can be found at https://github.com/GalSim-developers/GalSim.
Cosmological parameter estimation is entering a new era. Large collaborations need to coordinate high-stakes analyses using multiple methods; furthermore such analyses have grown in complexity due to ...sophisticated models of cosmology and systematic uncertainties. In this paper we argue that modularity is the key to addressing these challenges: calculations should be broken up into interchangeable modular units with inputs and outputs clearly defined. We present a new framework for cosmological parameter estimation, CosmoSIS, designed to connect together, share, and advance development of inference tools across the community. We describe the modules already available in CosmoSIS, including camb, Planck, cosmic shear calculations, and a suite of samplers. We illustrate it using demonstration code that you can run out-of-the-box with the installer available at http://bitbucket.org/joezuntz/cosmosis.
We present a sample of galaxies with the Dark Energy Survey (DES) photometry that replicates the properties of the BOSS CMASS sample. The CMASS galaxy sample has been well characterized by the Sloan ...Digital Sky Survey (SDSS) collaboration and was used to obtain the most powerful redshift-space galaxy clustering measurements to date. A joint analysis of redshift-space distortions (such as those probed by CMASS from SDSS) and a galaxy-galaxy lensing measurement for an equivalent sample from DES can provide powerful cosmological constraints. Unfortunately, the DES and SDSS-BOSS footprints have only minimal overlap, primarily on the celestial equator near the SDSS Stripe 82 region. Using this overlap, we build a robust Bayesian model to select CMASS-like galaxies in the remainder of the DES footprint. The newly defined DES-CMASS (DMASS) sample consists of 117,293 effective galaxies covering $1,244 {\rm deg}^2$. Through various validation tests, we show that the DMASS sample selected by this model matches well with the BOSS CMASS sample, specifically in the South Galactic cap (SGC) region that includes Stripe 82. Combining measurements of the angular correlation function and the clustering-z distribution of DMASS, we constrain the difference in mean galaxy bias and mean redshift between the BOSS CMASS and DMASS samples to be $\Delta b = 0.010^{+0.045}_{-0.052}$ and $\Delta z = \left( 3.46^{+5.48}_{-5.55} \right) \times 10^{-3}$ for the SGC portion of CMASS, and $\Delta b = 0.044^{+0.044}_{-0.043} $ and $\Delta z= ( 3.51^{+4.93}_{-5.91}) \times 10^{-3}$ for the full CMASS sample. These values indicate that the mean bias of galaxies and mean redshift in the DMASS sample is consistent with both CMASS samples within $1\sigma$.
We use data from the first-year observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature. We test our methodology in a sample of 1800 ...DES Y1-like mock catalogues. We use the pseudo-Cℓ method to estimate the APS and the mock catalogues to estimate its covariance matrix. We use templates to model the measured spectra and estimate template parameters firstly from the Cℓ’s of the mocks using two different methods, a maximum likelihood estimator and a Markov Chain Monte Carlo, finding consistent results with a good reduced χ2. Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. Finally, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalogue comprises galaxies within an effective area of 1318 deg2 and 0.6 < z < 1.0. We find that the DES Y1 data favour a model with BAO at the 2.6σ C.L. However, the goodness of fit is somewhat poor, with χ2/(d.o.f.) = 1.49. We identify a possible cause showing that using a theoretical covariance matrix obtained from Cℓ’s that are better adjusted to data results in an improved value of χ2/(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Lastly, our results correspond to a distance measurement of DA(zeff = 0.81)/rd = 10.65 ± 0.49, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.
ABSTRACT
We present and characterize the galaxy shape catalogue from the first 3 yr of Dark Energy Survey (DES) observations, over an effective area of 4143 deg2 of the southern sky. We describe our ...data analysis process and our self-calibrating shear measurement pipeline metacalibration, which builds and improves upon the pipeline used in the DES Year 1 analysis in several aspects. The DES Year 3 weak-lensing shape catalogue consists of 100 204 026 galaxies, measured in the riz bands, resulting in a weighted source number density of neff = 5.59 gal arcmin−2 and corresponding shape noise σe = 0.261. We perform a battery of internal null tests on the catalogue, including tests on systematics related to the point spread function (PSF) modelling, spurious catalogue B-mode signals, catalogue contamination, and galaxy properties.
Cosmic shear is a powerful probe of cosmological parameters, but its potential can be fully utilized only if galaxy shapes are measured with great accuracy. Two major effects have been identified ...which are likely to account for most of the bias seen for maximum likelihood methods in recent shear measurement challenges. Model bias occurs when the true galaxy shape is not well represented by the fitted model. Noise bias occurs due to the non-linear relationship between image pixels and galaxy shape. In this paper we investigate the potential interplay between these two effects when an imperfect model is used in the presence of high noise. We present analytical expressions for this bias, which depends on the residual difference between the model and real data. They can lead to biases not accounted for in previous calibration schemes. By measuring the model bias, noise bias and their interaction, we provide a complete statistical framework for measuring galaxy shapes with model fitting methods from GRavitational lEnsing Accuracy Testing (GREAT)-like images. We demonstrate the noise and model interaction bias using a simple toy model, which indicates that this effect can potentially be significant. Using real galaxy images from the Cosmological Evolution Survey (COSMOS) we quantify the strength of the model bias, noise bias and their interaction. We find that the interaction term is often of a similar size to the model bias term, and is smaller than the requirements of current and near future galaxy surveys.
We measure the weak lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey (DES). This pathfinder study is ...meant to (1) validate the Dark Energy Camera (DECam) imager for the task of measuring weak lensing shapes, and (2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, point spread function (PSF) modelling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well behaved, but the modelling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting Navarro-Frenk-White profiles to the clusters in this study, we determine weak lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1. (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.
We measure the weak-lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey. This pathfinder study is meant to ...1) validate the DECam imager for the task of measuring weak-lensing shapes, and 2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, PSF modeling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well-behaved, but the modeling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting NFW profiles to the clusters in this study, we determine weak-lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak-lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.
The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43 GHz and 94 GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization ...of the cosmic microwave background (CMB). QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, over 10,000 hr of data were collected, first with the 19 element 43 GHz array (3458 hr) and then with the 90 element 94 GHz array. Each array observes the same four fields, selected for low foregrounds, together covering Almost-Equal-To 1000 deg{sup 2}. This paper reports initial results from the 43 GHz receiver, which has an array sensitivity to CMB fluctuations of 69 {mu}K{radical}s. The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross-correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB, and EB power spectra in the multipole range l = 25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3{sigma} significance, the E-mode spectrum is consistent with the {Lambda}CDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r = 0.35{sup +1.06}{sub -0.87}. The combination of a new time-stream 'double-demodulation' technique, side-fed Dragonian optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r = 0.1.