We consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich (SZ) power spectrum and cluster counts analyses, highlighting subtle dependencies on the parameterisation and data combination. In $\Lambda$CDM, with an X-ray mass bias corresponding the expected hydrostatic mass bias, i.e., $b\simeq0.2$, our constraints from Planck SZ data are consistent with the latest results from SPT, DES-Y1 and KiDS+VIKING-450. In $\nu\Lambda$CDM, without prior information on $b$, our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint, as well as $(1-b)=0.64\pm0.04$~(68\%~CL). For forecasts, we find that a combination of a mock cosmic variance limited SZ power spectrum with primary CMB and BAO can improve the uncertainty on the total neutrino mass by 14\% with respect to CMB combined with BAO. This requires masking the heaviest clusters and probing the small-scale SZ power spectrum up to multipoles of $\ell_\mathrm{max}=10^4$. For a future competitive measurement of the total neutrino mass using CMB and the SZ power spectrum, but excluding BAO and lensing power spectrum, we find that a 1\% precision on the mass calibration is needed. Although SZ power spectrum-based measurements of the neutrino masses are challenging, we find that the SZ power spectrum can be used to tightly constrain intra-cluster medium properties.

We forecast the sensitivity of thirty-five different combinations of future Cosmic Microwave Background and Large Scale Structure data sets to cosmological parameters and to the total neutrino mass. We work under conservative assumptions accounting for uncertainties in the modelling of systematics. In particular, for galaxy redshift surveys, we remove the information coming from non-linear scales. We use Bayesian parameter extraction from mock likelihoods to avoid Fisher matrix uncertainties. Our grid of results allows for a direct comparison between the sensitivity of different data sets. We find that future surveys will measure the neutrino mass with high significance and will not be substantially affected by potential parameter degeneracies between neutrino masses, the density of relativistic relics, and a possible time-varying equation of state of Dark Energy.

Theoretical uncertainties on non-linear scales are among the main obstacles to exploit the sensitivity of forthcoming galaxy and hydrogen surveys like Euclid or the Square Kilometre Array (SKA). Here, we devise a new method to model the theoretical error that goes beyond the usual cut-off on small scales. The advantage of this more efficient implementation of the non-linear uncertainties is tested through a Markov-Chain-Monte-Carlo (MCMC) forecast of the sensitivity of Euclid and SKA to the parameters of the standard $\Lambda$CDM model, including massive neutrinos with total mass $M_\nu$, and to 3 extended scenarios, including 1) additional relativistic degrees of freedom ($\Lambda$CDM + $M_\nu$ + $N_\mathrm{eff}$), 2) a deviation from the cosmological constant ($\Lambda$CDM + $M_\nu$ + $w_0$), and 3) a time-varying dark energy equation of state parameter ($\Lambda$CDM + $M_\nu$ + $\left(w_0,w_a \right)$). We compare the sensitivity of 14 different combinations of cosmological probes and experimental configurations. For Euclid combined with Planck, assuming a plain cosmological constant, our method gives robust predictions for a high sensitivity to the primordial spectral index $n_{\rm s}$ ($\sigma(n_s)=0.00085$), the Hubble constant $H_0$ ($\sigma(H_0)=0.141 \, {\rm km/s/Mpc}$), the total neutrino mass $M_\nu$ ($\sigma(M_\nu)=0.020 \, {\rm eV}$). Assuming dynamical dark energy we get $\sigma(M_\nu)=0.030 \, {\rm eV}$ for the mass and $(\sigma(w_0), \sigma(w_a)) = (0.0214, 0.071)$ for the equation of state parameters. The predicted sensitivity to $M_\nu$ is mostly stable against the extensions of the cosmological model considered here. Interestingly, a significant improvement of the constraints on the extended model parameters is also obtained when combining Euclid with a low redshift HI intensity mapping survey by SKA1, demonstrating the importance of the synergy of Euclid and SKA.

In this work, we discuss the relation between the strength of the self-interaction of dark matter particles and the predicted properties of the inner density distributions of dark matter haloes. We present the results of $N$-body simulations for 28 haloes performed with the same initial conditions for cold dark matter and for self-interacting dark matter, for a range of cross-sections. We provide a simple phenomenological description of these results and compare with the semi-analytical model typically used in the literature. Using these results, we then predict how the inner dark matter surface density and core radius should depend on the self-interaction cross-section for observed haloes.

