The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish General Relativity from its potential modifications. We use the latest weak gravitational lensing dataset from the Kilo-Degree Survey, KiDS-1000, in conjunction with overlapping galaxy spectroscopic redshift surveys BOSS and 2dFLenS, to perform the most precise existing amplitude-ratio test. We measure the associated E_G statistic with 15-20% errors, in five dz = 0.1 tomographic redshift bins in the range 0.2 < z < 0.7, on projected scales up to 100 Mpc/h. The scale-independence and redshift-dependence of these measurements are consistent with the theoretical expectation of General Relativity in a Universe with matter density Omega_m = 0.27 +/- 0.04. We demonstrate that our results are robust against different analysis choices, including schemes for correcting the effects of source photometric redshift errors, and compare the performance of angular and projected galaxy-galaxy lensing statistics.

We present a combined tomographic weak gravitational lensing analysis of the Kilo Degree Survey (KV450) and the Dark Energy Survey (DES-Y1). We homogenize the analysis of these two public cosmic shear datasets by adopting consistent priors and modeling of nonlinear scales, and determine new redshift distributions for DES-Y1 based on deep public spectroscopic surveys. Adopting these revised redshifts results in a $0.8\sigma$ reduction in the DES-inferred value for $S_8$, which decreases to a $0.5\sigma$ reduction when including a systematic redshift calibration error model from mock DES data based on the MICE2 simulation. The combined KV450 + DES-Y1 constraint on $S_8 = 0.762^{+0.025}_{-0.024}$ is in tension with the Planck 2018 constraint from the cosmic microwave background at the level of $2.5\sigma$. This result highlights the importance of developing methods to provide accurate redshift calibration for current and future weak lensing surveys.

We present cosmological constraints from a joint cosmic shear analysis of the Kilo-Degree Survey (KV450) and the Dark Energy Survey (DES-Y1), conducted using Complete Orthogonal Sets of E/B-Integrals (COSEBIs). With COSEBIs we isolate any B-modes which have a non-cosmic shear origin and demonstrate the robustness of our cosmological E-mode analysis as no significant B-modes are detected. We highlight how COSEBIs are fairly insensitive to the amplitude of the non-linear matter power spectrum at high $k$-scales, mitigating the uncertain impact of baryon feedback in our analysis. COSEBIs, therefore, allow us to utilise additional small-scale information, improving the DES-Y1 joint constraints on $S_8=\sigma_8(\Omega_{\rm m}/0.3)^{0.5}$ and $\Omega_{\rm m}$ by $20\%$. Adopting a flat $\Lambda$CDM model we find $S_8=0.755^{+0.019}_{-0.021}$, which is in $3.2\sigma$ tension with the Planck Legacy analysis of the cosmic microwave background.

We explore the implications of a rapid appearance of dark energy between the redshifts ($z$) of one and two on the expansion rate and growth of perturbations. Using both Gaussian process regression and a parameteric model, we show that this is the preferred solution to the current set of low-redshift ($z<3$) distance measurements if $H_0=73~\rm km\,s^{-1}\,Mpc^{-1}$ to within 1\% and the high-redshift expansion history is unchanged from the $\Lambda$CDM inference by the Planck satellite. Dark energy was effectively non-existent around $z=2$, but its density is close to the $\Lambda$CDM model value today, with an equation of state greater than $-1$ at $z<0.5$. If sources of clustering other than matter are negligible, we show that this expansion history leads to slower growth of perturbations at $z<1$, compared to $\Lambda$CDM, that is measurable by upcoming surveys and can alleviate the $\sigma_8$ tension between the Planck CMB temperature and low-redshift probes of the large-scale structure.

