The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

It is not well understood whether the growth of early-type cluster galaxies proceeds inside-out, outside-in, or at the same pace at all radii. In this work we measured the galaxy size, defined by the radius including 80\% of the galaxy light, non-parametrically. We also determined a non-parametric estimate of galaxy light concentration, which measures the curvature of the surface brightness profile in the galaxy outskirts. We used an almost random sampling of a mass-limited sample formed by 128 morphologically early-type galaxies in clusters with $\log M/M_{\odot} \ga 10.7$ spanning the wide range $0.17<z<1.81$. From these data we derived the size-mass and concentration-mass relations, as well as their evolution. At 80\% light radius, early-type galaxies in clusters are about 2.7 times larger than at 50\% radius at all redshifts, and close to de Vaucouleurs profiles in the last 10 Gyr. While between $z=2$ and $z=0$ both half-light and 80\% light sizes increase by a factor of $1.7$, concentration stays constant within $2$\%, that is to say the size growth of early-type galaxies in cluster environments proceeds at the same pace at both radii. Existing physical explanations proposed in the literature are inconsistent with our results, demonstrating the need for dedicated numerical simulations to identify the physical mechanism affecting the galaxy structure.

The thermal Sunyaev-Zeldovich (SZ) effect and the X-ray emission offer separate and highly complementary probes of the thermodynamics of the intracluster medium. We present JoXSZ, the first publicly available code designed to jointly fit SZ and X-ray data coming from various instruments to derive the thermodynamic profiles of galaxy clusters. JoXSZ follows a fully Bayesian forward-modelling approach, accounts for the SZ calibration uncertainty and X-ray background level systematic. It improves upon most state-of-the-art, and not publicly available, analyses because it adopts the correct Poisson-Gauss expression for the joint likelihood, makes full use of the information contained in the observations, even in the case of missing values within the datasets, has a more inclusive error budget, and adopts a consistent temperature across the various parts of the code, allowing for differences between X-ray and SZ gas mass weighted temperatures when required by the user. JoXSZ accounts for beam smearing and data analysis transfer function, accounts for the temperature and metallicity dependencies of the SZ and X-ray conversion factors, adopts flexible parametrization for the thermodynamic profiles, and on user request allows either adopting or relaxing the assumption of hydrostatic equilibrium (HE). When HE holds, JoXSZ uses a physical (positive) prior on the radial derivative of the enclosed mass and derives the mass profile and overdensity radii $r_\Delta$. For these reasons, JoXSZ goes beyond simple SZ and electron density fits. We illustrate the use of JoXSZ by combining Chandra and NIKA data on the high-redshift cluster CL J1226.9+3332. The code is written in Python, it is fully documented and the users are free to customize their analysis in accordance with their needs and requirements. JoXSZ is publicly available on GitHub.

High-resolution observations of the thermal Sunyaev-Zeldovich (SZ) effect and of the X-ray emission of galaxy clusters are becoming more and more widespread, offering us an unique asset to the study of the thermodynamic properties of the intracluster medium. We present JoXSZ, a Bayesian forward-modelling Python code designed to jointly fit the SZ data and the three dimensional X-ray data cube. JoXSZ is able to derive the thermodynamic profiles of galaxy clusters for the first time making full and consistent use of all the information contained in the observations. JoXSZ will be publicly available on GitHub in the near future.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Galaxy cluster analyses based on high-resolution observations of the Sunyaev-Zeldovich (SZ) effect have become common in the last decade. We present PreProFit, the first publicly available code designed to fit the pressure profile of galaxy clusters from SZ data. PreProFit is based on a Bayesian forward-modelling approach, allows the analysis of data coming from different sources, adopts a flexible parametrization for the pressure profile, and fits the model to the data accounting for Abel integral, beam smearing, and transfer function filtering. PreProFit is computationally efficient, is extensively documented, has been released as an open source Python project, and was developed to be part of a joint analysis of X-ray and SZ data on galaxy clusters. PreProFit returns $\chi^2$, model parameters and uncertainties, marginal and joint probability contours, diagnostic plots, and surface brightness radial profiles. PreProFit also allows the use of analytic approximations for the beam and transfer functions useful for feasibility studies.

Our knowledge of the variety of galaxy clusters has been increasing in the last few years thanks to our progress in understanding the severity of selection effects on samples. To understand the reason for the observed variety, we study CL2015, a cluster easily missed in X-ray selected observational samples. Its core-excised X-ray luminosity is low for its mass M500, well below the mean relation for an X-ray selected sample, but only ~1.5 sigma below that derived for an X-ray unbiased sample. We derived thermodynamic profiles and hydrostatic masses with the acquired deep Swift X-ray data, and we used archival Einstein, Planck, and SDSS data to derive additional measurements, such as integrated Compton parameter, total mass, and stellar mass. The pressure and the electron density profiles of CL2015 are systematically outside the +/- 2 sigma range of the universal profiles; in particular the electron density profile is even lower than the one derived from Planck-selected clusters. CL2015 also turns out to be fairly different in the X-ray luminosity versus integrated pressure scaling compared to an X-ray selected sample, but it is a normal object in terms of stellar mass fraction. CL2015's hydrostatic mass profile, by itself or when is considered together with dynamical masses, shows that the cluster has an unusual low concentration and an unusual sparsity compared to clusters in X-ray selected samples. The different behavior of CL2015 is caused by its low concentration. When concentration differences are accounted for, the properties of CL2015 become consistent with comparison samples. CL2015 is perhaps the first known cluster with a remarkably low mass concentration for which high quality X-ray data exist. Objects similar to CL2015 fail to enter observational X-ray selected samples because of their low X-ray luminosity relative to their mass.

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 study the molecular gas properties of two star-forming galaxies separated by 6 kpc in the projected space and belonging to a galaxy cluster selected from the Irac Shallow Cluster Survey, at a redshift $z=1.2$, i.e., $\sim2$ Gyr after the cosmic star formation density peak. This work describes the first CO detection from $1<z<1.4$ star forming cluster galaxies with no reported clear evidence of AGN. We exploit observations taken with the NOEMA interferometer at $\sim3$ mm to detect CO(2-1) line emission from the two selected galaxies, unresolved by our observations. Based on the CO(2-1) spectrum we estimate a total molecular gas mass $M({\rm H_2})=(2.2^{+0.5}_{-0.4})\times10^{10}$ $M_\odot$ and dust mass $M_{\rm dust}<4.2\times10^8~M_\odot$ for the two blended sources. The two galaxies have similar stellar masses and a large relative velocity of $\sim$400 km/s estimated from the CO(2-1) line width. These findings tend to privilege a scenario where both sources contribute to the observed CO(2-1). By using the archival Spitzer MIPS flux at 24$\mu$m we estimate an ${\rm SFR(24\mu m)}=(28^{+12}_{-8})~M_\odot$/yr for each of the two galaxies. Assuming that the two sources equally contribute to the observe CO(2-1) our analysis yields a depletion time scale $\tau_{\rm dep}=(3.9^{+1.4}_{-1.8})\times10^8$ yr, and a molecular gas to stellar mass ratio $0.17\pm0.13$ for each of two sources, separately. Our results are in overall agreement with those of other distant cluster galaxies. The two target galaxies have molecular gas mass and depletion time that are marginally compatible with, but smaller than those of main sequence field galaxies, suggesting that the molecular gas has not been refueled enough. Higher resolution and higher frequency observations will enable us to spatially resolve the two sources and possibly distinguish between different gas processing mechanisms.

[abridged] This work aims to observationally investigate the history of size growth of early-type galaxies and how the growth depends on cosmic epoch and the mass of the halo in which they are embedded. We carried out a photometric and structural analysis in the rest-frame $V$ band of a mass-selected ($\log M/M_\odot >10.7$) sample of red-sequence early-type galaxies with spectroscopic/grism redshift in the general field up to $z=2$ to complement a previous work presenting an identical analysis but in halos 100 times more massive and 1000 times denser. We homogeneously derived sizes (effective radii) fully accounting for the multi-component nature of galaxies and the common presence of isophote twists and ellipticity gradients. By using these mass-selected samples, composed of 170 red-sequence early-type galaxies in the general field and 224 identically selected and analyzed in clusters, we isolate the effect on galaxy sizes of the halo in which galaxies are embedded and its dependence on epoch. We find that the $\log$ of the galaxy size at a fixed stellar mass, $\log M/M_\odot= 11$, has increased with epoch at a rate twice as fast in the field than in cluster in the last 10 Gyr ($0.26\pm0.03$ versus $0.13\pm0.02$ dex per unit redshift). Red-sequence early-type galaxies in the general field reached the size of their cousins in denser environment by $z=0.25\pm0.13$ in spite of being three times smaller at $z\sim2$. Data point toward a model where size growth is epoch-independent (i.e., $\partial \log r_e /\partial z = c$), but with a rate $c$ depending on environment, $\partial c /\partial \log M_{halo} \approx 0.05$. Environment determines the growth rate ($d \log r_e / dz$) at all redshifts, indicating an external origin for the galaxy growth without any clear epoch where it ceases to have an effect.

All estimates of cluster mass have some intrinsic scatter and perhaps some bias with true mass even in the absence of measurement errors for example caused by cluster triaxiality and large scale structure. Knowledge of the bias and scatter values is fundamental for both cluster cosmology and astrophysics. In this paper we show that the intrinsic scatter of a mass proxy can be constrained by measurements of the gas fraction because masses with higher values of intrinsic scatter with true mass produce more scattered gas fractions. Moreover, the relative bias of two mass estimates can be constrained by comparing the mean gas fraction at the same (nominal) cluster mass. Our observational study addresses the scatter between caustic (i.e., dynamically estimated) and true masses, and the relative bias of caustic and hydrostatic masses. For these purposes, we used the X-ray Unbiased Cluster Sample, a cluster sample selected independently from the intracluster medium content with reliable masses: 34 galaxy clusters in the nearby ($0.050<z<0.135$) Universe, mostly with $14<\log M_{500}/M_\odot \lesssim 14.5$, and with caustic masses. We found a 35\% scatter between caustic and true masses. Furthermore, we found that the relative bias between caustic and hydrostatic masses is small, $0.06\pm0.05$ dex, improving upon past measurements. The small scatter found confirms our previous measurements of a highly variable amount of feedback from cluster to cluster, which is the cause of the observed large variety of core-excised X-ray luminosities and gas masses.

We use a cluster sample selected independently of the intracluster medium content with reliable masses to measure the mean gas mass fraction and its scatter, the biases of the X-ray selection on gas mass fraction, and the covariance between the X-ray luminosity and gas mass. The sample is formed by 34 galaxy clusters in the nearby ($0.050<z<0.135$) Universe, mostly with $14<\log M_{500}/M_\odot \lesssim 14.5$, and with masses calculated with the caustic technique. First, we found that integrated gas density profiles have similar shapes, extending earlier results based on subpopulations of clusters such as those that are relaxed or X-ray bright for their mass. Second, the X-ray unbiased selection of our sample allows us to unveil a variegate population of clusters; the gas mass fraction shows a scatter of $0.17\pm0.04$ dex, possibly indicating a quite variable amount of feedback from cluster to cluster, which is larger than is found in previous samples targeting subpopulations of galaxy clusters, such as relaxed or X-ray bright clusters. The similarity of the gas density profiles induces an almost scatterless relation between X-ray luminosity, gas mass, and halo mass, and modulates selection effects in the halo gas mass fraction: gas-rich clusters are preferentially included in X-ray selected samples. The almost scatterless relation also fixes the relative scatters and slopes of the $L_X-M$ and $M_{gas}-M$ relations and makes core-excised X-ray luminosities and gas masses fully covariant. Therefore, cosmological or astrophysical studies involving X-ray or SZ selected samples need to account for both selection effects and covariance of the studied quantities with X-ray luminosity/SZ strength.

We carried out a photometric and structural analysis in the rest-frame $V$ band of a mass-selected ($\log M/M_\odot >10.7$) sample of red-sequence galaxies in 14 galaxy clusters, 6 of which are at $z>1.45$. To this end, we reduced/analyzed about 300 orbits of multicolor images taken with the Advanced Camera for Survey and the Wide Field Camera 3 on the Hubble Space Telescope. We uniformly morphologically classified galaxies from $z=0.023$ to $z=1.803$, and we homogeneously derived sizes (effective radii) for the entire sample. Furthermore, our size derivation allows, and therefore is not biased by, the presence of the usual variety of morphological structures seen in early-type galaxies, such as bulges, bars, disks, isophote twists, and ellipiticy gradients. By using such a mass-selected sample, composed of 244 red-sequence early-type galaxies, we find that the $\log$ of the galaxy size at a fixed stellar mass, $\log M/M_\odot= 11$ has increased with time at a rate of $0.023\pm0.002$ dex per Gyr over the last 10 Gyr, in marked contrast with the threefold increase found in the literature for galaxies in the general field over the same period. This suggests, at face value, that secular processes should be excluded as the primary drivers of size evolution because we observed an environmental environmental dependent size growth. Using spectroscopic ages of Coma early-type galaxies we also find that recently quenched early-type galaxies are a numerically minor population not different enough in size to alter the mean size at a given mass, which implies that the progenitor bias is minor, i.e., that the size evolution measured by selecting galaxies at the redshift of observation is indistinguishable from the one that compares ancestors and descendents.

We present a catalog of galaxy cluster masses derived by exploiting the tight correlation between mass and richness, i.e., a properly computed number of bright cluster galaxies. The richness definition adopted in this work is properly calibrated, shows a small scatter with mass, and has a known evolution, which means that we can estimate accurate ($0.16$ dex) masses more precisely than by adopting any other richness estimates or X-ray or SZ-based proxies based on survey data. We measured a few hundred galaxy clusters at $0.05<z<0.22$ in the low-extinction part of the Sloan Digital Sky Survey footprint that are in the 2015 catalog of Planck-detected clusters, that have a known X-ray emission, that are in the Abell catalog, or that are among the most most cited in the literature. Diagnostic plots and direct images of clusters are individually inspected and we improved cluster centers and, when needed, we revised redshifts. Whenever possible, we also checked for indications of contamination from other clusters on the line of sight, and found ten such cases. All this information, with the derived cluster mass values, are included in the distributed value-added cluster catalog of the 275 clusters with a derived mass larger than $10^{14}$ M$_{\odot}$. A web front-end is available at http://www.brera.mi.astro.it/$\sim$andreon/famous.html . Finally, in a technical appendix we illustrate with Planck clusters how to minimize the sensitivity of comparisons between masses listed in different catalogs to the specific overlapping of the considered subsamples, a problem recognized

