We present the first determination of the Hubble constant $H_0$ from strong lensing time delay data and type Ia supernova luminosity distances that is independent of the cosmological model. We also determine the spatial curvature model-independently. We assume that light propagation over long distances is described by the FLRW metric and geometrical optics holds, but make no assumption about the contents of the Universe or the theory of gravity on cosmological scales. We find $H_0=75.7^{+4.5}_{-4.4}$ km/s/Mpc and $\Omega_{K0}=0.12^{+0.27}_{-0.25}$. This is a 6\% determination of $H_0$. A weak prior from the cosmic microwave background on the distance to the last scattering surface improves this to $H_0=76.8^{+4.2}_{-3.8}$ km/s/Mpc and $\Omega_{K0}=0.18^{+0.25}_{-0.18}$. Assuming zero spatial curvature, we get $H_0=74.2^{+3.0}_{-2.9}$ km/s/Mpc, a precision of $4\%$. The measurements also provide a consistency test of the FLRW metric: we find no evidence against it.

Scattering events with compact objects are expected in the primordial black hole (PBH) cold dark matter (CDM) scenario due to close encounters between stars and PBH in the dense environments of dwarf spheroidals. We develop a Bayesian framework to search for correlations among Milky Way stellar trajectories and those of globular clusters and dwarf galaxies in the halo, and other nearby galaxies. We apply the method to a selection of hypervelocity stars (HVS) and globular clusters from Gaia DR2 catalog, and known nearby (mostly dwarf) galaxies with full phase-space and size measurements. We report positive evidence for trajectory intersection $\sim$20-40 Myr ago of up to 2 stars, depending on priors, with the Sagittarius dwarf Spheroidal (dSph) galaxy when assuming Marchetti et al. (2019) distance estimates. We verify that the result is compatible with their evolutionary status, setting a lower bound for the stellar age of $\sim$100 Myr. However, such scattering events are not confirmed when assuming Anders et al. (2019) distance estimates. We discuss shortcomings related to present data quality and future prospects for detection of HVS with the full Gaia catalog and Sagittarius dSph.

We compute the redshift-dependent angular bispectrum of galaxy number counts at tree-level, including nonlinear clustering bias and estimating numerically for the first time the effect of redshift space distortions (RSD). We show that for narrow redshift bins the amplitude of nonlinear RSD is comparable with the matter density perturbations. While our numerical results only include terms relevant on sub-horizon scales, the formalism can readily be extended to the full tree-level bispectrum. Our approach does not rely on the flat-sky approximation and it can be easily generalized to different sources by including the appropriate bias expansion. We test the accuracy of Limber approximation for different z-bins. We highlight the subtle but relevant differences in the angular bispectrum of galaxy number counts with respect to CMB, due to the different scale dependence of perturbations. Our formalism can also be directly applied to the angular HI intensity mapping bispectrum.

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.

We performed a detailed study of the evolution of the star formation rate (SFR) and stellar mass of the brightest group galaxies (BGGs) and their relative contribution to the total baryon budget within $R_{200}$ ($f^{BGG}_{b,200}$). The sample comprises 407 BGGs selected from X-ray galaxy groups ($M_{200}=10^{12.8}-10^{14} \;M_{\odot}$) out to $z\sim1.3$ identified in the COSMOS, XMM-LSS, and AEGIS fields. We find that BGGs constitute two distinct populations of quiescent and star-forming galaxies and their mean SFR is $\sim2$ dex higher than the median SFR at $ z<1.3 $. Both the mean and the median SFRs decline with time by $>2$ dex. The mean (median) of stellar mass of BGGs has grown by $0.3$ dex since $z=1.3$ to the present day. We show that up to $\sim45\% $ of the stellar mass growth in a star-forming BGG can be due to its star-formation activity. With respect to $f^{BGG}_{b,200}$, we find it to increase with decreasing redshift by $\sim0.35$ dex while decreasing with halo mass in a redshift dependent manner. We show that the slope of the relation between $f^{BGG}_{b,200}$ and halo mass increases negatively with decreasing redshift. This trend is driven by an insufficient star-formation in BGGs, compared to the halo growth rate. We separately show the BGGs with the 20\% highest $f^{BGG}_{b,200}$ are generally non-star-forming galaxies and grow in mass by processes not related to star formation (e.g., dry mergers and tidal striping). We present the $ M_\star-M_h $ and $ M_\star/M_h-M_h $ relations and compare them with semi-analytic model predictions and a number of results from the literature. We quantify the intrinsic scatter in stellar mass of BGGs at fixed halo mass ($\sigma_{log M_{\star}}$) and find that $\sigma_{log M_{\star}}$ increases from 0.3 dex at $ z\sim0.2 $ to 0.5 dex at $ z\sim1.0 $ due to the bimodal distribution of stellar mass.

