We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $\kappa_{\rm ext}$ at the position of B1608+656 is $\kappa_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

Massive neutrinos influence the background evolution of the Universe as well as the growth of structure. Being able to model this effect and constrain the sum of their masses is one of the key challenges in modern cosmology. Weak gravitational lensing by the large-scale structure has proven to be a powerful tool to achieve this goal, and its importance to precision cosmology is borne out in the scientific results of galaxy surveys such as the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES) and Hyper SuprimeCam (HSC). Weak-lensing cosmological constraints will also soon reach higher levels of precision with next-generation surveys like LSST, WFIRST and Euclid. In this context, we use the MassiveNus simulations to derive constraints on the sum of neutrino masses $M_{\nu}$, the present-day total matter density ${\Omega}_m$, and the primordial power spectrum normalization $A_s$ in a tomographic setting. Along with the lensing power spectrum, we use peak counts as higher-order statistics. We find that the combination of the two statistics does not add any relevant information over considering just the peaks alone. We also compare two different filtering techniques on noisy maps and show that a starlet (wavelet) filter produces tighter constraints than Gaussian filtering by 50%, 25%, and 38% on $M_{\nu}$, ${\Omega}_m$, and $A_s$, respectively. For the starlet case, we further find a minimum resolution that allows us to obtain constraints comparable to what is achieved with the full wavelet decomposition, and we show that the information contained in the coarse map cannot be neglected.