We aim to determine the intrinsic variety, at a given mass, of the properties of the intracluster medium in clusters of galaxies. This requires a cluster sample selected independently of the intracluster medium content for which reliable masses and subsequent X-ray data can be obtained. We present one such sample, consisting of 34 galaxy clusters selected independently of their X-ray properties in the nearby ($0.050<z<0.135$) Universe and mostly with $14<\log M_{500}/M_\odot \lesssim 14.5$, where masses are dynamically estimated. We collected the available X-ray observations from the archives and then observed the remaining clusters with the low-background Swift X-ray telescope, which is extremely useful for sampling a cluster population expected to have low surface brightness. We found that clusters display a large range (up to a factor 50) in X-ray luminosities within $r_{500}$ at a given mass, whether or not the central emission ($r<0.15 r_{500}$) is excised, unveiling a wider cluster population than seen in Sunayev-Zeldovich surveys or inferred from the population seen in X-ray surveys. The measured dispersion is $0.5$ dex in $L_X$ at a given mass.

Richness, i.e., the number of bright cluster galaxies, is known to correlate with the cluster mass, however, to exploit it as mass proxy we need a way to estimate the aperture in which galaxies should be counted that minimizes the scatter between mass and richness. In this work, using a sample of 39 clusters with accurate caustic masses at 0.1<z<0.22, we first show that the scatter between mass and richness derived from survey data is negligibly small, as small as best mass proxies. The scatter turns out to be smaller than in some previous works and has a 90% upper limit of 0.05 dex in mass. The current sample, adjoining 76 additional clusters analyzed in previous works, establishes an almost scatterless, minimally evolving (if at all), mass-richness scaling in the redshift range 0.03<z<0.55. We then exploit this negligible scatter to derive the reference aperture to be used to compute richness and to predict the mass of cluster samples. These predicted masses have a total 0.16 dex scatter with caustic mass, about half of which is not intrinsic to the proxy, but related to the noisiness of the caustic masses used for test proxy performances. These results make richness-based masses of best quality and available for large samples at a low observational cost.

One of the key ingredients of the Unified Model of Active Galactic Nuclei (AGN) is the presence of a torus-like optically thick medium composed by dust and gas around the putative supermassive black hole. However, the structure, size and composition of this circumnuclear medium are still matter of debate. To this end, the search for column density variations through X-ray monitoring on different timescales (months, weeks and few days) is fundamental to constrain size, kinematics and location of the X-ray absorber(s). Here we describe our project of mining the Swift-XRT archive to assemble a sample of AGN with extreme column density variability and determining the physical properties of the X-ray absorber(s). We also present the results obtained from a daily-weekly Swift-XRT follow-up monitoring recently performed on one of the most interesting new candidates for variability discovered so far, Mrk 915.

This work observationally addresses the relative distribution of total and optically luminous matter in galaxy clusters by computing the radial profile of the stellar-to-total mass ratio. We adopt state-of-the-art accurate lensing masses free from assumptions about the mass radial profile and we use extremely deep multicolor wide--field optical images to distinguish star formation from stellar mass, to properly calculate the mass in galaxies of low mass, those outside the red sequence, and to allow a contribution from galaxies of low mass that is clustercentric dependent. We pay special attention to issues and contributions that are usually underrated, yet are major sources of uncertainty, and we present an approach that allows us to account for all of them. Here we present the results for three very massive clusters at $z\sim0.45$, MACSJ1206.2-0847, MACSJ0329.6-0211, and RXJ1347.5-1145. We find that stellar mass and total matter are closely distributed on scales from about 150 kpc to 2.5 Mpc: the stellar-to-total mass ratio is radially constant. We find that the characteristic mass stays constant across clustercentric radii and clusters, but that the less-massive end of the galaxy mass function is dependent on the environment.

In light of the tension between cosmological parameters from Planck cosmic microwave background and galaxy clusters, we revised the Planck analysis of the $Y_{SZ}$-mass calibration to allow evolution to be determined by the data instead of being imposed as an external constraint. Our analysis uses the very same data and Malmquist bias corrections as used by the Planck team in order to emphasize that differences in the results come from differences in the assumptions. The evolution derived from 71 calibrating clusters, with $0.05<z<0.45$, is proportional to $E^{2.5\pm0.4}(z)$, so inconsistent with the self-similar evolution ($E^{2/3}$) assumed by previous analyses. When allowing for evolution, the slope of $Y_{SZ}$-mass relation turns out to be $1.51\pm0.07$, which is shallower by $4.8\sigma$ than the value derived when assuming self-similar evolution, introducing a mass-dependent bias. The non-self-similar evolution of $Y_{SZ}$ has to be accounted for in analyses aimed to establish the biases of Planck masses.

The quenching rate is known to depend on galaxy stellar mass and environment, however, possible dependences on the hosting halo properties, such as mass, richness, and dynamical status, are still debated. The determination of these dependences is hampered by systematics, induced by noisy estimates of cluster mass or by the lack of control on galaxy stellar mass, which may mask existing trends or introduce fake trends. We studied a sample of local clusters (20 with 0.02<z<0.1 and log(M200/Msun)>14), selected independent of the galaxy properties under study, having homogeneous optical photometry and X-ray estimated properties. Using those top quality measurements of cluster mass, hence of cluster scale, richness, iron abundance, and cooling time/presence of a cool-core, we study the simultaneous dependence of quenching on these cluster properties on galaxy stellar mass M and normalised cluster-centric distance r/r200. We found that the quenching rate can be completely described by two variables only, galaxy stellar mass and normalised cluster-centric distance, and is independent of halo properties (mass, richness, iron abundance, presence of a cool-core, and central cooling time). These halo properties change, in most cases, by less than 3% the probability that a galaxy is quenched, once the mass-size (M200-r200) scaling relation is accounted for through cluster-centric distance normalisation.

The metal content of clusters of galaxies and its relation to their stellar content is revisited making use of a cluster sample for which all four basic parameters are homogeneously measured within consistent radii, namely core-excised mass-weighted metallicity plus total, stellar and ICM masses. For clusters of total mass $M_{500} >$ $\sim 10^{14}$ $M_{\odot}$ nice agreement is found between their iron content and what expected from empirical supernova yields. For the same clusters, there also appears to be at least as much iron in the intracluster medium (ICM) as there is still locked into stars (i.e., the ICM/stars metal share is about unity). However, for more massive clusters the stellar mass fraction appears to drop substantially without being accompanied by a drop in the ICM metallicity, thus generating a major tension with the nucleosynthesis expectation and inflating the metal share to extremely high values (up to $\sim 6$). Various possible solutions of this conundrum are discussed, but are all considered either astrophysically implausible, or lacking an independent observational support. For this reason we still entertain the possibility that even some of the best cluster data may be faulty, though we are not able to identify any obvious bias. Finally, based on the stellar mass-metallicity relation for local galaxies we estimate the contribution of galaxies to the ICM enrichment as a function of their mass, concluding that even the most massive galaxies must have lost a major fraction of the metals they have produced.

We analysed the richness--mass scaling of 23 very massive clusters at $0.15<z<0.55$ with homogenously measured weak-lensing masses and richnesses within a fixed aperture of $0.5$ Mpc radius. We found that the richness--mass scaling is very tight (the scatter is $<0.09$ dex with 90 \% probability) and independent of cluster evolutionary status and morphology. This implies a close association between infall and evolution of dark matter and galaxies in the central region of clusters. We also found that the evolution of the richness-mass intercept is minor at most, and, given the minor mass evolution across the studied redshift range, the richness evolution of individual massive clusters also turns out to be very small. Finally, it was paramount to account for the cluster mass function and the selection function. Ignoring them would led to biases larger than the (otherwise quoted) errors. Our study benefits from: a) weak-lensing masses instead of proxy-based masses thereby removing the ambiguity between a real trend and one induced by an accounted evolution of the used mass proxy; b) the use of projected masses that simplify the statistical analysis thereby not requiring consideration of the unknown covariance induced by the cluster orientation/triaxiality; c) the use of aperture masses as they are free of the pseudo-evolution of mass definitions anchored to the evolving density of the Universe; d) a proper accounting of the sample selection function and of the Malmquist-like effect induced by the cluster mass function; e) cosmological simulations for the computation of the cluster mass function, its evolution, and the mass growth of each individual cluster.

We present Hubble Space Telescope imaging and grism spectroscopy in the field of the distant galaxy cluster JKCS 041 using the Wide Field Camera 3. We confirm that JKCS 041 is a rich cluster and derive a redshift z=1.80 via the spectroscopic identification of 19 member galaxies, of which 15 are quiescent. These are centered upon diffuse X-ray emission seen by the Chandra observatory. As JKCS 041 is the most distant known cluster with such a large, spectroscopically confirmed quiescent population, it provides a unique opportunity to study the effect of the environment on galaxy properties at early epochs. We construct high-quality composite spectra of the quiescent cluster members that reveal prominent Balmer and metallic absorption lines. Using these, we measure the mean stellar ages in two bins of stellar mass. The quiescent cluster members' ages agree remarkably closely with that inferred by Whitaker et al. for similarly selected samples in the field, supporting the idea that the cluster environment is more efficient at truncating star formation while not having a strong effect on the mean epoch of quenching. We find some evidence (90% confidence) for a lower fraction of disk-like quiescent systems in JKCS 041 compared to a sample of coeval field galaxies drawn from the CANDELS survey. Taking this into account, we do not detect a significant difference between the mass-radius relations of the quiescent JKCS 041 members and our z~1.8 field sample. Finally, we demonstrate how differences in the morphological mixture of quenched systems can complicate measures of the environmental dependence of size growth.

Using deep two-color near-infrared HST imaging and unbiased grism spectroscopy we present a detailed study of the z=1.803 JKCS041 cluster. Uniquely, for a high redshift cluster, we confirm a mass of $\log M=14.2$ in solar units using four different techniques based on the X-ray temperature, the X-ray luminosity and the cluster richness. JKCS041 is thus a progenitor of a local system like the Coma cluster. Our rich dataset and the abundant population of 14 spectroscopically-confirmed red sequence galaxies allows us to explore the past star formation history of this system in unprecedented detail. Remarkably, we find a prominent red sequence down to stellar masses as low as $\log M=9.8$, corresponding to a mass range of 2 dex. These quiescent galaxies are concentrated around the cluster center with a core radius of 330 kpc. Blue members are few and avoid the cluster center. In JKCS041 quenching was therefore largely completed by a look-back time of 10 Gyr and we can constrain the epoch at which this occurred via spectroscopic age-dating of the individual galaxies. Most galaxies were quenched about 1.1 Gyr prior to the epoch of observation. The less massive quiescent galaxies are somewhat younger, corresponding to a decrease in age of 650 Myr per mass dex, but the scatter in age at fixed mass is only 380 Myr (at $\log M=11$). The size-mass relation of quiescent galaxies in JKCS041 is consistent with that observed for local clusters within our uncertainties. Comparing our data on JKCS041 with 41 clusters at lower redshift, we find that the form of the mass function of red sequence galaxies has hardly evolved in the past 10 Gyr, both in terms of its faint end slope and characteristic mass. Despite observing JKCS041 soon after its quenching and the three fold expected increase in mass in the next 10 Gyr, it is already remarkably similar to present-day clusters.

Using Spitzer 3.6 micron data we derived the luminosity function and the mass function of galaxies in five z>1.4 clusters selected to have a firm intracluster medium detection. The five clusters differ in richness (ISCSJ1438.1+3414 and XMMXCSJ2215.9-1738 are twice as rich as ISCSJ1432.4+3250, IDCSJ1426.5+3508, and JKCS041) and morphological appareance. At the median redshift z=1.5, from the 150 member galaxies of the five clusters, we derived a characteristic magnitude of 16.92+/-0.13 in the [3.6] band and a characteristic mass of lgM*=11.30+/-0.05 Msol. We find that the characteristic luminosity and mass does not evolve between z=1 and 1.4<z<1.8, directly ruling out ongoing mass assembly between these epochs because massive galaxies are already present up to z=1.8. Lower-redshift build-up epochs have already been ruled out by previous works, leaving only z>1.8 as a possible epoch for the mass build up. However, the observed values of m^* at very high redshift are too bright for galaxies without any star formation immediately preceding the observed redshift and therefore imply a star formation episode not earlier than z_f=2.5. For the first time, mass/luminosity functions are able to robustly distinguish tiny differences between formation redshifts and to set {\it upper} limits to the epoch of the last star-formation episode.

Estimates of cosmological parameters using galaxy clusters have the scatter in the observable at a given mass as a fundamental parameter. This work computes the amplitude of the scatter for a newly introduced mass proxy, the product of the cluster total luminosity times the mass-to-light ratio, usually referred as stellar mass. The analysis of 12 galaxy clusters with excellent total masses shows a tight correlation between the stellar mass, or stellar fraction, and total mass within r500 with negligible intrinsic scatter: the 90% upper limit is 0.06 dex, the posterior mean is 0.027 dex. This scatter is similar to the one of best-determined mass proxies, such as Yx, i.e. the product of X-ray temperature and gas mass. The size of the cluster sample used to determine the intrinsic scatter is small, as in previous works proposing low-scatter proxies because very accurate masses are needed to infer very small values of intrinsic scatter. Three-quarters of the studied clusters have lgM <~14 Msol, which is advantageous from a cosmological perspective because these clusters are far more abundant than more massive clusters. At the difference of other mass proxies such as Yx, stellar mass can be determined with survey data up to at least z=0.9 using upcoming optical near-infrared surveys, such as DES and Euclid, or even with currently available surveys, covering however smaller solid angles. On the other end, the uncertainty about the predicted mass of a single cluster is large, 0.21 to 0.32 dex, depending on cluster richness. This is largely because the proxy itself has ~0.10 dex errors for clusters of lgM<~ 14 Msol mass.

