The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

Unbiased and precise mass calibration of galaxy clusters is crucial to fully exploit galaxy clusters as cosmological probes. Stacking of weak lensing signal allows us to measure observable-mass relations down to less massive halos halos without extrapolation. We propose a Bayesian inference method to constrain the intrinsic scatter of the mass proxy in stacked analyses. The scatter of the stacked data is rescaled with respect to the individual scatter based on the number of binned clusters. We apply this method to the galaxy clusters detected with the AMICO (Adaptive Matched Identifier of Clustered Objects) algorithm in the third data release of the Kilo-Degree Survey. The results confirm the optical richness as a low scatter mass proxy. Based on the optical richness and the calibrated weak lensing mass-richness relation, mass of individual objects down to ~10^13 solar masses can be estimated with a precision of ~20 per cent.

Relic galaxies are thought to be the progenitors of high-redshift red nuggets that for some reason missed the channels of size growth and evolved passively and undisturbed since the first star formation burst (at $z>2$). These local ultracompact old galaxies are unique laboratories for studying the star formation processes at high redshift and thus the early stage of galaxy formation scenarios. Counterintuitively, theoretical and observational studies indicate that relics are more common in denser environments, where merging events predominate. To verify this scenario, we compared the number counts of a sample of ultracompact massive galaxies (UCMGs) selected within the third data release of the Kilo Degree Survey, that is, systems with sizes $R_{\rm e} < 1.5 \, \rm kpc$ and stellar masses $M_{\rm \star} > 8 \times 10^{10}\, \rm M_{\odot}$, with the number counts of galaxies with the same masses but normal sizes in field and cluster environments. Based on their optical and near-infrared colors, these UCMGs are likely to be mainly old, and hence representative of the relic population. We find that both UCMGs and normal-size galaxies are more abundant in clusters and their relative fraction depends only mildly on the global environment, with denser environments penalizing the survival of relics. Hence, UCMGs (and likely relics overall) are not special because of the environment effect on their nurture, but rather they are just a product of the stochasticity of the merging processes regardless of the global environment in which they live.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

Since the discovery of the accelerated cosmic expansion, one of the most important tasks in observational cosmology is to determine the nature of the dark energy. We should build our understanding on a minimum of assumptions in order to avoid biases from assumed cosmological models. The two most important functions describing the evolution of the universe and its structures are the expansion function E(a) and the linear growth factor D_+(a). The expansion function has been determined in previous papers in a model-independent way using distance moduli to type-Ia supernovae and assuming only a metric theory of gravity, spatial isotropy and homogeneity. Here, we extend this analysis in three ways: (1) We extend the data sample by combining the Pantheon measurements of type-Ia supernovae with measurements of baryonic acoustic oscillations; (2) we substantially simplify and generalise our method for reconstructing the expansion function; and (3) we use the reconstructed expansion function to determine the linear growth factor of cosmic structures, equally independent of specific assumptions on an underlying cosmological model other than the usual spatial symmetries. We show that the result is quite insensitive to the initial conditions for solving the growth equation, leaving the present-day matter-density parameter {\Omega}_m0 as the only relevant parameter for an otherwise purely empirical and accurate determination of the growth factor.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

We present the mass calibration for galaxy clusters detected with the AMICO code in KiDS DR3 data. The cluster sample comprises $\sim$ 7000 objects and covers the redshift range 0.1 < $z$ < 0.6. We perform a weak lensing stacked analysis by binning the clusters according to redshift and two different mass proxies provided by AMICO, namely the amplitude $A$ (measure of galaxy abundance through an optimal filter) and the richness $\lambda^*$ (sum of membership probabilities in a consistent radial and magnitude range across redshift). For each bin, we model the data as a truncated NFW profile plus a 2-halo term, taking into account uncertainties related to concentration and miscentring. From the retrieved estimates of the mean halo masses, we construct the $A$-$M_{200}$ and the $\lambda^*$-$M_{200}$ relations. The relations extend over more than one order of magnitude in mass, down to $M_{200} \sim 2 (5) \times 10^{13} M_\odot/h$ at $z$ = 0.2 (0.5), with small evolution in redshift. The logarithmic slope is $\sim 2.0$ for the $A$-mass relation, and $\sim 1.7$ for the $\lambda^*$-mass relation, consistent with previous estimations on mock catalogues and coherent with the different nature of the two observables.