In cosmologies with massive neutrinos, the galaxy bias defined with respect to the total matter field (cold dark matter, baryons, and non-relativistic neutrinos) depends on the sum of the neutrino masses $M_{\nu}$, and becomes scale-dependent even on large scales. This effect has been usually neglected given the sensitivity of current surveys, but becomes a severe systematic for future surveys aiming to provide the first detection of non-zero $M_{\nu}$. The effect can be corrected for by defining the bias with respect to the density field of cold dark matter and baryons instead of the total matter field. In this work, we provide a simple prescription for correctly mitigating the neutrino-induced scale-dependent bias effect in a practical way. We clarify a number of subtleties regarding how to properly implement this correction in the presence of redshift-space distortions and non-linear evolution of perturbations. We perform a MCMC analysis on simulated galaxy clustering data that match the expected sensitivity of the \textit{Euclid} survey. We find that the neutrino-induced scale-dependent bias can lead to important shifts in both the inferred mean value of $M_{\nu}$, as well as its uncertainty. We show how these shifts propagate to other cosmological parameters correlated with $M_{\nu}$, such as the cold dark matter physical density $\Omega_{cdm} h^2$ and the scalar spectral index $n_s$. In conclusion, we find that correctly accounting for the neutrino-induced scale-dependent bias will be of crucial importance for future galaxy clustering analyses. We encourage the cosmology community to correctly account for this effect using the simple prescription we present in our work. The tools necessary to easily correct for the neutrino-induced scale-dependent bias will be made publicly available in an upcoming release of the Boltzmann solver \texttt{CLASS}.

MontePython is a parameter inference package for cosmology. We present the latest development of the code over the past couple of years. We explain, in particular, two new ingredients both contributing to improve the performance of Metropolis-Hastings sampling: an adaptation algorithm for the jumping factor, and a calculation of the inverse Fisher matrix, which can be used as a proposal density. We present several examples to show that these features speed up convergence and can save many hundreds of CPU-hours in the case of difficult runs, with a poor prior knowledge of the covariance matrix. We also summarise all the functionalities of MontePython in the current release, including new likelihoods and plotting options.

We consider a recently proposed model in which dark matter interacts with a thermal background of dark radiation. Dark radiation consists of relativistic degrees of freedom which allow larger values of the expansion rate of the universe today to be consistent with CMB data ($H_0$-problem). Scattering between dark matter and radiation suppresses the matter power spectrum at small scales and can explain the apparent discrepancies between $\Lambda$CDM predictions of the matter power spectrum and direct measurements of Large Scale Structure LSS ($\sigma_8$-problem). We go beyond previous work in two ways: 1. we enlarge the parameter space of our previous model and allow for an arbitrary fraction of the dark matter to be interacting and 2. we update the data sets used in our fits, most importantly we include LSS data with full $k$-dependence to explore the sensitivity of current data to the shape of the matter power spectrum. We find that LSS data prefer models with overall suppressed matter clustering due to dark matter - dark radiation interactions over $\Lambda$CDM at 3-4 $\sigma$. However recent weak lensing measurements of the power spectrum are not yet precise enough to clearly distinguish two limits of the model with different predicted shapes for the linear matter power spectrum. In two Appendices we give a derivation of the coupled dark matter and dark radiation perturbation equations from the Boltzmann equation in order to clarify a confusion in the recent literature, and we derive analytic approximations to the solutions of the perturbation equations in the two physically interesting limits of all dark matter weakly interacting or a small fraction of dark matter strongly interacting.

We perform dark-matter-only simulations of 28 relaxed massive cluster-sized haloes for Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, to study structural differences between the models at large radii, where the impact of baryonic physics is expected to be very limited. We find that the distributions for the radial profiles of the density, ellipsoidal axis ratios, and velocity anisotropies ($\beta$) of the haloes differ considerably between the models (at the $\sim1\sigma$ level), even at $\gtrsim10\%$ of the virial radius, if the self-scattering cross section is $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$. Direct comparison with observationally inferred density profiles disfavours SIDM for $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, but in an intermediate radial range ($\sim3\%$ of the virial radius), where the impact of baryonic physics is uncertain. At this level of the cross section, we find a narrower $\beta$ distribution in SIDM, clearly skewed towards isotropic orbits, with no SIDM (90\% of CDM) haloes having $\beta>0.12$ at $7\%$ of the virial radius. We estimate that with an observational sample of $\sim30$ ($\sim10^{15}$ M$_\odot$) relaxed clusters, $\beta$ can potentially be used to put competitive constraints on SIDM, once observational uncertainties improve by a factor of a few. We study the suppression of the memory of halo assembly history in SIDM clusters. For $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, we find that this happens only in the central halo regions ($\sim1/4$ of the scale radius of the halo), and only for haloes that assembled their mass within this region earlier than a formation redshift $z_f\sim2$. Otherwise, the memory of assembly remains and is reflected in ways similar to CDM, albeit with weaker trends.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