This review article considers some of the most common methods used in astronomy for regressing one quantity against another in order to estimate the model parameters or to predict an observationally expensive quantity using trends between object values. These methods have to tackle some of the awkward features prevalent in astronomical data, namely heteroscedastic (point-dependent) errors, intrinsic scatter, non-ignorable data collection and selection effects, data structure and non-uniform population (often called Malmquist bias), non-Gaussian data, outliers and mixtures of regressions. We outline how least square fits, weighted least squares methods, Maximum Likelihood, survival analysis, and Bayesian methods have been applied in the astrophysics literature when one or more of these features is present. In particular we concentrate on errors-in-variables regression and we advocate Bayesian techniques.

In this paper we describe the first data release of the the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is a ~12degree^2 survey in the near-infrared Z,Y,J,H and K_s bands, specifically designed to enable the evolution of galaxies and large structures to be traced as a function of both epoch and environment from the present day out to z=4, and active galactic nuclei (AGN) and the most massive galaxies up to and into the epoch of reionization. With its depth and area, VIDEO will be able to fully explore the period in the Universe where AGN and starburst activity were at their peak and the first galaxy clusters were beginning to virialize. VIDEO therefore offers a unique data set with which to investigate the interplay between AGN, starbursts and environment, and the role of feedback at a time when it was potentially most crucial. We provide data over the VIDEO-XMM3 tile, which also covers the Canada-France-Hawaii-Telescope Legacy Survey Deep-1 field (CFHTLS-D1). The released VIDEO data reach a 5-sigma AB-magnitude depth of Z=25.7, Y=24.5, J=24.4, H=24.1 and K_s=23.8 in 2 arcsec diameter apertures (the full depth of Y=24.6 will be reached within the full integration time in future releases). The data are compared to previous surveys over this field and we find good astrometric agreement with the Two-Micron All Sky Survey, and source counts in agreement with the recently released UltraVISTA survey data. The addition of the VIDEO data to the CFHTLS-D1 optical data increases the accuracy of photometric redshifts and significantly reduces the fraction of catastrophic outliers over the redshift range 0<z<1 from 5.8 to 3.1 per cent in the absence of an i-band luminosity prior. (Truncated Abstract)

This work attains a threefold objective: first, we derived the richness-mass scaling in the local Universe from data of 53 clusters with individual measurements of mass. We found a 0.46+-0.12 slope and a 0.25+-0.03 dex scatter measuring richness with a previously developed method. Second, we showed on a real sample of 250 0.06<z<0.9 clusters, most of which are at z<0.3, with spectroscopic redshift that the colour of the red sequence allows us to measure the clusters' redshift to better than Delta z=0.02. Third, we computed the predicted prior of the richness-mass scaling to forecast the capabilities of future wide-field-area surveys of galaxy clusters to constrain cosmological parameters. We computed the uncertainty and the covariance matrix of the (evolving) richness-mass scaling of a PanStarrs 1+Euclid-like survey accounting for a large suite of sources of errors. We find that the richness-mass scaling parameters, which are the input ingredients of cosmological forecasts using cluster counts, can be determined 10^5 times better than estimated in previous works that did not use weak-lensing mass estimates. The better knowledge of the scaling parameters likely has a strong impact on the relative importance of the different probes used to constrain cosmological parameters. Richness-mass scaling parameters were recovered, but only if the cluster mass function and the weak-lensing redshift-dependent selection function were accounted for in the fitting of the mass-richness scaling. This emphasizes the limitations of often adopted simplifying assumptions, such as having a mass-complete redshift-independent sample. The fitting code used for computing the predicted prior, including the treatment of the mass function and of the weak-lensing selection function, is provided in the appendix. [Abridged]

This work measures the evolution of the iron content in galaxy clusters by a rigorous analysis of the data of 130 clusters at 0.1<z<1.3. This task is made difficult by a) the low signal-to-noise ratio of abundance measurements and the upper limits, b) possible selection effects, c) boundaries in the parameter space, d) non-Gaussian errors, e) the intrinsic variety of the objects studied, and f) abundance systematics. We introduce a Bayesian model to address all these issues at the same time, thus allowing cross-talk (covariance). On simulated data, the Bayesian fit recovers the input enrichment history, unlike in standard analysis. After accounting for a possible dependence on X-ray temperature, for metal abundance systematics, and for the intrinsic variety of studied objects, we found that the present-day metal content is not reached either at high or at low redshifts, but gradually over time: iron abundance increases by a factor 1.5 in the 7 Gyr sampled by the data. Therefore, feedback in metal abundance does not end at high redshift. Evolution is established with a moderate amount of evidence, 19 to 1 odds against faster or slower metal enrichment histories. We quantify, for the first time, the intrinsic spread in metal abundance, 18+/-3 %, after correcting for the effect of evolution, X-ray temperature, and metal abundance systematics. Finally, we also present an analytic approximation of the X-ray temperature and metal abundance likelihood functions, which are useful for other regression fitting involving these parameters. The data for the 130 clusters and code used for the stochastic computation are provided with the paper.

This short note points out two of the incongruences that I find in the Loredo (2012) comments on Andreon (2012), i.e. on my chapter written for the book "Astrostatistical Challenges for the New Astronomy". First, I find illogic the Loredo decision of putting my chapter among those presenting simple models, because one of the models illustrated in my chapter is qualified by him as "impressing for his complexity". Second, Loredo criticizes my chapter at one location confusing it with another paper by another author, because my chapter do not touch the subject mentioned by Loredo (2012) critics, the comparison between Bayesian and frequentist fitting models.

Galaxy mass and environment are known to play a key role in galaxy evolution: looking at galaxy colors at different redshifts, fixed galaxy mass and environment, offers a powerful diagnosis to disentangle the role of each. In this work, we study the simulateneous dependence of the fraction of blue galaxies fblue on secular evolution, environment and galaxy mass with a well-controlled cluster sample. We are thus able to study the evolution and respective role of the cessation of star formation history (SFH) in clusters due to galaxy mass ("mass quenching") or to environment ("environmental quenching"). We define an homogenous X-ray selected cluster sample (25 clusters with 0 < z < 1 and one cluster at z \sim 2.2), having similar masses and well-defined sizes. Using multicolor photometry and a large spectroscopic sample to calibrate photometric redshifts, we carefully estimate fblue for each cluster at different galaxy mass and cluster-centric distance bins. We then fit with a simple model the dependence of fblue on redshift (z), environment (r/r200) and galaxy mass (M). fblue increases with cluster-centric distance with a slope $1.2^{+0.4}_{-0.3}$, decreases with galaxy mass with a slope $-3.8^{+0.6}_{-0.5}$, and increases with redshift with a slope $3.2^{+0.7}_{-0.5}$. The data also require for the first time a differential evolution with galaxy mass of fblue with redshift, with lower mass galaxies evolving slower by a factor $-4.1^{+1.1}_{-0.9}$. Our study shows that the processes responsible for the cessation of star formation in clusters are effective at all epochs (z<2.2), and more effective in denser environments and for more massive galaxies. We found that the mass and environmental quenchings are separable, that environmental quenching does not change with epoch, and that mass quenching is a dynamical process, i.e. its evolutionary rate is mass-dependent. [Abridged]

Current analysis of astronomical data are confronted with the daunting task of modeling the awkward features of astronomical data, among which heteroscedastic (point-dependent) errors, intrinsic scatter, non-ignorable data collection (selection effects), data structure, non-uniform populations (often called Malmquist bias), non-Gaussian data, and upper/lower limits. This chapter shows, by examples, how modeling all these features using Bayesian methods. In short, one just need to formalize, using maths, the logical link between the involved quantities, how the data arise and what we already known on the quantities we want to study. The posterior probability distribution summarizes what we known on the studied quantities after the data, and we should not be afraid about their actual numerical computation, because it is left to (special) Monte Carlo programs such as JAGS. As examples, we show how to predict the mass of a new object disposing of a calibrating sample, how to constraint cosmological parameters from supernovae data and how to check if the fitted data are in tension with the adopted fitting model. Examples are given with their coding. These examples can be easily used as template for completely different analysis, on totally unrelated astronomical objects, requiring to model the same awkward data features.

The dependence of galaxy star formation activity on environment - especially in clusters - at high redshift is still poorly understood, as illustrated by the still limited number of z>1.4 clusters on the one hand, and by the still debated star formation-density relation at high redshift on the other hand. The zphot~2.2 JKCS041 cluster allows to probe such environmental dependence of star formation activity at an unprecedented combination of redshifts and environments. Its study permits to enlarge the knowledge of high redshift clusters and to put strong leverage on observational constraints for galaxy evolution models. We analyze deep u\astg'r'i'z'JHKs images from the CFHTLS/WIRDS surveys, which cover JKCS041 cluster field. We first estimate photometric redshifts based on multi-wavelength photometry. We then lead a careful analysis to test the presence of a Butcher-Oemler effect. We work on galaxies within 2\timesr200 and with masses >1.34\times10^11 Msun, and use two comparison clusters at z=0 and z=1 of similar mass. We estimate the radial profiles of the fraction of blue galaxies, taking into account the star aging with decreasing redshift. After confirming the high redshift nature of JKCS041, we find no evidence for a Butcher-Oemler effect between z~2.2 and z~0 for galaxies more massive than 1.34\times10^11 Msun. In the cluster center, a change greater than \Deltafblue/\Deltaz=0.16 between z~0 and z~2.2 would be easily detected. We also find that JKCS041 shows a consistent and systematic increase of the fraction of star-forming galaxies with cluster-centric distance, hence with decreasing density, for both a M>1.34\times10^11 Msun selected sample and a lower mass sample. In particular, very few (less than 15%) star-forming galaxies are found within r200/2 among high mass (M>1.34\times10^11 Msun) galaxies. Our results show that the present-day star formation-density relation is already in place at z~2.2.

[Abridged] We study the X-ray properties of a color-selected sample of clusters at 0.1<z<0.3, to quantify the real aboundance of the population of X-ray dark or underluminous clusters and at the same time the spurious detection contamination level of color-selected cluster catalogs. A careful selection allowed us to have an X-ray-unbiased sample of 33 clusters to measure the Lx-richness relation. Swift 1.4 Ms X-ray observations show that at least 89% of the color-detected clusters are real objects with a potential well deep enough to heat and retain an intracluster medium. The percentage rises to 94 % when one includes the single spectroscopically confirmed color-selected cluster whose X-ray emission is not secured. Looking at our results from the opposite perspective, the percentage of X-ray dark clusters among color-selected clusters is very low: at most about 11 per cent (at 90 % confidence). Supplementing our data with those from literature, we conclude that X-ray- and color- cluster surveys sample the same population and consequently that in this regard we can safely use clusters selected with any of the two methods for cosmological purposes. This is an essential and promising piece of information for upcoming surveys in both the optical/IR (DES, EUCLID) and X-ray (eRosita). Richness correlates with X-ray luminosity with a large scatter, 0.51+/-0.08 (0.44+/-0.07) dex in lgLx at a given richness, when Lx is measured in a 500 (1070) kpc aperture. We release data and software to estimate the X-ray flux, or its upper limit, of a source with over-Poisson background fluctuations (found in this work to be ~20% on cluster angular scales) and to fit X-ray luminosity vs richness if there is an intrinsic scatter. These Bayesian applications rigorously account for boundaries (e.g. the X-ray luminosity and the richness cannot be negative).

Above redshift z ~ 1.4, we known little or nothing about the stellar ages of red galaxies in clusters, yet at these high redshifts important changes are predicted by current renditions of galaxy formation models embedded in the standard hierarchical paradigm of structure formation. Red-sequence galaxies in the cluster JKCS041 at z=2.2 show a tight distribution in colour that indicates a star formation history that is highly synchronized across galaxies. Specifically, we measure a spread in stellar age of 160+/-30 Myr, in marked disagreement with the current understanding of how massive red galaxies form in clusters, i.e. if they are produced somewhat stochastically in merging episodes that sometimes involve gas, hence star formation. The existence of a tight distribution in colour when the universe was at one quarter of its current age implies that mechanisms that have not yet been implemented in current galaxy formation scenarios long ago began to shape the star formation history of red cluster galaxies.

The evolution of the properties of the hot gas that fills the potential well of galaxy clusters is poorly known, since models are unable to give robust predictions and observations lack a sufficient redshift leverage and are affected by selection effects. Here, with just two high redshift, z approx 1.8, clusters avoiding selection biases, we obtain a significant extension of the redshift range and we begin to constrain the possible evolution of the X-ray luminosity vs temperature relation. The two clusters, JKC041 at z=2.2 and ISCSJ1438+3414 at z=1.41, are respectively the most distant cluster overall, and the second most distant that can be used for studying scaling relations. Their location in the X-ray luminosity vs temperature plane, with an X-ray luminosity 5 times lower than expected, suggests at the 95 % confidence that the evolution of the intracluster medium has not been self-similar in the last three quarters of the Universe age. Our conclusion is reinforced by data on a third, X-ray selected, high redshift cluster, too faint for its temperature when compared to a sample of similarly selected objects. Our data suggest that non-gravitational effects, such as the baryon physics, influence the evolution of galaxy cluster. Precise knowledge of evolution is central for using galaxy clusters as cosmological probes in planned X-ray surveys such as WFXT or JDEM.

New deep z'-J data readly show a narrow red sequence co-centered with, and similary concentrated to, the extended X-ray emission of the cluster of galaxies JKCS041. The JKCS041 red sequence is 0.32+/-0.06 mag redder in z'-J than the red sequence of the zspec=1.62 IRC0218A cluster, putting JKCS041 at z<<1.62. The colour difference of the two red sequences gives a red-sequence based redshift of z=2.20+/-0.11 for JKCS041, where the uncertainty accounts for uncertainties in stellar synthesis population models, in photometric calibration and in the red sequence colour of both JKCS041 and IRC0218A clusters.