We reanalyse the anisotropic galaxy clustering measurement from the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that using the full shape information provides cosmological constraints that are comparable to other low-redshift probes. We find $\Omega_\mathrm{m} = 0.317^{+0.015}_{-0.019}$, $\sigma_8 = 0.710\pm 0.049$, and $h = 0.704\pm 0.024$ for flat $\Lambda$CDM cosmologies using uninformative priors on $\Omega_\mathrm{c}h^2$, $100\theta_\mathrm{MC}$, $\ln 10^{10} A_{s}$, and $n_{s}$, and a prior on $\Omega_\mathrm{b}h^2$ that is much wider than current constraints. We quantify the agreement between the Planck 2018 constraints from the cosmic microwave background and BOSS, finding the two data sets to be consistent within a flat $\Lambda$CDM cosmology using the Bayes factor as well as the prior-insensitive suspiciousness statistic. Combining two low-redshift probes, we jointly analyse the clustering of BOSS galaxies with weak lensing measurements from the Kilo-Degree Survey (KV450). The combination of BOSS and KV450 improves the measurement by up to 45%, constraining $\sigma_8 = 0.702\pm 0.029$ and $S_8 = \sigma_8\sqrt{\Omega_\mathrm{m}/0.3} = 0.728\pm 0.026$. Over the full 5D parameter space, the odds in favour of single cosmology describing galaxy clustering, lensing, and the cosmic microwave background are $7\pm2$. The suspiciousness statistic signals a $2.1\pm0.3\sigma$ tension between the combined low-redshift probes and measurements from the cosmic microwave background.

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $\xi \ge 5$ and an additional 204 systems at $4<\xi<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $\xi>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($\lambda-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $\sigma$ level---with the relations intersecting at $\lambda=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

3D data compression techniques can be used to determine the natural basis of radial eigenmodes that encode the maximum amount of information in a tomographic large-scale structure survey. We explore the potential of the Karhunen-Lo\`eve decomposition in reducing the dimensionality of the data vector for cosmic shear measurements, and apply it to the final data from the \cfh survey. We find that practically all of the cosmological information can be encoded in one single radial eigenmode, from which we are able to reproduce compatible constraints with those found in the fiducial tomographic analysis (done with 7 redshift bins) with a factor of ~30 fewer datapoints. This simplifies the problem of computing the two-point function covariance matrix from mock catalogues by the same factor, or by a factor of ~800 for an analytical covariance. The resulting set of radial eigenfunctions is close to ell-independent, and therefore they can be used as redshift-dependent galaxy weights. This simplifies the application of the Karhunen-Lo\`eve decomposition to real-space and Fourier-space data, and allows one to explore the effective radial window function of the principal eigenmodes as well as the associated shear maps in order to identify potential systematics. We also apply the method to extended parameter spaces and verify that additional information may be gained by including a second mode to break parameter degeneracies. The data and analysis code are publicly available at https://github.com/emiliobellini/kl_sample.

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018.

The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

We directly constrain the non-linear alignment (NLA) model of intrinsic galaxy alignments, analysing the most representative and complete flux-limited sample of spectroscopic galaxies available for cosmic shear surveys. We measure the projected galaxy position-intrinsic shear correlations and the projected galaxy clustering signal using high-resolution imaging from the Kilo Degree Survey (KiDS) overlapping with the GAMA spectroscopic survey, and data from the Sloan Digital Sky Survey. Separating samples by colour, we make no significant detection of blue galaxy alignments, constraining the blue galaxy NLA amplitude $A_{\textrm{IA}}^{\textrm{B}}=0.21^{+0.37}_{-0.36}$ to be consistent with zero. We make robust detections ($\sim9\sigma$) for red galaxies, with $A_{\textrm{IA}}^{\textrm{R}}=3.18^{+0.47}_{-0.46}$, corresponding to a net radial alignment with the galaxy density field, and we find no evidence for any scaling of alignments with galaxy luminosity. We provide informative priors for current and future weak lensing surveys, an improvement over de facto wide priors that allow for unrealistic levels of intrinsic alignment contamination. For a colour-split cosmic shear analysis of the final KiDS survey area, we forecast that our priors will improve the constraining power on $S_{8}$ and the dark energy equation of state $w_{0}$, by up to $62\%$ and $51\%$, respectively. Our results indicate, however, that the modelling of red/blue-split galaxy alignments may be insufficient to describe samples with variable central/satellite galaxy fractions.

We present three tiers of Bayesian consistency tests for the general case of $correlated$ datasets. Building on duplicates of the model parameters assigned to each dataset, these tests range from Bayesian evidence ratios as a global summary statistic, to posterior distributions of model parameter differences, to consistency tests in the data domain derived from posterior predictive distributions. For each test we motivate meaningful threshold criteria for the internal consistency of datasets. Without loss of generality we focus on mutually exclusive, correlated subsets of the same dataset in this work. As an application, we revisit the consistency analysis of the two-point weak lensing shear correlation functions measured from KiDS-450 data. We split this dataset according to large vs. small angular scales, tomographic redshift bin combinations, and estimator type. We do not find any evidence for significant internal tension in the KiDS-450 data, with significances below $3\, \sigma$ in all cases. Software and data used in this analysis can be found at http://kids.strw.leidenuniv.nl/sciencedata.php

We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $\Lambda$CDM Universe and assuming GR, we find $\Omega_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}.

We present a weak gravitational lensing analysis of 815 square degree of $i$-band imaging from the Kilo-Degree Survey (KiDS-$i$-800). In contrast to the deep $r$-band observations, which take priority during excellent seeing conditions and form the primary KiDS dataset (KiDS-$r$-450), the complementary yet shallower KiDS-$i$-800 spans a wide range of observing conditions. The overlapping KiDS-$i$-800 and KiDS-$r$-450 imaging therefore provides a unique opportunity to assess the robustness of weak lensing measurements. In our analysis, we introduce two new `null' tests. The `nulled' two-point shear correlation function uses a matched catalogue to show that the calibrated KiDS-$i$-800 and KiDS-$r$-450 shear measurements agree at the level of $1 \pm 4$\%. We use five galaxy lens samples to determine a `nulled' galaxy-galaxy lensing signal from the full KiDS-$i$-800 and KiDS-$r$-450 surveys and find that the measurements agree to $7 \pm 5$\% when the KiDS-$i$-800 source redshift distribution is calibrated using either spectroscopic redshifts, or the 30-band photometric redshifts from the COSMOS survey.

We present results from an ongoing investigation using the Kilo Degree Survey (KiDS) on the VLT Survey Telescope (VST) to provide a census of ultra-compact massive galaxies (UCMGs), defined as galaxies with stellar masses $M_{\rm \star} > 8 \times 10^{10} \rm M_{\odot}$ and effective radii $R_{\rm e} < 1.5\,\rm kpc$. UCMGs, which are expected to have undergone very few merger events, provide a unique view on the accretion history of the most massive galaxies in the Universe. Over an effective sky area of nearly 330 square degrees, we select UCMG candidates from KiDS multi-colour images, which provide high quality structural parameters, photometric redshifts and stellar masses. Our sample of $\sim 1000$ photometrically selected UCMGs at $z < 0.5$ represents the largest sample of UCMG candidates assembled to date over the largest sky area. In this paper we present the first effort to obtain their redshifts using different facilities, starting with first results for 28 candidates with redshifts $z < 0.5$, obtained at NTT and TNG telescopes. We confirmed, as bona fide UCMGs, 19 out of the 28 candidates with new redshifts. A further 46 UCMG candidates are confirmed with literature spectroscopic redshifts (35 at $z < 0.5$), bringing the final cumulative sample of spectroscopically-confirmed lower-z UCMGs to 54 galaxies, which is the largest sample at redshifts below $0.5$. We use these spectroscopic redshifts to quantify systematic errors in our photometric selection, and use these to correct our UCMG number counts. We finally compare the results to independent datasets and simulations.

In the past few years, several independent collaborations have presented cosmological constraints from tomographic cosmic shear analyses. These analyses differ in many aspects: the datasets, the shear and photometric redshift estimation algorithms, the theory model assumptions, and the inference pipelines. To assess the robustness of the existing cosmic shear results, we present in this paper a unified analysis of four of the recent cosmic shear surveys: the Deep Lens Survey (DLS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Science Verification data from the Dark Energy Survey (DES-SV), and the 450 deg$^{2}$ release of the Kilo-Degree Survey (KiDS-450). By using a unified pipeline, we show how the cosmological constraints are sensitive to the various details of the pipeline. We identify several analysis choices that can shift the cosmological constraints by a significant fraction of the uncertainties. For our fiducial analysis choice, considering a Gaussian covariance, conservative scale cuts, assuming no baryonic feedback contamination, identical cosmological parameter priors and intrinsic alignment treatments, we find the constraints (mean, 16% and 84% confidence intervals) on the parameter $S_{8}\equiv \sigma_{8}(\Omega_{\rm m}/0.3)^{0.5}$ to be $S_{8}=0.94_{-0.045}^{+0.046}$ (DLS), $0.66_{-0.071}^{+0.070}$ (CFHTLenS), $0.84_{-0.061}^{+0.062}$ (DES-SV) and $0.76_{-0.049}^{+0.048}$ (KiDS-450). From the goodness-of-fit and the Bayesian evidence ratio, we determine that amongst the four surveys, the two more recent surveys, DES-SV and KiDS-450, have acceptable goodness-of-fit and are consistent with each other. The combined constraints are $S_{8}=0.79^{+0.042}_{-0.041}$, which is in good agreement with the first year of DES cosmic shear results and recent CMB constraints from the Planck satellite.

We model the expansion history of the Universe as a Gaussian Process and find constraints on the dark energy density and its low-redshift evolution using distances inferred from the Luminous Red Galaxy (LRG) and Lyman-alpha (Ly$\alpha$) datasets of the Baryon Oscillation Spectroscopic Survey, supernova data from the Joint Light-curve Analysis (JLA) sample, Cosmic Microwave Background (CMB) data from the Planck satellite, and local measurement of the Hubble parameter from the Hubble Space Telescope ($\mathsf H0$). Our analysis shows that the CMB, LRG, Ly$\alpha$, and JLA data are consistent with each other and with a $\Lambda$CDM cosmology, but the ${\mathsf H0}$ data is inconsistent at moderate significance. Including the presence of dark radiation does not alleviate the ${\mathsf H0}$ tension in our analysis. While some of these results have been noted previously, the strength here lies in that we do not assume a particular cosmological model. We calculate the growth of the gravitational potential in General Relativity corresponding to these general expansion histories and show that they are well-approximated by $\Omega_{\rm m}^{0.55}$ given the current precision. We assess the prospects for upcoming surveys to measure deviations from $\Lambda$CDM using this model-independent approach.

We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very complete spectroscopic data overlapping with KiDS are available, the performance of the ML photo-zs clearly surpasses that of BPZ, currently the primary photo-z method for KiDS. Using the Galaxy And Mass Assembly (GAMA) spectroscopic survey as calibration, we furthermore study how photo-zs improve for bright sources when photometric parameters additional to magnitudes are included in the photo-z derivation, as well as when VIKING and WISE infrared bands are added. While the fiducial four-band ugri setup gives a photo-z bias $\delta z=-2e-4$ and scatter $\sigma_z<0.022$ at mean z = 0.23, combining magnitudes, colours, and galaxy sizes reduces the scatter by ~7% and the bias by an order of magnitude. Once the ugri and IR magnitudes are joined into 12-band photometry spanning up to 12 $\mu$, the scatter decreases by more than 10% over the fiducial case. Finally, using the 12 bands together with optical colours and linear sizes gives $\delta z<4e-5$ and $\sigma_z<0.019$. This paper also serves as a reference for two public photo-z catalogues accompanying KiDS DR3, both obtained using the ANNz2 code. The first one, of general purpose, includes all the 39 million KiDS sources with four-band ugri measurements in DR3. The second dataset, optimized for low-redshift studies such as galaxy-galaxy lensing, is limited to r<20, and provides photo-zs of much better quality than in the full-depth case thanks to incorporating optical magnitudes, colours, and sizes in the GAMA-calibrated photo-z derivation.

We present cosmological parameter constraints from a joint analysis of three cosmological probes: the tomographic cosmic shear signal in $\sim$450 deg$^2$ of data from the Kilo Degree Survey (KiDS), the galaxy-matter cross-correlation signal of galaxies from the Galaxies And Mass Assembly (GAMA) survey determined with KiDS weak lensing, and the angular correlation function of the same GAMA galaxies. We use fast power spectrum estimators that are based on simple integrals over the real-space correlation functions, and show that they are practically unbiased over relevant angular frequency ranges. We test our full pipeline on numerical simulations that are tailored to KiDS and retrieve the input cosmology. By fitting different combinations of power spectra, we demonstrate that the three probes are internally consistent. For all probes combined, we obtain $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}=0.800_{-0.027}^{+0.029}$, consistent with Planck and the fiducial KiDS-450 cosmic shear correlation function results. Marginalising over wide priors on the mean of the tomographic redshift distributions yields consistent results for $S_8$ with an increase of $28\%$ in the error. The combination of probes results in a $26\%$ reduction in uncertainties of $S_8$ over using the cosmic shear power spectra alone. The main gain from these additional probes comes through their constraining power on nuisance parameters, such as the galaxy intrinsic alignment amplitude or potential shifts in the redshift distributions, which are up to a factor of two better constrained compared to using cosmic shear alone, demonstrating the value of large-scale structure probe combination.

We compute the spherical-sky weak-lensing power spectrum of the shear and convergence. We discuss various approximations, such as flat-sky, and first- and second- order Limber equations for the projection. We find that the impact of adopting these approximations is negligible when constraining cosmological parameters from current weak lensing surveys. This is demonstrated using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that the reported tension with Planck Cosmic Microwave Background (CMB) temperature anisotropy results cannot be alleviated. For future large-scale surveys with unprecedented precision, we show that the spherical second-order Limber approximation will provide sufficient accuracy. In this case, the cosmic-shear power spectrum is shown to be in agreement with the full projection at the sub-percent level for l > 3, with the corresponding errors an order of magnitude below cosmic variance for all l. When computing the two-point shear correlation function, we show that the flat-sky fast Hankel transformation results in errors below two percent compared to the full spherical transformation. In the spirit of reproducible research, our numerical implementation of all approximations and the full projection are publicly available within the package nicaea at http://www.cosmostat.org/software/nicaea.

We perform a combined analysis of cosmic shear tomography, galaxy-galaxy lensing tomography, and redshift-space multipole power spectra (monopole and quadrupole) using 450 deg$^2$ of imaging data by the Kilo Degree Survey (KiDS) overlapping with two spectroscopic surveys: the 2-degree Field Lensing Survey (2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict the galaxy-galaxy lensing and multipole power spectrum measurements to the overlapping regions with KiDS, and self-consistently compute the full covariance between the different observables using a large suite of $N$-body simulations. We methodically analyze different combinations of the observables, finding that galaxy-galaxy lensing measurements are particularly useful in improving the constraint on the intrinsic alignment amplitude (by 30%, positive at $3.5\sigma$ in the fiducial data analysis), while the multipole power spectra are useful in tightening the constraints along the lensing degeneracy direction (e.g. factor of two stronger matter density constraint in the fiducial analysis). The fully combined constraint on $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.742 \pm 0.035$, which is an improvement by 20% compared to KiDS alone, corresponds to a $2.6\sigma$ discordance with Planck, and is not significantly affected by fitting to a more conservative set of scales. Given the tightening of the parameter space, we are unable to resolve the discordance with an extended cosmology that is simultaneously favored in a model selection sense, including the sum of neutrino masses, curvature, evolving dark energy, and modified gravity. The complementarity of our observables allows for constraints on modified gravity degrees of freedom that are not simultaneously bounded with either probe alone, and up to a factor of three improvement in the $S_8$ constraint in the extended cosmology compared to KiDS alone.

We present measurements of the weak gravitational lensing shear power spectrum based on $450$ sq. deg. of imaging data from the Kilo Degree Survey. We employ a quadratic estimator in two and three redshift bins and extract band powers of redshift auto-correlation and cross-correlation spectra in the multipole range $76 \leq \ell \leq 1310$. The cosmological interpretation of the measured shear power spectra is performed in a Bayesian framework assuming a $\Lambda$CDM model with spatially flat geometry, while accounting for small residual uncertainties in the shear calibration and redshift distributions as well as marginalising over intrinsic alignments, baryon feedback and an excess-noise power model. Moreover, massive neutrinos are included in the modelling. The cosmological main result is expressed in terms of the parameter combination $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ yielding $S_8 = \ 0.651 \pm 0.058$ (3 z-bins), confirming the recently reported tension in this parameter with constraints from Planck at $3.2\sigma$ (3 z-bins). We cross-check the results of the 3 z-bin analysis with the weaker constraints from the 2 z-bin analysis and find them to be consistent. The high-level data products of this analysis, such as the band power measurements, covariance matrices, redshift distributions, and likelihood evaluation chains are available at http://kids.strw.leidenuniv.nl/

We test extensions to the standard cosmological model with weak gravitational lensing tomography using 450 deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). In these extended cosmologies, which include massive neutrinos, nonzero curvature, evolving dark energy, modified gravity, and running of the scalar spectral index, we also examine the discordance between KiDS and cosmic microwave background measurements from Planck. The discordance between the two datasets is largely unaffected by a more conservative treatment of the lensing systematics and the removal of angular scales most sensitive to nonlinear physics. The only extended cosmology that simultaneously alleviates the discordance with Planck and is at least moderately favored by the data includes evolving dark energy with a time-dependent equation of state (in the form of the $w_0-w_a$ parameterization). In this model, the respective $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ constraints agree at the $1\sigma$ level, and there is `substantial concordance' between the KiDS and Planck datasets when accounting for the full parameter space. Moreover, the Planck constraint on the Hubble constant is wider than in LCDM and in agreement with the Riess et al. (2016) direct measurement of $H_0$. The dark energy model is moderately favored as compared to LCDM when combining the KiDS and Planck measurements, and remains moderately favored after including an informative prior on the Hubble constant. In both of these scenarios, marginalized constraints in the $w_0-w_a$ plane are discrepant with a cosmological constant at the $3\sigma$ level. Moreover, KiDS constrains the sum of neutrino masses to 4.0 eV (95% CL), finds no preference for time or scale dependent modifications to the metric potentials, and is consistent with flatness and no running of the spectral index. The analysis code is public at https://github.com/sjoudaki/kids450

We present the tomographic cross-correlation between galaxy lensing measured in the Kilo Degree Survey (KiDS-450) with overlapping lensing measurements of the cosmic microwave background (CMB), as detected by Planck 2015. We compare our joint probe measurement to the theoretical expectation for a flat $\Lambda$CDM cosmology, assuming the best-fitting cosmological parameters from the KiDS-450 cosmic shear and Planck CMB analyses. We find that our results are consistent within $1\sigma$ with the KiDS-450 cosmology, with an amplitude re-scaling parameter $A_{\rm KiDS} = 0.86 \pm 0.19$. Adopting a Planck cosmology, we find our results are consistent within $2\sigma$, with $A_{\it Planck} = 0.68 \pm 0.15$. We show that the agreement is improved in both cases when the contamination to the signal by intrinsic galaxy alignments is accounted for, increasing $A$ by $\sim 0.1$. This is the first tomographic analysis of the galaxy lensing -- CMB lensing cross-correlation signal, and is based on five photometric redshift bins. We use this measurement as an independent validation of the multiplicative shear calibration and of the calibrated source redshift distribution at high redshifts. We find that constraints on these two quantities are strongly correlated when obtained from this technique, which should therefore not be considered as a stand-alone competitive calibration tool.

Superluminous supernovae are beginning to be discovered at redshifts as early as the epoch of reionization. A number of candidate mechanisms is reviewed, together with the discovery programs.

We present a new training set for estimating empirical photometric redshifts of galaxies, which was created as part of the 2dFLenS project. This training set is located in a 700 sq deg area of the KiDS South field and is randomly selected and nearly complete at r<19.5. We investigate the photometric redshift performance obtained with ugriz photometry from VST-ATLAS and W1/W2 from WISE, based on several empirical and template methods. The best redshift errors are obtained with kernel-density estimation, as are the lowest biases, which are consistent with zero within statistical noise. The 68th percentiles of the redshift scatter for magnitude-limited samples at r<(15.5, 17.5, 19.5) are (0.014, 0.017, 0.028). In this magnitude range, there are no known ambiguities in the colour-redshift map, consistent with a small rate of redshift outliers. In the fainter regime, the KDE method produces p(z) estimates per galaxy that represent unbiased and accurate redshift frequency expectations. The p(z) sum over any subsample is consistent with the true redshift frequency plus Poisson noise. Further improvements in redshift precision at r<20 would mostly be expected from filter sets with narrower passbands to increase the sensitivity of colours to small changes in redshift.

We develop a statistical estimator to infer the redshift probability distribution of a photometric sample of galaxies from its angular cross-correlation in redshift bins with an overlapping spectroscopic sample. This estimator is a minimum variance weighted quadratic function of the data: a quadratic estimator. This extends and modifies the methodology presented by McQuinn & White (2013). The derived source redshift distribution is degenerate with the source galaxy bias, which must be constrained via additional assumptions. We apply this estimator to constrain source galaxy redshift distributions in the Kilo-Degree imaging survey through cross-correlation with the spectroscopic 2-degree Field Lensing Survey, presenting results first as a binned step-wise distribution in the range z < 0.8, and then building a continuous distribution using a Gaussian process model. We demonstrate the robustness of our methodology using mock catalogues constructed from N-body simulations, and comparisons with other techniques for inferring the redshift distribution.

We investigate the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyze these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the nonlinear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favor of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck datasets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another, and that the level of concordance between the CFHTLenS and Planck datasets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favored by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited

We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

We present the Red-sequence Cluster Lensing Survey (RCSLenS), an application of the methods developed for the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) to the ~785deg$^2$, multi-band imaging data of the Red-sequence Cluster Survey 2 (RCS2). This project represents the largest public, sub-arcsecond seeing, multi-band survey to date that is suited for weak gravitational lensing measurements. With a careful assessment of systematic errors in shape measurements and photometric redshifts we extend the use of this data set to allow cross-correlation analyses between weak lensing observables and other data sets. We describe the imaging data, the data reduction, masking, multi-colour photometry, photometric redshifts, shape measurements, tests for systematic errors, and a blinding scheme to allow for more objective measurements. In total we analyse 761 pointings with r-band coverage, which constitutes our lensing sample. Residual large-scale B-mode systematics prevent the use of this shear catalogue for cosmic shear science. The effective number density of lensing sources over an unmasked area of 571.7deg$^2$ and down to a magnitude limit of r~24.5 is 8.1 galaxies per arcmin$^2$ (weighted: 5.5 arcmin$^{-2}$) distributed over 14 patches on the sky. Photometric redshifts based on 4-band griz data are available for 513 pointings covering an unmasked area of 383.5 deg$^2$ We present weak lensing mass reconstructions of some example clusters as well as the full survey representing the largest areas that have been mapped in this way. All our data products are publicly available through CADC at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/community/rcslens/query.html in a format very similar to the CFHTLenS data release.

