We describe S2LET, a fast and robust implementation of the scale-discretised wavelet transform on the sphere. Wavelets are constructed through a tiling of the harmonic line and can be used to probe ...spatially localised, scale-dependent features of signals on the sphere. The reconstruction of a signal from its wavelets coefficients is made exact here through the use of a sampling theorem on the sphere. Moreover, a multiresolution algorithm is presented to capture all information of each wavelet scale in the minimal number of samples on the sphere. In addition S2LET supports the HEALPix pixelisation scheme, in which case the transform is not exact but nevertheless achieves good numerical accuracy. The core routines of S2LET are written in C and have interfaces in Matlab, IDL and Java. Real signals can be written to and read from FITS files and plotted as Mollweide projections. The S2LET code is made publicly available, is extensively documented, and ships with several examples in the four languages supported. At present the code is restricted to axisymmetric wavelets but will be extended to directional, steerable wavelets in a future release.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Abstract
We compare two rotationally invariant decomposition techniques on linear polarization data: the spin-2 spherical harmonic decomposition in two opposite parities, the E- and B-mode, and the ...multiscale analysis of the gradient of linear polarization,
$|\nabla \boldsymbol {P}|$
. We demonstrate that both decompositions have similar properties in the image domain and the spatial frequency domain. They can be used as complementary tools for turbulence analysis of interstellar magnetic fields in order to develop a better understanding of the origin of energy sources for the turbulence, the origin of peculiar magnetic field structures and their underlying physics. We also introduce a new quantity |∇EB| based on the E- and B-modes and we show that in the intermediate- and small-scale limit
$|\nabla EB| \simeq |\nabla \boldsymbol {P}|$
. Analysis of the 2.3 GHz S-band Polarization All Sky Survey shows many extended coherent filament-like features appearing as ‘double jumps’ in the
$|\nabla \boldsymbol {P}|$
map that are correlated with negative and positive filaments of B-type polarization. These local asymmetries between the two polarization types, E and B, of the non-thermal Galactic synchrotron emission have an influence on the E- and B-mode power spectra analyses. The wavelet-based formalism of the polarization gradient analysis allows us to locate the position of E- or B-mode features responsible for the local asymmetries between the two polarization types. In analysed subregions, the perturbations of the magnetic field are trigged by star clusters associated with H ii regions, the Orion–Eridanus superbubble and the North Polar Spur at low Galactic latitude.
It has been claimed recently that massive sterile neutrinos could bring about a new concordance between observations of the cosmic microwave background, the large-scale structure of the Universe, and ...local measurements of the Hubble constant, H(0). We demonstrate that this apparent concordance results from combining data sets which are in significant tension, even within this extended model, possibly indicating remaining systematic biases in the measurements. We further show that this tension remains when the cosmological model is further extended to include significant tensor modes, as suggested by the recent BICEP2 results. Using the Bayesian evidence, we show that the cold dark matter model with a cosmological constant is strongly favored over its neutrino extensions by various combinations of data sets. Robust data combinations yield stringent limits of ∑m(ν) ≲ 0.3 eV and m(ν,sterile)(eff) ≲ 0.3 eV at 95% C.L. for the sum of active and sterile neutrinos, respectively.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
Exact Wavelets on the Ball Leistedt, B.; McEwen, J. D.
IEEE transactions on signal processing,
12/2012, Volume:
60, Issue:
12
Journal Article
Peer reviewed
Open access
We develop an exact wavelet transform on the three-dimensional ball (i.e. on the solid sphere), which we name the flaglet transform. For this purpose we first construct an exact transform on the ...radial half-line using damped Laguerre polynomials and develop a corresponding quadrature rule. Combined with the spherical harmonic transform, this approach leads to a sampling theorem on the ball and a novel three-dimensional decomposition which we call the Fourier-Laguerre transform. We relate this new transform to the well-known Fourier-Bessel decomposition and show that band-limitedness in the Fourier-Laguerre basis is a sufficient condition to compute the Fourier-Bessel decomposition exactly. We then construct the flaglet transform on the ball through a harmonic tiling, which is exact thanks to the exactness of the Fourier-Laguerre transform (from which the name flaglets is coined). The corresponding wavelet kernels are well localised in real and Fourier-Laguerre spaces and their angular aperture is invariant under radial translation. We introduce a multiresolution algorithm to perform the flaglet transform rapidly, while capturing all information at each wavelet scale in the minimal number of samples on the ball. Our implementation of these new tools achieves floating-point precision and is made publicly available. We perform numerical experiments demonstrating the speed and accuracy of these libraries and illustrate their capabilities on a simple denoising example.
Abstract
We describe the Dark Energy Survey (DES) photometric data set assembled from the first three years of science operations to support DES Year 3 cosmologic analyses, and provide usage notes ...aimed at the broad astrophysics community.
Y3
GOLD
improves on previous releases from DES,
Y1
GOLD
, and Data Release 1 (DES DR1), presenting an expanded and curated data set that incorporates algorithmic developments in image detrending and processing, photometric calibration, and object classification.
Y3
GOLD
comprises nearly 5000 deg
2
of
grizY
imaging in the south Galactic cap, including nearly 390 million objects, with depth reaching a signal-to-noise ratio ∼10 for extended objects up to
i
AB
∼ 23.0, and top-of-the-atmosphere photometric uniformity <3 mmag. Compared to DR1, photometric residuals with respect to Gaia are reduced by 50%, and per-object chromatic corrections are introduced.
Y3
GOLD
augments DES DR1 with simultaneous fits to multi-epoch photometry for more robust galactic color measurements and corresponding photometric redshift estimates.
Y3
GOLD
features improved morphological star–galaxy classification with efficiency >98% and purity >99% for galaxies with 19 <
i
AB
< 22.5. Additionally, it includes per-object quality information, and accompanying maps of the footprint coverage, masked regions, imaging depth, survey conditions, and astrophysical foregrounds that are used to select the cosmologic analysis samples.
We describe the creation, content, and validation of the Dark Energy Survey (DES) internal year-one cosmology data set, Y1A1 GOLD, in support of upcoming cosmological analyses. The Y1A1 GOLD data set ...is assembled from multiple epochs of DES imaging and consists of calibrated photometric zero-points, object catalogs, and ancillary data products-e.g., maps of survey depth and observing conditions, star-galaxy classification, and photometric redshift estimates-that are necessary for accurate cosmological analyses. The Y1A1 GOLD wide-area object catalog consists of million objects detected in co-added images covering in the DES grizY filters. The 10 limiting magnitude for galaxies is , , , , and . Photometric calibration of Y1A1 GOLD was performed by combining nightly zero-point solutions with stellar locus regression, and the absolute calibration accuracy is better than 2% over the survey area. DES Y1A1 GOLD is the largest photometric data set at the achieved depth to date, enabling precise measurements of cosmic acceleration at z 1.
We introduce redMaGiC, an automated algorithm for selecting luminous red galaxies (LRGs). The algorithm was specifically developed to minimize photometric redshift uncertainties in photometric ...large-scale structure studies. redMaGiC achieves this by self-training the colour cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photo-zs are very nearly as accurate as the best machine learning-based methods, yet they require minimal spectroscopic training, do not suffer from extrapolation biases, and are very nearly Gaussian. We apply our algorithm to Dark Energy Survey (DES) Science Verification (SV) data to produce a redMaGiC catalogue sampling the redshift range z ∈ 0.2, 0.8. Our fiducial sample has a comoving space density of 10−3 (h
−1 Mpc)−3, and a median photo-z bias (z
spec − z
photo) and scatter (σ
z
/(1 + z)) of 0.005 and 0.017, respectively. The corresponding 5σ outlier fraction is 1.4 per cent. We also test our algorithm with Sloan Digital Sky Survey Data Release 8 and Stripe 82 data, and discuss how spectroscopic training can be used to control photo-z biases at the 0.1 per cent level.
We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based ...photometric redshift methods-annz2, bpz calibrated against BCC-Ufig simulations, skynet, and tpz-are analyzed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-z's. From the galaxies in the DES SV shear catalogue, which have mean redshift 0.72 + or - 0.01 over the range 0.3 < z< 1.3, we construct three tomographic bins with means of z= {0.45,0.67,1.00}. These bins each have systematic uncertainties delta sub(z)<, ~ 0.05 in the mean of the fiducial skynet photo-z n(z). We propagate the errors in the redshift distributions through to their impact on cosmological parameters estimated with cosmic shear, and find that they cause shifts in the value of sigma sub(8) of approximately 3%. This shift is within the one sigma statistical errors on sigma sub(8) for the DES SV shear catalogue. We further study the potential impact of systematic differences on the critical surface density, capital sigma sub(crit), finding levels of bias safely less than the statistical power of DES SV data. We recommend a final Gaussian prior for the photo-z bias in the mean of n(z) of width 0.05 for each of the three tomographic bins, and show that this is a sufficient bias model for the corresponding cosmology analysis.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with ...the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range 0.2 < z
phot < 1.2, a cross-correlation signal is detected at 6σ and 4σ with SPT and Planck , respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant (>2σ) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the z < 1.2 universe is 0.73 ± 0.16 times as large as predicted in the Λ cold dark matter Planck cosmology, a 1.7σ deviation.
Abstract
We construct the largest curved-sky galaxy weak lensing mass map to date from the DES first-year (DES Y1) data. The map, about 10 times larger than the previous work, is constructed over a ...contiguous ≈1500 deg2, covering a comoving volume of ≈10 Gpc3. The effects of masking, sampling, and noise are tested using simulations. We generate weak lensing maps from two DES Y1 shear catalogues, MetaCalibration and Im3shape, with sources at redshift 0.2 < z < 1.3, and in each of four bins in this range. In the highest signal-to-noise map, the ratio between the mean signal to noise in the E-mode map and the B-mode map is ∼1.5 (∼2) when smoothed with a Gaussian filter of σG = 30 (80) arcmin. The second and third moments of the convergence κ in the maps are in agreement with simulations. We also find no significant correlation of κ with maps of potential systematic contaminants. Finally, we demonstrate two applications of the mass maps: (1) cross-correlation
with different foreground tracers of mass and (2) exploration of the largest peaks and voids in the maps.