This paper aims at robustly determining the redshift of the cluster of galaxies JKCS041 and at putting constraints on the formation epoch of the color-magnitude sequence in two very high redshift clusters. New deep z'-J data show a clear narrow red sequence that is co-centered with, and similarly concentrated on, the extended X-ray emission of the cluster of galaxies JKCS041. The JKCS041 red sequence is 0.32+/-0.06 mag redder in z'-J than the red sequence of the zspec=1.62 IRC0218A cluster, putting JKCS041 at z>>1.62 and ruling out z<~1.49 the latter claimed by a recent paper. The color difference of the two red sequences gives a red-sequence-based redshift of z=2.20+/-0.11 for JKCS041, where the uncertainty accounts for uncertainties in stellar synthesis population models, in photometric calibration, and in the red sequence color of both JKCS041 and IRC0218A clusters. We do not observe any sign of truncation of the red sequence for both clusters down to J=23 mag (1.0e+11 solar masses), which suggests that it is already in place in clusters rich and massive enough to heat and retain hot gas at these high redshifts.

We present gas constraints from Sunyaev-Zel'dovich (SZ) effect measurements in a sample of eleven X-ray and infrared (IR) selected galaxy clusters at z >=1, using data from the Sunyaev-Zel'dovich Array (SZA). The cylindrically integrated Compton-y parameter, Y , is calculated by fitting the data to a two-parameter gas pressure profile. Where possible, we also determine the temperature of the hot intra-cluster plasma from Chandra and XMM-Newton data, and constrain the gas mass within the same aperture (r_2500 ) as Y . The SZ effect is detected in the clusters for which the X-ray data indicate gas masses above ~ 10^13 Msun, including XMMU J2235-2557 at redshift z = 1.39, which to date is one of the most distant clusters detected using the SZ effect. None of the IR-selected targets are detected by the SZA measurements, indicating low gas masses for these objects. For these and the four other undetected clusters, we quote upper limits on Y and Mgas_SZ , with the latter derived from scaling relations calibrated with lower redshift clusters. We compare the constraints on Y and X-ray derived gas mass Mgas_X-ray to self-similar scaling relations between these observables determined from observations of lower redshift clusters, finding consistency given the measurement error.

[Abridged] The analysis of a sample of 52 clusters with precise and hypothesis-parsimonious measurements of mass shows that low mass clusters and groups are not simple scaled-down versions of their massive cousins in terms of stellar content: lighter clusters have more stars per unit cluster mass. The same analysis also shows that the stellar content of clusters and groups displays an intrinsic spread at a given cluster mass, i.e. clusters are not similar each other in the amount of stars they contain, not even at a fixed cluster mass. The stellar mass fraction depends on halo mass with (logarithmic) slope -0.55+/-0.08 and with 0.15+/-0.02 dex of intrinsic scatter at a fixed cluster mass. The intrinsic scatter at a fixed cluster mass we determine for gas mass fractions is smaller, 0.06+/-0.01 dex. The intrinsic scatter in both the stellar and gas mass fractions is a distinctive signature that the regions from which clusters and groups collected matter, a few tens of Mpc, are yet not representative, in terms of gas and baryon content, of the mean matter content of the Universe. The observed stellar mass fraction values are in marked disagreement with gasdynamics simulations with cooling and star formation of clusters and groups. We found the the baryon (gas+stellar) fraction is fairly constant for clusters and groups with 13.7<lg(mass)<15.0 solar masses and it is offset from the WMAP-derived value by about 6 sigmas. The offset could be related to the possible non universality of the baryon fraction pointed out by our measurements of the intrinsic scatter. Our analysis is the first that does not assume that clusters are identically equal at a given halo mass and it is also more accurate in many aspects. The data and code used for the stochastic computation are distributed with the paper.

We use a sample of 53 galaxy clusters at 0.03 < z < 0.1 with available masses derived from the caustic technique and with velocity dispersions computed using 208 galaxies on average per cluster, in order to investigate the scaling between richness, mass and velocity dispersion. A tight scaling between richness and mass is found, with an intrinsic scatter of only 0.19 dex in mass and with a slope one, i.e. clusters which have twice as many galaxies are twice as massive. When richness is measured without any knowledge of the cluster mass or linked parameters (such as r200), it can predict mass with an uncertainty of 0.29+/-0.01 dex. As a mass proxy, richness competes favourably with both direct measurements of mass given by the caustic method, which has typically 0.14 dex errors (vs 0.29) and X-ray luminosity, which offers a similar 0.30 dex uncertainty. The similar performances of X-ray luminosity and richness in predicting cluster masses has been confirmed using cluster masses derived from velocity dispersion fixed by numerical simulations. These results suggest that cluster masses can be reliably estimated from simple galaxy counts, at least at the redshift and masses explored in this work. This has important applications in the estimation of cosmological parameters from optical cluster surveys, because in current surveys clusters detected in the optical range outnumber, by at least one order of magnitude, those detected in X-ray. Our analysis is robust from astrophysical and statistical perspectives. The data and code used for the stochastic computation is distributed with the paper. [Abridged]

[Abridged] We report the discovery of JKCS041, a massive near-infrared selected cluster of galaxies z=1.9. The cluster was originally discovered using a modified red-sequence method and was also detected in follow-up Chandra data as extended X-ray source. Optical and near-infrared imaging data alone allow us to show that the detection of JKCS041 is as secure, even in absence of the X-ray data. We investigate the possibility that JKCS041 is not a discrete galaxy cluster at z=1.9, and find other explanations to be unlikely. The X-ray detection and statistical arguments rule out the hypothesis that JKCS041 is actually a blend of groups along the line of sight, and we find that the X-ray emitting gas is too hot and dense to be a filament projected along the line of sight. The cluster has an X-ray core radius of 36.6 arcsec (about 300 kpc), an X-ray temperature of 7.4 keV, a bolometric X-ray luminosity within R500 of 7.6 10^44 erg/s, an estimated mass of M500=2.9 10^14 Msol, the latter derived under the usual (and strong) assumptions. The cluster is composed of 16.4 galaxies within 1.5 arcmin (750 kpc) brighter than K~20.7 mag. The high redshift of JKCS041 is determined from the detection colour, from the detection of the cluster in a galaxy sample formed by z_phot>1.6 galaxies and from a photometric redshift based on 11-band spectral energy distribution fitting. By means of the latter we find the cluster redshift to be 1.84<z<2.12 at 68 % confidence. Therefore, JKCS041 is a cluster of galaxies at z=1.9 with deep potential well, making it the most distant cluster with extended X-ray emission known.