Supernovae (SNe) could be powerful probes of the properties of stars and galaxies at high redshifts in future surveys. Wide fields and longer exposure times are required to offset diminishing star formation rates and lower fluxes to detect useful numbers of events at high redshift. In principle, the Large Synoptic Survey Telescope (LSST) could discover large numbers of early SNe because of its wide fields but only at lower redshifts because of its AB mag limit of ~ 24. But gravitational lensing by galaxy clusters and massive galaxies could boost flux from ancient SNe and allow LSST to detect them at earlier times. Here, we calculate detection rates for lensed SNe at z ~ 5 - 7 for LSST. We find that the LSST Wide Deep Fast survey could detect up to 120 lensed Population (Pop) I and II SNe but no lensed Pop III SNe. Deep-drilling programs in a single 10 square degree FoV could detect Pop I and II core-collapse SNe and Pop III pair-instability SNe at AB magnitudes of 27 - 28 and 26, respectively. An alternative deep survey over 80 nights with a one-year cadence could find ~ 8 Pop III SNe.

We present the latest version of the ray-tracing simulation code Skylens, which can be used to develop image simulations that reproduce strong lensing observations by any mass distribution with a high level of realism. Improvements of the code with respect to previous versions include the implementation of the multi-lens plane formalism, the use of denoised source galaxies from the Hubble eXtreme Deep Field, and the simulation of substructures in lensed arcs and images, based on a morphological analysis of bright nearby galaxies. Skylens can simulate observations with virtually any telescope. We present examples of space- and ground-based observations of a galaxy cluster through the Wide Field Channel on the Advanced Camera for Surveys of the Hubble Space Telescope, the Near Infrared Camera of the James Webb Space Telescope, the Wide Field Imager of the Wide Field Infrared Survey Telescope, the Hyper Suprime Camera of the Subaru telescope, and the Visible Imaging Channel of the Euclid space mission.

We present the first catalogue of galaxy cluster candidates derived from the third data release of the Kilo Degree Survey (KiDS-DR3). The sample of clusters has been produced using the Adaptive Matched Identifier of Clustered Objects (AMICO) algorithm. In this analysis AMICO takes advantage of the luminosity and spatial distribution of galaxies only, not considering colours. In this way, we prevent any selection effect related to the presence or absence of the red-sequence in the clusters. The catalogue contains 7988 candidate galaxy clusters in the redshift range 0.1<z<0.8 down to S/N>3.5 with a purity approaching 95% over the entire redshift range. In addition to the catalogue of galaxy clusters we also provide a catalogue of galaxies with their probabilistic association to galaxy clusters. We quantify the sample purity, completeness and the uncertainties of the detection properties, such as richness, redshift, and position, by means of mock galaxy catalogues derived directly from the data. This preserves their statistical properties including photo-z uncertainties, unknown absorption across the survey, missing data, spatial correlation of galaxies and galaxy clusters. Being based on the real data, such mock catalogues do not have to rely on the assumptions on which numerical simulations and semi-analytic models are based on. This paper is the first of a series of papers in which we discuss the details and physical properties of the sample presented in this work.

We present EasyCritics, an algorithm to detect strongly-lensing groups and clusters in wide-field surveys without relying on a direct recognition of arcs. EasyCritics assumes that light traces mass in order to predict the most likely locations of critical curves from the observed fluxes of luminous red early-type galaxies in the line of sight. The positions, redshifts and fluxes of these galaxies constrain the idealized gravitational lensing potential as a function of source redshift up to five free parameters, which are calibrated on few known lenses. From the lensing potential, EasyCritics derives the critical curves for a given, representative source redshift. The code is highly parallelized, uses fast Fourier methods and, optionally, GPU acceleration in order to process large datasets efficiently. The search of a $\smash{1 \, \mathrm{deg}^2}$ field of view requires less than 1 minute on a modern quad-core CPU, when using a pixel resolution of $0.25''/\mathrm{px}$. In this first part of a paper series on EasyCritics, we describe the main underlying concepts and present a first demonstration on data from the Canada-France-Hawaii-Telescope Lensing Survey. We show that EasyCritics is able to identify known group- and cluster-scale lenses, including a cluster with two giant arc candidates that were previously missed by automated arc detectors.