We present the 2-degree Field Lensing Survey (2dFLenS), a new galaxy redshift survey performed at the Anglo-Australian Telescope. 2dFLenS is the first wide-area spectroscopic survey specifically targeting the area mapped by deep-imaging gravitational lensing fields, in this case the Kilo-Degree Survey. 2dFLenS obtained 70,079 redshifts in the range z < 0.9 over an area of 731 sq deg, and is designed to extend the datasets available for testing gravitational physics and promote the development of relevant algorithms for joint imaging and spectroscopic analysis. The redshift sample consists first of 40,531 Luminous Red Galaxies (LRGs), which enable analyses of galaxy-galaxy lensing, redshift-space distortion, and the overlapping source redshift distribution by cross-correlation. An additional 28,269 redshifts form a magnitude-limited (r < 19.5) nearly-complete sub-sample, allowing direct source classification and photometric-redshift calibration. In this paper, we describe the motivation, target selection, spectroscopic observations, and clustering analysis of 2dFLenS. We use power spectrum multipole measurements to fit the redshift-space distortion parameter of the LRG sample in two redshift ranges 0.15 < z < 0.43 and 0.43 < z < 0.7 as beta = 0.49 +/- 0.15 and beta = 0.26 +/- 0.09, respectively. These values are consistent with those obtained from LRGs in the Baryon Oscillation Spectroscopic Survey. 2dFLenS data products will be released via our website http://2dflens.swin.edu.au.