We evaluate the effect of structure formation on the average expansion rate with a statistical treatment where density peaks and troughs are modelled as homogeneous ellipsoids. This extends earlier work that used spherical regions. We find that the shear and the presence of filamentary and planar structures have only a small impact on the results. The expansion rate times the age of the universe $Ht$ increases from 2/3 to 0.83 at late times, in order of magnitude agreement with observations, although the change is slower and takes longer than in the real universe. We discuss shortcomings that have to be addressed for this and similar statistical models in the literature to develop into realistic quantitative treatment of backreaction.

Cosmological models with Galileon gravity are an alternative to the standard $\Lambda {\rm CDM}$ paradigm with testable predictions at the level of its self-accelerating solutions for the expansion history, as well as large-scale structure formation. Here, we place constraints on the full parameter space of these models using data from the cosmic microwave background (CMB) (including lensing), baryonic acoustic oscillations (BAO) and the Integrated Sachs-Wolfe (ISW) effect. We pay special attention to the ISW effect for which we use the cross-spectra, $C_\ell^{\rm T g}$, of CMB temperature maps and foreground galaxies from the WISE survey. The sign of $C_\ell^{\rm T g}$ is set by the time evolution of the lensing potential in the redshift range of the galaxy sample: it is positive if the potential decays (like in $\Lambda {\rm CDM}$), negative if it deepens. We constrain three subsets of Galileon gravity separately known as the Cubic, Quartic and Quintic Galileons. The cubic Galileon model predicts a negative $C_\ell^{\rm T g}$ and exhibits a $7.8\sigma$ tension with the data, which effectively rules it out. For the quartic and quintic models the ISW data also rule out a significant region of the parameter space but permit regions where the goodness-of-fit is comparable to $\Lambda {\rm CDM}$. The data prefers a non zero sum of the neutrino masses ($\sum m_\nu\approx 0.5$eV) with $ \sim \! 5\sigma$ significance in these models. The best-fitting models have values of $H_0$ consistent with local determinations, thereby avoiding the tension that exists in $\Lambda {\rm CDM}$. We also identify and discuss a $\sim \! 2\sigma$ tension that Galileon gravity exhibits with recent BAO measurements. Our analysis shows overall that Galileon cosmologies cannot be ruled out by current data but future lensing, BAO and ISW data hold strong potential to do so.

If the FRW metric is a good approximation on large scales, then the distance and the expansion rate, as well different notions of distance, satisfy certain consistency conditions. We fit the JLA SNIa distance data to determine the expected amplitude of the violation of these conditions if accelerated expansion is due to backreaction. Adding cosmic clock and BAO expansion rate data, we also model-independently determine the current observational limits on such violation. We find that the predicted maximum backreaction amplitude $|k_H|\lesssim1$ (95% C.I.) is of the same order as the current observational constraints $|k_H|\lesssim1$, the precise numbers depending on the adopted fitting method (polynomials or splines) and stellar population evolution model. We also find that constraints on the value of $H_0$ determined from expansion rate data are sensitive to the stellar evolution model. We forecast constraints from projected LSST+Euclid-like SNIa plus Euclid galaxy differential age data. We find improvement by factor of 6 for the backreaction case and 3 for the model-independent case, probing an interesting region of possible signatures.