We report the results of $EasyCritics$, a fully automated algorithm for the efficient search of strong-lensing (SL) regions in wide-field surveys, applied to the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). By using only the photometric information of the brightest elliptical galaxies distributed over a wide redshift range ($\smash{0.2 \lesssim z \lesssim 0.9}$) and without requiring the identification of arcs, our algorithm produces lensing potential models and catalogs of critical curves of the entire survey area. We explore several parameter set configurations in order to test the efficiency of our approach. In a specific configuration, $EasyCritics$ generates only $\sim1200$ possibly super-critical regions in the CFHTLS area, drastically reducing the effective area for inspection from $154$ sq. deg to $\sim0.623$ sq. deg, $i.e.$ by more than two orders of magnitude. Among the pre-selected SL regions, we identify 32 of the 44 previously known lenses on the group and cluster scale, and discover 9 new promising lens candidates. The detection rate can be easily improved to $\sim82\%$ by a simple modification in the parameter set, but at the expense of increasing the total number of possible SL candidates. Note that $EasyCritics$ is fully complementary to other arc-finders since we characterize lenses instead of directly identifying arcs. Although future comparisons against numerical simulations are required for fully assessing the efficiency of $EasyCritics$, the algorithm seems very promising for upcoming surveys covering $\smash{10^{4}}$ sq. deg, such as the $Euclid$ mission and $LSST$, where the pre-selection of candidates for any kind of SL analysis will be indispensable due to the expected enormous data volume.

We present AMICO (Adaptive Matched Identifier of Clustered Objects), a new algorithm for the detection of galaxy clusters in photometric surveys. AMICO is based on the Optimal Filtering technique, which allows to maximise the signal-to-noise ratio of the clusters. In this work we focus on the new iterative approach to the extraction of cluster candidates from the map produced by the filter. In particular, we provide a definition of membership probability for the galaxies close to any cluster candidate, which allows us to remove its imprint from the map, allowing the detection of smaller structures. As demonstrated in our tests, this method allows the deblending of close-by and aligned structures in more than $50\%$ of the cases for objects at radial distance equal to $0.5 \times R_{200}$ or redshift distance equal to $2 \times \sigma_z$, being $\sigma_z$ the typical uncertainty of photometric redshifts. Running AMICO on mocks derived from N-body simulations and semi-analytical modelling of the galaxy evolution, we obtain a consistent mass-amplitude relation through the redshift range $0.3 < z < 1$, with a logarithmic slope $\sim 0.55$ and a logarithmic scatter $\sim 0.14$. The fraction of false detections is steeply decreasing with S/N, and negligible at S/N > 5.

In this paper, we present the tools used to search for galaxy clusters in the Kilo Degree Survey (KiDS), and our first results. The cluster detection is based on an implementation of the optimal filtering technique that enables us to identify clusters as over-densities in the distribution of galaxies using their positions on the sky, magnitudes, and photometric redshifts. The contamination and completeness of the cluster catalog are derived using mock catalogs based on the data themselves. The optimal signal to noise threshold for the cluster detection is obtained by randomizing the galaxy positions and selecting the value that produces a contamination of less than 20%. Starting from a subset of clusters detected with high significance at low redshifts, we shift them to higher redshifts to estimate the completeness as a function of redshift: the average completeness is ~ 85%. An estimate of the mass of the clusters is derived using the richness as a proxy. We obtained 1858 candidate clusters with redshift 0 < z_c < 0.7 and mass 13.5 < log(M500/Msun) < 15 in an area of 114 sq. degrees (KiDS ESO-DR2). A comparison with publicly available Sloan Digital Sky Survey (SDSS)-based cluster catalogs shows that we match more than 50% of the clusters (77% in the case of the redMaPPer catalog). We also cross-matched our cluster catalog with the Abell clusters, and clusters found by XMM and in the Planck-SZ survey; however, only a small number of them lie inside the KiDS area currently available.

According to the theory of general relativity, masses deflect light in a way similar to convex glass lenses. This gravitational lensing effect is astigmatic, giving rise to image distortions. These distortions allow to quantify cosmic structures statistically on a broad range of scales, and to map the spatial distribution of dark and visible matter. We summarise the theory of weak gravitational lensing and review applications to galaxies, galaxy clusters and larger-scale structures in the Universe.

We present a method to model optical images of galaxies using Expectation Maximization Principal Components Analysis (EMPCA). The method relies on the data alone and does not assume any pre-established model or fitting formula. It preserves the statistical properties of the sample, minimizing possible biases. The precision of the reconstructions appears to be suited for photometric, morphological and weak lensing analysis, as well as the realization of mock astronomical images. Here, we put some emphasis on the latter because weak gravitational lensing is entering a new phase in which systematics are becoming the major source of uncertainty. Accurate simulations are necessary to perform a reliable calibration of the ellipticity measurements on which the final bias depends. As a test case, we process $7038$ galaxies observed with the ACS/WFC stacked images of the Hubble eXtreme Deep Field (XDF) and measure the accuracy of the reconstructions in terms of their moments of brightness, which turn out to be comparable to what can be achieved with well-established weak-lensing algorithms.