AIMS: We present a sample of candidate quasars selected using the KX-technique. The data cover 0.68 deg^2 of the X-ray Multi-Mirror (XMM) Large-Scale Structure (LSS) survey area where overlapping multi-wavelength imaging data permits an investigation of the physical nature of selected sources. METHODS: The KX method identifies quasars on the basis of their optical (R and z') to near-infrared (Ks) photometry and point-like morphology. We combine these data with optical (u*,g'r',i',z') and mid-infrared (3.6-24 micron) wavebands to reconstruct the spectral energy distributions (SEDs) of candidate quasars. RESULTS: Of 93 sources selected as candidate quasars by the KX method, 25 are classified as quasars by the subsequent SED analysis. Spectroscopic observations are available for 12/25 of these sources and confirm the quasar hypothesis in each case. Even more, 90% of the SED-classified quasars show X-ray emission, a property not shared by any of the false candidates in the KX-selected sample. Applying a photometric redshift analysis to the sources without spectroscopy indicates that the 25 sources classified as quasars occupy the interval 0.7 < z < 2.5. The remaining 68/93 sources are classified as stars and unresolved galaxies.

[Abridged]We present the morphological analysis based on HST-NICMOS observations in the F160W filter of a sample of 32 early-type galaxies (ETGs) at 1<z<2 with spectroscopic confirmation of their redshift and spectral type. We find that ETGs at are composed of two distinct populations, an older population (oETGs) and a younger population (yETGs) whose mean ages differ by about 1.5-2 Gyr. Young ETGs are not denser than local ones since they follow the size-mass relation of local ETGs and luminosity evolution brings them onto the local Kormendy and size-luminosity relations. On the constrary, old ETGs do not follow the size-mass relation of local ETGs and luminosity evolution does not account for the discrepancy they show with respect to the local size-luminosity and Kormendy relations. An increase of their effective radius R_e by a factor 2.5-3 (a density decrease by a factor 15-30) from z~1.5-2 to z~0 is required to bring these galaxies onto the local scaling relations. The different behaviour shown by the two populations implies different formation and evolution scenarios. The older population must have formed at higher-z in a sort of dissipative gas-rich collapse ables to produce remnants which at z~2 are old and compact. The size evolution they must experience from z~2 to z~0 has to leave unchanged their mass to not exceed the local number of high-mass ETGs. Thus, major merging cannot fit this requirement. The younger population of ETGs can be formed later through subsequent episodes of merging which increased progressively their size and assembled their mass down to z~2. At z<2 they evolve purely in luminosity since episodes of major merging would bring them far from the local scaling relations.

This paper presents morphological type, membership, and U-V color for a sample of galaxies in the Coma cluster direction, complete down to M_B=-15.00 mag and extending down to M_B=-14.25 mag. We have examined 1155 objects from the GMP 1983 catalog on B and V images of the CFH12K camera, and obtained the Hubble type in most cases. Coma cluster membership for 473 galaxies was derived using morphology, apparent size, and surface brightness, and, afterward, redshift. The comparison among morphology- and redshift- memberships and among luminosity functions derived from this morphologically-selected sample, or by using statistical members or spectroscopic members, all show that the morphological membership provided here can be trusted. For the first time, the morphological classification of Coma galaxies reaches magnitudes that are faint enough to observe the whole magnitude range of the giant types, E, S0, and spiral stages. The data presented in this paper makes our sample the richest environment where membership and morphology for complete samples down to faint magnitudes M_B~-15 mag are available, thereby enlarging the baseline of environmental studies.

We measure the relative evolution of the number of bright and faint (as faint as 0.05 L*) red galaxies in a sample of 28 clusters, of which 16 are at 0.50<= z<=1.27, all observed through a pair of filters bracketing the 4000 Angstrom break rest-frame. The abundance of red galaxies, relative to bright ones, is constant over all the studied redshift range, 0<z<1.3, and rules out a differential evolution between bright and faint red galaxies as large as claimed in some past works. Faint red galaxies are largely assembled and in place at z=1.3 and their deficit does not depend on cluster mass, parametrized by velocity dispersion or X-ray luminosity. Our analysis, with respect to previous one, samples a wider redshift range, minimizes systematics and put a more attention to statistical issues, keeping at the same time a large number of clusters.

We present an analysis of deep XMM-Newton and Chandra observations of the z=1.05 galaxy cluster XLSSJ022403.9-041328 (hereafter XLSSC 029), detected in the XMM-Newton large scale structure survey. Density and temperature profiles of the X-ray emitting gas were used to perform a hydrostatic mass analysis of the system. This allowed us to measure the total mass and gas fraction in the cluster and define overdensity radii R500 and R2500. The global properties of XLSSC 029 were measured within these radii and compared with those of the local population. The gas mass fraction was found to be consistent with local clusters. The mean metal abundance was 0.18 +0.17 -0.15 Zsol, with the cluster core regions excluded, consistent with the predicted and observed evolution. The properties of XLSSC 029 were then used to investigate the position of the cluster on the M-kT, YX-M, and LX-M scaling relations. In all cases the observed properties of XLSSC 029 agreed well with the simple self-similar evolution of the scaling relations. This is the first test of the evolution of these relations at z > 1 and supports the use of the scaling relations in cosmological studies with distant galaxy clusters.

We present the morphological analysis based on HST-NICMOS observations in the F160W filter of a sample of 30 early-type galaxies spectroscopically confirmed at 1.2<z<2. We derive the effective radius R_e and the mean surface brightness mu_e of galaxies in the rest-frame R-band. We find that early-types at z~1.5 are characterized by a surface brightness (SB) much higher then their local counterparts with comparable R_e. In particular, we find that the mean SB of these early-types should get fainter by ~2.5 mag from z~1.5 to z=0 to match the SB of the local early-types with comparable R_e. This evolution exceeds by a factor two the luminosity evolution expected for early-types in this redshift range and more than a factor three the one derived from the observed luminosity function of galaxies. Consequently, an evolution of the effective radius R_e from the epoch of their formation towards z=0 has to be invoked and the hypothesis of fixed size rejected.

We present deep I and z' imaging of the colour-selected cluster RzCS 052 and study the color-magnitude relation of this cluster, its scatter, the morphological distribution on the red sequence, the luminosity and stellar mass functions of red galaxies and the cluster blue fraction. We find that the stellar populations of early type galaxies in this cluster are uniformly old and that their luminosity function does not show any sign of evolution other than the passive evolution of their stellar populations. We rule out a significant contribution from mergers in the buildup of the red sequence of RzCS 052. The cluster has a large (~30%) blue fraction and and we infer that the evolution of the blue galaxies is faster than an exponentially declining star formation model and that these objects have probably experienced starburst episodes. Mergers are unlikely to be the driver of the observed colour evolution, because of the measured constancy of the mass function, as derived from near-infrared photometry of 32 clusters, including RzCS 052, presented in a related paper. Mechanisms with clustercentric radial dependent efficiencies are disfavored as well, because of the observed constant blue fraction with clustercentric distance.

We report on the discovery of the z=1.016 cluster RzCS 052 using a modified red sequence method, followup spectroscopy and X-ray imaging. This cluster has a velocity dispersion of 710+-150 km/s, a virial mass of 4.0e14 Msol (based on 21 spectroscopically confirmed members) and an X-ray luminosity of (0.68+- 0.47)e44 ergs/s in the [1-4] keV band. This optically selected cluster appears to be of richness class 3 and to follow the known L_X-sigma_v relation for high redshift X-ray selected clusters. Using these data, we find that the halo occupation number for this cluster is only marginally consistent with what expected assuming a self-similar evolution of cluster scaling relations, suggesting perhaps a break of them at z~1. We also rule out a strong galaxy merging activity between z=1 and today. Finally, we present a Bayesian approach to measuring cluster velocity dispersions and X-ray luminosities in the presence of a background: we critically reanalyze recent claims for X-ray underluminous clusters using these techniques and find that the clusters can be accommodated within the existing L_X -sigma_v relation.

We present a sample of 29 galaxy clusters from the XMM-LSS survey over an area of some 5deg2 out to a redshift of z=1.05. The sample clusters, which represent about half of the X-ray clusters identified in the region, follow well defined X-ray selection criteria and are all spectroscopically confirmed. For all clusters, we provide X-ray luminosities and temperatures as well as masses. The cluster distribution peaks around z=0.3 and T =1.5 keV, half of the objects being groups with a temperature below 2 keV. Our L-T(z) relation points toward self-similar evolution, but does not exclude other physically plausible models. Assuming that cluster scaling laws follow self-similar evolution, our number density estimates up to z=1 are compatible with the predictions of the concordance cosmology and with the findings of previous ROSAT surveys. Our well monitored selection function allowed us to demonstrate that the inclusion of selection effects is essential for the correct determination of the evolution of the L-T relation, which may explain the contradictory results from previous studies. Extensive simulations show that extending the survey area to 10deg2 has the potential to exclude the non-evolution hypothesis, but that constraints on more refined ICM models will probably be limited by the large intrinsic dispersion of the L-T relation. We further demonstrate that increasing the dispersion in the scaling laws increases the number of detectable clusters, hence generating further degeneracy [in addition to sigma8, Omega_m, L(M,z) and T(M,z)] in the cosmological interpretation of the cluster number counts. We provide useful empirical formulae for the cluster mass-flux and mass-count-rate relations as well as a comparison between the XMM-LSS mass sensitivity and that of forthcoming SZ surveys.

Following the presentation of the XMM-LSS X-ray source detection package by Pacaud et al., we provide the source lists for the first 5.5 surveyed square degrees. The catalogues pertain to the [0.5-2] and [2-10] keV bands and contain in total 3385 point-like or extended sources above a detection likelihood of 15 in either band. The agreement with deep logN-logS is excellent. The main parameters considered are position, countrate, source extent with associated likelihood values. A set of additional quantities such as astrometric corrections and fluxes are further calculated while errors on the position and countrate are deduced from simulations. We describe the construction of the band-merged catalogue allowing rapid sub-sample selection and easy cross-correlation with external multi-wavelength catalogues. A small optical CFHTLS multi-band subset of objects is associated wich each source along with an X-ray/optical overlay. We make the full X-ray images available in FITS format. The data are available at CDS and, in a more extended form, at the Milan XMM-LSS database.

Our goal is to probe the populations of obscured and unobscured AGN investigating their optical-IR and X-ray properties as a function of X-ray flux, luminosity and redshift within a hard X-ray selected sample of 136 X-ray sources in the XMM Medium Deep Survey (XMDS) with wide multiwavelength coverage. The XMDS area is covered with optical photometry from the VVDS and CFHTLS surveys and infrared Spitzer data. Based on the X-ray luminosity and X-ray to optical ratio, 132 sources are likely AGN, of which 122 have unambiguous optical - IR identification. The observed optical and IR spectral energy distributions of sources are fitted with AGN/galaxy templates in order to classify them and compute photometric redshifts. 70% of the AGN are fitted by a type 2 AGN or a star forming galaxy template and are grouped together in a single class of ``optically obscured'' AGN. They have ``red'' optical colors and generally show significant X-ray absorption from X-ray spectra or hardness ratios (N$_H > 10^{22}$ cm$^{-2}$). Sources with SEDs typical of type 1 AGN have ``blue'' optical colors and exhibit X-ray absorption in about 30% of cases. We performed a stacking analysis for obscured and type 1 AGN. The stacked X-ray spectrum of obscured AGN is flatter than that of type 1 AGN and has an average spectral slope of Gamma = 1.6. The subsample of objects fitted by a galaxy template has an even harder stacked spectrum, with Gamma = 1.2 - 1.3. The obscured fraction is larger at lower fluxes, lower redshifts and lower luminosities. X-ray absorption is less common than ``optical'' obscuration and its incidence is nearly constant with redshift and luminosity. This implies that X-ray absorption is not necessarily related to optical obscuration.

We present the spectral energy distributions (SEDs) of a hard X-ray selected sample. The sample contains 136 sources with F(2-10 keV)>10^-14 erg/cm^2/s and 132 are AGNs. The sources are detected in a 1 square degree area of the XMM-Newton-Medium Deep Survey where optical data from the VVDS, CFHTLS surveys, and infrared data from the SWIRE survey are available. Based on a SED fitting technique we derive photometric redshifts with sigma(1+z)=0.11 and 6% of outliers and identify AGN signatures in 83% of the objects. This fraction is higher than derived when a spectroscopic classification is available. The remaining 17+9-6% of AGNs shows star-forming galaxy SEDs (SF class). The sources with AGN signatures are divided in two classes, AGN1 (33+6-1%) and AGN2 (50+6-11). The AGN1 and AGN2 classes include sources whose SEDs are fitted by type 1 and type 2 AGN templates, respectively. On average, AGN1s show soft X-ray spectra, consistent with being unabsorbed, while AGN2s and SFs show hard X-ray spectra, consistent with being absorbed. The analysis of the average SEDs as a function of X-ray luminosity shows a reddening of the IR SEDs, consistent with a decreasing contribution from the host galaxy at higher luminosities. The AGNs in the SF classes are likely obscured in the mid-infrared, as suggested by their low L(3-20micron)/Lcorr(0.5-10 keV) ratios. We confirm the previously found correlation for AGNs between the radio luminosity and the X-ray and the mid-infrared luminosities. The X-ray-radio correlation can be used to identify heavily absorbed AGNs. However, the estimated radio fluxes for the missing AGN population responsible for the bulk of the background at E>10 keV are too faint to be detected even in the deepest current radio surveys.

We have measured the 3.6 micron luminosity and mass evolution of about 1000 galaxies in 32 clusters at 0.2<z<1.25 with a special attention to methodological issues, as emphasized in this proceeding contribution. We find that the luminosity of our galaxies evolves as an old and passively evolving population formed at high redshift without any need for additional redshift-dependent evolution. Models with a prolonged stellar mass growth are rejected by the data with high confidence. The data also reject models in which the age of the stars is the same at all redshifts. Similarly, the characteristic stellar mass evolves, in the last two thirds of the universe age, as expected for a stellar population formed at high redshift. Together with the old age of stellar populations derived from fundamental plane studies, our data seems to suggest that massive early-type cluster galaxies have been completely assembled at high redshift, and not only that their stars are old. The quality of the data allows us to derive the LF and mass evolution homogeneously over the whole redshift range, using a single estimator. The Schechter function describes the galaxy luminosity function well. The characteristic luminosity at z=0.5 is found to be 16.30 mag, with an uncertainty of 10 per cent.

We present details of the discovery of XLSSJ022303.0-043622, a z=1.2 cluster of galaxies. This cluster was identified from its X-ray properties and selected as a z>1 candidate from its optical/near-IR characteristics in the XMM Large-Scale Structure Survey (XMM-LSS). It is the most distant system discovered in the survey to date. We present ground-based optical and near IR observations of the system carried out as part of the XMM-LSS survey. The cluster has a bolometric X-ray luminosity of 1.1 +/- 0.7 x 10^44 erg/s, fainter than most other known z>1 X-ray selected clusters. In the optical it has a remarkably compact core, with at least a dozen galaxies inside a 125 kpc radius circle centred on the X-ray position. Most of the galaxies within the core, and those spectroscopically confirmed to be cluster members, have stellar masses similar to those of massive cluster galaxies at low redshift. They have colours comparable to those of galaxies in other z>1 clusters, consistent with showing little sign of strong ongoing star formation. The bulk of the star formation within the galaxies appears to have ceased at least 1.5 Gyr before the observed epoch. Our results are consistent with massive cluster galaxies forming at z>1 and passively evolving thereafter. We also show that the system is straightforwardly identified in Spitzer/IRAC 3.6 and 4.5 micron data obtained by the SWIRE survey emphasising the power and utility of joint XMM and Spitzer searches for the most distant clusters.

We present the XMM-LSS cluster catalogue corresponding to the CFHTLS D1 area. The list contains 13 spectroscopically confirmed, X-ray selected galaxy clusters over 0.8 deg2 to a redshift of unity and so constitutes the highest density sample of clusters to date. Cluster X-ray bolometric luminosities range from 0.03 to 5x10^{44} erg/s. In this study, we describe our catalogue construction procedure: from the detection of X-ray cluster candidates to the compilation of a spectroscopically confirmed cluster sample with an explicit selection function. The procedure further provides basic X-ray products such as cluster temperature, flux and luminosity. We detected slightly more clusters with a (0.5-2.0 keV) X-ray fluxes of >2x10^{-14} erg/s/cm^{-2} than we expected based on expectations from deep ROSAT surveys. We also present the Luminosity-Temperature relation for our 9 brightest objects possessing a reliable temperature determination. The slope is in good agreement with the local relation, yet compatible with a luminosity enhancement for the 0.15 < z< 0.35 objects having 1 < T < 2 keV, a population that the XMM-LSS is identifying systematically for the first time. The present study permits the compilation of cluster samples from XMM images whose selection biases are understood. This allows, in addition to studies of large-scale structure, the systematic investigation of cluster scaling law evolution, especially for low mass X-ray groups which constitute the bulk of our observed cluster population. All cluster ancillary data (images, profiles, spectra) are made available in electronic form via the XMM-LSS cluster database.

The colour and luminosity distributions of red galaxies in the cluster Abell 1185 (z=0.0325) were studied down to M*+8 in the B, V and R bands. The colour-magnitude (hereafter CM) relation is linear without evidence for a significant bending down to absolute magnitudes which are seldom probed in literature (M_R=-12.5 mag). The CM relation is thin (+/-0.04 mag) and its thickness is quite independent from the magnitude. The luminosity function of red galaxies in Abell 1185 is adequately described by a Schechter function, with a characteristic magnitude and a faint end slope that also well describe the LF of red galaxies in other clusters. There is no passband dependency of the LF shape other than an obvious M* shift due to the colour of the considered population. Finally, we conclude that, based on colours and luminosity, red galaxies form an homogeneous population over four decades in stellar mass, providing a second evidence against faint red galaxies being a recent cluster population.

Using one of the deepest datasets available, we determine that the red sequence of the massive cluster MS1054-0321 at z=0.831 is well populated at all studied magnitudes, showing no deficit of faint (down to M^*+3.5) red galaxies: the faint end of the colour-magnitude relation is neither empty nor underpopulated. The effect is quantified by the computation of the luminosity function (LF) of red galaxies. We found a flat slope, showing that the abundance of red galaxies is similar at faint and at intermediate magnitudes. Comparison with present-day and z~0.4 LFs suggests that the slope of the LF is not changed, within the errors, between z=0.831 and z=0. Therefore, the analysis of the LF shows no evidence for a decreasing (with magnitude or redshift) number of faint red galaxies. The presence of faint red galaxies in high redshift clusters disfavours scenarios where the evolution of red galaxies is mass-dependent, because the mass dependency should differentially depauperate the red sequence, while the MS1054-0321 colour-magnitude relation is populated as in nearby clusters and as in z~0.4 clusters. The presence of abundant faint red galaxies in the high redshift cluster MS1054-0321 restricts the room for allocating descendants of Butcher-Oemler galaxies, because they should change the faint end slope of the LF of red galaxies, while instead the same faint end slopes are observed in MS1054-0321, at z~0 and at z~0.4. In the rich MS1054-0321 cluster, the colour-magnitude relation seems to be fully in place at z=0.831 and therefore red galaxies of all magnitudes were wholly assembled at higher redshift.

We have measured the 3.6 micron luminosity evolution of about 1000 galaxies in 32 clusters at 0.2<z<1.25, without any a priori assumption about luminosity evolution, i.e. in a logically rigorous way. We find that the luminosity of our galaxies evolves as an old and passively evolving population formed at high redshift without any need for additional redshift-dependent evolution. Models with a prolonged stellar mass growth are rejected by the data with high confidence. The data also reject models in which the age of the stars is the same at all redshifts. Similarly, the characteristic stellar mass evolves, in the last two thirds of the universe age, as expected for a stellar population formed at high redshift. Together with the old age of stellar populations derived from fundamental plane studies, our data seems to suggest that early-type cluster galaxies have been completely assembled at high redshift, and not only that their stars are old. The quality of the data allows us to derive the LF and mass evolution homogeneously over the whole redshift range, using a single estimator. The Schechter function describes the galaxy luminosity function well. The characteristic luminosity at z=0.5 is is found to be 16.30 mag, with an uncertainty of 10 per cent.

The present paper focuses on the much debated Butcher-Oemler effect: the increase with redshift of the fraction of blue galaxies in clusters. Considering a representative cluster sample made of seven group/clusters at z~0.35, we have measured the blue fraction from the cluster core to the cluster outskirts and the field mainly using wide field CTIO images. This sample represents a random selection of a volume complete x-ray selected cluster sample, selected so that there is no physical connection with the studied quantity (blue fraction), to minimize observational biases. In order to statistically assess the significance of the Butcher-Oemler effect, we introduce the tools of Bayesian inference. Furthermore, we modified the blue fraction definition in order to take into account the reduced age of the universe at higher redshifts, because we should no longer attempt to reject an unphysical universe in which the age of the Universe does depend on redshift, whereas the age of its content does not. We measured the blue fraction from the cluster center to the field and we find that the cluster affects the properties of the galaxies up to two virial radii at z~0.35. Data suggest that during the last 3 Gyrs no evolution of the blue fraction, from the cluster core to the field value, is seen beyond the one needed to account for the varying age with redshift of the Universe and of its content. The agreement of the radial profiles of the blue fraction at z=0 and z~0.35 implies that the pattern infall did not change over the last 3 Gyr, or, at least, its variation has no observational effect on the studied quantity.

These short notes collect some thoughts about the way Butcher-Oemler-like and evolutionary studies are performed. We emphasize the shortcomings of several overlooked technical ingredients in the above type of studies, with the hope of making a useful cookbook for future works. We also briefly report new results of an ongoing Butcher-Oemler study at intermediate redshift.

We present X-ray and optical spectroscopic observations of twelve galaxy groups and clusters identified within the XMM Large-Scale Structure (LSS) survey. Groups and clusters are selected as extended X-ray sources from a 3.5 deg2 XMM image mosaic above a flux limit 8e-15 ergs/s/cm2 in the [0.5-2] keV energy band. Deep BVRI images and multi-object spectroscopy confirm each source as a galaxy concentration located within the redshift interval 0.29<z<0.56. We combine line-of-sight velocity dispersions with the X-ray properties of each structure computed from a two-dimensional surface brightness model and a single temperature fit to the XMM spectral data. The resulting distribution of X-ray luminosity, temperature and velocity dispersion indicate that the XMM-LSS survey is detecting low-mass clusters and galaxy groups to redshifts z < 0.6. Confirmed systems display little or no evidence for X-ray luminosity evolution at a given X-ray temperature compared to lower redshift X-ray group and cluster samples. A more complete understanding of these trends will be possible with the compilation of a statistically complete sample of galaxy groups and clusters anticipated within the continuing XMM-LSS survey.

In this paper we present a rigorous derivation of the luminosity function (LF) in presence of a background. Our approach is free from the logical contradictions of assigning negative values to positively defined quantities and avoid the use of incorrect estimates for the 68 % confidence interval (error bar). It accounts for Poisson fluctuations ignored in previous approaches and does not requires binning of the data. The method is extensible to more complex situations, does not require the existence of an environment--independent LF, and clarifies issues common to field LF derivations. We apply the method to two clusters of galaxies at intermediate redshift (z~0.3) with among the deepest and widest K_s observations ever taken. Finally, we point out short-comings of flip--flopping magnitudes.

We present results of an initial search for clusters of galaxies at z~1 and above, using data from 2.9 square degrees of XMM-Newton images. By selecting weak potentially extended X-ray sources with faint or no identifications in deep, ground-based optical imaging, we have constructed a starting sample of 19 high redshift cluster candidates. Near-IR and R,z' imaging of these fields identified nine of them as high redshift systems. Six of these were confirmed spectroscopically, three at z~1.0 and the other three in 0.8<z<0.92 range. The remaining three systems have solid photometric evidence to be at z_phot ~ 0.8, 1.0 and 1.3. The present sample significantly increases the number of such clusters. The measured density of z>~1 clusters, after discarding "low" redshift systems at z <~ 0.92 is about 1.7 /deg2 (with 68 % confidence interval equal to [1.0,2.9]) for f_X >~2.5 10^{-15} ergs/cm2/s ([0.5-2] keV) and this is a lower limit, having screened not all potential z~1 candidate clusters. Coordinates, x-ray measures and evidence for nine x-ray selected high redshift clusters is given.

We have designed a medium deep large area X-ray survey with XMM - the XMM Large Scale Structure survey, XMM-LSS - with the scope of extending the cosmological tests attempted using ROSAT cluster samples to two redshift bins between 0<z<1 while maintaining the precision of earlier studies. Two main goals have constrained the survey design: the evolutionary study of the cluster-cluster correlation function and of the cluster number density. The results are promising and, so far, in accordance with our predictions as to the survey sensitivity and cluster number density. The feasibility of the programme is demonstrated and further X-ray coverage is awaited in order to proceed with a truly significant statistical analysis. (Abridged)

An impressive amount of different astrophysical data converges towards the picture of a spatially flat universe undergoing a today phase of accelerated expansion. The nature of the dark energy dominating the energy content of the universe is still unknown and a lot of different scenarios are viable candidates to explain cosmic acceleration. Most of the methods employed to test these cosmological models are essentially based on distance measurements to a particular class of objects. A different method, based on the lookback time to galaxy clusters and the age of the universe, is used here. In particular, we constrain the characterizing parameters of three classes of dark energy cosmological models to see whether they are in agreement with this kind of data, based on time measurements rather than distance observations.

We present a combined study of the colour-magnitude relation, colour distribution and luminosity function of a sample of 24 clusters at redshifts 0.3<z<1. The sample is largely composed of X-ray selected/detected clusters. Most of the clusters at redshifts z<0.6 display X-ray luminosity or richness typical of poor clusters or groups, rather than the more typical, massive clusters studied in literature at redshifts z>~0.3. All our clusters, including groups, display a colour-magnitude relation consistent with a passively evolving stellar population formed at a redshift z_f>~2, in accordance with observed galaxy populations in more massive clusters studied at comparable redshifts. Colours and luminosity functions show that the cluster galaxy population is consistent with the presence of at least two components: old systems formed at high redshift that have evolved passively from that epoch, together with a galaxy population displaying more recent star formation. The former population forms at ~2< z_f <~5, the latter at redshifts z<1. A model in which stars do not evolve is clearly rejected by both by the colour of reddest galaxies and by the characteristic luminosity m* measures. All clusters (with one possible exception) are detected independently by an almost three dimensional optical search employing sky position and colour - this despite the primary X-ray selection and low X-ray flux/optical richness displayed by the majority of the sample.

We present first results on the use of a colour space filter for detecting galaxy clusters at cosmological redshifts in the XMM Large Scale Structure Survey. All clusters studied, but one, are successfully colour--detected in spite of their large redshift (0.3<z<1.0), X-ray selection and intrinsic low richness (R=0 or below). We experimentally show that the cluster redshift can be derived, with good accuracy, from the colour of the red sequence, at 0.06<z<1 using a sample of about 160 clusters.

The XMM--LSS project is detecting distant clusters of low mass, quite comparable in mass to the ones in the local universe. This allows a direct comparison of galaxy properties at different redshifts in ``similar'' clusters. We present here first results on the evolution of the reddest galaxies in ~25 clusters/groups at ~0.3< z <~1.0 and for the whole galaxy population in the same clusters. The emerging picture from the current study is that the counterparts of present day clusters tend to show two or more distinct populations : a relatively old (z_f>2-5) population evolving passively together with a younger population, ostensibly responsible for the apparent brightening of the characteristic magnitudes, m*.

We present five newly found galaxy clusters at z>0.6 from the XMM Large-Scale Structure Survey (XMM-LSS). All five objects are extended X-ray sources in the XMM images. For three of them we have sufficient spectroscopically confirmed member galaxies that an estimate of the velocity dispersion is possible: XLSSC 001 at z=0.613 and \sigma_V=867^{+80}_{-60} km/s, XLSSC 002 at z=0.772 and \sigma_V=524^{+267}_{-116} km/s and XLSSC 003 at z=0.839 and \sigma_V=780^{+137}_{-75} km/s. These three clusters have X-ray bolometric luminosities L_X \sim 1-3 \times 10^{44} erg/s and temperatures 2-4 keV, and consequently are less massive than previously known clusters at similar redshifts, but nevertheless they follow the low redshift scaling relations between L_X, T and \sigma_V, within the limits of the measurement errors. One of the clusters, XLSSC 004, is detected independently as an overdensity of galaxies of a colour R-z'=1.4 that matches the redshift of the central galaxy z=0.87, although it cannot unambiguously be confirmed by the spectroscopic observations alone. The highest redshift candidate cluster pertaining to this paper, XLSSC 005, is most likely a double cluster complex at a redshift around unity, associated with an extended X-ray source with probable substructure.

We have observed three clusters at z~0.7, of richness comparable to the low redshift sample of Butcher & Oemler (BO), and determined their fraction of blue galaxies. When adopting the standard error definition, two clusters have a low blue fraction for their redshifts, whereas the fraction of the third one is compatible with the expected value. A detailed analysis of previous BO-like studies that adopted different definitions of the blue fraction shows that the modified definitions are affected by contaminating signals: colour segregation in clusters affects blue fractions derived in fixed metric apertures, differential evolution of early and late type spirals potentially affects blue fractions derived with a non standard choice of the colour cut, the younger age of the Universe at high redshift affects blue fractions computed with a colour cut taken relatively to a fixed non evolving colour. Adopting these definitions we find largely varying blue fractions. This thorough analysis of the drawbacks of the different possible definitions of the blue fraction should allow future studies to perform measures in the same scale. Finally, if one adopts a more refined error analysis to deal with BO and our data, a constant blue fraction with redshift cannot be excluded, showing that the BO effect is still far from being detected beyond doubt.

The cluster and field luminosity functions (LFs) determined on large homogeneous samples (N>2200 galaxies each) are almost indistinguishable, down to M^*+4 in the r and i filters, hence suggesting that the effect of the cluster environment on the galaxy properties does not affect the galaxy luminosity function in red bands. The similarity of the red band LFs in different environments suggests that the galaxy mass function is preserved during the galaxy infall in the cluster. By analyzing a large sample of galaxies in clusters, ideal from many points of view (multicolor data, large size, many clusters, metric magnitudes) we found that luminosity evolution is required by the data but the latter do not unambiguously derive its flavour if a differential luminosity evolution between bright and faint galaxies is allowed. We show that the LF parameters (slope, characteristic magnitude and their evolution) and errors depend on assumptions in a way seldom recognized in literature. We also point out logical inconsistencies between hypothesis assumed in deriving literature LF and presented results, suggesting caution in interpreting similar published results.

We studied the homogeneity, across clusters, of the color of the red sequence (the intercept of the color--magnitude relation) of 158 clusters and groups detected in the Early Data Release of the Sloan Digital Sky Survey (EDR-SDSS) in the redshift range 0.06<z<0.34. We found a high degree of homogeneity: the color of the red sequence shows an intrinsic scatter of 0.02 mag across clusters, suggesting that either galaxies on the red sequence formed a long time ago (z>2) or else their star formation is universally delayed with preservation of a small spread in age formation. The latter possibility is ruled out by the mere existence of galaxies at high redshift. While the old age of ellipticals was already been claimed for a small heterogeneous collection of clusters, most of which are rich ones, we found that it holds for ten to one hundred large sample, representative of all clusters and groups detected on the EDR-SDSS. Hence we suggest the possible universality of the color of the galaxies on the red sequence. The observed red sequence color does not depend on cluster/group richness at a level of 0.02 mag, while a ~0.23 mag effect is expected according to the hierarchical prediction. Therefore, the stellar population of red sequence galaxies is similar in clusters and groups, in spite of different halo histories. Finally, since the observed rest--frame color of the red sequence does not depend on environment and redshift, it can be used as a distance indicator, with an error sigma_z=0.018. [Abridged]

We (Pierre et al. 2003) have designed a medium deep large area X-ray survey with XMM - the XMM Large Scale Structure survey, XMM-LSS - with the scope of extending the cosmological tests attempted using ROSAT cluster samples to two redshift bins between 0<z<1 while maintaining the precision of earlier studies. The optimal survey design was found to be an 8x8 deg2 area, paved with 10 ks XMM pointing separated by 20 arcmin. This area is the target of several complementary surveys, from the ultraviolet to the radio wavelengths, allowing, beside cosmological studies, detailed studies of the objects included in our area.

We studied the scatter across clusters of the color of the red sequence in a representative and large sample of clusters (more than 200) detected on the Early Data Release of the Sloan Digital Sky Survey (EDR--SDSS) in the redshift range 0.06<z<0.34. We found an extreme degree of homogeneity in the color of the red sequence (the intrinsic scatter is about 0.02 to 0.03 mag) suggesting that either galaxies on the red sequence formed a long time ago (z>2) or else their star formation is universally delayed with preservation of a small spread in age formation. The latter possibility is ruled out by the mere existence of galaxies at high redshift. While the old age of ellipticals was already been claimed for a small heterogeneous collection of clusters, most of which are rich ones, we found that it holds for a ten to one hundred larger sample, representative of all clusters and groups detected on the EDR--SDSS. Hence, we claim the possible universality of the color of the galaxies on the red sequence. Furthermore, the sample includes a large number of very poor clusters (also called groups), not studied in previous works, for which the hierarchical and monolithic scenarios of elliptical formation predict different colors for the brightest ellipticals. The observed red sequence color does not depend on cluster/group richness at a level of 0.02 mag, while a 0.23 mag effect is expected according to the hierarchical prediction. Therefore, the stellar population of red sequence galaxies is similar in clusters and groups, in spite of different halo histories. Finally, we propose the use of color of the red sequence as distance indicator. (Abridged)

The multiplicity function (MF) of groups and clusters of galaxies is determined using galaxy catalogues extracted from a set of Digitized Palomar Sky Survey (DPOSS) plates. The two different types of structures (of low and high richness) were identified using two different algorithms: a modified version of the van Albada method for groups, and a peak finding algorithm for larger structures. In a 300 sq. deg. area up to z<0.2, we find 2944 groups and 179 clusters. Our MF covers a wide range of richnesses, from 2 to 200, and the two MF's derived by the two algorithms match smoothly without the need for additional conditions or normalisations. The resulting multiplicity function, of slope alpha= -2.08 +- 0.07, strongly resembles a power law.

The surface brightness distribution (SBD) function describes the number density of galaxies as measured against their central surface brightness. Because detecting galaxies with low central surface brightnesses is both time-consuming and complicated, determining the shape of this distribution function can be difficult. In a recent paper Cross, et al. suggested a bell-shaped SBD disk-galaxy function which peaks near the canonical Freeman value of 21.7 and then falls off significantly by 23.5 B mag arcsec-2. This is in contradiction to previous studies which have typically found flat (slope=0) SBD functions out to 24 - 25 B mag arcsec^-2 (the survey limits). Here we take advantage of a recent surface-brightness limited survey by Andreon & Cuillandre which reaches considerably fainter magnitudes than the Cross, et.al sample (M_B reaches fainter than -12 for Andreon & Cuillandre while the Cross, et.al sample is limited to M_B < -16) to re-evaluate both the SBD function as found by their data and the SBD for a wide variety of galaxy surveys, including the Cross, et al. data. The result is a SBD function with a flat slope out through the survey limits of 24.5 B mag arcsec^-2, with high confidence limits.

We present a neural network based approach to the determination of photometric redshift. The method was tested on the Sloan Digital Sky Survey Early Data Release (SDSS-EDR) reaching an accuracy comparable and, in some cases, better than SED template fitting techniques. Different neural networks architecture have been tested and the combination of a Multi Layer Perceptron with 1 hidden layer (22 neurons) operated in a Bayesian framework, with a Self Organizing Map used to estimate the accuracy of the results, turned out to be the most effective. In the best experiment, the implemented network reached an accuracy of 0.020 (interquartile error) in the range 0<zphot<0.3, and of 0.022 in the range 0<zphot<0.5.

Using new observations of the galaxy cluster A118 at intermediate redshift (z=0.31) in the Ks band, we were able to detect the cluster from the center to half the Abell radius (1.5 Mpc, Ho=50 km/s/Mpc) and possibly to 2.0 Mpc. The analysis of both the spatial distribution of galaxies of various luminosities and of the luminosity function (LF) of galaxies in different cluster locations strongly confirms and extends to larger clustercentric radii the luminosity segregation found in a previous analysis of this cluster restricted to a smaller cluster area: there is an excess of bright galaxies in the cluster core (inside 250 Kpc) or a deficit of dwarfs in the remain part of the cluster. Outside the cluster core and as far as 1.5 or even 2 Mpc, the giant-to-dwarf ratio is constant. Because of the luminosity segregation, the LF of the AC118 shows a larger number of bright galaxies per unit dwarf in the core than in other cluster locations. All non-core LFs, computed at several cluster locations, are compatible each other. These results hold both including or excluding the galaxies located in an overdensity found in the far South of AC118 and in the second clump in galaxy density at the cluster North-West. Since the near-infrared emission is a good tracer of the stellar mass, we interpret the segregation found as a mass segregation.

[Abridged] We determined the luminosity function (LF) of the Coma cluster galaxies selected by luminosity, and the LF bi-variate in central brightness. The Coma cluster and control fields were imaged using the CFH12K (42x28 arcmin) and UH8K (28x28 arcmin) wide-field cameras at the CFHT. Selected HST images were used for testing. Quantities were derived from measurements in at least two colors, which have the following features: (1) Galaxies as faint as three times the luminosity of the brightest globular clusters are in the completeness region of our data. (2) We have a complete census (in the explored region) of low surface brightness galaxies with central surface brightness galaxies almost as low as the faintest so far cataloged ones. (3) The explored area is among the largest ever sampled with CCDs at comparable depth for any cluster of galaxies. (4) The error budget includes all sources of errors known to date. Using HST images we also discovered that blends of globular clusters, not resolved in individual components due to seeing, look like dwarf galaxies when observed from the ground and are numerous and bright. The derived Coma LF is relatively steep (alpha=-1.4) over the 11 magnitudes sampled, but the slope and shape depend on color. A large population of faint low surface brightness galaxies was discovered, representing the largest contributor (in number) to the LF at faint magnitudes. We found a clear progression for a faintening of the LF from high surface brightness galaxies (mu~20 mag/arcsec2) to galaxies of very faint central brightnesses (mu~24.5 mag/arcsec2), and some evidence for a steepening. Compact galaxies, usually classified as stars and therefore not included in the LF, are found to be a minor population in Coma.

By using high-resolution and deep Ks band observations of early-type galaxies of the nearby Universe and of a cluster at z=0.3 we show that the two luminosity functions (LFs) of the local universe derived from 2MASS data miss a fair fraction of the flux of the galaxies (more than 20 to 30%) and a whole population of galaxies of central brightness fainter than the isophote used for detection, but bright enough to be included in the published LFs. In particular, the fraction of lost flux increases as the galaxy surface brightness become fainter. Therefore, the so far derived LF slopes and characteristic luminosity as well as luminosity density are underestimated. Other published near-infrared LFs miss flux in general, including the LF of the distant field computed in a 3 arcsec aperture.

By means of CCD photometry in three bands (Gunn g, r, i) we investigate the existence of 12 candidate clusters extracted via a model independent peak finding algorithm (\cite{memsait}) from DPOSS data. The derived color-magnitude diagrams allow us to confirm the physical nature of 9 of the cluster candidates, and to estimate their photometric redshifts. Of the other candidates, one is a fortuitous detection of a true cluster at z~0.4, one is a false detection and the last is undecidable on the basis of the available data. The accuracy of the photometric redshifts is tested on an additional sample of 8 clusters with known spectroscopic redshifts. Photometric redshifts turn out to be accurate within z~0.01 (interquartile range).

The composite galaxy luminosity function (hereafter LF) of 39 Abell clusters of galaxies is derived by computing the statistical excess of galaxy counts in the cluster direction with respect to control fields. Due to the wide field coverage of the digitised POSS-II plates, we can measure field counts around each cluster in a fully homogeneous way. Furthermore, the availability of virtually unlimited sky coverage allows us to directly compute the LF errors without having to rely on the estimated variance of the background. The wide field coverage also allows us to derive the LF of the whole cluster, including galaxies located in the cluster outskirts. The global composite LF has a slope alpha ~ -1.1+/-0.2 with minor variations from blue to red filters, and M* ~ -21.7,-22.2,-22.4 mag (H_0=50 km/s/Mpc) in g, r and i filters, respectively. These results are in quite good agreement with several previous determinations and in particular with the LF determined for the inner region of a largely overlapping set of clusters, but derived making use of a completely different method for background subtraction. The similarity of the two LFs suggests the existence of minor differences between the LF in the cluster outskirts and in the central region, or a negligible contribution of galaxies in the cluster outskirts to the global LF.

By means of wide-field (42x28 arcmin) CCD images of the Coma cluster of galaxies and of control fields, taken at the Canada-France-Hawaii telescope with the CFH12K and UH8K cameras, we determine the luminosity function bi-variate in central brightness. We found a clear progression for a steeping and weakening of the FL going from high surface brightness galaxies (mu~20$ mag/arcsec2) to galaxies of very-faint central brightnesses (mu 24.5 mag/arcsec2). Compact galaxies, usually rejected in the star/galaxy classification, are found to be a minor population in Coma. A huge population of low surface brightness galaxies red and faint are discovered, representing the largest contributor to the LF at faint magnitudes. Among them, there could be the remnant of high redshift starbursting galaxies responsible for the Butcher-Oemler effect.

We present the composite luminosity function (hereafter LF) of galaxies for 81 Abell clusters studied in our survey of the Northern Hemisphere, using DPOSS data processed by the CRoNaRio collaboration. The derived LF is very accurate due to the use of homogeneous data both for the clusters and the control fields and to the local estimate of the background, which takes into account the presence of large-scale structures and of foreground clusters and groups. The global composite LF is quite flat down to $M^*+5$ has a slope $\alpha\sim-1.0\pm0.2$ with minor variations from blue to red filters, and $M^*\sim-21.8,-22.0,-22.3$ mag ($H_0=50$ km s$^{-1}$ Mpc$^{-1}$) in the $g, r$ and $i$ filters, respectively (errors are detailed in the text). We find a significant difference between rich and poor clusters thus arguing in favour of a dependence of the LF on the properties of the environment.

[Abridged] The Coma cluster luminosity function (LF) from ultraviolet (2000 AA) to the near-infrared (H band) is summarized. In the UV the LF is very steep, much steeper than in the optical. The steep Coma UV LF implies that faint and bright galaxies give similar contributions to the total UV flux and to the total metal production rate. The Coma UV LF is dominated by star forming galaxies, not by massive and large galaxies. The optical Coma LF is relatively steep (alpha=-1.4) over the 11 magnitudes sampled. We found a clear steeping of the FL going from B to R bands, indicative of the presence of a large number of red dwarfs, as faint as three bright globular clusters. Furthermore, using Hubble Space Telescope images, we discover that blends of globular clusters, not resolved in individual components due to seeing, look like dwarf galaxies when observed from the ground and are numerous and bright. The existence of these fake extended sources increases the steepness of the LF at faint magnitudes, if not deal on. The near-infrared LF was computed on a near-infrared selected sample of galaxies which photometry is complete down to the typical dwarf (M^*+5) luminosity. The Coma LF can be described by a Schechter function with intermediate slope (alpha~-1.3), plus a dip at M_H~-22 mag. The shape of the Coma LF in H band is quite similar to the one found in the B band. The similarity of the LF in the optical and H bands implies that in the central region of Coma there is no new population of galaxies which is too faint to be observed in the optical band (because dust enshrouded, for instance), down to the magnitudes of dwarfs.

The deep near-infrared luminosity function of AC118, a cluster of galaxies at z=0.3, is presented. AC118 is a bimodal cluster, as evidenced both by our near-infrared images of lensed galaxies, by public X-ray Rosat images and by the spatial distribution of bright galaxies. Taking advantage of the extension and depth of our data, which sample an almost unexplored region in the depth vs. observed area diagram, we derive the luminosity function (LF), down to the dwarf regime (M*+5), computed in several cluster portions. The overall LF, computed on a 2.66 Mpc2 areas (H_0=50 km/s/Mpc), has an intermediate slope (alpha=-1.2). However, the LF parameters depend on the surveyed cluster region: the central concentration has 2.6^{+5.1}_{-1.7} times more bright galaxies and 5.3^{+7.2}_{-2.3} times less dwarfs per typical galaxy than the outer region, which includes galaxies at an average projected distance of ~580 kpc (errors are quoted at the 99.9 % confidence level). The LF in the secondary AC118 clump is intermediate between the central and outer one. In other words, the near-infrared AC118 LF steepens going from high to low density regions. At an average clustercentric distance of ~580 kpc, the AC118 LF is statistically indistinguishable from the LF of field galaxies at similar redshift, thus suggesting that the hostile cluster environment plays a minor role in shaping the LF at large clustercentric distances, while it strongly affects the LF at higher galaxy density.

We present the composite luminosity function (hereafter LF) of galaxies for 80 Abell clusters studied in our survey of the Northern Hemisphere, using DPOSS data in the framework of the CRoNaRio collaboration. Our determination of the LF has been computed with very high accuracy thanks to the use of homogeneous data both for the clusters and the control fields and to a local estimate of the background, which takes into account the presence of large-scale structures and of foreground clusters and groups. The global composite LF is quite flat down to M^+5 and it has a slope alpha=-1.0 +/- 0.2 with minor variations from blue to red filters, and M^* = -21.9,-22.0,-22.3 mag (H_0=50 km s^-1 Mpc^-1) in g, r and i filters, respectively (errors are detailed in the text). We find a significant difference between rich and poor clusters supporting the existence of an LF's dependence on the environment.

We discuss the implementation and validation of a procedure aimed to detect groups and clusters of galaxies in DPOSS catalogs.

[Abriged] Astronomical Wide Field Imaging performed with new large format CCD detectors poses data reduction problems of unprecedented scale which are difficult to deal with traditional interactive tools. We present here NExt (Neural Extractor): a new Neural Network (NN) based package capable to detect objects and to perform both deblending and star/galaxy classification in an automatic way. Traditionally, in astronomical images, objects are first discriminated from the noisy background by searching for sets of connected pixels having brightnesses above a given threshold and then they are classified as stars or as galaxies through diagnostic diagrams having variables choosen accordingly to the astronomer's taste and experience. In the extraction step, assuming that images are well sampled, NExt requires only the simplest a priori definition of "what an object is" (id est, it keeps all structures composed by more than one pixels) and performs the detection via an unsupervised NN approaching detection as a clustering problem which has been thoroughly studied in the artificial intelligence literature. In order to obtain an objective and reliable classification, instead of using an arbitrarily defined set of features, we use a NN to select the most significant features among the large number of measured ones, and then we use their selected features to perform the classification task. In order to optimise the performances of the system we implemented and tested several different models of NN. The comparison of the NExt performances with those of the best detection and classification package known to the authors (SExtractor) shows that NExt is at least as effective as the best traditional packages.

We exploited the large areal coverage offered by the DPOSS to analyze the outermost regions of the galactic globular cluster M 92. Two independent photometric reduction programs (SKICAT and DAOPHOT) were used to construct a color-magnitude diagram and a surface density profile for this cluster, based on J- and F-band DPOSS plates. A strong similarity has been found in the performance of the two programs in the low--crowded outermost cluster regions. After removing the background contribution, we obtained the cluster outer surface density profile down to a surface brightness magnitude of mu_V ~ 31 mag/sq.arcsec and matched it with the inner profile of Trager et al. 1995. The profile shapes match very well: since our data are uncalibrated, the shift in magnitude between the profiles has been also used to calibrate our profile. The analysis shows that the cluster has an extra tidal halo extending out to ~ 30 arcmin from the cluster center at a 3 sigma level over the background noise. This halo is revealed to be almost circular.

We present the near-infrared H band luminosity function (hereafter LF) of the Coma cluster of galaxies. It is the deepest ever computed in the near-infrared, for any type of environment, extending over 7 magnitudes, down to ~M_H*+6. The LF was computed on a near-infrared selected sample of galaxies which photometry, complete down to the typical dwarf luminosity, is presented in a companion paper. The Coma LF can be described by a Schechter function with intermediate slope (alpha~-1.3), plus a dip at M_H~-22 mag. The shape of the Coma LF in H band is quite similar to the one found in the B band and, with less confidence, to the R band LF as well. The similarity of the LF in the optical and H bands implies that in the central region of Coma there is no new population of galaxies which is too faint to be observed in the optical band (because dust enshrouded, for instance), down to the magnitudes of dwarfs. The exponential cut of the LF at the bright end is in good agreement with the one derived from shallower near-infrared samples of galaxies, both in clusters and in the field. This fact is suggestive of a similarity of the tip of the mass function of galaxies, irrespective of the environment where they are found. The dip at M_H~-22 mag is instead unique among all the so far measured near-infrared LF, although several published observations are not deep enough or spanning a suitable wide field to distinctly detect this feature. The faint end of the LF, reaching M_H~-19 mag (roughly M_B~-15), is steep, but less than previously suggested from shallower near-infrared observations of an adjacent region in the Coma cluster. The differences between our measured LF and that measured previously in other regions suggests a dependency on environment of the faint end of the mass function (below M*+2.5).

We present extensive and accurate photometry in the near-infrared H band of a complete sample of objects in an area of about 400 arcmin2 toward the Coma cluster of galaxies. The sample, including about 300 objects, is complete down to H~17 mag, the exact value depending on the type of magnitude (isophotal, aperture, Kron) and the particular region studied. This is six magnitudes below the characteristic magnitude of galaxies, well into the dwarfs' regime at the distance of the Coma cluster. For each object (star or galaxy) we provide aperture magnitudes computed within five different apertures, the magnitude within the 22 mag arcsec2 isophote, the Kron magnitude and radius, magnitude errors, as well as the coordinates, the isophotal area, and a stellarity index. Photometric errors are 0.2 mag at the completness limit. This sample is meant to be the zero-redshift reference for evolutionary studies of galaxies.

Low Surface Brightness Galaxies (LSBGs) today constitute the less known fraction of galaxies: their number density and their physical properties (e.g. luminosity, colour, dynamics) are still largely uncertain. This is mainly due to the difficoulty in finding them (since their presence goes often unnoticed when astronomical images are analyzed with standard methods) and in measuring their properties. We present here the first results of a novel application to DPOSS plates, a project within the CRoNaRio collaboration. The aim is to build a large catalog of LSBGs, with much better understanding of completness and biases. First results are quite encouraging: in the same regions covered by existing catalogues, done on the same plates, we find several new LSBGs in addition to the known ones.

We present the composite luminosity function of galaxies for 24 Abell clusters studied in our survey of the Northern Hemisphere, using DPOSS data in the framework of the CRoNaRio collaboration. Our determination of the luminosity function has been computed with very high accuracy thanks to 1) the use of homogeneous data both in the cluster and in the control field, 2) a local estimate of the background, which takes into account the presence of large-scale structures (which tipically make the background contribution larger) and of foreground clusters and groups; 3) the inclusion in the error calculation of the variance of background counts, which is tipically larger than Poissonian fluctuation. We plan at least a tenfold increase of the number of clusters.