We present an optimized variant of the halo model, designed to produce accurate matter power spectra well into the non-linear regime for a wide range of cosmological models. To do this, we introduce physically motivated free parameters into the halo-model formalism and fit these to data from high-resolution $N$-body simulations. For a variety of $\Lambda$CDM and $w$CDM models the halo-model power is accurate to $\simeq 5$ per cent for $k\leq 10h\,\mathrm{Mpc}$ and $z\leq 2$. An advantage of our new halo model is that it can be adapted to account for the effects of baryonic feedback on the power spectrum. We demonstrate this by fitting the halo model to power spectra from the OWLS hydrodynamical simulation suite via parameters that govern halo internal structure. We are able to fit all feedback models investigated at the 5 per cent level using only two free parameters, and we place limits on the range of these halo parameters for feedback models investigated by the OWLS simulations. Accurate predictions to high-$k$ are vital for weak lensing surveys, and these halo parameters could be considered nuisance parameters to marginalize over in future analyses to mitigate uncertainty regarding the details of feedback. Finally, we investigate how lensing observables predicted by our model compare to those from simulations and from HALOFIT for a range of $k$-cuts and feedback models and quantify the angular scales at which these effects become important. Code to calculate power spectra from the model presented in this paper can be found at https://github.com/alexander-mead/hmcode