The statistics of peaks in weak gravitational lensing maps is a promising technique to constrain cosmological parameters in present and future surveys. Here we investigate its power when using general extreme value statistics which is very sensitive to the exponential tail of the halo mass function. To this end, we use an analytic method to quantify the number of weak lensing peaks caused by galaxy clusters, large-scale structures and observational noise. Doing so, we further improve the method in the regime of high signal-to-noise ratios dominated by non-linear structures by accounting for the embedding of those counts into the surrounding shear caused by large scale structures. We derive the extreme value and order statistics for both over-densities (positive peaks) and under-densities (negative peaks) and provide an optimized criterion to split a wide field survey into sub-fields in order to sample the distribution of extreme values such that the expected objects causing the largest signals are mostly due to galaxy clusters. We find good agreement of our model predictions with a ray-tracing $N$-body simulation. For a Euclid-like survey, we find tight constraints on $\sigma_8$ and $\Omega_\text{m}$ with relative uncertainties of $\sim 10^{-3}$. In contrast, the equation of state parameter $w_0$ can be constrained only with a $10\%$ level, and $w_\text{a}$ is out of reach even if we include redshift information.

We present the first reported case of the simultaneous metallicity determination of a gamma-ray burst (GRB) host galaxy, from both afterglow absorption lines as well as strong emission-line diagnostics. Using spectroscopic and imaging observations of the afterglow and host of the long-duration Swift GRB121024A at z = 2.30, we give one of the most complete views of a GRB host/environment to date. We observe a strong damped Ly-alpha absorber (DLA) with a hydrogen column density of log N(HI) = 21.88 +/- 0.10, H2 absorption in the Lyman-Werner bands (molecular fraction of log(f)~ -1.4; fourth solid detection of molecular hydrogen in a GRB-DLA), the nebular emission lines H-alpha, H-beta, [O II], [O III] and [N II], as well as metal absorption lines. We find a GRB host galaxy that is highly star-forming (SFR ~ 40 solar masses/yr ), with a dust-corrected metallicity along the line of sight of [Zn/H]corr = -0.6 +/- 0.2 ([O/H] ~ -0.3 from emission lines), and a depletion factor [Zn/Fe] = 0.85 +/- 0.04. The molecular gas is separated by 400 km/s (and 1-3 kpc) from the gas that is photoexcited by the GRB. This implies a fairly massive host, in agreement with the derived stellar mass of log(M/M_solar ) = 9.9+/- 0.2. We dissect the host galaxy by characterising its molecular component, the excited gas, and the line-emitting star-forming regions. The extinction curve for the line of sight is found to be unusually flat (Rv ~15). We discuss the possibility of an anomalous grain size distributions. We furthermore discuss the different metallicity determinations from both absorption and emission lines, which gives consistent results for the line of sight to GRB 121024A.

We present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (\Sigma)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

Strong gravitational lensing provides fundamental insights into the understanding of the dark matter distribution in massive galaxies, galaxy clusters and the background cosmology. Despite their importance, the number of gravitational arcs discovered so far is small. The urge for more complete, large samples and unbiased methods of selecting candidates is rising. A number of methods for the automatic detection of arcs have been proposed in the literature, but large amounts of spurious detections retrieved by these methods forces observers to visually inspect thousands of candidates per square degree in order to clean the samples. This approach is largely subjective and requires a huge amount of eye-ball checking, especially considering the actual and upcoming wide field surveys, which will cover thousands of square degrees. In this paper we study the statistical properties of colours of gravitational arcs detected in the 37 deg^2 of the CARS survey. We have found that most of them lie in a relatively small region of the (g'-r',r'-i') colour-colour diagram. To explain this property, we provide a model which includes the lensing optical depth expected in a LCDM cosmology that, in combination with the sources' redshift distribution of a given survey, in our case CARS, peaks for sources at redshift z~1. By further modelling the colours derived from the SED of the galaxies dominating the population at that redshift, the model well reproduces the observed colours. By taking advantage of the colour selection suggested by both data and model, we show that this multi-band filtering returns a sample 83% complete and a contamination reduced by a factor of ~6.5 with respect to the single-band arcfinder sample. New arc candidates are also proposed.