We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray-tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass $M_{\nu} =$ 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of $M_{\nu}$. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($\Lambda$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $\Lambda$CDM from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis.

Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak counts and Minkowski functionals are combined into a joint feature vector and compared to the descriptors and statistics that are common to the field of digital image processing. To learn new features directly from the data we use a Convolutional Neural Network (CNN). For the mapping between feature vectors and the predictions of their underlying model, we implement two different classifiers; one based on a nearest-neighbour search and one that is based on a fully connected neural network. We find that the neural network provides a much more robust classification than the nearest-neighbour approach and that the CNN provides the most discriminating representation of the data. It achieves the cleanest separation between the different models and the highest classification success rate of 59% for a single source redshift. Once we perform a tomographic CNN analysis, the total classification accuracy increases significantly to 76% with no observational degeneracies remaining. Visualising the filter responses of the CNN at different network depths provides us with the unique opportunity to learn from very complex models and to understand better why they perform so well.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

General relativity (GR) has been well tested up to solar system scales, but it is much less certain that standard gravity remains an accurate description on the largest, that is, cosmological, scales. Many extensions to GR have been studied that are not yet ruled out by the data, including by that of the recent direct gravitational wave detections. Degeneracies among the standard model ($\Lambda$CDM) and modified gravity (MG) models, as well as among different MG parameters, must be addressed in order to best exploit information from current and future surveys and to unveil the nature of dark energy. We propose various higher-order statistics in the weak-lensing signal as a new set of observables able to break degeneracies between massive neutrinos and MG parameters. We have tested our methodology on so-called $f(R)$ models, which constitute a class of viable models that can explain the accelerated universal expansion by a modification of the fundamental gravitational interaction. We have explored a range of these models that still fit current observations at the background and linear level, and we show using numerical simulations that certain models which include massive neutrinos are able to mimic $\Lambda$CDM in terms of the 3D power spectrum of matter density fluctuations. We find that depending on the redshift and angular scale of observation, non-Gaussian information accessed by higher-order weak-lensing statistics can be used to break the degeneracy between $f(R)$ models and $\Lambda$CDM. In particular, peak counts computed in aperture mass maps outperform third- and fourth-order moments.

Merging galaxy clusters present a unique opportunity to study the properties of dark matter in an astrophysical context. These are rare and extreme cosmic events in which the bulk of the baryonic matter becomes displaced from the dark matter halos of the colliding subclusters. Since all mass bends light, weak gravitational lensing is a primary tool to study the total mass distribution in such systems. Combined with X-ray and optical analyses, mass maps of cluster mergers reconstructed from weak-lensing observations have been used to constrain the self-interaction cross-section of dark matter. The dynamically complex Abell 520 (A520) cluster is an exceptional case, even among merging systems: multi-wavelength observations have revealed a surprising high mass-to-light concentration of dark mass, the interpretation of which is difficult under the standard assumption of effectively collisionless dark matter. We revisit A520 using a new sparsity-based mass-mapping algorithm to independently assess the presence of the puzzling dark core. We obtain high-resolution mass reconstructions from two separate galaxy shape catalogs derived from Hubble Space Telescope observations of the system. Our mass maps agree well overall with the results of previous studies, but we find important differences. In particular, although we are able to identify the dark core at a certain level in both data sets, it is at much lower significance than has been reported before using the same data. As we cannot confirm the detection in our analysis, we do not consider A520 as posing a significant challenge to the collisionless dark matter scenario.

Peak statistics in weak lensing maps access the non-Gaussian information contained in the large-scale distribution of matter in the Universe. They are therefore a promising complement to two-point and higher-order statistics to constrain our cosmological models. To prepare for the high-precision data of next-generation surveys, we assess the constraining power of peak counts in a simulated Euclid-like survey on the cosmological parameters $\Omega_\mathrm{m}$, $\sigma_8$, and $w_0^\mathrm{de}$. In particular, we study how the Camelus model--a fast stochastic algorithm for predicting peaks--can be applied to such large surveys. We measure the peak count abundance in a mock shear catalogue of ~5,000 sq. deg. using a multiscale mass map filtering technique. We then constrain the parameters of the mock survey using Camelus combined with approximate Bayesian computation (ABC). We find that peak statistics yield a tight but significantly biased constraint in the $\sigma_8$-$\Omega_\mathrm{m}$ plane, indicating the need to better understand and control the model's systematics. We calibrate the model to remove the bias and compare results to those from the two-point correlation functions (2PCF) measured on the same field. In this case, we find the derived parameter $\Sigma_8=\sigma_8(\Omega_\mathrm{m}/0.27)^\alpha=0.76_{-0.03}^{+0.02}$ with $\alpha=0.65$ for peaks, while for 2PCF the value is $\Sigma_8=0.76_{-0.01}^{+0.02}$ with $\alpha=0.70$. We therefore see comparable constraining power between the two probes, and the offset of their $\sigma_8$-$\Omega_\mathrm{m}$ degeneracy directions suggests that a combined analysis would yield tighter constraints than either measure alone. As expected, $w_0^\mathrm{de}$ cannot be well constrained without a tomographic analysis, but its degeneracy directions with the other two varied parameters are still clear for both peaks and 2PCF. (abridged)

We study effects of inhomogeneities on distance measures in an exact relativistic Swiss-cheese model of the universe, focusing on the distance modulus. The model has LCDM background dynamics, and the `holes' are non-symmetric structures described by the Szekeres metric. The Szekeres exact solution of Einstein's equations, which is inhomogeneous and anisotropic, allows us to capture potentially relevant effects on light propagation due to nontrivial evolution of structures in an exact framework. Light beams traversing a single Szekeres structure in different ways can experience either magnification or demagnification, depending on the particular path, and we explore a small additional effect due time evolution of the structures. We study the probability distributions of $\Delta\mu=\mu_{LCDM}-\mu_{SC}$ for sources at different redshifts in various Swiss-cheese constructions, where the light beams travel through a large number of randomly oriented holes with random impact parameters. We find for $\Delta\mu$ the dispersions $0.004\le \sigma_{\Delta\mu} \le 0.008$ mag for sources with redshifts $1.0\le z \le 1.5$, which are smaller than the intrinsic dispersion of, for example, magnitudes of type Ia supernovae. The shapes of the distributions we obtain for our Swiss-cheese constructions are peculiar in the sense that they are not consistently skewed toward the demagnification side, as they are in analyses of lensing in cosmological simulations, reflecting a limitation of these constructions. This could be the result of requiring the continuity of Einstein's equations throughout the overall spacetime patchwork, which imposes the condition that compensating overdense shells must accompany the underdense void regions in the holes. The possibility to explore other uses of these constructions that could circumvent this limitation and lead to different statistics remains open. (Abridged)

The study of relativistic, higher order and nonlinear effects has become necessary in recent years in the pursuit of precision cosmology. We develop and apply here a framework to study gravitational lensing in exact models in general relativity that are not restricted to homogeneity and isotropy, and where full nonlinearity and relativistic effects are included. We apply the framework to a specific, anisotropic galaxy cluster model which is based on a modified NFW halo density profile and described by the Szekeres metric. We examine the effects of increasing levels of anisotropy in the galaxy cluster on lensing observables like the convergence and shear for various lensing geometries, finding a strong nonlinear response in both the convergence and shear for rays passing through anisotropic regions of the cluster. Deviation from the expected values in a spherically symmetric structure are asymmetric with respect to path direction and thus will persist as a statistical effect when averaged over some ensemble of such clusters. The resulting relative difference in various geometries can be as large as approximately 2%, 8%, and 24% in the measure of convergence for levels of anisotropy of 5%, 10%, and 15%, respectively, as a fraction of total cluster mass. For the total magnitude of shear, the relative difference can grow near the center of the structure to be as large as 15%, 32%, and 44% for the same levels of anisotropy, averaged over the two extreme geometries. The convergence is impacted most strongly for rays which pass in directions along the axis of maximum dipole anisotropy in the structure, while the shear is most strongly impacted for rays which pass in directions orthogonal to this axis, as expected. These effects due to anisotropic structures will affect lensing measurements and must be fully examined in an era of precision cosmology.

We study the effects and implications of anisotropies at the scale of galaxy clusters by building an exact general relativistic model of a cluster using the inhomogeneous and anisotropic Szekeres metric. The model is built from a modified Navarro-Frenk-White (NFW) density profile. We compare this to a corresponding spherically symmetric structure in the Lemaitre-Tolman (LT) model and quantify the impact of introducing varying levels of anisotropy. We examine two physical measures of gravitational infall -- the growth rate of density and the velocity of the source dust in the model. We introduce a generalization of the LT dust velocity profile for the Szekeres metric and demonstrate its consistency with the growth rate of density. We find that the growth rate of density in one substructure increases by 0.5%, 1.5%, and 3.75% for 5%, 10%, and 15% levels of introduced anisotropy, which is measured as the fractional displaced mass relative to the spherically symmetric case. The infall velocity of the dust is found to increase by 2.5, 10, and 20 km/s (0.5%, 2%, and 4.5%), respectively, for the same three levels of anisotropy. This response to the anisotropy in a structure is found to be strongly nonlinear with respect to the strength of anisotropy. These relative velocities correspond to an equivalent increase in the total mass of the spherically symmetric structure of 1%, 3.8%, and 8.4%, indicating that not accounting for the presence of anisotropic mass distributions in cluster models can strongly bias the determination of physical properties like the total mass.

Probes of cosmic expansion constitute the main basis for arguments to support or refute a possible apparent acceleration due to different expansion rates in the universe as described by inhomogeneous cosmological models. We present in this Letter a separate argument based on results from an analysis of the growth rate of large-scale structure in the universe as modeled by the inhomogeneous cosmological models of Szekeres. We use the models with no assumptions of spherical or axial symmetries. We find that while the Szekeres models can fit very well the observed expansion history without a $\Lambda$, they fail to produce the observed late-time suppression in the growth unless $\Lambda$ is added to the dynamics. A simultaneous fit to the supernova and growth factor data shows that the cold dark matter model with a cosmological constant ($\Lambda$CDM) provides consistency with the data at a confidence level of 99.65% while the Szekeres model without $\Lambda$ achieves only a 60.46% level. When the data sets are considered separately, the Szekeres with no $\Lambda$ fits the supernova data as well as the $\Lambda$CDM does, but provides a very poor fit to the growth data with only 31.31% consistency level compared to 99.99% for the $\Lambda$CDM. This absence of late-time growth suppression in inhomogeneous models without a $\Lambda$ is consolidated by a physical explanation.

We use the Szekeres inhomogeneous cosmological models to study the growth of large-scale structure in the universe including nonzero spatial curvature and a cosmological constant. In particular, we use the Goode and Wainwright formulation, as in this form the models can be considered to represent exact nonlinear perturbations of an averaged background. We identify a density contrast in both classes I and II of the models, for which we derive growth evolution equations. By including Lambda, the time evolution of the density contrast as well as kinematic quantities can be tracked through the matter- and Lambda-dominated cosmic eras up to the present and into the future. In various models of class I and class II, the growth rate is found to be stronger than that of the LCDM cosmology, and it is suppressed at later times due to the presence of Lambda. We find that there are Szekeres models able to provide a growth history similar to that of LCDM while requiring less matter content and nonzero curvature, which speaks to the importance of including the effects of large-scale inhomogeneities in analyzing the growth of large-scale structure. Using data for the growth factor f from redshift space distortions and the Lyman-alpha forest, we obtain best fit parameters for class II models and compare their ability to match observations with LCDM. We find that there is negligible difference between best fit Szekeres models with no priors and those for LCDM, both including and excluding Lyman-alpha data. We also find that the growth index gamma parametrization cannot be applied in a simple way to the growth in Szekeres models, so a direct comparison of the function f to the data is performed. We conclude that the Szekeres models can provide an exact framework for the analysis of large-scale growth data that includes inhomogeneities and allows for different interpretations of observations. (abridged)

This study belongs to a series devoted to using Szekeres inhomogeneous models to develop a theoretical framework where observations can be investigated with a wider range of possible interpretations. We look here into the growth of large-scale structure in the models. The Szekeres models are exact solutions to Einstein's equations that were originally derived with no symmetries. We use a formulation of the models that is due to Goode and Wainwright, who considered the models as exact perturbations of an FLRW background. Using the Raychaudhuri equation, we write for the two classes of the models, exact growth equations in terms of the under/overdensity and measurable cosmological parameters. The new equations in the overdensity split into two informative parts. The first part, while exact, is identical to the growth equation in the usual linearly perturbed FLRW models, while the second part constitutes exact non-linear perturbations. We integrate numerically the full exact growth rate equations for the flat and curved cases. We find that for the matter-dominated era, the Szekeres growth rate is up to a factor of three to five stronger than the usual linearly perturbed FLRW cases, reflecting the effect of exact Szekeres non-linear perturbations. The growth is also stronger than that of the non-linear spherical collapse model, and the difference between the two increases with time. This highlights the distinction when we use general inhomogeneous models where shear and a tidal gravitational field are present and contribute to the gravitational clustering. Additionally, it is worth observing that the enhancement of the growth found in the Szekeres models during the matter-dominated era could suggest a substitute to the argument that dark matter is needed when using FLRW models to explain the enhanced growth and resulting large-scale structures that we observe today (abridged)

The peculiar velocity of the Local Group of galaxies manifested in the Cosmic Microwave Background dipole is found to decompose into three dominant components. The three components are clearly separated because they arise on distinct spatial scales and are fortuitously almost orthogonal in their influences. The nearest, which is distinguished by a velocity discontinuity at ~7 Mpc, arises from the evacuation of the Local Void. We lie in the Local Sheet that bounds the void. Random motions within the Local Sheet are small. Our Galaxy participates in the bulk motion of the Local Sheet away from the Local Void. The component of our motion on an intermediate scale is attributed to the Virgo Cluster and its surroundings, 17 Mpc away. The third and largest component is an attraction on scales larger than 3000 km/s and centered near the direction of the Centaurus Cluster. The amplitudes of the three components are 259, 185, and 455 km/s, respectively, adding collectively to 631 km/s in the reference frame of the Local Sheet. Taking the nearby influences into account causes the residual attributed to large scales to align with observed concentrations of distant galaxies and reduces somewhat the amplitude of motion attributed to their pull. On small scales, in addition to the motion of our Local Sheet away from the Local Void, the nearest adjacent filament, the Leo Spur, is seen to be moving in a direction that will lead to convergence with our filament. Finally, a good distance to an isolated galaxy within the Local Void reveals that this dwarf system has a motion of at least 230 km/s away from the void center. Given the velocities expected from gravitational instability theory in the standard cosmological paradigm, the distance to the center of the Local Void must be at least 23 Mpc from our position. The Local Void is large!

Upcoming surveys for galaxy clusters using the Sunyaev-Zel'dovich effect are potentially sensitive enough to create a peculiar velocity catalog. The statistics of these peculiar velocities are sensitive to cosmological parameters. We develop a method to explore parameter space using N-body simulations in order to quantify dark matter halo velocity statistics which will be useful for cluster peculiar velocity observations. We show that mass selection bias from a kinetic Sunyaev-Zel'dovich velocity catalog forecasts rms peculiar velocities with a much more complicated Omega_m dependency than suggested by linear theory. In addition, we show that both two-point functions for velocities disagree with linear theory predictions out to separations of approx. 40 Mpc/h. A pedagogical appendix is included developing linear theory notation with respect to the two-point peculiar velocity functions.

COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter operating at a wavelength of 1 cm. The entire instrument was built specifically for CMB polarization studies. We report here on observations of February 2001 - April 2001 using this system. We set an upper limit on E-mode polarized anisotropies of 33.5 uK (95% confidence limit) in the l-range 200-600.

COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter. The entire instrument was built specifically for CMB polarization studies. Careful attention was given to receiver and optics design, stability of the pointing platform, avoidance of systematic offsets, and development of data analysis techniques. Here we describe the experiment, its strengths and weaknesses, and the various things we have learned that may benefit future efforts to measure the polarization of the CMB.

We present an analytic method for rapidly forecasting the accuracy of gravitational potential reconstruction possible from measurement of radial peculiar velocities of every galaxy cluster with M > M_th in solid angle \theta^2 and over redshift range z_min < z < z_max. These radial velocities can be determined from measurement of the kinetic and thermal Sunyaev-Zeldovich effects. For a shallow survey with 0.2 < z < 0.4, coincident with the SDSS photometric survey, one mode of the gravitational potential (on length scales > 60 Mpc) can be reconstructed for every ~8 cluster velocity determinations. Deeper surveys require measurement of more clusters per S/N > 1 mode. Accuracy is limited by the ``undersampling noise'' due to our non-observation of the large fraction of mass that is not in galaxy clusters. Determining the gravitational potential will allow for detailed study of the relationship between galaxies and their surrounding large-scale density fields over a wide range of redshifts, and test the gravitational instability paradigm on very large scales. Observation of weak lensing by large-scale structure provides complementary information since lensing is sensitive to the tangential modes that do not affect the velocity.

Galaxy cluster peculiar velocities can be inferred from high-sensitivity, high-resolution multiple-frequency observations in the 30 to 400 GHz range. While galaxy cluster counts and power spectra are sensitive to the growth factor, peculiar velocities are sensitive to the time-derivative of the growth factor and are hence complementary. Using linear perturbation theory, we forecast constraints on Omega(matter), the Hubble constant today, and the dark-energy equation of state parameter, given 820 densely sampled cluster locations (near redshift 1) from a Lambda-CDM N-body simulation and 820 sparsely sampled cluster locations in a broader redshift range.

We have developed a fast, accurate and generally applicable method for inferring the power spectrum and its uncertainties from maps of the cosmic microwave background (CMB) in the presence of inhomogeneous and correlated noise. For maps with 10 to 100 thousand pixels, we apply an exact power spectrum estimation algorithm to submaps of the data at various resolutions, and then combine the results in an optimal manner. To analyze larger maps efficiently one must resort to sub-optimal combinations in which cross-map power spectrum error correlations are only calculated approximately. We expect such approximations to work well in general, and in particular for the megapixel maps to come from the next generation of satellite missions.