The projected multiplicity function of galaxies gives (per square degree) the density of galaxy aggregates formed by N members. We use the CRoNaRio matched catalogues to derive the low N (N<15) multiplicity function from D-POSS data. The van Albada (1982) algorithm was implemented and used to identify candidate groups in the CRoNaRio catalogues. In absence of redshift surveys complete down to the magnitude limit of the DPOSS material and covering a wide enough area, the performances of the algorithm were tested on realistically simulated catalogues. The application to a set of 13 CRoNaRio catalogues allowed us to derive a multiplicity function as accurate as those available in literature obtained from redshift surveys.

We present preliminary results of a simple peak finding algorithm applied to catalogues of galaxies, extracted from the Second Palomar Sky Survey in the framework of the CRoNaRio project. All previously known Abell and Zwicky clusters in a test region of 5x5 sq. deg. are recovered and new candidate clusters are also detected. This algorithm represents an alternative way of searching for galaxy clusters with respect to that implemented and tested at Caltech on the same type of data (Gal et al. 1998).

The UV (2000 A) luminosity function (hereafter UV LF) of Coma cluster galaxies, based on more than 120 members, is computed as the statistical difference between counts in the Coma direction and in the field. Our UV LF is an up-date of a preliminary constrain on the UV LF previously computed without the essential background counts. The UV LF is well described by a power law with slope alpha~0.46, or equivalently, by a Schechter function with M* much brighter than the brightest cluster galaxy and with a slope alpha_S~-2.0 or larger. In spite of what happens in the optical band, low luminosity galaxies give a large contribution to the integral luminosity, and by inference, to the total metal production rate. Galaxies blue in UV-b and/or b-r dominate the Coma cluster UV LF, both in number and luminosity. The major source of error in the estimate of the UV LF cames from the background determination in the Coma direction, which is still uncertain, even though constrained at high and low amplitudes by redshift surveys covering the studied field.