We perform a model independent reconstruction of the cosmic expansion rate based on type Ia supernova data. Using the Union 2.1 data set, we show that the Hubble parameter behaviour allowed by the data without making any hypothesis about cosmological model or underlying gravity theory is consistent with a flat LCDM universe having H_0 = 70.43 +- 0.33 and Omega_m=0.297 +- 0.020, weakly dependent on the choice of initial scatter matrix. This is in closer agreement with the recently released Planck results (H_0 = 67.3 +- 1.2, Omega_m = 0.314 +- 0.020) than other standard analyses based on type Ia supernova data. We argue this might be an indication that, in order to tackle subtle deviations from the standard cosmological model present in type Ia supernova data, it is mandatory to go beyond parametrized approaches.

According to the current standard model of Cosmology, matter in the Universe arranges itself along a network of filamentary structure. These filaments connect the main nodes of this so-called 'Cosmic Web', which are clusters of galaxies. Although its large-scale distribution is clearly characterized by numerical simulations, constraining the dark matter content of the cosmic web in reality turns out to be difficult. The natural method of choice is gravitational lensing. However, the direct detection and mapping of the elusive filament signal is challenging and in this work we present two methods,specifically tailored to achieve this task. A linear matched filter aims at the detection of the smooth mass component of filaments and is optimized to perform a shear decomposition that follows the anisotropic component of the lensing signal. Filaments clearly inherit this property due to their morphology. At the same time, the contamination arising from the central massive cluster is controlled in a natural way. The filament 1 {\sigma} detection is of about {\kappa} ~ 0.01-0.005 according to the filter's template width and length, enabling the detection of structures out of reach with other approaches. The second, complementary method seeks to detect the clumpy component of filaments. The detection is determined by the number density of sub-clump identifications in an area enclosing the potential filament, as it was found within the observed field with the filter approach. We test both methods against Mock observations based on realistic N-Body simulations of filamentary structure and prove the feasibility of detecting filaments with ground-based data.

In this work we investigate the problem concerning the presence of additional bodies gravitationally bounded with the WASP-3 system. We present eight new transits of this planet and analyse all the photometric and radial velocity data published so far. We did not observe significant periodicities in the Fourier spectrum of the observed minus calculated (O-C) transit timing and radial velocity diagrams (the highest peak having false-alarm probabilities of 56 per cent and 31 per cent, respectively) or long-term trends. Combining all the available information, we conclude that the radial velocity and transit timing techniques exclude, at 99 per cent confidence limit, any perturber more massive than M \gtrsim 100 M_Earth with periods up to 10 times the period of the inner planet. We also investigate the possible presence of an exomoon on this system and determined that considering the scatter of the O-C transit timing residuals a coplanar exomoon would likely produce detectable transits. This hypothesis is however apparently ruled out by observations conducted by other researchers. In case the orbit of the moon is not coplanar the accuracy of our transit timing and transit duration measurements prevents any significant statement. Interestingly, on the basis of our reanalysis of SOPHIE data we noted that WASP-3 passed from a less active (log R'_hk=-4.95) to a more active (log R'_hk=-4.8) state during the 3 yr monitoring period spanned by the observations. Despite no clear spot crossing has been reported for this system, this analysis claims for a more intensive monitoring of the activity level of this star in order to understand its impact on photometric and radial velocity measurements.

We studied the effect of primordial non-Gaussianity with varied bispectrum shapes on the number counts of signal-to-noise peaks in wide field cosmic shear maps. The two cosmological contributions to this particular weak lensing statistic, namely the chance projection of Large Scale Structure and the occurrence of real, cluster-sized dark matter halos, have been modeled semi-analytically, thus allowing to easily introduce the effect of non-Gaussian initial conditions. We performed a Fisher matrix analysis by taking into account the full covariance of the peak counts in order to forecast the joint constraints on the level of primordial non-Gaussianity and the amplitude of the matter power spectrum that are expected by future wide field imaging surveys. We find that positive-skewed non-Gaussianity increases the number counts of cosmic shear peaks, more so at high signal-to-noise values, where the signal is mostly dominated by massive clusters as expected. The increment is at the level of ~1 for f_NL=10 and ~10 for f_NL=100 for a local shape of the primordial bispectrum, while different bispectrum shapes give generically a smaller effect. For a future survey on the model of the proposed ESA space mission Euclid and by avoiding the strong assumption of being capable to distinguish the weak lensing signal of galaxy clusters from chance projection of Large Scale Structures we forecasted a 1-sigma error on the level of non-Gaussianity of ~30-40 for the local and equilateral models, and of ~100-200 for the less explored enfolded and orthogonal bispectrum shapes.