The unknown nature of dark energy motivates continued cosmological tests of large-scale gravitational physics. We present a new consistency check based on the relative amplitude of non-relativistic galaxy peculiar motions, measured via redshift-space distortion, and the relativistic deflection of light by those same galaxies traced by galaxy-galaxy lensing. We take advantage of the latest generation of deep, overlapping imaging and spectroscopic datasets, combining the Red Cluster Sequence Lensing Survey (RCSLenS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the WiggleZ Dark Energy Survey and the Baryon Oscillation Spectroscopic Survey (BOSS). We quantify the results using the "gravitational slip" statistic E_G, which we estimate as 0.48 +/- 0.10 at z=0.32 and 0.30 +/- 0.07 at z=0.57, the latter constituting the highest redshift at which this quantity has been determined. These measurements are consistent with the predictions of General Relativity, for a perturbed Friedmann-Robertson-Walker metric in a Universe dominated by a cosmological constant, which are E_G = 0.41 and 0.36 at these respective redshifts. The combination of redshift-space distortion and gravitational lensing data from current and future galaxy surveys will offer increasingly stringent tests of fundamental cosmology.

We present measurements of both scale- and time-dependent deviations from the standard gravitational field equations. These late-time modifications are introduced separately for relativistic and non-relativistic particles, by way of the parameters $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ using two bins in both scale and time, with transition wavenumber $0.01$ Mpc$^{-1}$ and redshift 1. We emphasize the use of two dynamical probes to constrain this set of parameters, galaxy power spectrum multipoles and the direct peculiar velocity power spectrum, which probe fluctuations on different scales. The multipole measurements are derived from the WiggleZ and BOSS Data Release 11 CMASS galaxy redshift surveys and the velocity power spectrum is measured from the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with additional cosmological probes including baryon acoustic oscillations, Type Ia SNe, the cosmic microwave background (CMB), lensing of the CMB, and the temperature--galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we find the inferred best-fit parameter values of $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ to be consistent with the standard model at the $95\%$ confidence level. Furthermore, accounting for the Alcock-Paczynski effect, we perform joint fits for the expansion history and growth index gamma; we measure $\gamma = 0.665 \pm 0.0669$ ($68\%$ C.L) for a fixed expansion history, and $\gamma = 0.73^{+0.08}_{-0.10}$ ($68\%$ C.L) when the expansion history is allowed to deviate from $\Lambda$CDM. With a fixed expansion history the inferred value is consistent with GR at the $95\%$ C.L; alternatively, a $2\sigma$ tension is observed when the expansion history is not fixed, this tension is worsened by the combination of growth and SNe data.

We present novel statistical tools to cross-correlate frequency cleaned thermal Sunyaev-Zel'dovich (tSZ) maps and tomographic weak lensing (wl) convergence maps. Moving beyond the lowest order cross-correlation, we introduce a hierarchy of mixed higher-order statistics, the cumulants and cumulant correlators, to analyze non-Gaussianity in real space, as well as corresponding polyspectra in the harmonic domain. Using these moments, we derive analytical expressions for the joint two-point probability distribution function (2PDF) for smoothed tSZ (y_s) and convergence (\kappa_s) maps. The presence of tomographic information allows us to study the evolution of higher order {\em mixed} tSZ-weak lensing statistics with redshift. We express the joint PDFs p_{\kappa y}(\kappa_s,y_s) in terms of individual one-point PDFs (p_{\kappa}(\kappa_s), p_y(y_s)) and the relevant bias functions (b_{\kappa}(\kappa_s), b_y(y_s)). Analytical results for two different regimes are presented that correspond to the small and large angular smoothing scales. Results are also obtained for corresponding {\em hot spots} in the tSZ and convergence maps. In addition to results based on hierarchical techniques and perturbative methods, we present results of calculations based on the lognormal approximation. The analytical expressions derived here are generic and applicable to cross-correlation studies of arbitrary tracers of large scale structure including e.g. that of tSZ and soft X-ray background.

From a combination of probes including the cosmic microwave background (WMAP7+SPT), Hubble constant (HST), baryon acoustic oscillations (SDSS+2dFGRS), and supernova distances (Union2), we have explored the extent to which the constraints on the effective number of neutrinos and sum of neutrino masses are affected by our ignorance of other cosmological parameters, including the curvature of the universe, running of the spectral index, primordial helium abundance, evolving late-time dark energy, and early dark energy. In a combined analysis of the effective number of neutrinos and sum of neutrino masses, we find mild (2.2 sigma) preference for additional light degrees of freedom. However, the effective number of neutrinos is consistent with the canonical expectation of 3 massive neutrinos and no extra relativistic species to within 1 sigma when allowing for evolving dark energy and relaxing the strong inflation prior on the curvature and running. The agreement improves with the possibility of an early dark energy component, itself constrained to be less than 5% of the critical density (95% CL) in our expanded parameter space. In extensions of the standard cosmological model, the derived amplitude of linear matter fluctuations sigma_8 is found to closely agree with low-redshift cluster abundance measurements. The sum of neutrino masses is robust to assumptions of the effective number of neutrinos, late-time dark energy, curvature, and running at the level of 1.2 eV (95% CL). The upper bound degrades to 2.0 eV (95% CL) when further including the early dark energy density and primordial helium abundance as additional free parameters. Even in extended cosmological parameter spaces, Planck alone could determine the possible existence of extra relativistic species at 4 sigma confidence and constrain the sum of neutrino masses to 0.2 eV (68% CL).

We find that the viability of a cosmological model that incorporates 2 sterile neutrinos with masses around 1 eV each, as favored by global neutrino oscillation analyses including short baseline results, is significantly dependent on the choice of datasets included in the analysis and the ability to control the systematic uncertainties associated with these datasets. Our analysis includes a variety of cosmological probes including the cosmic microwave background (WMAP7+SPT), Hubble constant (HST), galaxy power spectrum (SDSS-DR7), and supernova distances (SDSS and Union2 compilations). In the joint observational analysis, our sterile neutrino model is equally favored as a LCDM model when using the MLCS light curve fitter for the supernova measurements, and strongly disfavored by the data at \Delta\chi^2 ~ 18 when using the SALT2 fitter. When excluding the supernova measurements, the sterile neutrino model is disfavored by the other datasets at \Delta\chi^2 ~ 12, and at best becomes mildly disfavored at \Delta\chi^2 ~ 3 when allowing for curvature, evolving dark energy, additional relativistic species, running of the spectral index, and freedom in the primordial helium abundance. No single additional parameter accounts for most of this effect. Therefore, if laboratory experiments continue to favor a scenario with roughly eV mass sterile neutrinos, and if this becomes decisively disfavored by cosmology, then a more exotic cosmological model than explored here may become necessary.

At high angular frequencies, beyond the damping tail of the primary power spectrum, the dominant contribution to the power spectrum of cosmic microwave background (CMB) temperature fluctuations is the thermal Sunyaev-Zel'dovich (tSZ) effect. We investigate various important statistical properties of the Sunyaev-Zel'dovich maps, using well-motivated models for dark matter clustering to construct statistical descriptions of the tSZ effect to all orders enabling us to determine the entire probability distribution function (PDF). Any generic deterministic biasing scheme can be incorporated in our analysis and the effects of projection, biasing and the underlying density distribution can be analysed separately and transparently in this approach. We introduce the cumulant correlators as tools to analyse tSZ catalogs and relate them to corresponding statistical descriptors of the underlying density distribution. The statistics of hot spots in frequency-cleaned tSZ maps are also developed in a self-consistent way to an arbitrary order, to obtain results complementary to those found using the halo model. We also consider different beam sizes, to check the extent to which the PDF can be extracted from various observational configurations. The formalism is presented with two specific models for underlying matter clustering: (1) the hierarchical ansatz; and (2) the lognormal distribution. We find both models to be in very good agreement with the simulation results, though the lognormal model has an edge over the hierarchical model.

We forecast the expected cosmological constraints from a combination of probes of both the universal expansion rate and matter perturbation growth, in the form of weak lensing tomography, galaxy tomography, supernovae, and the cosmic microwave background incorporating all cross-correlations between the observables for an extensive cosmological parameter set. We allow for non-zero curvature and parameterize our ignorance of the early universe by allowing for a non-negligible fraction of dark energy (DE) at high redshifts. We find that early DE density can be constrained to 0.2% of the critical density of the universe with Planck combined with a ground-based LSST-like survey, while curvature can be constrained to 0.06%. However, these additional degrees of freedom degrade our ability to measure late-time dark energy and the sum of neutrino masses. We find that the combination of cosmological probes can break degeneracies and constrain the sum of neutrino masses to 0.04 eV, present DE density also to 0.2% of the critical density, and the equation of state to 0.01 - roughly a factor of two degradation in the constraints overall compared to the case without allowing for early DE. The constraints for a space-based mission are similar. Even a modest 1% dark energy fraction of the critical density at high redshift, if not accounted for in future analyses, biases the cosmological parameters by up to 2 sigma. Our analysis suggests that throwing out nonlinear scales (multipoles > 1000) may not result in significant degradation in future parameter measurements when multiple cosmological probes are combined. We find that including cross-correlations between the different probes can result in improved constraints by up to a factor of 2 for the sum of neutrino masses and early dark energy density.

Primordial non-Gaussianity is a crucial test of inflationary cosmology. We consider the impact of non-Gaussianity on the ionization power spectrum from 21 cm emission during the epoch of reionization. We focus on the power spectrum on large scales at redshifts of 7 to 8 and explore the expected constraint on the local non-Gaussianity parameter f_NL for current and next-generation 21 cm experiments. We show that experiments such as SKA and MWA could measure f_NL values of order 10. This can be improved by an order of magnitude with a fast-Fourier transform telescope like Omniscope.

The Sunyaev-Zel'dovich (SZ) effect has a distinct spectral signature that allows its separation from fluctuations in the cosmic microwave background (CMB) and foregrounds. Using CMB anisotropies measured in Wilkinson Microwave Anisotropy Probe's five-year maps, we constrain the SZ fluctuations at large, degree angular scales corresponding to multipoles in the range from 10 to 400. We provide upper bounds on SZ fluctuations at multipoles greater than 50, and find evidence for a hemispherically asymmetric signal at ten degrees angular scales. The amplitude of the detected signal cannot be easily explained with the allowed number density and temperature of electrons in the Galactic halo. We have failed to explain the excess signal as a residual from known Galactic foregrounds or instrumental uncertainties such as 1/f-noise.

At high angular frequencies the thermal Sunyaev-Zel'dovich (tSZ) effect constitutes the dominant signal in the CMB sky. The tSZ effect is caused by large scale pressure fluctuations in the baryonic distribution in the Universe so its statistical properties provide estimates of corresponding properties of the projected 3D pressure fluctuations. It's power spectrum is a sensitive probe of the density fluctuations, and the bispectrum can be used to separate the bias associated with pressure. The bispectrum is often probed with a one-point real-space analogue, the skewness. In addition to the skewness the morphological properties, as probed by the well known Minkowski Functionals (MFs), also require the generalized one-point statistics, which at the lowest order are identical to the skewness parameters. The concept of generalized skewness parameters can be extended to define a set of three associated generalized skew-spectra. We use these skew-spectra to probe the morphology of the tSZ sky or the y-sky. We show how these power spectra can be recovered from the data in the presence of arbitrary mask and noise templates using the well known Pseudo-Cl (PCL) approach for arbitrary beam shape. We also employ an approach based on the halo model to compute the tSZ bispectrum. The bispectrum from each of these models is then used to construct the generalized skew-spectra. We consider the performance of an all-sky survey with Planck-type noise and compare the results against a noise-free ideal experiment using a range of smoothing angles. We find that the skew-spectra can be estimated with very high signal-to-noise ratio from future frequency cleaned tSZ maps that will be available from experiments such as Planck. This will allow their mode by mode estimation for a wide range of angular frequencies and will help us to differentiate them from various other sources of non-Gaussianity.

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We consider corrections to the underlying cosmology due to secondary contributions from weak gravitational lensing, the integrated Sachs-Wolfe effect, and the Sunyaev-Zel'dovich effect contained in the trispectrum. We incorporate these additional contributions to the covariance of a binned angular power spectrum of temperature anisotropies in the analysis of current and prospective data sets. Although recent experiments such as ACBAR and CBI are not particularly sensitive to these additional non-Gaussian effects, the interpretation of Planck and CMBPol anisotropy spectra will require an accounting of non-Gaussian covariance leading to a degradation in cosmological parameter estimates by up to 20% and 30%, respectively.

We examine the influence of percent-level dark energy corrections to the nonlinear matter power spectrum on constraints of the dark energy equation of state from future weak lensing probes. We explicitly show that a poor approximation (off by > 10%) to the nonlinear corrections causes a > 1 sigma bias on the determination of the dark energy equation of state. Future weak lensing surveys must therefore incorporate dark energy modifications to the nonlinear matter power spectrum accurate to the percent-level, to avoid introducing significant bias in their measurements. For the WMAP5 cosmology, the more accurate power spectrum is more sensitive to dark energy properties, resulting in a factor of two improvement in dark energy equation of state constraints. We explore the complementary constraints on dark energy from future weak lensing and supernova surveys. A space-based, JDEM-like survey measures the equation of state in five independent redshift bins to ~10%, while this improves to ~5% for a wide-field ground-based survey like LSST. These constraints are contingent upon our ability to control weak lensing systematic uncertainties to the sub-percent level.

Our ignorance of the dark energy is generally described by a two-parameter equation of state. In these approaches a particular {\it ad hoc} functional form is assumed, and only two independent parameters are incorporated. We propose a model-independent, multi-parameter approach to fitting the dark energy, and show that next-generation surveys will constrain the equation of state in three or more independent redshift bins to better than 10%. Future knowledge of the dark energy will surpass two numbers (e.g., [$w_0$,$w_1$] or [$w_0$,$w_a$]), and we propose a more flexible approach to the analysis of present and future data.