We demonstrate the importance of including the lensing contribution in galaxy clustering analyses with large galaxy redshift surveys. It is well known that radial cross-correlations between different redshift bins of galaxy surveys are dominated by lensing. But we show here that also neglecting lensing in the autocorrelations within one bin severely biases cosmological parameter estimation with redshift surveys. It leads to significant shifts for several cosmological parameters, most notably the scalar spectral index and the neutrino mass scale. Especially the latter parameter is one of the main targets of future galaxy surveys.

We address the impact of consistent modifications of gravity on the largest observable scales, focusing on relativistic effects in galaxy number counts and the cross-correlation between the matter large scale structure (LSS) distribution and the cosmic microwave background (CMB). Our analysis applies to a very broad class of general scalar-tensor theories encoded in the Horndeski Lagrangian and is fully consistent on linear scales, retaining the full dynamics of the scalar field and not assuming quasi-static evolution. As particular examples we consider self-accelerating Covariant Galileons, Brans-Dicke theory and parameterizations based on the effective field theory of dark energy, using the \hiclass\, code to address the impact of these models on relativistic corrections to LSS observables. We find that especially effects which involve integrals along the line of sight (lensing convergence, time delay and the integrated Sachs-Wolfe effect -- ISW) can be considerably modified, and even lead to $\mathcal{O}(1000\%)$ deviations from General Relativity in the case of the ISW effect for Galileon models, for which standard probes such as the growth function only vary by $\mathcal{O}(10\%)$. These effects become dominant when correlating galaxy number counts at different redshifts and can lead to $\sim 50\%$ deviations in the total signal that might be observable by future LSS surveys. Because of their integrated nature, these deep-redshift cross-correlations are sensitive to modifications of gravity even when probing eras much before dark energy domination. We further isolate the ISW effect using the cross-correlation between LSS and CMB temperature anisotropies and use current data to further constrain Horndeski models (abridged).

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We modified the CLASS code in order to include relativistic galaxy number counts in spatially curved geometries; we present the formalism and study the effect of relativistic corrections on spatial curvature. The new version of the code is now publicly available. Using a Fisher matrix analysis, we investigate how measurements of the spatial curvature parameter $\Omega_K$ with future galaxy surveys are affected by relativistic effects, which influence observations of the large scale galaxy distribution. These effects include contributions from cosmic magnification, Doppler terms and terms involving the gravitational potential. As an application, we consider angle and redshift dependent power spectra, which are especially well suited for model independent cosmological constraints. We compute our results for a representative deep, wide and spectroscopic survey, and our results show the impact of relativistic corrections on the spatial curvature parameter estimation. We show that constraints on the curvature parameter may be strongly biased if, in particular, cosmic magnification is not included in the analysis. Other relativistic effects turn out to be subdominant in the studied configuration. We analyze how the shift in the estimated best-fit value for the curvature and other cosmological parameters depends on the magnification bias parameter, and find that significant biases are to be expected if this term is not properly considered in the analysis.

We estimate the sensitivity of future galaxy surveys to cosmological parameters, using the redshift dependent angular power spectra of galaxy number counts, $C_\ell(z_1,z_2)$, calculated with all relativistic corrections at first order in perturbation theory. We pay special attention to the redshift dependence of the non-linearity scale and present Fisher matrix forecasts for Euclid-like and DES-like galaxy surveys. We compare the standard $P(k)$ analysis with the new $C_\ell(z_1,z_2)$ method. We show that for surveys with photometric redshifts the new analysis performs significantly better than the $P(k)$ analysis. For spectroscopic redshifts, however, the large number of redshift bins which would be needed to fully profit from the redshift information, is severely limited by shot noise. We also identify surveys which can measure the lensing contribution and we study the monopole, $C_0(z_1,z_2)$.

We present some accurate and efficient computations of large scale structure observables, obtained with a modified version of the CLASS code which is made publicly available. This code includes all relativistic corrections and computes both the power spectrum Cl(z1,z2) and the corresponding correlation function xi(theta,z1,z2) in linear perturbation theory. For Gaussian initial perturbations, these quantities contain the full information encoded in the large scale matter distribution at the level of linear perturbation theory. We illustrate the usefulness of our code for cosmological parameter estimation through a few simple examples.