We define an optimal basis system into which cosmological observables can be decomposed. The basis system can be optimised for a specific cosmological model or for an ensemble of models, even if based on drastically different physical assumptions. The projection coefficients derived from this basis system, the so-called features, provide a common parameterisation for studying and comparing different cosmological models independently of their physical construction. They can be used to directly compare different cosmologies and study their degeneracies in terms of a simple metric separation. This is a very convenient approach, since only very few realisations have to be computed, in contrast to Markov-Chain Monte Carlo methods. Finally, the proposed basis system can be applied to reconstruct the Hubble expansion rate from supernova luminosity distance data with the advantage of being sensitive to possible unexpected features in the data set. We test the method both on mock catalogues and on the SuperNova Legacy Survey data set.

Galaxy clusters are usually detected in blind optical surveys via suitable filtering methods. We present an optimal matched filter which maximizes their signal-to-noise ratio by taking advantage of the knowledge we have of their intrinsic physical properties and of the data noise properties. In this paper we restrict our application to galaxy magnitudes, positions and photometric redshifts if available, and we also apply the filter separately to weak lensing data. The method is suitable to be naturally extended to a multi-band approach which could include not only additional optical bands but also observables with different nature such as X-rays. For each detection, the filter provides its significance, an estimate for the richness and for the redshift even if photo-z are not given. The provided analytical error estimate is tested against numerical simulations. We finally apply our method to the COSMOS field and compare the results with previous cluster detections obtained with different methods. Our catalogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including both optical and weak-lensing information.

Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal SZ signal, and compare them to X-ray cluster detections. Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the simulation volume to identify true and spurious detections. Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations with the characteristics of ACT detect clusters with masses M>6-9e13 M_o/h, quite independent of redshift, reach 50% completeness at ~1e14 M_o/h and 100% completeness at ~2e14 M_o/h. Samples are contaminated by only a few per cent spurious detections. This is broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and X-ray maps confirms that the filters introduce at most a negligible bias.

We construct two linear filtering techniques based on weak gravitational lensing to constrain the inner slope $\alpha$ of the density profile of dark-matter halos. Both methods combine all available information into an estimate of this single number. Under idealised assumptions, $\alpha$ is constrained to ~15% if the halo concentration c is known, and to < 30% if not. We argue that the inevitable degeneracy between density-profile slope and halo concentration cannot be lifted under realistic conditions, and show by means of Fisher-matrix methods which linear combination of $\alpha$ and c is best constrained by our filtering of lensing data. This defines a new parameter, called P1, which is then constrained to ~15% for a single massive halo. If the signals of many halos can be stacked, their density profiles should thus be well constrained by the linear filters proposed here with the advantagebe insensitive to the cluster substructures.

We develop and apply an analytic method to predict peak counts in weak-lensing surveys. It is based on the theory of Gaussian random fields and suitable to quantify the level of spurious detections caused by chance projections of large-scale structures as well as the shape and shot noise contributed by the background galaxies. We compare our method to peak counts obtained from numerical ray-tracing simulations and find good agreement at the expected level. The number of peak detections depends substantially on the shape and size of the filter applied to the gravitational shear field. Our main results are that weak-lensing peak counts are dominated by spurious detections up to signal-to-noise ratios of 3--5 and that most filters yield only a few detections per square degree above this level, while a filter optimised for suppressing large-scale structure noise returns up to an order of magnitude more.

We seek to find a shapelet-based scheme for deconvolving galaxy images from the PSF which leads to unbiased shear measurements. Based on the analytic formulation of convolution in shapelet space, we construct a procedure to recover the unconvolved shapelet coefficients under the assumption that the PSF is perfectly known. Using specific simulations, we test this approach and compare it to other published approaches. We show that convolution in shapelet space leads to a shapelet model of order $n_{max}^h = n_{max}^g + n_{max}^f$ with $n_{max}^f$ and $n_{max}^g$ being the maximum orders of the intrinsic galaxy and the PSF models, respectively. Deconvolution is hence a transformation which maps a certain number of convolved coefficients onto a generally smaller number of deconvolved coefficients. By inferring the latter number from data, we construct the maximum-likelihood solution for this transformation and obtain unbiased shear estimates with a remarkable amount of noise reduction compared to established approaches. This finding is particularly valid for complicated PSF models and low $S/N$ images, which renders our approach suitable for typical weak-lensing conditions.

