The intrinsic alignments of galaxies are recognized as a contaminant to weak gravitational lensing measurements. In this work, we study the alignment of galaxy shapes and spins at low redshift (z ∼ ...0.5) in Horizon-AGN, an adaptive-mesh-refinement hydrodynamical cosmological simulation box of 100 h
− 1 Mpc a side with AGN feedback implementation. We find that spheroidal galaxies in the simulation show a tendency to be aligned radially towards overdensities in the dark matter density field and other spheroidals. This trend is in agreement with observations, but the amplitude of the signal depends strongly on how shapes are measured and how galaxies are selected in the simulation. Disc galaxies show a tendency to be oriented tangentially around spheroidals in three dimensions. While this signal seems suppressed in projection, this does not guarantee that disc alignments can be safely ignored in future weak lensing surveys. The shape alignments of luminous galaxies in Horizon-AGN are in agreement with observations and other simulation works, but we find less alignment for lower luminosity populations. We also characterize the systematics of galaxy shapes in the simulation and show that they can be safely neglected when measuring the correlation of the density field and galaxy ellipticities.
We present candl , an automatically differentiable python likelihood for analysing cosmic microwave background power spectrum measurements. candl is powered by JAX, which makes it fast and easy to ...calculate derivatives of the likelihood. This facilitates, for example, robust Fisher matrices without finite-difference methods. We show the benefits of candl through a series of example calculations, covering forecasting, robustness tests, and gradient-based Markov chain Monte Carlo sampling. These also include optimising the band power bin width to minimise parameter errors of a realistic mock data set. Moreover, we calculate the correlation of parameter constraints from correlated and partially overlapping subsets of the SPT-3G 2018 TT/TE/EE data release. In a traditional analysis framework, these tasks are slow and require careful fine-tuning to obtain stable results. As such, a fully differentiable pipeline allows for a higher level of scrutiny; we argue that this is the paradigm shift required to leverage incoming data from ground-based experiments, which will significantly improve the cosmological parameter constraints from the Planck mission. candl comes with the latest primary and lensing power spectrum data from the South Pole Telescope and Atacama Cosmology Telescope collaborations and will be used as part of the upcoming SPT-3G TT/TE/EE and ϕϕ data releases. Along with the core code, we release a series of auxiliary tools, which simplify common analysis tasks and interface the likelihood with other cosmological software. candl is pip-installable and publicly available on Github.
The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z = 1.2 ...synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation horizon-agn. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field ‘filamentary’ direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalogue are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to ∼10 h
− 1 Mpc comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate-mass galaxies on a similar scale at a level of 10−4 which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step towards the use of very realistic catalogue of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.
A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more ...than 150 000 galaxies per time step in the redshift range 1.2 < z < 1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3 × 1010 M⊙. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.
Intrinsic galaxy shape and angular momentum alignments can arise in cosmological large-scale structure due to tidal interactions or galaxy formation processes. Cosmological hydrodynamical simulations ...have recently come of age as a tool to study these alignments and their contamination to weak gravitational lensing. We probe the redshift and luminosity evolution of intrinsic alignments in Horizon-AGN between z = 0 and 3 for galaxies with an r-band absolute magnitude of Mr
≤ −20. Alignments transition from being radial at low redshifts and high luminosities, dominated by the contribution of ellipticals, to being tangential at high redshift and low luminosities, where discs dominate the signal. This cannot be explained by the evolution of the fraction of ellipticals and discs alone: intrinsic evolution in the amplitude of alignments is necessary. The alignment amplitude of elliptical galaxies alone is smaller in amplitude by a factor of ≃2, but has similar luminosity and redshift evolution as in current observations and in the non-linear tidal alignment model at projected separations of ≳1 Mpc. Alignments of discs are null in projection and consistent with current low-redshift observations. The combination of the two populations yields an overall amplitude a factor of ≃4 lower than observed alignments of luminous red galaxies with a steeper luminosity dependence. The restriction on accurate galaxy shapes implies that the galaxy population in the simulation is complete only to Mr
≤ −20. Higher resolution simulations will be necessary to avoid extrapolation of the intrinsic alignment predictions to the range of luminosities probed by future surveys.
We have performed a 70 billion dark-matter particles N-body simulation in a 2 h-1 Gpc periodic box, using the concordance, cosmological model as favored by the latest WMAP3 results. We have computed ...a full-sky convergence map with a resolution of $\Delta \theta \simeq 0.74$ arcmin2, spanning 4 orders of magnitude in angular dynamical range. Using various high-order statistics on a realistic cut sky, we have characterized the transition from the linear to the nonlinear regime at $\ell \simeq 1000$ and shown that realistic galactic masking affects high-order moments only below $\ell < 200$. Each domain (Gaussian and non-Gaussian) spans 2 decades in angular scale. This map is therefore an ideal tool for testing map-making algorithms on the sphere. As a first step in addressing the full map reconstruction problem, we have benchmarked in this paper two denoising methods: 1) Wiener filtering applied to the Spherical Harmonics decomposition of the map and 2) a new method, called MRLens, based on the modification of the Maximum Entropy Method on a Wavelet decomposition. While the latter is optimal on large spatial scales, where the signal is Gaussian, MRLens outperforms the Wiener method on small spatial scales, where the signal is highly non-Gaussian. The simulated full-sky convergence map is freely available to the community to help the development of new map-making algorithms dedicated to the next generation of weak-lensing surveys.
We present the first measurements of the weak gravitational lensing signal induced by the large-scale mass distribution in the universe from data obtained as part of the ongoing Canada-France-Hawaii ...Telescope Legacy Survey (CFHTLS). The data used in this analysis are from the Wide Synoptic Survey, which aims to image 6170 deg super(2) in five filters. We have analyzed an effective area of 622 deg super(2) (31 pointings) of i data spread over two of the three survey fields. These data are of excellent quality, and the results bode well for the remainder of the survey: we do not detect a significant "B" mode, suggesting that residual systematics are negligible at the current level of accuracy. Assuming a cold dark matter model and marginalizing over the Hubble parameter h e 0.6, 0.8, the source redshift distribution, and systematics, we constrain s sub(8), the amplitude of the matter power spectrum. At a fiducial matter density sub(m) = 0.3 we find s sub(8) = 0.85 c 0.06. This estimate is in excellent agreement with previous studies. A combination of our results with those from the Deep component of the CFHTLS enables us to place a constraint on a constant equation of state for the dark energy, based on cosmic shear data alone. We find that w sub(0) < -0.8 at 68% confidence.
We present a comprehensive analysis of weak gravitational lensing by large-scale structure in the Hubble Space Telescope Cosmic Evolution Survey (COSMOS), in which we combine space-based galaxy shape ...measurements with ground-based photometric redshifts to study the redshift dependence of the lensing signal and constrain cosmological parameters. After applying our weak lensing-optimized data reduction, principal-component interpolation for the spatially, and temporally varying ACS point-spread function, and improved modelling of charge-transfer inefficiency, we measured a lensing signal that is consistent with pure gravitational modes and no significant shape systematics. We carefully estimated the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation, including the full non-Gaussian sampling variance. We tested our lensing pipeline on simulated space-based data, recalibrated non-linear power spectrum corrections using the ray-tracing analysis, employed photometric redshift information to reduce potential contamination by intrinsic galaxy alignments, and marginalized over systematic uncertainties. We find that the weak lensing signal scales with redshift as expected from general relativity for a concordance ΛCDM cosmology, including the full cross-correlations between different redshift bins. Assuming a flat ΛCDM cosmology, we measure $\sigma_8$($\Omega_\mathrm{m}$/0.3$)^{0.51}$ = 0.75±0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints $\Omega_\mathrm{m}$ = $0.266^{+0.025}_{-0.023}$, $\sigma_8$ = $0.802^{+0.028}_{-0.029}$ (all 68.3% conf.). Dropping the assumption of flatness and using priors from the HST Key Project and Big-Bang nucleosynthesis only, we find a negative deceleration parameter q0 at 94.3% confidence from the tomographic lensing analysis, providing independent evidence of the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w ∈ -2,0, we obtain w <-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration of the usefulness of tomographic weak lensing measurements from space.
Aims.We present an exploration of weak lensing by large-scale structure in the linear regime, using the third-year (T0003) CFHTLS Wide data release. Our results place tight constraints on the scaling ...of the amplitude of the matter power spectrum $\sigma_8$ with the matter density $\Omega_{\rm m}$. Methods.Spanning 57 square degrees to $i^{\prime}_{AB} = 24.5$ over three independent fields, the unprecedented contiguous area of this survey permits high signal-to-noise measurements of two-point shear statistics from 1 arcmin to 4 degrees. Understanding systematic errors in our analysis is vital in interpreting the results. We therefore demonstrate the percent-level accuracy of our method using STEP simulations, an E/B-mode decomposition of the data, and the star-galaxy cross correlation function. We also present a thorough analysis of the galaxy redshift distribution using redshift data from the CFHTLS T0003 Deep fields that probe the same spatial regions as the Wide fields. Results.We find $\sigma_8(\Omega_{\rm m} / 0.25)^{0.64} = 0.785$ ± 0.043 using the aperture-mass statistic for the full range of angular scales for an assumed flat cosmology, in excellent agreement with WMAP3 constraints. The largest physical scale probed by our analysis is 85 Mpc, assuming a mean redshift of lenses of 0.5 and a ΛCDM cosmology. This allows for the first time to constrain cosmology using only cosmic shear measurements in the linear regime. Using only angular scales $\theta> 85$ arcmin, we find $\sigma_8(\Omega_{\rm m} / 0.25)_{{\rm lin}}^{0.53} = 0.837$ ± 0.084, which agree with the results from our full analysis. Combining our results with data from WMAP3, we find $\Omega_{\rm m}=0.248$ ± 0.019 and $\sigma_8=0.771$ ± 0.029.
We propose a novel approximation to the low-ℓ joint likelihood of the angular spectrum Cℓ of masked cosmic microwave background temperature maps which is both very accurate and very fast to evaluate. ...We show that, for a flat prior, the posterior distribution of each Cℓ closely follows an inverse gamma distribution even with partial sky coverage and that the posterior correlation is weak enough that a copula approximation to the joint likelihood is quite accurate. In this paper, the quantities needed to build such a copula approximation (inverse gamma parameters at each angular frequency and a correlation matrix) are computed from an exploration of the posterior using adaptive importance sampling. The accuracy of the proposed approximation is assessed using statistical criteria as well as a mock cosmological parameter fit. When applied to the Wilkinson Microwave Anisotropy Probe 5 data set, the copula approximation yields cosmological parameter estimates at the same level of accuracy as the best current techniques.