One of the most relevant problems in the extraction of scientifically useful information from wide field astronomical images (both photographic plates and CCD frames) is the recognition of the objects against a noisy background and their classification in unresolved (star-like) and resolved (galaxies) sources. In this paper we present a neural network based method capable to perform both tasks and discuss in detail the performance of object detection in a representative celestial field. The performance of our method is compared to that of other methodologies often used within the astronomical community.

Using PSPC {\it Rosat} data, we measure x-ray surface brightness profiles, size and luminosity of the Butcher-Oemler (BO) sample of clusters of galaxies. The cluster x-ray size, as measured by the Petrosian r_{\eta=2} radius, does not change with redshift and is independent from x-ray luminosity. On the other hand, the x-ray luminosity increases with redshift. Considering that fair samples show no-evolution, or negative luminosity evolution, we conclude that the BO sample is not formed from the same class of objects observed at different look-back times. This is in conflict with the usual interpretation of the Butcher-Oemler as an evolutionary (or redshift-dependent) effect, based on the assumption that we are comparing the same class of objects at different redshifts. Other trends present in the BO sample reflect selection criteria rather than differences in look-back time, as independently confirmed by the fact that trends loose strength when we enlarge the sample with x-ray selected sample of clusters. The variety of optical sizes and shapes of the clusters in the Butcher-Oemler sample, and the Malmquist-like bias, are the reasons for these selection effects that mimic the trends usually interpreted as changes due to evolution.