We study the reliability of dark-matter halo detections with three different linear filters applied to weak-lensing data. We use ray-tracing in the multiple lens-plane approximation through a large cosmological simulation to construct realizations of cosmic lensing by large-scale structures between redshifts zero and two. We apply the filters mentioned above to detect peaks in the weak-lensing signal and compare them with the true population of dark matter halos present in the simulation. We confirm the stability and performance of a filter optimized for suppressing the contamination by large-scale structure. It allows the reliable detection of dark-matter halos with masses above a few times 1e13 M_sun/h with a fraction of spurious detections below ~10%. For sources at redshift two, 50% of the halos more massive than ~7e13 M_sun/h are detected, and completeness is reached at ~2e14 M_sun/h.

Observations of the Cosmic Microwave Background (CMB) have revealed an unexpected quadrupole-octopole alignment along a preferred axis pointing toward the Virgo cluster. We here investigate whether this feature can be explained in the framework of the concordance model by secondary anisotropies produced by the non-linear evolution of the gravitational potential, the so-called Rees-Sciama (RS) effect. We focus on the effect caused by the local superclusters, which we calculate using a constrained high-resolution hydrodynamical simulation, based on the IRAS 1.2-Jy all-sky galaxy redshift survey, which reproduces the main structures of our Universe out to a distance of 110 Mpc from our Galaxy. The resulting RS effect peaks at low multipoles and has a minimum/maximum amplitude of -6.6\mu K 1.9\mu K. Even though its quadrupole is well aligned with the one measured for the CMB, its amplitude is not sufficient to explain the observed magnitude of the quadrupole/octopole alignment. In addition, we analyze the WMAP-3 data with a linear matched filter in an attempt to determine an upper limit for the RS signal amplitude on large scales. We found that it is possible to infer a weak upper limit of 30\mu K for its maximum amplitude.

We consider the role of the galactic kinetic Sunyaev Zeldovich (SZ) effect as a CMB foreground. While the galactic thermal Sunyaev Zeldovich effect has previously been studied and discarded as a potential CMB foreground, we find that the kinetic SZ effect is dominant in the galactic case. We analyse the detectability of the kinetic SZ effect by means of an optimally matched filter technique applied to a simulation of an ideal observation. We obtain no detection, getting a S/N ratio of 0.1, thereby demonstrating that the kinetic SZ effect can also safely be ignored as a CMB foreground. However we provide maps of the expected signal for inclusion in future high precision data processing. Furthermore, we rule out the significant contamination of the polarised CMB signal by second scattering of galactic kinetic Sunyaev-Zeldovich photons, since we show that the scattering of the CMB quadrupole photons by galactic electrons is a stronger effect than the Sunyaev Zeldovich second scattering, and has already been shown to produce no significant polarised contamination. We confirm the latter assessment also by means of an optimally matched filter.

The polarization sensitivity of the upcoming millimetric observatories will open new possibilities for studying the properties of galaxy clusters and for using them as powerful cosmological probes. For this reason it is necessary to investigate in detail the characteristics of the polarization signals produced by their highly ionized intra-cluster medium (ICM). This work is focussed on the polarization effect induced by the ICM bulk motions, the so-called kpSZ signal, which has an amplitude proportional to the optical depth and to the square of the tangential velocity. In particular we study how this polarization signal is affected by the internal dynamics of galaxy clusters and what is its dependence on the physical modelling adopted to describe the baryonic component. This is done by producing realistic kpSZ maps starting from the outputs of two different sets of high-resolution hydrodynamical N-body simulations. The first set (17 objects) follows only non-radiative hydrodynamics, while for each of 9 objects of the second set we implement four different kinds of physical processes. Our results shows that the kpSZ signal turns out to be a very sensitive probe of the dynamical status of galaxy clusters. We find that major merger events can amplify the signal up to one order of magnitude with respect to relaxed clusters, reaching amplitude up to about 100 nuK. This result implies that the internal ICM dynamics must be taken into account when evaluating this signal because simplicistic models, based on spherical rigid bodies, may provide wrong estimates. Finally we find that the dependence on the physical modelling of the baryonic component is relevant only in the very inner regions of clusters.