We discuss the dominant terms of the relativistic galaxy number counts to second order in cosmological perturbation theory on sub-Hubble scales and on intermediate to large redshifts. In particular, we determine their contribution to the bispectrum. In addition to the terms already known from Newtonian second order perturbation theory, we find that there are a series of additional `lensing-like' terms which contribute to the bispectrum. We derive analytical expressions for the full leading order bispectrum and we evaluate it numerically for different configurations, indicating how they can be measured with upcoming surveys. In particular, the new `lensing-like' terms are not negligible for wide redshift bins and even dominate the bispectrum at well separated redshifts. This offers us the possibility to measure them in future surveys.

We investigate how the lensing potential can be measured tomographically with future galaxy surveys using their number counts. Such a measurement is an independent test of the standard $\Lambda$CDM framework and can be used to discern modified theories of gravity. We perform a Fisher matrix forecast based on galaxy angular-redshift power spectra, assuming specifications consistent with future photometric Euclid-like surveys and spectroscopic SKA-like surveys. For the Euclid-like survey we derive a fitting formula for the magnification bias. Our analysis suggests that the cross correlation between different redshift bins is very sensitive to the lensing potential such that the survey can measure the amplitude of the lensing potential at the same level of precision as other standard $\Lambda$CDM cosmological parameters.

We determine the number counts to second order in cosmological perturbation theory in the Poisson gauge and allowing for anisotropic stress. The calculation is performed using an innovative approach based on the recently proposed "geodesic light-cone" gauge. This allows us to determine the number counts in a purely geometric way, without using Einstein's equation. The result is valid for general dark energy models and (most) modified gravity models. We then evaluate numerically the relevant contributions to the number counts bispectrum. In particular we consider the terms involving the density, redshift space distortion and lensing.

We investigate the cosmological dependence and the constraining power of large-scale galaxy correlations, including all redshift-distortions, wide-angle, lensing and gravitational potential effects on linear scales. We analyze the cosmological information present in the lensing convergence and in the gravitational potential terms describing the so-called "relativistic effects," and we find that, while smaller than the information contained in intrinsic galaxy clustering, it is not negligible. We investigate how neglecting them does bias cosmological measurements performed by future spectroscopic and photometric large-scale surveys such as SKA and Euclid. We perform a Fisher analysis using the CLASS code, modified to include scale-dependent galaxy bias and redshift-dependent magnification and evolution bias. Our results show that neglecting relativistic terms introduces an error in the forecasted precision in measuring cosmological parameters of the order of a few tens of percent, in particular when measuring the matter content of the Universe and primordial non-Gaussianity parameters. Therefore, we argue that radial correlations and integrated relativistic terms need to be taken into account when forecasting the constraining power of future large-scale number counts of galaxy surveys.

We argue that from observations alone, only the transverse power spectrum $C_\ell(z_1,z_2)$ and the corresponding correlation function $\xi(\theta,z_1,z_2)$ can be measured and that these contain the full three dimensional information. We determine the two point galaxy correlation function at linear order in perturbation theory. Redshift space distortions are taken into account for arbitrary angular and redshift separations. We discuss the shape of the longitudinal and the transversal correlation functions which are very different from each other and from their real space counterpart. We then go on and suggest how to measure both, the Hubble parameter, $H(z)$, and the angular diameter distance, $D_A(z)$, separately from these correlation functions and perform an Alcock-Paczynski test.

The fitting formula for the location of the first acoustic peak in the matter power spectrum is revised. We discuss the physics that leads to baryon acoustic oscillations: the recombination history, the tight coupling approximation and the velocity overshoot effect. A new fitting formula is proposed, which is in accordance within 5% with numerical results for a suitable range of cosmological parameters, whereas previous results yield deviations of up to 20%. The crucial improvement turns out to be the accuracy of the recombination history.