We constructed the composite Luminosity Function (LF) of cluster galaxies in the g,r and i bands from the photometry of a mixed (Abell and X-ray selected) sample of the cores of 65 clusters, ranging in redshift from 0.05 to 0.25. The composite LF has been obtained from complete samples of ~2200 galaxies in the magnitude range -23<M<-17.5 (-18 in i). Cluster membership has been determined on the basis of color-color plots for each cluster and the resulting outlier counts have been checked against field counts in the same bands. We find that the galaxy density of the environment determines the shape of the LF, in the sense that bright galaxies are more abundant in dense clusters.

Contradictory results on the evolution of lenticular galaxies (S0s) at intermediate redshift have been recently published. A careful analysis of the used methodology shows that in order to gain insights on the evolution of the morphological types it is necessary: (1) to use quantitative morphological types, (2) to propagate morphological errors in the quantity used to measure the evolution of the S0s and (3) to adopt the same morphological classes at high as at low redshift, testing the correspondence with standard types. Some published analysis fail to satisfy at least one of the above three conditions, thus undermining the claim for a young age of S0s of intermediate redshift clusters.

We find that the bright part (M_B<-18.5 mag) of luminosity function of each Hubble type has the same shape in three poor clusters (Virgo, Fornax and Centaurus) and in two rich clusters (Coma and Cl0939+4713). The fact that these type-dependent luminosity functions are invariant in shape from poor to rich environments gives support to the hypothesis that they may be universal. We also present in tabular form improved type-dependent luminosity functions, based on a much larger sample than previous works.

This paper have been withdraw by the autors, because of a too early submission.

Two contradictory results on the evolution of SO galaxies now exist in the recent literature; either S0s are old ($z_{formation}>2$) and are evoloving passively or most of them formed at $z<0.5$, as implied by the deficit of S0s in intermediate redshift ($z\sim0.5$) clusters. The resolution of this controversy may be that the apparent deficit of S0s has been derived from a quantity -- the E to S0 ratio -- which is prone to morphological classification errors. Once all sources of error are taken into account, the E to S0 ratios of clusters at different redshifts are fully compatible, and no additional creation of S0s at $z<0.5$ is required by the data. Furthermore, there is no deficit at all of S0s in the intermediate redshift cluster for which we have morphological types of very high quality, and thus derive an E to S0 ratio with a small error.

In Coma and two other clusters of galaxies the primary morphological segregation is the one with respect to a privileged direction. Only when this segregation is not considered it appears that the morphological types are segregated in density or in clustercentric distance.

We present an optical and near-infrared survey of galaxies in nearby clusters aimed at determining fundamental quantities of galaxies, such as multivariate luminosity function and color distribution for each Hubble type. The main characteristics of our survey are completeness in absolute magnitude, wide wavelength coverage and faint limiting magnitudes.

The CRONA project is a joint enterprise between Caltech and the astronomical observatory of Rome and Naples aimed to the photometric calibration of the Second Palomar Sky Survey (POSS-II) plates and to the construction of a complete photometric and astrometric catalogue of all objects present on the POSS-II. The calibrated final catalogues of objects will be made available to the community. We present here a comparison between the performances of the SKICAT and DAOPHOT packages applied to POSS-II images of the globular cluster M92.

We present the results of an isophotal shape analysis of galaxies in the Coma and Perseus clusters. These data, together with those of two previous papers, provide two complete samples of galaxies with reliable Hubble types in rich clusters: 1) all galaxies brighter than m_b = 16.5 falling within one degree (=2.3 Mpc) from the center of the Coma cluster (187 galaxies), 2) all galaxies brighter than m_{Zwicky}=15.7 in a region of 5 deg 3' times 5 deg 27' around the center of the Perseus cluster (139 galaxies). These two complete samples cover 5 orders of magnitude in galaxy density and span areas of 91 and 17 Mpc^2, clustercentric radii up to 2.3 and 6.4 Mpc, for Perseus and Coma respectively. They will be used in subsequent papers to study the dependence of galaxy types on cluster environment and as reference samples in comparisons with distant clusters.

The Gunn g, r and i CCD images of a representative sample of 17 Brightest Cluster Galaxies (BCM) have been analyzed in order to derive surface brightness and color profiles, together with geometrical parameters like eccentricity and position angle. The sample includes both X-ray and optically selected clusters, ranging in redshift from z=0.049 to z=0.191. We find that BCMs are substantially well described by de Vaucouleurs' law out to radii of $\sim 60-80$ kpc, and that color gradients are generally absent. Only in two cases we find a surface brightness excess with respect to the $r^{1/4}$ law, which for A150 is coupled with a change in the $g-r$ color. The rest frame colors of BCMs do not show any intrinsic dispersion. By parametrizing the environment with the local galaxy number density, we find that it is correlated with the BCM extension, i.e. BCMs with larger effective radii are found in denser environments.

We have performed an isophotal analysis of galaxies in the distant (z=0.4) cluster Cl0939+4713 (Abell 851), using post-refurbished Hubble Space Telescope images. Morphological type estimates for the galaxies are given. A rigorous comparison of the properties of the types shows that early-type galaxies in Cl0939+4713 are, within the present statistical or systematic errors, indistinguishable from their counterparts in Coma in all their studied properties, namely ellipticity profile, slope of the color-magnitude relation from the near-ultraviolet to the near-infrared colors, mean surface brightness and luminosity function in the restframe photographic J band, homogeneity in color (around the color-magnitude relation). Furthermore ellipticals and lenticulars are separately homogeneous in their mean surface brightness, and have similar relative mean surface brightnesses in both clusters. Spirals are overabundant in Cl0939+4713 with respect to Coma, but twice less than previously estimated, and are more similar to field spirals than to cluster spirals. We suggest that the differences in the photometric properties of the spirals in the two clusters arise from differences in cluster gas density distribution which ultimately bring star formation in the spirals to a stop in Coma, but not in Cl0939+4713. The morphological types are segregated along a privileged direction that coincides with the position angle of the major axis of the outer X-ray isophotes in the cluster, just like in Coma and Perseus.

The traditional method of morphological classification, by visual inspection of images of uniform quality and by reference to standards for each type, is critically examined. The rate of agreement among traditional morphologists on the morphological type of galaxies is estimated from published classification works, and is estimated at about 20 %, when galaxies are classified into three bins (E, S0, S+Irr). The advantages of the quantitative method of structural classification for classifying galaxies in clusters are outlined. This method is based on the isophotal analysis of galaxy images, and on the examination of quantitative structural parameters derived from this analysis, such as the profiles of luminosity, ellipticity and deviations from ellipticity of the galaxy. The structural and traditional methods are compared on a complete sample of 190 galaxies in the Coma cluster. The morphological types derived by both methods agree to within 15 or 20 %, the same rate as among traditional morphologists alone, thus showing that our morphological classes do coincide with the traditional ones. The galaxies with discrepant types are mostly faint (mag_B > 16.0), or have features typical of spirals, but which have not been detected, noticed or taken into account by traditional morphologists. The rate of agreement is also good for galaxies in a distant cluster (z ~0.4). The structural method, gives highly reproducible results, never reached by traditional estimations of the morphological type. Thus, the morphological types obtained with this method should be preferred. The advantages of the method are further demonstrated by new results on the properties of galaxies in clusters.

We present the results of an isophotal shape analysis of three samples of galaxies in the Coma cluster. Quantitative morphology, together with structural and photometric parameters, is given for each galaxy. Special emphasis has been placed on the detailed classification of early-type galaxies. The three samples are: i) a sample of 97 early-type galaxies brighter than m$_B$ = 17.00 falling within one degree from the center of the Coma cluster; these galaxies were observed with CCD cameras, mostly in good to excellent resolution conditions; ii) a magnitude complete sample of 107 galaxies of all morphological types down to m$_B$ = 17.00 falling in a circular region of 50 arcmin diameter, slightly offcentered to the North-West of the cluster center; the images for this and the next sample come from digitized photographic plates; iii) a complete comparison sample of 26 galaxies of all morphological types down to m$_R$ = 16.05 (or m$_B$ $\simeq$ 17.5), also in a region of 50 arcmin diameter, but centered 2.6 degrees West of the cluster center. The reliability of our morphological classifications and structural parameters of galaxies, down to the adopted magnitude limits, is assessed by comparing the results on those galaxies for which we had images taken with different instrumentation and/or seeing conditions, and by comparing our results with similar data from other observers.

Intrinsic color profiles of a sample of nine Brightest Cluster Galaxies (BCGs) are recovered from the observed color profiles by evaluating spurious gradients introduced by errors in the determination of the sky levels and by different seeing conditions between the observations. Isophote shapes and surface brightness profiles are presented for the four newly observed BCGs. Three out of nine BCGs show color gradients of the order of 0.10 mag per decade in radius. Five BCGs do not possess any color gradient larger than 0.01 mag per decade in radius. We do not see any correlation between the presence (or the sign) of a color gradient and the BCG morphology (slope and shape of its surface brightness profile), or the cluster richness. This argues against a strong and recent influence of the environment on the BCG morphology. The slope of the BCG surface brightness profile is correlated to the cluster richness, posing a constraint on the initial conditions, or on the physical mechanism that is responsible for the present BCG morphology. Finally, only one BCG can be a cD despite visual classification as such of all studied BCGs.