The mass function of galaxy clusters is an important cosmological probe. Differences in the selection method could potentially lead to biases when determining the mass function. From the optical and X-ray data of the XMM-Newton Follow-Up Survey, we obtained a sample of galaxy cluster candidates using weak gravitational lensing, the optical Postman matched filter method, and a search for extended X-ray sources. We developed our weak-lensing search criteria by testing the performance of the aperture mass statistic on realistic ray-tracing simulations matching our survey parameters and by comparing two filter functions. We find that the dominant noise source for our survey is shape noise at almost all significance levels and that spurious cluster detections due to projections of large-scale structures are negligible, except possibly for highly significantly detected peaks. Our full cluster catalog has 155 cluster candidates, 116 found with the Postman matched filter, 59 extended X-ray sources, and 31 shear selected potential clusters. Most of these cluster candidates were not previously known. The present catalog will be a solid foundation for studying possible selection effects in either method.

We apply the linear filter for the weak-lensing signal of dark-matter halos developed in Maturi et al. (2005) to the cosmic-shear data extracted from the Garching-Bonn-Deep-Survey (GaBoDS). We wish to search for dark-matter halos through weak-lensing signatures which are significantly above the random and systematic noise level caused by intervening large-scale structures. We employ a linear matched filter which maximises the signal-to-noise ratio by minimising the number of spurious detections caused by the superposition of large-scale structures (LSS). This is achieved by suppressing those spatial frequencies dominated by the LSS contamination. We confirm the improved stability and reliability of the detections achieved with our new filter compared to the commonly-used aperture mass (Schneider, 1996; Schneider et al., 1998) and to the aperture mass based on the shear profile expected for NFW haloes (see e.g. Schirmer et al., 2004; Hennawi & Spergel, 2005). Schirmer et al.~(2006) achieved results comparable to our filter, but probably only because of the low average redshift of the background sources in GaBoDS, which keeps the LSS contamination low. For deeper data, the difference will be more important, as shown by Maturi et al. (2005). We detect fourteen halos on about eighteen square degrees selected from the survey. Five are known clusters, two are associated with over-densities of galaxies visible in the GaBoDS image, and seven have no known optical or X-ray counterparts.

We discuss how to use the Rees-Sciama (RS) effect associated with merging clusters of galaxies to measure their kinematic properties. In a previous work (Rubino-Martin et al. 2004), the morphology and symmetries of the effect were examined by means of a simplified model. Here, we use realistic N-body simulations to better describe the effect, and to confirm that the signal has a characteristic quadrupole structure. From the amplitude of the signal obtained, we conclude that it is necessary to combine several cluster mergers in order to achieve a detection. Using the extended Press-Schechter formalism, we characterized the expected distribution of the parameters describing the mergers, and we used these results to generate realistic mock catalogues of cluster mergers. To optimize the extraction of the RS signal, we developed an extension of the spatial filtering method described in Haehnelt & Tegmark (1996). This extended filter has a general definition, so it can be applied in many other fields, such as gravitational lensing of the CMB or lensing of background galaxies. It has been applied to our mock catalogues, and we show that with the announced sensitivities of future experiments like the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) or the Atacama Large Millimeter Array (ALMA), a detection of the signal will be possible if we consider of the order of 1,000 cluster mergers.

We construct a linear filter optimised for detecting dark-matter halos in weak-lensing data. The filter assumes a mean radial profile of the halo shear pattern and modifies that shape by the noise power spectrum. Aiming at separating dark-matter halos from spurious peaks caused by large-scale structure lensing, we model the noise as being composed of weak lensing by large-scale structures and Poisson noise from random galaxy positions and intrinsic ellipticities. Optimal filtering against the noise requires the optimal filter scale to be smaller than typical halo sizes. Although a perfect separation of halos from spurious large-scale structure peaks is strictly impossible, we use numerical simulations to demonstrate that our filter produces substantially more sensitive, reliable and stable results than the conventionally used aperture-mass statistic.

We adapt a non-linear filter proposed by Hu 2001 for detecting lensing of the CMB by large-scale structures to recover surface-density profiles of galaxy clusters from their localised, weak gravitational lensing effect on CMB fields. Shifting the band-pass of the filter to smaller scales, and normalising it such as to reproduce the convergence rather than the deflection angle, we find that the mean density profile of a sample of 100 clusters can be recovered to better than 10% from well within the scale radius to almost the virial radius. The kinetic Sunyaev-Zel'dovich effect is shown to be a negligible source of error. We test the filter applying it to data simulated using the characteristics of the Atacama Cosmology Telescope (ACT), showing that it will be possible to recover mean cluster profiles outside a radius of 1' corresponding to ACT's angular resolution.