The present day universe consists of galaxies, galaxy clusters, one-dimensional filaments and two-dimensional sheets or pancakes, all of which combine to form the cosmic web. The so called "Zeldovich pancakes", are very difficult to observe, because their overdensity is only slightly greater than the average density of the universe. Falco et al (2014) presented a method to identify Zeldovich pancakes in observational data, and these were used as a tool for estimating the mass of galaxy clusters. Here we expand and refine that observational detection method. We study two pancakes on scales of 10 Mpc, identified from spectroscopically observed galaxies near the Coma cluster, and compare with twenty numerical pancakes. We find that the observed structures have velocity dispersion about 100 km/sec, which is relatively low compared to typical groups and filaments. These velocity dispersions are consistent with those found for the numerical pancakes. We also confirm that the identified structures are in fact two-dimensional structures. Finally, we estimate the stellar to total mass of the observational pancakes to be $2 \times 10^{-4}$, within one order of magnitude, which is smaller than that of clusters of galaxies.

The interplay between cosmology and earth based experiments is crucial in order to pin down neutrino physics. Indeed cosmology can provide very tight, yet model dependent, constraints on some neutrino properties. Here we focus on the neutrino mass sum, reviewing the up to date current bounds and showing the results of our forecast of the sensitivity of future experiments. Finally, we discuss the case for sterile neutrinos, explaining how non standard sterile neutrino self-interactions can reconcile the oscillation anomalies with cosmology.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

Future Cosmic Microwave Background experiments together with upcoming galaxy and 21-cm surveys will provide extremely accurate measurements of different cosmological observables located at different epochs of the cosmic history. The new data will be able to constrain the neutrino mass sum with the best precision ever. In order to exploit the complementarity of the different redshift probes, a deep understanding of the physical effects driving the impact of massive neutrinos on CMB and large scale structures is required. The goal of this work is to describe these effects, assuming a summed neutrino mass close to its minimum allowed value. We find that parameter degeneracies can be removed by appropriate combinations, leading to robust and model independent constraints. A joint forecast of the sensitivity of Euclid and DESI surveys together with a CORE-like CMB experiment leads to a $1\sigma$ uncertainty of $14$~meV on the summed neutrino mass. However this particular combination gives rise to a peculiar degeneracy between $M_\nu$ and the optical depth at reionization. Independent constraints from 21-cm surveys can break this degeneracy and decrease the $1\sigma$ uncertainty down to $12$~meV.

The mass measurement of galaxy clusters is an important tool for the determination of cosmological parameters describing the matter and energy content of the Universe. However, the standard methods rely on various assumptions about the shape or the level of equilibrium of the cluster. We present a novel method of measuring cluster masses. It is complementary to most of the other methods, since it only uses kinematical information from outside the virialized cluster. Our method identifies objects, as galaxy sheets or filaments, in the cluster outer region, and infers the cluster mass by modeling how the massive cluster perturbs the motion of the structures from the Hubble flow. At the same time, this technique allows to constrain the three-dimensional orientation of the detected structures with a good accuracy. We use a cosmological numerical simulation to test the method. We then apply the method to the Coma cluster, where we find two galaxy sheets, and measure the mass of Coma to be Mvir=(9.2\pm2.4)10^{14} Msol, in good agreement with previous measurements obtained with the standard methods.

We report the discovery of a unique gravitational lens system, SDSSJ2222+2745, producing five spectroscopically confirmed images of a z_s=2.82 quasar lensed by a foreground galaxy cluster at z_l=0.49. We also present photometric and spectroscopic evidence for a sixth lensed image of the same quasar. The maximum separation between the quasar images is 15.1". Both the large image separations and the high image multiplicity of the lensed quasar are in themselves exceptionally rare, and observing the combination of these two factors is an exceptionally unlikely occurrence in present datasets. This is only the third known case of a quasar lensed by a cluster, and the only one with six images. The lens system was discovered in the course of the Sloan Giant Arcs Survey, in which we identify candidate lenses in the Sloan Digital Sky Survey and target these for follow up and verification with the 2.56m Nordic Optical Telescope. Multi-band photometry obtained over multiple epochs from September 2011 to September 2012 reveal significant variability at the ~10-30% level in some of the quasar images, indicating that measurements of the relative time delay between quasar images will be feasible. In this lens system we also identify a bright (g = 21.5) giant arc corresponding to a strongly lensed background galaxy at z_s=2.30. We fit parametric models of the lens system, constrained by the redshift and positions of the quasar images and the redshift and position of the giant arc. The predicted time delays between different pairs of quasar images range from ~100 days to ~6 years.