We measure the central stellar velocity dispersion function for quiescent galaxies in a set of nine northern clusters in the redshift range $0.18 < z < 0.29$ and with strong lensing arcs in Hubble Space Telescope images. The velocity dispersion function links galaxies directly to their dark matter halos. From dense SDSS and MMT/Hectospec spectroscopy we identify $\sim 380$ spectroscopic members in each cluster. We derive physical properties of cluster members including redshift, $D_{n}4000$, and central stellar velocity dispersion and we include a table of these measurements for 3419 cluster members. We construct the velocity dispersion functions for quiescent galaxies with $D_{n}4000 > 1.5$ and within $R_{200}$. The cluster velocity dispersion functions all show excesses at $\sigma \gtrsim 250 km s^{-1}$ compared to the field velocity dispersion function. The velocity dispersion function slope at large velocity dispersion ($\sigma > 160 km s^{-1}$) is steeper for more massive clusters, consistent with the trend observed for cluster luminosity functions. The spatial distribution of galaxies with large velocity dispersion at radii larger than $R_{200}$ further underscores the probable major role of dry mergers in the growth of massive cluster galaxies during cluster assembly.

We derive the distribution of the phantom dark matter in the eight classical dwarf galaxies surrounding the Milky Way, under the assumption that MOND is the correct theory of gravity. According to their observed shape, we model the dwarves as axi-symmetric systems, rather than spherical systems, as usually assumed. In addition, as required by the assumption of the MOND framework, we realistically include the external gravitational field of the Milky Way and of the large-scale structure beyond the Local Group. For the dwarves where the external field dominates over the internal gravitational field, the phantom dark matter has, from the star distribution, an offset of ~0.1-0.2 kpc, depending on the mass-to-light ratio adopted. This offset is a substantial fraction of the dwarf half-mass radius. For Sculptor and Fornax, where the internal and external gravitational fields are comparable, the phantom dark matter distribution appears disturbed with spikes at the locations where the two fields cancel each other; these features have little connection with the distribution of the stars within the dwarves. Finally, we find that the external field due to the large-scale structure beyond the Local Group has a very minor effect. The features of the phantom dark matter we find represent a genuine prediction of MOND and could thus falsify this theory of gravity, in the version we adopt here, if they were not observationally confirmed.

We use a new spherical accretion recipe tested on N-body simulations to measure the observed mass accretion rate (MAR) of 129 clusters in the CIRS and HeCS. The observed clusters cover the redshift range $0.01 < z < 0.30$ and the mass range $\sim 10^{14}-10^{15}h^{-1}\text{M}_\odot$. Based on three-dimensional mass profiles of simulated clusters reaching beyond the virial radius, our recipe returns MARs that agree with MARs based on merger trees. We adopt this recipe to estimate the MAR of real clusters based on measurement of the mass profile out to $\sim 3R_{200}$. We use the caustic method to measure the mass profiles to these large radii. We demonstrate the solidity of our estimates by applying the same approach to a set of mock redshift surveys of a sample of 2000 simulated clusters with median mass $M_{200}= 10^{14} {\,h^{-1}\rm{M_{\odot}}}$, and a sample of 50 simulated clusters with median mass $M_{200}= 10^{15} {\,h^{-1}\rm{M_{\odot}}}$: the median MARs based on the caustic mass profiles of the simulated clusters are unbiased and agree within $17\%$ with the median MARs based on the real mass profile of the clusters. The MAR of the CIRS and HeCS clusters increases with the mass and the redshift of the accreting cluster, in excellent agreement with the growth of clusters in the $\Lambda$CDM model.

We test if Refracted Gravity (RG) can describe the dynamics of disk galaxies without resorting to dark matter. RG is a classical theory of gravity where the standard Poisson equation is modified by the gravitational permittivity, $\epsilon$, a universal monotonic function of the local mass density. We use the rotation curves and the vertical velocity dispersions of 30 galaxies in the DiskMass Survey (DMS) to determine $\epsilon$. RG describes the kinematic profiles with mass-to-light ratios in agreement with SPS models, and disk thicknesses in agreement with observations, once observational biases are considered. Our results rely on setting the three free parameters of RG for each galaxy. However, we show that the differences of these parameters from galaxy to galaxy could be ascribed to statistical fluctuations. We adopt an approximate method to find a single set of parameters that might properly describe the kinematics of the entire sample and suggest the universality of $\epsilon$. We finally show that the RG models of the individual rotation curves can only partly describe the radial acceleration relation (RAR). Evidently, the RG models underestimate the observed accelerations by 0.1-0.3 dex at low Newtonian accelerations. Another problem is the strong correlations, at largely more than 5$\sigma$, between the residuals of the RAR models and three radially-dependent properties of galaxies, whereas the DMS data show considerably less significant correlations, at more than 4$\sigma$, for only two of them. These correlations might originate the non-null intrinsic scatter of the RG models, at odds with the observed intrinsic scatter of galaxy samples, different from DMS, which is consistent with 0. Further studies are required to assess if these discrepancies in the RAR originate from the DMS sample, which might not be ideal for deriving the RAR, or if they are genuine failures of RG.

We show how Conformal Gravity (CG) has to satisfy a fine-tuning condition to describe the rotation curves of disk galaxies without the aid of dark matter. Interpreting CG as a gauge natural theory yields conservation laws and their associated superpotentials without ambiguities. We consider the light deflection of a point-like lens and impose that the two Schwarzschild-like metrics with and without the lens are identical at infinite distances from the lens. The energy conservation law implies that the parameter $\gamma$ in the linear term of the metric has to vanish, otherwise the two metrics are physically inaccessible from each other. This linear term is responsible to mimic the role of dark matter in disk galaxies and gravitational lensing systems. Our analysis shows that removing the need of dark matter with CG thus relies on a fine-tuning condition on $\gamma$. We also illustrate why the results of previous investigations of gravitational lensing in CG largely disagree. These discrepancies derive from the erroneous use of the deflection angle definition adopted in General Relativity, where the vacuum solution is asymptotically flat, unlike CG. In addition, the lens mass is identified with various combinations of the metric parameters. However, these identifications are arbitrary, because the mass is not a conformally invariant quantity, unlike the conserved charge associated to the energy conservation law. Based on this conservation law and by removing the fine-tuning condition on $\gamma$, i.e. by setting $\gamma=0$, the energy difference between the metric with the point-like lens and the metric without it defines a conformally invariant quantity that can in principle be used for (1) a proper derivation of light deflection in CG, and (2) the identification of the lens mass with a function of the parameters $\beta$ and $k$ of the Schwarzschild-like metric.

We explore connections between brightest cluster galaxies (BCGs) and their host clusters. We first construct a HeCS-omnibus cluster sample including 227 galaxy clusters within $0.02 < z < 0.30$; the total number of spectroscopic members from MMT/Hectospec and SDSS observations is 52325. Taking advantage of the large spectroscopic sample, we compute physical properties of the clusters including the dynamical mass and cluster velocity dispersion ($\sigma_{cl}$). We also measure the central stellar velocity dispersion of the BCGs ($\sigma_{*,BCGs}$) to examine the relation between BCG velocity dispersion and cluster velocity dispersion for the first time. The observed relation between BCG velocity dispersion and the cluster velocity dispersion is remarkably tight. Interestingly, the $\sigma_{*, BCG} / \sigma_{cl}$ ratio decreases as a function of $\sigma_{cl}$ unlike the prediction from the numerical simulation of Dolag et al. (2010). The trend in $\sigma_{*, BCG} / \sigma_{cl}$ suggests that the BCG formation is more efficient in lower mass halos.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

We explore the structure of galaxy cluster Abell 2029 and its surroundings based on intensive spectroscopy along with X-ray and weak lensing observations. The redshift survey includes 4376 galaxies (1215 spectroscopic cluster members) within 40 arcmin of the cluster center; the redshifts are included here. Two subsystems, A2033 and a Southern Infalling Group (SIG) appear in the infall region based on the spectroscopy as well as on the weak lensing and X-ray maps. The complete redshift survey of A2029 also identifies at least 12 foreground and background systems (10 are extended X-ray sources) in the A2029 field; we include a census of their properties. The X-ray luminosities $L_{X}$ - velocity dispersions ($\sigma_{cl}$) scaling relations for A2029, A2033, SIG, and the foreground/background systems are consistent with the known cluster scaling relations. The combined spectroscopy, weak lensing, and X-ray observations provide a robust measure of the masses of A2029, A2033, and SIG. The total mass of the infalling groups (A2033 and SIG) is $\sim 60\%$ of the M200 of the primary cluster, A2029. Simple dynamical consid- erations suggest that A2029 will accrete these subsystems in next few Gyr. In agreement with simulations and with other clusters observed in a similar redshift range, the total mass in the A2029 infall region is comparable with the A2029 M200 and will mostly be accreted in the long-term future.

To investigate the dynamics of the galaxy cluster A2142, we compile an extended catalog of 2239 spectroscopic redshifts of sources, including newly measured 237 redshifts, within 30 arcmin from the cluster center. With the $\sigma$-plateau algorithm from the caustic method, we identify 868 members and a number of substructures in the galaxy distribution both in the outskirts, out to $\sim$3.5 Mpc from the cluster center, and in the central region. In the outskirts, one substructure overlaps a falling clump of gas previously identified in the X-ray band. These substructures suggests the presence of multiple minor mergers, which are responsible for the complex dynamics of A2142, and the absence of recent or ongoing major mergers. We show that the distribution of the galaxies in the cluster core and in several substructures are consistent with the mass distribution inferred from the weak lensing signal. Moreover, we use spatially-resolved X-ray spectroscopy to measure the redshift of different regions of the intracluster medium within $\sim$3 arcmin from the cluster center. We find a ring of gas near the two X-ray cold fronts identified in previous analyses and measure a velocity of this ring of $810 \pm 330 {\rm km s^{-1}}$ larger than the cluster mean velocity. Our analysis suggests the presence of another ring surrounding the core, whose velocity is $660 \pm 300 {\rm km s^{-1}}$ larger than the cluster velocity. These X-ray features are not associated to any optical substructures, and support the core-sloshing scenario suggested in previous work.

We develop the Blooming Tree Algorithm, a new technique that uses spectroscopic redshift data alone to identify the substructures and the surrounding groups of galaxy clusters, along with their member galaxies. Based on the estimated binding energy of galaxy pairs, the algorithm builds a binary tree that hierarchically arranges all the galaxies in the field of view. The algorithm searches for buds, corresponding to gravitational potential minima on the binary tree branches; for each bud, the algorithm combines the number of galaxies, their velocity dispersion and their average pairwise distance into a parameter that discriminates between the buds that do not correspond to any substructure or group, and thus eventually die, and the buds that correspond to substructures and groups, and thus bloom into the identified structures. We test our new algorithm with a sample of 300 mock redshift surveys of clusters in different dynamical states; the clusters are extracted from a large cosmological $N$-body simulation of a $\Lambda$CDM model. We limit our analysis to substructures and surrounding groups identified in the simulation with mass larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within 6 $h^{-1}$~Mpc from the cluster center, the technique recovers $ 80$\% of the real substructures and $ 60$\% of the surrounding groups; in $57$\% of the identified structures, at least 60\% of the member galaxies of the substructures and groups belong to the same real structure. These results improve by roughly a factor of two the performance of the best substructure identification algorithm currently available, the $\sigma$ plateau algorithm, and suggest that our Blooming Tree Algorithm can be an invaluable tool for detecting substructures of galaxy clusters and investigating their complex dynamics.

We apply a friends-of-friends algorithm to the HectoMAP redshift survey and cross-identify associated X-ray emission in the ROSAT All-Sky Survey data (RASS). The resulting flux limited catalog of X-ray cluster survey is complete to a limiting flux of $\sim3 \times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ and includes 15 clusters (7 newly discovered) with redshift $z \leq 0.4$. HectoMAP is a dense survey ($\sim1200$ galaxies deg$^{-2}$) that provides $\sim50$ members (median) in each X-ray cluster. We provide redshifts for the 1036 cluster members. Subaru/Hyper Suprime-Cam imaging covers three of the X-ray systems and confirms that they are impressive clusters. The HectoMAP X-ray clusters have an $L_{X} - {\sigma}_{cl}$ scaling relation similar to that of known massive X-ray clusters. The HectoMAP X-ray cluster sample predicts $\sim 12000 \pm3000$ detectable X-ray clusters in the RASS to the limiting flux, comparable with previous estimates.

We use the dense HectoMAP redshift survey to explore the properties of 104 redMaPPer cluster candidates. The redMaPPer systems in HectoMAP cover the full range of richness and redshift (0.08 $< z <$ 0.60). Fifteen systems included in the Subaru/Hyper Suprime-Cam public data release are bona fide clusters. The median number of spectroscopic members per cluster is $\sim20$. We include redshifts of 3547 member candidates listed in the redMaPPer catalog whether they are cluster members or not. We evaluate the redMaPPer membership probability spectroscopically. The scaled richness ({\lambda}rich/S) provided by redMaPPer correlates tightly with the spectroscopic richness regardless of the cluster redshift and appears to be a better mass proxy than the original richness, {\lambda}rich. The purity (number of real systems) in redMaPPer exceeds 90% even at the lowest richness; however, there is some incompleteness. Five massive galaxy clusters (M $\gtrsim 2 \times 10^{13}$ M$_{\odot}$) associated with X-ray emission in the HectoMAP region are missing from the catalog.

We use dense redshift surveys to explore the properties of galaxy clusters selected from the redMaPPer catalog of overdensities of red galaxies. Our new survey, HeCS-red (Hectospec Cluster Survey of red-sequence selected clusters), includes 10,589 new or remeasured redshifts from MMT/Hectospec observations of redMaPPer clusters at redshifts $z$=0.08-0.25 with large estimated richnesses (richness estimate $\lambda >64$). Our spectra confirm that each of these candidate clusters corresponds to an overdensity in redshift space. The redMaPPer photometric redshifts have a slight bias towards higher redshifts. We measure the scaling relation between velocity dispersion $\sigma_p$ and redMaPPer richness estimates $\lambda$. The observed relation shows intrinsic scatter of 24\% in velocity dispersion at fixed richness, and a range of a factor of two in measured $\sigma_p$ at fixed richness. We extend our analysis to HeCS-red-ext, a sample that includes several clusters selected by X-ray flux or SZ signal. The heterogeneous sample of 121 clusters in HeCS-red-ext shows similar intrinsic scatter, but the range of $\sigma_p$ at fixed richness increases to a factor of three. We evaluate the membership probability estimates $P_{mem}$ for individual galaxies provided by redMaPPer. The spectroscopic membership fraction is larger than $P_{mem}$ for $0.05\leq P_{mem}\leq 0.7$; conversely, it is smaller than $P_{mem}$ at $P_{mem}\geq 0.8$. We compare spectroscopic richness estimates to redMaPPer richness estimates and find good agreement on average, but a range of a factor of two in spectroscopic richness at fixed redMaPPer richness. Overall, within the high-richness and low-redshift cut of our sample, spectroscopically estimated parameters such as velocity dispersion correlate well with photometric richness estimates, although the relations contain substantial scatter.

Compact radio galaxies with a GHz-peaked spectrum (GPS) and/or compact-symmetric-object (CSO) morphology (GPS/CSOs) are increasingly detected in the X-ray domain. Their radio and X-ray emissions are affected by significant absorption. However, the locations of the X-ray and radio absorbers are still debated. We investigated the relationship between the column densities of the total ($N_{\mathrm{H}}$) and neutral ($N_{\mathrm{HI}}$) hydrogen to statistically constrain the picture. We compiled a sample of GPS/CSOs including both literature data and new radio data that we acquired with the Westerbork Synthesis Radio Telescope for sources whose X-ray emission was either established or under investigation. In this sample, we compared the X-ray and radio hydrogen column densities, and found that $N_{\mathrm{H}}$ and $N_{\mathrm{HI}}$ display a significant positive correlation, with $N_{\mathrm{HI}} \propto N_{\mathrm{H}}^b$, where $b=0.47$ and $b=0.35$, depending on the subsample. The $N_{\mathrm{H}}$ - $N_{\mathrm{HI}}$ correlation suggests that the X-ray and radio absorbers are either co-spatial or different components of a continuous structure. The correlation displays a large intrinsic spread that we suggest to originate from fluctuations, around a mean value, of the ratio between the spin temperature and the covering factor of the radio absorber, $T_{\rm s}/C_{\rm f}$.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Galaxies in clusters are strongly affected by their environment. They evolve according to several physical mechanisms that are active in clusters. Their efficiency can strongly depend on the orbital configuration of the galaxies. Our aim is to analyse the orbits of the galaxies in the cluster Abell 85, based on the study of the galaxy velocity anisotropy parameter. We have solved the Jeans equation under the assumption that the galaxies in A85 are collisionless objects, within the spherically symmetric gravitational potential of the virialized cluster. The mass of the cluster was estimated with X-ray and caustic analyses. We find that the anisotropy profile of the full galaxy population in A85 is an increasing monotonic function of the distance from the cluster centre: on average, galaxies in the central region (r/r200 < 0.3) are on isotropic orbits, while galaxies in the outer regions are on radial orbits. We also find that the orbital properties of the galaxies strongly depend on their stellar colour. In particular, blue galaxies are on less radial orbits than red galaxies. The different families of cluster galaxies considered here have the pseudo phase-space density profiles Q(r) and Qr(r) consistent with the profiles expected in virialized dark matter halos in $N$-body simulations. This result suggests that the galaxies in A85 have reached dynamical equilibrium within the cluster potential. Our results indicate that the origin of the blue and red colour of the different galaxy populations is the different orbital shape rather than the accretion time.

We carried out the deep spectroscopic observations of the nearby cluster A2151 with AF2/WYFFOS@WHT. The caustic technique enables us to identify 360 members brighter than $M_r = -16$ and within 1.3$R_{200}$. We separated the members into subsamples according to photometrical and dynamical properties such as colour, local environment and infall time. The completeness of the catalogue and our large sample allow us to analyse the velocity dispersion and the luminosity functions of the identified populations. We found evidence of a cluster still in its collapsing phase. The LF of the red population of A2151 shows a deficit of dwarf red galaxies. Moreover, the normalized LFs of the red and blue populations of A2151 are comparable to the red and blue LFs of the field, even if the blue galaxies start dominating one magnitude fainter and the red LF is well represented by a single Schechter function rather than a double Schechter function. We discuss how the evolution of cluster galaxies depends on their mass: bright and intermediate galaxies are mainly affected by dynamical friction and internal/mass quenching, while the evolution of dwarfs is driven by environmental processes which need time and a hostile cluster environment to remove the gas reservoirs and halt the star formation.

Based on an extensive redshift survey for galaxy cluster Abell 2029 and Coma, we measure the luminosity functions (LFs), stellar mass functions (SMFs) for the entire cluster member galaxies. Most importantly, we measure the velocity dispersion functions (VDFs) for quiescent members. The MMT/Hectospec redshift survey for galaxies in A2029 identifies 982 spectroscopic members; for 838 members we derive the central velocity dispersion from the spectroscopy. Coma is the only other cluster surveyed as densely. The LFs, SMFs and VDFs for A2029 and Coma are essentially identical. The SMFs of the clusters are consistent with simulations. The A2029 and Coma VDFs for quiescent galaxies have a significantly steeper slope than those of field galaxies for velocity dispersion $\lesssim 100$ km s$^{-1}$. The cluster VDFs also exceed the field at velocity dispersion $\gtrsim 250$ km s$^{-1}$. The differences between cluster and field VDFs are potentially important tests of simulations and of the formation of structure in the universe.

We discuss the phenomenology of gravitational lensing in the purely metric $f\left(\chi\right)$ gravity, an $f(R)$ gravity where the action of the gravitational field depends on the source mass. We focus on the strong lensing regime in galaxy-galaxy lens systems and in clusters of galaxies. Using an approximate metric solution accurate to second order of the velocity field $v/c$, we show how, in the $f\left(\chi\right)=\chi^{3/2}$ gravity, the same light deflection can be produced by point-like lenses with masses smaller than in General Relativity; this mass difference increases with increasing impact parameter and decreasing lens mass. However, for sufficiently massive point-like lenses and small impact parameters, $f\left(\chi\right)=\chi^{3/2}$ and GR yield indistinguishable light deflection angles: this regime occurs both in observed galaxy-galaxy lens systems and in the central regions of galaxy clusters. In the former systems, the GR and $f\left(\chi\right)$ masses are compatible with the mass of standard stellar populations and little or no dark matter, whereas, on the scales of the core of galaxy clusters, the presence of substantial dark matter is required both in General Relativity, and in our approximate $f\left(\chi\right)=\chi^{3/2}$ point-like lens solution. We thus conclude that our approximate metric solution of $f\left(\chi\right)=\chi^{3/2}$ is unable to describe the observed phenomenology of the strong lensing regime without the aid of dark matter.

For the first time, we explore the dynamics of the central region of a galaxy cluster within $r_{500}\sim 600h^{-1}$~kpc from its center by combining optical and X-ray spectroscopy. We use (1) the caustic technique that identifies the cluster substructures and their galaxy members with optical spectroscopic data, and (2) the X-ray redshift fitting procedure that estimates the redshift distribution of the intracluster medium (ICM). We use the spatial and redshift distributions of the galaxies and of the X-ray emitting gas to associate the optical substructures to the X-ray regions. When we apply this approach to Abell 85 (A85), a complex dynamical structure of A85 emerges from our analysis: a galaxy group, with redshift $z=0.0509 \pm 0.0021$ is passing through the cluster center along the line of sight dragging part of the ICM present in the cluster core; two additional groups, at redshift $z=0.0547 \pm 0.0022$ and $z=0.0570 \pm 0.0020$, are going through the cluster in opposite directions, almost perpendicularly to the line of sight, and have substantially perturbed the dynamics of the ICM. An additional group in the outskirts of A85, at redshift $z=0.0561 \pm 0.0023$, is associated to a secondary peak of the X-ray emission, at redshift $z=0.0583^{+0.0039}_{-0.0047}$. Although our analysis and results on A85 need to be confirmed by high-resolution spectroscopy, they demonstrate how our new approach can be a powerful tool to constrain the formation history of galaxy clusters by unveiling their central and surrounding structures.

We compare X-ray and caustic mass profiles for a sample of 16 massive galaxy clusters. We assume hydrostatic equilibrium in interpreting the X-ray data, and use large samples of cluster members with redshifts as a basis for applying the caustic technique. The hydrostatic and caustic masses agree to better than $\approx20\%$ on average across the radial range covered by both techniques $(\sim[0.2-1.25]R_{500})$. The mass profiles were measured independently and do not assume a common functional form. Previous studies suggest that, at $R_{500}$, the hydrostatic and caustic masses are biased low and high respectively. We find that the ratio of hydrostatic to caustic mass at $R_{500}$ is $1.20^{+0.13}_{-0.11}$; thus it is larger than 0.9 at $\approx3\sigma$ and the combination of under- and over-estimation of the mass by these two techniques is $\approx10\%$ at most. There is no indication of any dependence of the mass ratio on the X-ray morphology of the clusters, indicating that the hydrostatic masses are not strongly systematically affected by the dynamical state of the clusters. Overall, our results favour a small value of the so-called hydrostatic bias due to non-thermal pressure sources.

The Andromeda galaxy is observed to have a system of two large dwarf ellipticals and ~13 smaller satellite galaxies that are currently co-rotating in a thin plane, in addition to 2 counter-rotating satellite galaxies. We explored the consistency of those observations with a scenario where the majority of the co-rotating satellite galaxies originated from a subhalo group, where NGC 205 was the host and the satellite galaxies occupied dark matter sub-subhalos. We ran N-body simulations of a close encounter between NGC 205 and M31. In the simulations, NGC 205 was surrounded by massless particles to statistically sample the distribution of the sub-subhalos expected in a subhalo that has a mass similar to NGC 205. We made Monte Carlo samplings and found that, using a set of reference parameters, the probability of producing a thinner distribution of sub-subhalos than the observed NGC 205 + 15 smaller satellites (thus including the 2 counter-rotators, but excluding M32) increased from <1e-8 for the initial distribution to ~0.01 at pericentre. The probability of the simulated sub-subhalos occupying the locations of the observed co-rotating satellites in the line of sight velocity versus projected on-sky distance plane is at most 0.002 for 11 out of 13 satellites. Increasing the mass of M31 and the extent of the initial distribution of sub-subhalos gives a maximum probability of 0.004 for all 13 co-rotating satellites, but the probability of producing the thinness would drop to ~ 0.001.

We present the study of the colour-magnitude diagram of the cluster Abell 2151 (A2151), with a particular focus on the low-mass end. The deep spectroscopy with AF2/WYFFOS@WHT and the caustic method enable us to obtain 360 members within 1.3 R_200 and absolute magnitude M_r < M*_r+6. This nearby cluster shows a well defined red sequence up to M_r ~ - 18.5; at fainter magnitudes only 36% of the galaxies lie on the extrapolation of the red sequence. We compare the red sequences of A2151 and Abell 85, which is another nearby cluster with similar spectroscopic data, but with different mass and dynamical state. Both clusters show similar red sequences at the bright end (M_r < -19.5), whereas large differences appear at the faint end. This result suggests that the reddening of bright galaxies is independent of environment, unlike the dwarf population (M_r > -18.0).

We explore the possibility of measuring the mass accretion rate (MAR) of galaxy clusters from their mass profiles beyond the virial radius $R_{200}$. We derive the accretion rate from the mass of a spherical shell whose inner radius is $2R_{200}$, whose thickness changes with redshift, and whose infall velocity is assumed to be equal to the mean infall velocity of the spherical shells of dark matter halos extracted from $N$-body simulations. This approximation is rather crude in hierarchical clustering scenarios where both smooth accretion and aggregation of smaller dark matter halos contribute to the mass accretion of clusters.Nevertheless, in the redshift range $z=[0,2]$, our prescription returns an average MAR within $20-40 \%$ of the average rate derived from the merger trees of dark matter halos extracted from $N$-body simulations. The MAR of galaxy clusters has been the topic of numerous detailed numerical and theoretical investigations, but so far it has remained inaccessible to measurements in the real universe. Since the measurement of the mass profile of clusters beyond their virial radius can be performed with the caustic technique applied to dense redshift surveys of the cluster outer regions, our result suggests that measuring the mean MAR of a sample of galaxy clusters is actually feasible. We thus provide a new potential observational test of the cosmological and structure formation models.

We investigate the proof of concept and the implications of \textit{refracted gravity}, a novel modified gravity aimed to solve the discrepancy between the luminous and the dynamical mass of cosmic structures without resorting to dark matter. Inspired by the behavior of electric fields in matter, refracted gravity introduces a gravitational permittivity that depends on the local mass density and modifies the standard Poisson equation. The resulting gravitational field can become more intense than the Newtonian field and can mimic the presence of dark matter. We show that the refracted gravitational field correctly describes (1) the rotation curves and the Tully-Fisher relation of disk galaxies; and (2) the observed temperature profile of the X-ray gas of galaxy clusters. According to these promising results, we conclude that refracted gravity deserves further investigation.

We present a new deep spectroscopic catalogue for Abell 85, within 3.0 $\times$ 2.6 Mpc$^2$ and down to $M_{r} \sim M_{r}^* +6$. Using the Visible Multi-Object Spectrograph at the Very Large Telescope (VIMOS@VLT) and the AutoFiber 2 at the William Herschel Telescope (AF2@WHT), we obtained almost 1,430 new redshifts for galaxies with $m_r \leq 21$ mag and $\langle \mu_{e,r} \rangle \leq 24$ mag arcsec$^{-2}$. These redshifts, together with SDSS-DR6 and NED spectroscopic information, result in 460 confirmed cluster members. This dataset allows the study of the luminosity function (LF) of the cluster galaxies covering three orders of magnitudes in luminosities. The total and radial LFs are best modelled by a double Schechter function. The normalized LFs show that their bright ($M_{r} \leq -21.5$) and faint ($M_{r}\geq -18.0$) ends are independent of clustercentric distance and similar to the field LFs unlike the intermediate luminosity range ($-21.5 \leq M_{r} \leq -18.0$). Similar results are found for the LFs of the dominant types of galaxies: red, passive, virialized and early-infall members. On the contrary, the LFs of blue, star forming, non-virialized and recent-infall galaxies are well described by a single Schechter function. These populations contribute to a small fraction of the galaxy density in the innermost cluster region. However, in the outskirts of the cluster, they have similar densities to red, passive, virialized and early-infall members at the LF faint end. These results confirm a clear dependence of the colour and star formation of Abell 85 members in the cluster centric distance.

HectoMAP is a dense redshift survey of red galaxies covering a 53 $deg^{2}$ strip of the northern sky. HectoMAP is 97\% complete for galaxies with $r<20.5$, $(g-r)>1.0$, and $(r-i)>0.5$. The survey enables tests of the physical properties of large-scale structure at intermediate redshift against cosmological models. We use the Horizon Run 4, one of the densest and largest cosmological simulations based on the standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) model, to compare the physical properties of observed large-scale structures with simulated ones in a volume-limited sample covering 8$\times10^6$ $h^{-3}$ Mpc$^3$ in the redshift range $0.22<z<0.44$. We apply the same criteria to the observations and simulations to identify over- and under-dense large-scale features of the galaxy distribution. The richness and size distributions of observed over-dense structures agree well with the simulated ones. Observations and simulations also agree for the volume and size distributions of under-dense structures, voids. The properties of the largest over-dense structure and the largest void in HectoMAP are well within the distributions for the largest structures drawn from 300 Horizon Run 4 mock surveys. Overall the size, richness and volume distributions of observed large-scale structures in the redshift range $0.22<z<0.44$ are remarkably consistent with predictions of the standard $\Lambda$CDM model.

Dynamical analysis of compact groups provides important tests of models of compact group formation and evolution. By compiling 2066 redshifts from FLWO/FAST, from the literature, and from SDSS DR12 in the fields of compact groups in \citet{McC09}, we construct the largest sample of compact groups with complete spectroscopic redshifts in the redshift range $0.01 < z < 0.22$. This large redshift sample shows that the interloper fraction in the \citet{McC09} compact group candidates is $\sim 42\%$. A secure sample of 332 compact groups includes 192 groups with four or more member galaxies and 140 groups with three members. The fraction of early-type galaxies in these compact groups is 62\%, slightly higher than for the original Hickson compact groups. The velocity dispersions of early- and late-type galaxies in compact groups change little with groupcentric radius; the radii sampled are less than $100 ~h^{-1}$ kpc, smaller than the radii typically sampled by members of massive clusters of galaxies. The physical properties of our sample compact groups include size, number density, velocity dispersion, and local environment; these properties slightly differ from those derived for the original Hickson compact groups and for the DPOSS II compact groups. Differences result from subtle differences in the way the group candidates were originally selected. The space density of the compact groups changes little with redshift over the range covered by this sample. The approximate constancy of the space density for this sample is a potential constraint on the evolution of compact groups on a few Gigayear timescale.

With the Westerbork Synthesis Radio Telescope, we performed HI observations of a sample of known X-ray emitting Gigahertz-peaked-spectrum galaxies with compact-symmetric-object morphology (GPS/CSOs) that lacked an HI absorption detection. We combined radio and X-ray data of the full sample of X-ray emitting GPS/CSOs and found a significant, positive correlation between the column densities of the total and neutral hydrogen ($N_{\rm H}$ and $N_{\rm HI}$, respectively). Using a Bayesian approach, we simultaneously quantified the parameters of the $N_{\rm H} - N_{\rm HI}$ relation and the intrinsic spread of the data set. For a specific subset of our sample, we found $N_{\rm H} \propto N_{\rm HI}^b$, with $b=0.93^{+0.49}_{-0.33}$, and $\sigma_{int} (N_{\rm H})= 1.27^{+1.30}_{-0.40}$. The $N_{\rm H} - N_{\rm HI}$ correlation suggests a connection between the physical properties of the radio and X-ray absorbing gas.

We investigate the power of the caustic technique for identifying substructures of galaxy clusters from optical redshift data alone. The caustic technique is designed to estimate the mass profile of galaxy clusters to radii well beyond the virial radius, where dynamical equilibrium does not hold. Two by-products of this technique are the identification of the cluster members and the identification of the cluster substructures. We test the caustic technique as a substructure detector on two samples of 150 mock redshift surveys of clusters; the clusters are extracted from a large cosmological $N$-body simulation of a $\Lambda$CDM model and have masses of $M_{200} \sim 10^{14} h^{-1} M_{\odot}$ and $M_{200} \sim 10^{15} h^{-1} M_{\odot}$ in the two samples. We limit our analysis to substructures identified in the simulation with masses larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within $3R_{200}$, (1) the caustic technique recovers $\sim 30-50$\% of the real substructures, and (2) $\sim 15-20$\% of the substructures identified by the caustic technique correspond to real substructures of the central cluster, the remaining fraction being low-mass substructures, groups or substructures of clusters in the surrounding region, or chance alignments of unrelated galaxies. These encouraging results show that the caustic technique is a promising approach for investigating the complex dynamics of galaxy clusters.

We estimate cluster masses and velocity dispersions for 123 clusters from optical spectroscopy to compare the Sunyaev-Zeldovich (SZ) mass proxy and dynamical masses. Our new survey, HeCS-SZ (Hectospec Cluster Survey of SZ-selected clusters), includes 7,721 new or remeasured redshifts from MMT/Hectospec observations of 24 SZ-selected clusters at redshifts $z$=0.05-0.20 and not in previous surveys. We supplement the Hectospec data with spectra from the Sloan Digital Sky Survey (SDSS) and cluster data from the Cluster Infall Regions in SDSS (CIRS) project and the Hectospec Cluster Survey (HeCS), our Hectospec survey of clusters selected by X-ray flux. We measure the scaling relation between velocity dispersion and SZ mass estimates from the integrated Compton parameter for an SZ complete sample of 83 clusters. The observed relation agrees very well with a simple virial scaling from mass (based on SZ) to velocity dispersion. The SZ mass estimates (calibrated with hydrostatic X-ray mass estimates) are not significantly biased. Further, the velocity dispersion of cluster galaxies is consistent with the expected velocity dispersion of dark matter particles, indicating that galaxies are good dynamical tracers (i.e., velocity bias is small). Significant mass bias in SZ mass estimates could relieve tension between cosmological results from Planck SZ cluster counts and Planck CMB data. However, the excellent agreement between our measured velocity dispersions and those predicted from a virial scaling relation suggests that any SZ mass bias is too small to reconcile SZ and CMB results. In principle, SZ mass bias and velocity bias of galaxies could conspire to yield good agreement, but the required velocity bias is $\sigma_{galaxy}\approx 0.77\sigma_{DM}$, outside the range of plausible models of velocity bias in the literature.

This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk's mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.

The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.

We suggest how we can use the mass profile of galaxy clusters beyond their virial radius to measure their mass accretion rate, a key prediction of structure formation models. The mass profile can be estimated by applying the caustic technique to dense redshift surveys of clusters and their outskirts, where dynamical equilibrium does not necessarily hold. An additional probe of the mass growth of clusters is their mass fraction in substructures. We show that the caustic technique, that identifies cluster substructures as a by-product, returns catalogs of substructures with mass larger than a few $10^{13}h^{-1}M_\odot$ that are between 60% and 80% complete, depending on the density of the redshift survey.

We study the mid-infrared (MIR) properties of the galaxies in the A2199 supercluster at z = 0.03 to understand the star formation activity of galaxy groups and clusters in the supercluster environment. Using the Wide-field Infrared Survey Explorer data, we find no dependence of mass-normalized integrated SFRs of galaxy groups/clusters on their virial masses. We classify the supercluster galaxies into three classes in the MIR color-luminosity diagram: MIR blue cloud (massive, quiescent and mostly early-type), MIR star-forming sequence (mostly late-type), and MIR green valley galaxies. These MIR green valley galaxies are distinguishable from the optical green valley galaxies, in the sense that they belong to the optical red sequence. We find that the fraction of each MIR class does not depend on virial mass of each group/cluster. We compare the cumulative distributions of surface galaxy number density and cluster/group-centric distance for the three MIR classes. MIR green valley galaxies show the distribution between MIR blue cloud and MIR SF sequence galaxies. However, if we fix galaxy morphology, early- and late-type MIR green valley galaxies show different distributions. Our results suggest a possible evolutionary scenario of these galaxies: 1) Late-type MIR SF sequence galaxies -> 2) Late-type MIR green valley galaxies -> 3) Early-type MIR green valley galaxies -> 4) Early-type MIR blue cloud galaxies. In this sequence, star formation of galaxies is quenched before the galaxies enter the MIR green valley, and then morphological transformation occurs in the MIR green valley.

We use dense redshift surveys of nine galaxy clusters at $z\sim0.2$ to compare the galaxy distribution in each system with the projected matter distribution from weak lensing. By combining 2087 new MMT/Hectospec redshifts and the data in the literature, we construct spectroscopic samples within the region of weak-lensing maps of high (70--89%) and uniform completeness. With these dense redshift surveys, we construct galaxy number density maps using several galaxy subsamples. The shape of the main cluster concentration in the weak-lensing maps is similar to the global morphology of the number density maps based on cluster members alone, mainly dominated by red members. We cross correlate the galaxy number density maps with the weak-lensing maps. The cross correlation signal when we include foreground and background galaxies at 0.5$z_{\rm cl}<z<2z_{\rm cl}$ is $10-23$% larger than for cluster members alone at the cluster virial radius. The excess can be as high as 30% depending on the cluster. Cross correlating the galaxy number density and weak-lensing maps suggests that superimposed structures close to the cluster in redshift space contribute more significantly to the excess cross correlation signal than unrelated large-scale structure along the line of sight. Interestingly, the weak-lensing mass profiles are not well constrained for the clusters with the largest cross correlation signal excesses ($>$20% for A383, A689 and A750). The fractional excess in the cross correlation signal including foreground and background structures could be a useful proxy for assessing the reliability of weak-lensing cluster mass estimates.

We present a spectrophotometric analysis of the galaxy pop. in the area of the merging cluster Abell 3921 at redshift 0.093. We investigate the impact of the complex cluster environment on galaxy properties such as morphology or star formation rate. We combine multi-object spectroscopy from the 2dF spectrograph with optical imaging taken with the ESO WFI. We carry out a redshift analysis and determine cluster velocity dispersions using biweight statistics. Applying a Dressler-Shectman (DS-)test we seek evidence for cluster substructure. Cluster and field galaxies are investigated with respect to [OII] and H{\alpha} equivalent width, SFR and morphological descriptors such as concentration index and Gini coefficient. We study these cluster galaxy properties as a function of clustercentric distance and investigate the spatial distribution of various galaxy types. Applying the DS-test we find a 3rd component (A3921-C) in addition to the two main subclusters (A3921-A and A3921-B) already known. The re-determined mass ratio between the main components A and B is approx. 2:1. Similar to prev. studies of galaxy clusters, we find that a large fraction of the disk galaxies close to the cluster core show no detectable star formation. These are likely systems that are quenched due to ram pressure stripping. We also find quenched spirals at rather large distances of 3 to 4 Mpc from the cluster core. A3921-C might be a group of galaxies falling onto the main cluster components. We speculate that the unexpected population of quenched spirals at large clustercentric radii in A3921-A and A3921-B might be an effect of the ongoing cluster merger: shocks in the ICM might give raise to enhanced ram pressure stripping and at least in part be the cause for the quenching of star formation. These quenched spirals might be an interm. stage in the morphological transformation of field spirals into cluster S0s.

We present a new deep determination of the spectroscopic LF within the virial radius of the nearby and massive Abell\,85 (A85) cluster down to the dwarf regime (M* + 6) using VLT/VIMOS spectra for $\sim 2000$ galaxies with m$_r \leq 21$ mag and $\langle \mu_{e,r} \rangle \leq 24$ mag arcsec$^{-2}$. The resulting LF from 438 cluster members is best modelled by a double Schechter function due to the presence of a statistically significant upturn at the faint-end. The amplitude of this upturn ($\alpha_{f} = -1.58^{+0.19}_{-0.15}$), however, is much smaller than that of the SDSS composite photometric cluster LF by Popesso et al. 2006, $\alpha_{f} \sim$ -2. The faint-end slope of the LF in A85 is consistent, within the uncertainties, with that of the field. The red galaxy population dominates the LF at low luminosities, and is the main responsible for the upturn. The fact that the slopes of the spectroscopic LFs in the field and in a cluster as massive as A85 are similar suggests that the cluster environment does not play a major role in determining the abundance of low-mass galaxies.

In this article we continue to test cosmological models centred on Modified Newtonian Dynamics (MOND) with light sterile neutrinos, which could in principle be a way to solve the fine-tuning problems of the standard model on galaxy scales while preserving successful predictions on larger scales. Due to previous failures of the simple MOND cosmological model, here we test a speculative model where the modified gravitational field is produced only by the baryons and the sterile neutrinos produce a purely Newtonian field (hence Modified Baryonic Dynamics). We use two component cosmological simulations to separate the baryonic N-body particles from the sterile neutrino ones. The premise is to attenuate the over-production of massive galaxy cluster halos which were prevalent in the original MOND plus light sterile neutrinos scenario. Theoretical issues with such a formulation notwithstanding, the Modified Baryonic Dynamics model fails to produce the correct amplitude for the galaxy cluster mass function for any reasonable value of the primordial power spectrum normalisation.

The classical dwarf spheroidals (dSphs) provide a critical test for Modified Newtonian Dynamics (MOND) because they are observable satellite galactic systems with low internal accelerations and low, but periodically varying, external acceleration. This varying external gravitational field is not commonly found acting on systems with low internal acceleration. Using Jeans modelling, Carina in particular has been demonstrated to require a V-band mass-to-light ratio greater than 5, which is the nominal upper limit for an ancient stellar population. We run MOND N-body simulations of a Carina-like dSph orbiting the Milky Way to test if dSphs in MOND are stable to tidal forces over the Hubble time and if those same tidal forces artificially inflate their velocity dispersions and therefore their apparent mass-to-light ratio. We run many simulations with various initial total masses for Carina, and Galactocentric orbits (consistent with proper motions), and compare the simulation line of sight velocity dispersions (losVDs) with the observed losVDs of Walker et al. (2007). We find that the dSphs are stable, but that the tidal forces are not conducive to artificially inflating the losVDs. Furthermore, the range of mass-to-light ratios that best reproduces the observed line of sight velocity dispersions of Carina is 5.3 to 5.7 and circular orbits are preferred to plunging orbits. Therefore, some tension still exists between the required mass-to-light ratio for the Carina dSph in MOND and those expected from stellar population synthesis models. It remains to be seen whether a careful treatment of the binary population or triaxiality might reduce this tension.

In this letter we examine the evolution of the radial metallicity gradient induced by secular processes, in the disk of an $N$-body Milky Way-like galaxy. We assign a [Fe/H] value to each particle of the simulation according to an initial, cosmologically motivated, radial chemical distribution and let the disk dynamically evolve for 6 Gyr. This direct approach allows us to take into account only the effects of dynamical evolution and to gauge how and to what extent they affect the initial chemical conditions. The initial [Fe/H] distribution increases with R in the inner disk up to R ~ 10 kpc and decreases for larger R. We find that the initial chemical profile does not undergo major transformations after 6 Gyr of dynamical evolution. The final radial chemical gradients predicted by the model in the solar neighborhood are positive and of the same order of those recently observed in the Milky Way thick disk. We conclude that: 1) the spatial chemical imprint at the time of disk formation is not washed out by secular dynamical processes, and 2) the observed radial gradient may be the dynamical relic of a thick disk originated from a stellar population showing a positive chemical radial gradient in the inner regions.

Abell 383 is a famous rich cluster (z = 0.1887) imaged extensively as a basis for intensive strong and weak lensing studies. Nonetheless there are few spectroscopic observations. We enable dynamical analyses by measuring 2360 new redshifts for galaxies with r$_{petro} \leq 20.5$ and within 50$^\prime$ of the BCG (Brightest Cluster Galaxy: R.A.$_{2000} = 42.014125^\circ$, Decl$_{2000} = -03.529228^\circ$). We apply the caustic technique to identify 275 cluster members within 7$h^{-1}$ Mpc of the hierarchical cluster center. The BCG lies within $-11 \pm 110$ km s$^{-1}$ and 21 $\pm 56 h^{-1}$ kpc of the hierarchical cluster center; the velocity dispersion profile of the BCG appears to be an extension of the velocity dispersion profile based on cluster members. The distribution of cluster members on the sky corresponds impressively with the weak lensing contours of Okabe et al. (2010) especially when the impact of foreground and background structure is included. The values of R$_{200}$ = $1.22\pm 0.01 h^{-1}$ Mpc and M$_{200}$ = $(5.07 \pm 0.09)\times 10^{14} h^{-1}$ M$_\odot$ obtained by application of the caustic technique agree well with recent completely independent lensing measures. The caustic estimate extends direct measurement of the cluster mass profile to a radius of $\sim 5 h^{-1}$ Mpc.

This work is focused on the study of the distribution in the CPPS of passive and star-forming galaxy populations and also, the intense and quiescent star-forming populations for a set of 9 nearby z<0.05 galaxy clusters. Furthermore, we compare the CPPS distribution of the passive galaxy population with the accreted halo population of a set of 28 simulated clusters and the star-forming population with the non-accreted population. We consider various cluster accretion epochs and accretion radii where it is assumed that star formation in galaxies becomes quenched, in order to segregate the accreted population from the non-accreted population. Just applying this segregation in simulations, we get a qualitative agreement between the CPPS distributions of the passive and the accreted populations and also between the star-forming and the non-accreted populations. The uncertainty in cluster centering strongly affects the pronounced cuspy profiles of the projected density and also, it can explain the main difference (i.e. inner slope) between the CPPS distribution of passive and accreted populations. The CPPS density of star-forming galaxies and the intensity of ram-pressure stripping present an opposite trend throughout the CPPS. This implies that ram-pressure stripping significantly contributes to modulate the observed CPPS distribution of star-forming galaxies inside and around cluster virial regions. The significant fraction of star-forming galaxies at the projected center of clusters are mainly those galaxies with low l-o-s velocities and they can be mainly identified as those galaxies with a remaining star formation activity (quiescent star-forming galaxies) inside the physical virial region or, in a lower degree, as galaxy interlopers i.e. outside the physical virial region.(abridged)

We use our Modified Newtonian Dynamics (MOND) cosmological particle-mesh N-body code to investigate the feasibility of structure formation in a framework involving MOND and light sterile neutrinos in the mass range 11 - 300 eV: always assuming that \Omega_{\nu_s}=0.225 for H_o=72 \kms Mpc^{-1}. We run a suite of simulations with variants on the expansion history, cosmological variation of the MOND acceleration constant, different normalisations of the power spectrum of the initial perturbations and interpolating functions. Using various box sizes, but typically with ones of length 256 Mpc/h, we compare our simulated halo mass functions with observed cluster mass functions and show that (i) the sterile neutrino mass must be larger than 30 eV to account for the low mass (M_{200}<10^{14.6} solar masses) clusters of galaxies in MOND and (ii) regardless of sterile neutrino mass or any of the variations we mentioned above, it is not possible to form the correct number of high mass (M_{200}>10^{15.1} solar masses) clusters of galaxies: there is always a considerable over production. This means that the ansatz of considering the weak-field limit of MOND together with a component of light sterile neutrinos to form structure from z ~ 200 fails. If MOND is the correct description of weak-field gravitational dynamics, it could mean that subtle effects of the additional fields in covariant theories of MOND render the ansatz inaccurate, or that the gravity generated by light sterile neutrinos (or by similar hot dark matter particles) is different from that generated by the baryons.

The caustic technique measures the mass of galaxy clusters in both their virial and infall regions and, as a byproduct, yields the list of cluster galaxy members. Here we use 100 galaxy clusters with mass M200>=1E14 Msun/h extracted from a cosmological N-body simulation of a LambdaCDM universe to test the ability of the caustic technique to identify the cluster galaxy members. We identify the true three-dimensional members as the gravitationally bound galaxies. The caustic technique uses the caustic location in the redshift diagram to separate the cluster members from the interlopers. We apply the technique to mock catalogues containing 1000 galaxies in the field of view of 12 Mpc/h on a side at the cluster location. On average, this sample size roughly corresponds to 180 real galaxy members within 3r200, similar to recent redshift surveys of cluster regions. The caustic technique yields a completeness, the fraction of identified true members, fc=0.95 (+- 0.03) within 3r200. The contamination increases from fi=0.020 (+0.046;-0.015) at r200 to fi=0.08 (+0.11;-0.05) at 3r200. No other technique for the identification of the members of a galaxy cluster provides such large completeness and small contamination at these large radii. The caustic technique assumes spherical symmetry and the asphericity of the cluster is responsible for most of the spread of the completeness and the contamination. By applying the technique to an approximately spherical system obtained by stacking the individual clusters, the spreads decrease by at least a factor of two. We finally estimate the cluster mass within 3r200 after removing the interlopers: for individual clusters, the mass estimated with the virial theorem is unbiased and within 30 per cent of the actual mass; this spread decreases to less than 10 per cent for the spherically symmetric stacked cluster.

The infall regions of galaxy clusters represent the largest gravitationally bound structures in a $\Lambda$CDM universe. Measuring cluster mass profiles into the infall regions provides an estimate of the ultimate mass of these haloes. We use the caustic technique to measure cluster mass profiles from galaxy redshifts obtained with the Hectospec Cluster Survey (HeCS), an extensive spectroscopic survey of galaxy clusters with MMT/Hectospec. We survey 58 clusters selected by X-ray flux at 0.1$<$$z$$<$0.3. The survey includes 21,314 unique MMT/Hectospec redshifts for individual galaxies; 10,275 of these galaxies are cluster members. For each cluster we acquired high signal-to-noise spectra for $\sim 200$ cluster members and a comparable number of foreground/background galaxies. The cluster members trace out infall patterns around the clusters. The members define a very narrow red sequence. The velocity dispersions decline with radius; we demonstrate that the determination of the velocity dispersion is insensitive to the inclusion of bluer members (a small fraction of the cluster population). We apply the caustic technique to define membership and estimate the mass profiles to large radii. The ultimate halo mass of clusters (the mass that remains bound in the far future of a $\Lambda$CDM universe) is on average (1.99$\pm$0.11)$M_{200}$, a new observational cosmological test in essential agreement with simulations. Summed profiles binned in $M_{200}$ and in $L_X$ demonstrate that the predicted NFW form of the density profile is a remarkably good representation of the data in agreement with weak lensing results extending to large radius. The concentration of these summed profiles is also consistent with theoretical predictions.

Weak gravitational lensing by galaxy clusters on faint higher redshift galaxies has been traditionally used to study the cluster mass distribution and as a tool to identify clusters as peaks in the shear maps. However, it becomes soon clear that peaks statistics can also be used as a way to constrain the underlying cosmological model due to its dependence on both the cosmic expansion rate and the growth rate of structures. This feature makes peak statistics particularly interesting from the point of view of discriminating between General Relativity and modified gravity. Here we consider a general class of $f(R)$ theories and compute the observable mass function based on the aperture mass statistics. We complement our theoretical analysis with a Fisher matrix forecast of the constraints that an Euclid\,-\,like survey can impose on the $f(R)$ model parameters. We show that peak statistics alone can in principle discriminate between General Relativity and $f(R)$ models and strongly constrain the $f(R)$ parameters that are sensitive to the non-linear growth of structure. However, further analysis is needed in order to include possible selection function in the peaks redshift determination.

The bulge dominated galaxy NGC 7814 provides one of the strongest dynamical tests possible for Modified Newtonian Dynamics (MOND). Spitzer 3.6 micron photometry fixes the bulge parameterisation and strongly constrains the properties of the sub-dominant stellar disk. Furthermore, the distance is known to better than 5 percent, virtually eliminating it as a free parameter. The rotation curve is easily measured, since the H I (and stellar) disks are edge on, and both the receding and approaching sides agree very well. In this paper we explore the agreement between the model and observed rotation curves in MOND given that the only two free parameters available are the mass-to-light ratios of the bulge and disk. We use a grid based MOND Poisson solver that accurately solves for the MOND gravity and produces our model rotation curves from a given mass distribution. The input to the Poisson solver is a 3D distribution of N particles which is generated from modelling the observed distribution of stars and gas in the galaxy. By ensuring a superior fit to the radial surface brightness profile than previous works, by virtue of a double Sersic fit to the bulge, we were able to produce excellent fits to the rotation curve with typical values for both mass-to-light ratios. We conclude that the model rotation curve of a mass distribution in MOND is extremely sensitive to the bulge-disk decomposition and even slight deviation from the observed mass distribution can produce large differences in the model rotation curve.

Cluster mass profiles are tests of models of structure formation. Only two current observational methods of determining the mass profile, gravitational lensing and the caustic technique, are independent of the assumption of dynamical equilibrium. Both techniques enable determination of the extended mass profile at radii beyond the virial radius. For 19 clusters, we compare the mass profile based on the caustic technique with weak lensing measurements taken from the literature. This comparison offers a test of systematic issues in both techniques. Around the virial radius, the two methods of mass estimation agree to within about 30%, consistent with the expected errors in the individual techniques. At small radii, the caustic technique overestimates the mass as expected from numerical simulations. The ratio between the lensing profile and the caustic mass profile at these radii suggests that the weak lensing profiles are a good representation of the true mass profile. At radii larger than the virial radius, the lensing mass profile exceeds the caustic mass profile possibly as a result of contamination of the lensing profile by large-scale structures within the lensing kernel. We highlight the case of the closely neighboring clusters MS0906+11 and A750 to illustrate the potential seriousness of contamination of the the weak lensing signal by unrelated structures.

General Relativity is able to describe the dynamics of galaxies and larger cosmic structures only if most of the matter in the Universe is dark, namely it does not emit any electromagnetic radiation. Intriguingly, on the scale of galaxies, there is strong observational evidence that the presence of dark matter appears to be necessary only when the gravitational field inferred from the distribution of the luminous matter falls below an acceleration of the order of 10^(-10) m/s^2. In the standard model, which combines Newtonian gravity with dark matter, the origin of this acceleration scale is challenging and remains unsolved. On the contrary, the full set of observations can be neatly described, and were partly predicted, by a modification of Newtonian dynamics, dubbed MOND, that does not resort to the existence of dark matter. On the scale of galaxy clusters and beyond, however, MOND is not as successful as on the scale of galaxies, and the existence of some dark matter appears unavoidable. A model combining MOND with hot dark matter made of sterile neutrinos seems to be able to describe most of the astrophysical phenomenology, from the power spectrum of the cosmic microwave background anisotropies to the dynamics of dwarf galaxies. Whether there exists a yet unknown covariant theory that contains General Relativity and Newtonian gravity in the weak field limit, and MOND as the ultra-weak field limit is still an open question.

We use Wide-field Infrared Survey Explorer (WISE) data covering the entire region (~130 deg^2) of the A2199 supercluster at z=0.03 to study the mid-infrared (MIR) properties of supercluster galaxies. We identify a `MIR star-forming sequence' in the WISE [3.4]-[12] color-12 \mu m luminosity diagram, consisting of late-type, star-forming galaxies. At a fixed star formation rate (SFR), the MIR-detected galaxies at 22 \mu m or 12 \mu m tend to be more metal rich and to have higher surface brightness than those without MIR detection. Using these MIR-detected galaxies, we construct the IR luminosity function (LF) and investigate its environmental dependence. Both total IR (TIR) and 12 \mu m LFs are dominated by late-type, star-forming galaxies. The contribution of active galactic nuclei (AGN)-host galaxies increases with both TIR and 12 \mu m luminosities. The contribution of early-type galaxies to the 12 \mu m LFs increases with decreasing luminosity. The faint-end slope of the TIR LFs does not change with environment, but the change of faint-end slope in the 12 \mu m LFs with the environment is significant: there is a steeper faint-end slope in the cluster core than in the cluster outskirts. This steepening results primarily from the increasing contribution of early-type galaxies toward the cluster. These galaxies are passively evolving, and contain old stellar populations with weak MIR emission from the circumstellar dust around asymptotic giant branch stars.

SHELS (Smithsonian Hectospec Lensing Survey) is a dense redshift survey covering a 4 square degree region to a limiting R = 20.6. In the construction of the galaxy catalog and in the acquisition of spectroscopic targets, we paid careful attention to the survey completeness for lower surface brightness dwarf galaxies. Thus, although the survey covers a small area, it is a robust basis for computation of the slope of the faint end of the galaxy luminosity function to a limiting M_R = -13.3 + 5logh. We calculate the faint end slope in the R-band for the subset of SHELS galaxies with redshif ts in the range 0.02 <= z < 0.1, SHELS_{0.1}. This sample contains 532 galaxies with R< 20.6 and with a median surface brightness within the half light radius of SB_{50,R} = 21.82 mag arcsec^{-2}. We used this sample to make one of the few direct measurements of the dependence of the faint end of the galaxy luminosity function on surface brightness. For the sample as a whole the faint end slope, alpha = -1.31 +/- 0.04, is consistent with both the Blanton et al. (2005b) analysis of the SDSS and the Liu et al. (2008) analysis of the COSMOS field. This consistency is impressive given the very different approaches of th ese three surveys. A magnitude limited sample of 135 galaxies with optical spectroscopic reds hifts with mean half-light surface brightness, SB_{50,R} >= 22.5 mag arcsec^{-2} is unique to SHELS_{0.1}. The faint end slope is alpha_{22.5} = -1.52+/- 0.16. SHELS_{0.1} shows that lower surface brightness objects dominate the faint end slope of the l uminosity function in the field, underscoring the importance of surface brightness limits in evaluating measurements of the faint end slope and its evolution.

The observed accelerated cosmic expansion can be a signature of fourth\,-\,order gravity theories, where the acceleration of the Universe is a consequence of departures from Einstein General Relativity, rather than the sign of the existence of a fluid with negative pressure. In the fourth\,-\,order gravity theories, the gravity Lagrangian is described by an analytic function $f(R)$ of the scalar curvature $R$ subject to the demanding conditions that no detectable deviations from standard GR is observed on the Solar System scale. Here we consider two classes of $f(R)$ theories able to pass Solar System tests and investigate their viability on cosmological scales. To this end, we fit the theories to a large dataset including the combined Hubble diagram of Type Ia Supernovae and Gamma Ray Bursts, the Hubble parameter $H(z)$ data from passively evolving red galaxies, Baryon Acoustic Oscillations extracted from the seventh data release of the Sloan Digital Sky Survey (SDSS) and the distance priors from the Wilkinson Microwave Anisotropy Probe seven years (WMAP7) data. We find that both classes of $f(R)$ fit very well this large dataset with the present\,-\,day values of the matter density, Hubble constant and deceleration parameter in agreement with previous estimates; however, the strong degeneracy among the $f(R)$ parameters prevents us from strongly constraining their values. We also derive the growth factor $g = d\ln{\delta}/d\ln{a}$, with $\delta = \delta \rho_M/\rho_M$ the matter density perturbation, and show that it can still be well approximated by $g(z) \propto \Omega_M(z)^{\gamma}$. We finally constrain $\gamma$ (on some representative scales) and investigate its redshift dependence to see whether future data can discriminate between these classes of $f(R)$ theories and standard dark energy models.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

LambdaCDM, for the currently preferred cosmological density Omega_0 and cosmological constant Omega_Lambda, predicts that the Universe expansion decelerates from early times to redshift z~0.9 and accelerates at later times. On the contrary, the cosmological model based on conformal gravity predicts that the cosmic expansion has always been accelerating. To distinguish between these two very different cosmologies, we resort to gamma-ray bursts (GRBs), which have been suggested to probe the Universe expansion history at z>1, where identified type Ia supernovae (SNe) are rare. We use the full Bayesian approach to infer the cosmological parameters and the additional parameters required to describe the GRB data available in the literature. For the first time, we use GRBs as cosmological probes without any prior information from other data. In addition, when we combine the GRB samples with SNe, our approach neatly avoids all the inconsistencies of most numerous previous methods that are plagued by the so-called circularity problem. In fact, when analyzed properly, current data are consistent with distance moduli of GRBs and SNe that can respectively be, in a variant of conformal gravity, ~15 and ~3 magnitudes fainter than in LambdaCDM. Our results indicate that the currently available SN and GRB samples are accommodated equally well by both LambdaCDM and conformal gravity and do not exclude a continuous accelerated expansion. We conclude that GRBs are currently far from being effective cosmological probes, as they are unable to distinguish between these two very different expansion histories.

Redshift surveys are a powerful tool of modern cosmology. We discuss two aspects of their power to map the distribution of mass and light in the universe: (1) measuring the mass distribution extending into the infall regions of rich clusters and (2) applying deep redshift surveys to the selection of clusters of galaxies and to the identification of very large structures (Great Walls). We preview the HectoMAP project, a redshift survey with median redshift z = 0.34 covering 50 square degrees to r= 21. We emphasize the importance and power of spectroscopy for exploring and understanding the nature and evolution of structure in the universe.

The highly organised distribution of satellite galaxies surrounding the Milky Way is a serious challenge to the concordance cosmological model. Perhaps the only remaining solution, in this framework, is that the dwarf satellite galaxies fall into the Milky Way's potential along one or two filaments, which may or may not plausibly reproduce the observed distribution. Here we test this scenario by making use of the proper motions of the Fornax, Sculptor, Ursa Minor and Carina dwarf spheroidals, and trace their orbits back through several variations of the Milky Way's potential and account for dynamical friction. The key parameters are the proper motions and total masses of the dwarf galaxies. Using a simple model we find no tenable set of parameters that can allow Fornax to be consistent with filamentary infall, mainly because the 1 sigma error on its proper motion is relatively small. The other three must walk a tightrope between requiring a small pericentre (less than 20 kpc) to lose enough orbital energy to dynamical friction and avoiding being tidally disrupted. We then employed a more realistic model with host halo mass accretion, and found that the four dwarf galaxies must have fallen in at least 5 Gyrs ago. This time interval is longer than organised distribution is expected to last before being erased by the randomisation of the satellite orbits.

We present the cosmic shear signal predicted by two viable cosmological models in the framework of modified-action f(R) theories. We use f(R) models where the current accelerated expansion of the Universe is a direct consequence of the modified gravitational Lagrangian rather than Dark Energy (DE), either in the form of vacuum energy/cosmological constant or of a dynamical scalar field (e.g. quintessence). We choose Starobinsky's (St) and Hu & Sawicki's (HS) f(R) models, which are carefully designed to pass the Solar System gravity tests. In order to further support - or rule out - f(R) theories as alternative candidates to the DE hypothesis, we exploit the power of weak gravitational lensing, specifically of cosmic shear. We calculate the tomographic shear matrix as it would be measured by the upcoming ESA Cosmic Vision Euclid satellite. We find that in the St model the cosmic shear signal is almost completely degenerate with LCDM, but it is easily distinguishable in the HS model. Moreover, we compute the corresponding Fisher matrix for both the St and HS models, thus obtaining forecasts for their cosmological parameters. Finally, we show that the Bayes factor for cosmic shear will definitely favour the HS model over LCDM if Euclid measures a value larger than ~0.02 for the extra HS parameter n_HS.

We use extreme value statistics to assess the significance of two of the most dramatic structures in the local Universe: the Shapley supercluster and the Sloan Great Wall. If we assume that Shapley (volume ~ 1.2 x 10^5 (Mpc/h)^3) evolved from an overdense region in the initial Gaussian fluctuation field, with currently popular choices for the background cosmological model and the shape and amplitude sigma8 of the initial power spectrum, we estimate that the total mass of the system is within 20 percent of 1.8 x 10^16 Msun/h. Extreme value statistics show that the existence of this massive concentration is not unexpected if the initial fluctuation field was Gaussian, provided there are no other similar objects within a sphere of radius 200 Mpc/h centred on our Galaxy. However, a similar analysis of the Sloan Great Wall, a more distant (z ~ 0.08) and extended concentration of structures (volume ~ 7.2 x 10^5 (Mpc/h)^3) suggests that it is more unusual. We estimate its total mass to be within 20 percent of 1.2 x 10^17 Msun/h; even if it is the densest such object of its volume within z=0.2, its existence is difficult to reconcile with Gaussian initial conditions if sigma8 < 0.9. This tension can be alleviated if this structure is the densest within the Hubble volume. Finally, we show how extreme value statistics can be used to address the likelihood that an object like Shapley exists in the same volume which contains the Great Wall, finding, again, that Shapley is not particularly unusual. It is straightforward to incorporate other models of the initial fluctuation field into our formalism.

We present a new Particle-Mesh cosmological N-body code for accurately solving the modified Poisson equation of the Quasi Linear formulation of MOND. We generate initial conditions for the Angus (2009) cosmological model, which is identical to LCDM except that the cold dark matter is switched for a single species of thermal sterile neutrinos. We set the initial conditions at z=250 for a (512 Mpc/h)^3 box with 256^3 particles and we evolve them down to z=0. We clearly demonstrate the necessity of MOND for developing the large scale structure in a hot dark matter cosmology and contradict the naive expectation that MOND cannot form galaxy clusters. We find that the correct order of magnitude of X-ray clusters (with T_X > 4.5 keV) can be formed, but that we overpredict the number of very rich clusters and seriously underpredict the number of lower mass clusters. The latter is a shortcoming of the resolution of our simulations, whereas we suggest that the over production of very rich clusters might be prevented by incorporating a MOND acceleration constant that varies with redshift and an expansion history that cannot be described by the usual Friedmann models. We present evidence that suggests the density profiles of our simulated clusters are compatible with those of observed X-ray clusters in MOND. It remains to be seen if the low mass end of the cluster mass function can be reproduced and if the high densities of dark matter in the central 20 kpc of groups and clusters of galaxies, measured in the MOND framework, can be achieved. As a last test, we computed the relative velocity between pairs of halos within 10 Mpc and find that pairs with velocities larger than 3000 km/s like the bullet cluster, can form without difficulty.

Flexion is the weak lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. The flexion signal-to-noise ratio can be an order of magnitude larger than that of shear. For the first time, we show how this makes flexion an invaluable tool for discriminating among alternative cosmological models. We analyse a scalar field model of unified dark matter and dark energy, a brane-world cosmology and two f(R) modified-action theories. We show that these models can be distinguished from LCDM at several standard deviations by measuring the power spectrum of cosmic flexion.

We present parameter estimation forecasts for future 3D cosmic shear surveys for a class of Unified Dark Matter (UDM) models, where a single scalar field mimics both Dark Matter (DM) and Dark Energy (DE). These models have the advantage that they can describe the dynamics of the Universe with a single matter component providing an explanation for structure formation and cosmic acceleration. A crucial feature of the class of UDM models we use in this work is characterised by a parameter, c_inf (in units of the speed of light c=1), that is the value of the sound speed at late times, and on which structure formation depends. We demonstrate that the properties of the DM-like behaviour of the scalar field can be estimated with very high precision with large-scale, fully 3D weak lensing surveys. We found that 3D weak lensing significantly constrains c_inf, and we find minimal errors 0.00003, for the fiducial value c_inf=0.001, and 0.000026, for c_inf=0.012. Moreover, we compute the Bayesian evidence for UDM models over the LCDM model as a function of c_inf. For this purpose, we can consider the LCDM model as a UDM model with c_inf=0. We find that the expected evidence clearly shows that the survey data would unquestionably favour UDM models over the LCDM model, for the values c_inf>0.001.

A phenomenological extension of the well-known brane-world cosmology of Dvali, Gabadadze and Porrati (eDGP) has recently been proposed. In this model, a cosmological-constant-like term is explicitly present as a non-vanishing tension sigma on the brane, and an extra parameter alpha tunes the cross-over scale r_c, the scale at which higher dimensional gravity effects become non negligible. Since the Hubble parameter in this cosmology reproduces the same LCDM expansion history, we study how upcoming weak lensing surveys, such as Euclid and DES (Dark Energy Survey), can confirm or rule out this class of models. We perform Markov Chain Monte Carlo simulations to determine the parameters of the model, using Type Ia Supernov\ae, H(z) data, Gamma Ray Bursts and Baryon Acoustic Oscillations. We also fit the power spectrum of the temperature anisotropies of the Cosmic Microwave Background to obtain the correct normalisation for the density perturbation power spectrum. Then, we compute the matter and the cosmic shear power spectra, both in the linear and non-linear regimes. The latter is calculated with the two different approaches of Hu and Sawicki (2007) (HS) and Khoury and Wyman (2009) (KW). With the eDGP parameters coming from the Markov Chains, KW reproduces the LCDM matter power spectrum at both linear and non-linear scales and the LCDM and eDGP shear signals are degenerate. This result does not hold with the HS prescription: Euclid can distinguish the eDGP model from LCDM because their expected power spectra roughly differ by the 3sigma uncertainty in the angular scale range 700<l<3000; on the contrary, the two models differ at most by the 1sigma uncertainty over the range 500<l<3000 in the DES experiment and they are virtually indistinguishable.

Modified Newtonian Dynamics (MOND) has been shown to be able to fit spiral galaxy rotation curves as well as giving a theoretical foundation for empirically determined scaling relations, such as the Tully - Fisher law, without the need for a dark matter halo. As a complementary analysis, one should investigate whether MOND can also reproduce the dynamics of early - type galaxies (ETGs) without dark matter. As a first step, we here show that MOND can indeed fit the observed central velocity dispersion $\sigma_0$ of a large sample of ETGs assuming a simple MOND interpolating functions and constant anisotropy. We also show that, under some assumptions on the luminosity dependence of the Sersic n parameter and the stellar M/L ratio, MOND predicts a fundamental plane for ETGs : a log - linear relation among the effective radius $R_{eff}$, $\sigma_0$ and the mean effective intensity $\langle I_e \rangle$. However, we predict a tilt between the observed and the MOND fundamental planes.

The caustic technique uses galaxy redshifts alone to measure the escape velocity and mass profiles of galaxy clusters to clustrocentric distances well beyond the virial radius, where dynamical equilibrium does not necessarily hold. We provide a detailed description of this technique and analyse its possible systematic errors. We apply the caustic technique to clusters with mass M_200>=10^{14}h^{-1} M_sun extracted from a cosmological hydrodynamic simulation of a LambdaCDM universe. With a few tens of redshifts per squared comoving megaparsec within the cluster, the caustic technique, on average, recovers the profile of the escape velocity from the cluster with better than 10 percent accuracy up to r~4 r_200. The caustic technique also recovers the mass profile with better than 10 percent accuracy in the range (0.6-4) r_200, but it overestimates the mass up to 70 percent at smaller radii. This overestimate is a consequence of neglecting the radial dependence of the filling function F_beta(r). The 1-sigma uncertainty on individual escape velocity profiles increases from ~20 to ~50 percent when the radius increases from r~0.1 r_200 to ~4 r_200. Individual mass profiles have 1-sigma uncertainty between 40 and 80 percent within the radial range (0.6-4) r_200. We show that the amplitude of these uncertainties is completely due to the assumption of spherical symmetry, which is difficult to drop. Alternatively, we can apply the technique to synthetic clusters obtained by stacking individual clusters: in this case, the 1-sigma uncertainty on the escape velocity profile is smaller than 20 percent out to 4 r_200. The caustic technique thus provides reliable average profiles which extend to regions difficult or impossible to probe with other techniques.

Using the caustic method, we identify the member stars of five dwarf spheroidal (dSph) galaxies of the Milky Way, the smallest dark matter (DM) dominated systems in the Universe. After our interloper rejection, we compute line-of-sight velocity dispersion profiles that are substantially smoother than previous results. Moreover, two dSphs have line-of-sight velocity dispersions 20% smaller than previous calculations suggested. Our Jeans modelling confirms that the DM content interior to 300pc is roughly constant with satellite luminosity. Finally, if we assume that MOND provides the true law of gravity, our identification of interlopers implies that four dSphs have mass-to-light ratios in agreement with stellar population synthesis models, whereas Carina still has a mass-to-light ratio a factor of two too large and remains a problem for MOND.

We present the first comparison of virial masses of galaxy clusters with their Sunyaev-Zel'dovich Effect (SZE) signals. We study 15 clusters from the Hectospec Cluster Survey (HeCS) with MMT/Hectospec spectroscopy and published SZE signals. We measure virial masses of these clusters from an average of 90 member redshifts inside the radius $r_{100}$. The virial masses of the clusters are strongly correlated with their SZE signals (at the 99% confidence level using a Spearman rank-sum test). This correlation suggests that $Y_{SZ}$ can be used as a measure of virial mass. Simulations predict a powerlaw scaling of $Y_{SZ}\propto M_{200}^\alpha$ with $\alpha\approx$1.6. Observationally, we find $\alpha$=1.11$\pm$0.16, significantly shallower (given the formal uncertainty) than the theoretical prediction. However, the selection function of our sample is unknown and a bias against less massive clusters cannot be excluded (such a selection bias could artificially flatten the slope). Moreover, our sample indicates that the relation between velocity dispersion (or virial mass estimate) and SZE signal has significant intrinsic scatter, comparable to the range of our current sample. More detailed studies of scaling relations are therefore needed to derive a robust determination of the relation between cluster mass and SZE.

In a dynamical-radiative model we recently developed to describe the physics of compact, GHz-Peaked-Spectrum (GPS) sources, the relativistic jets propagate across the inner, kpc-sized region of the host galaxy, while the electron population of the expanding lobes evolves and emits synchrotron and inverse-Compton (IC) radiation. Interstellar-medium gas clouds engulfed by the expanding lobes, and photoionized by the active nucleus, are responsible for the radio spectral turnover through free-free absorption (FFA) of the synchrotron photons. The model provides a description of the evolution of the spectral energy distribution (SED) of GPS sources with their expansion, predicting significant and complex high-energy emission, from the X-ray to the gamma-ray frequency domain. Here, we test this model with the broad-band SEDs of a sample of eleven X-ray emitting GPS galaxies with Compact-Symmetric-Object (CSO) morphology, and show that: (i) the shape of the radio continuum at frequencies lower than the spectral turnover is indeed well accounted for by the FFA mechanism; (ii) the observed X-ray spectra can be interpreted as non-thermal radiation produced via IC scattering of the local radiation fields off the lobe particles, providing a viable alternative to the thermal, accretion-disk dominated scenario. We also show that the relation between the hydrogen column densities derived from the X-ray (N_H) and radio (N_HI) data of the sources is suggestive of a positive correlation, which, if confirmed by future observations, would provide further support to our scenario of high-energy emitting lobes.

We present a new multi-phase sub-resolution model for star formation and feedback in SPH numerical simulations of galaxy formation. Our model, called MUPPI (MUlti-Phase Particle Integrator), describes each gas particle as a multi-phase system, with cold and hot gas phases, coexisting in pressure equilibrium, and a stellar component. Cooling of the hot tenuous gas phase feeds the cold gas phase. Stars are formed out of molecular gas with a given efficiency, which scales with the dynamical time of the cold phase. Our prescription for star formation is not based on imposing the Schmidt-Kennicutt relation, which is instead naturally produced by MUPPI. Energy from supernova explosions is deposited partly into the hot phase of the gas particles, and partly to that of neighboring particles. Mass and energy flows among the different phases of each particle are described by a set of ordinary differential equations which we explicitly integrate for each gas particle, instead of relying on equilibrium solutions. This system of equations also includes the response of the multi-phase structure to energy changes associated to the thermodynamics of the gas. We apply our model to two isolated disk galaxy simulations and two spherical cooling flows. MUPPI is able to reproduce the Schmidt-Kennicutt relation for disc galaxies. It also reproduces the basic properties of the inter-stellar medium in disc galaxies, the surface densities of cold and molecular gas, of stars and of star formation rate, the vertical velocity dispersion of cold clouds and the flows connected to the galactic fountains. Quite remarkably, MUPPI also provides efficient stellar feedback without the need to include a scheme of kinetic energy feedback. [abridged]

We use the Fifth Data Release of the Sloan Digital Sky Survey to study X-ray-selected galaxy groups and compare their properties to clusters. We search for infall patterns around the groups and use these to measure group mass profiles to large radii. In previous work, we analyzed infall patterns for an X-ray-selected sample of 72 clusters from the ROSAT All-Sky Survey. Here, we extend this approach to a sample of systems with smaller X-ray fluxes selected from the 400 Square Degree serendipitous survey of clusters and groups in ROSAT pointed observations. We identify 16 groups with SDSS DR5 spectroscopy, search for infall patterns, and compute mass profiles out to 2-6 Mpc from the group centers with the caustic technique. No other mass estimation methods are currently available at such large radii for these low-mass groups, because the virial estimate requires dynamical equilibrium and the gravitational lensing signal is too weak. Despite the small masses of these groups, most display recognizable infall patterns. We use caustic and virial mass estimates to measure the scaling relations between different observables, extending these relations to smaller fluxes and luminosities than many previous surveys. Close inspection reveals that three of the groups are subclusters in the outskirts of larger clusters. A fourth group is apparently undergoing a group-group merger. These four merging groups represent the most extreme outliers in the scaling relations. Excluding these groups, we find $L_X\propto\sigma_p^{3.4\pm1.6}$, consistent with previous determinations for both clusters and groups. Understanding cluster and group scaling relations is crucial for measuring cosmological parameters from clusters.

We present a spectroscopic sample of 910 distant halo stars from the Hypervelocity Star survey from which we derive the velocity dispersion profile of the Milky Way halo. The sample is a mix of 74% evolved horizontal branch stars and 26% blue stragglers. We estimate distances to the stars using observed colors, metallicities, and stellar evolution tracks. Our sample contains twice as many objects with R>50 kpc as previous surveys. We compute the velocity dispersion profile in two ways: with a parametric method based on a Milky Way potential model, and with a non-parametric method based on the caustic technique originally developed to measure galaxy cluster mass profiles. The resulting velocity dispersion profiles are remarkably consistent with those found by two independent surveys based on other stellar populations: the Milky Way halo exhibits a mean decline in radial velocity dispersion of -0.38+-0.12 km/s/kpc over 15<R<75 kpc. This measurement is a useful basis for calculating the total mass and mass distribution of the Milky Way halo.

(Abridged) Modified Newtonian Dynamics (MOND) can fit a broad range of galaxy kinematic data, but struggles with clusters of galaxies. MONDian clusters need dark matter, and here we test the 11 eV sterile neutrino - used to fit the first three acoustic peaks of the cosmic microwave background - by investigating their equilibrium distributions in 30 groups and clusters over a wide range of temperatures. Knowing both the sterile neutrino density and velocity dispersion tells us the Tremaine-Gunn phase-space limit at all radii. We find that all 30 systems serendipitously reach the Tremaine-Gunn limit by the centre, which means a portion of the dynamical mass must always be covered by the brightest cluster galaxy. Interestingly, the typical fitted K-band mass-to-light ratio is unity and at most 1.2, which is very consistent - although leaving no margin for error - with stellar population synthesis models. Amidst the sample there are several special cases including the Coma cluster (for which dark matter was first proposed), NGC 720 (where geometrical evidence for dark matter was found) and the bullet cluster (where dark matter - of some kind - in clusters was directly proven to exist). We demonstrate that 11 eV sterile neutrinos are unlikely to influence spiral galaxy rotation curves, as they don't influence even some very massive early-types (NGC 4125 and NGC 6482). Finally, we conclude that it is intriguing that the minimum mass of sterile neutrino particle that can match the cosmic microwave background is identical to the minimum mass found here to be consistent with equilibrium configurations of MONDian clusters of galaxies.

Fourth - order gravity theories have received much interest in recent years thanks to their ability to provide an accelerated cosmic expansion in a matter only universe. In these theories, the Lagrangian density of the gravitational field has the form R + f(R), and the explicit choice of the arbitrary function f(R) must meet the local tests of gravity and the constraints from the primordial abundance of the light elements. Two popular classes of f(R) models, which are expected to fulfill all the above requirements, have recently been proposed. However, neither of these models has ever been quantitatively tested against the available astrophysical data. Here, by combining Type Ia Supernovae and Gamma Ray Bursts, we investigate the ability of these models to reproduce the observed Hubble diagram over the redshift range (0, 7). We find that both models fit very well this dataset with the present day values of the matter density and deceleration parameters which agree with previous estimates. However, the strong degeneracy among the f(R) parameters prevents us from putting strong constraints on the values of these parameters; nevertheless, we can identify the regions of the parameter space that should, in principle, be carefully explored with future data and dynamical probes in order to discriminate among f(R) theories and standard dark energy models.

Weak gravitational lensing is a powerful tool for studying both the geometry and the dynamics of the Universe. Its power spectrum contains information on the sources emitting photons and on the large--scale structures that these hotons cross. We calculate the weak lensing cosmic convergence and shear power spectra, in linear theory and Limber's approximation, for two different classes of cosmological models: the standard \LambdaCDM and Unified Dark Matter (UDM) models. The latter models attempt to solve the problems of the dark matter in the dynamics of galaxies and galaxy clusters and of the late-time acceleration of the Universe expansion by introducing a scalar field that mimics both dark matter and dark energy. A crucial feature of the UDM models is the speed of sound c_\infty, that is the value of the sound speed at late times, on which structure formation depends. In this paper, we provide the predictions of the UDM models for the weak lensing signal,with various values of c_\infty. We consider both the Cosmic Microwave Background and background galaxies at different redshifts as sources for the weak lensing power spectra. We find that the power spectra in UDM models are more sensitive to the variations of c_\infty for sources located at low redshifts. Moreover, we find that for c_\infty>10^{-3} (in units of the speed of light), the UDM weak lensing convergence power spectrum C^{\kappa\kappa}(l) for background galaxies is strongly suppressed with respect to the \LambdaCDM spectrum, particularly at small angular scales l\gtrsim100.

We re-investigate the old problem of the survival of the five globular clusters orbiting the Fornax dwarf galaxy in both standard and modified Newtonian dynamics. For the first time in the history of the topic, we use accurate mass models for the Fornax dwarf, obtained through Jeans modelling of the recently published line of sight velocity dispersion data, and we are also not resigned to circular orbits for the globular clusters. Previously conceived problems stem from fixing the starting distances of the globulars to be less than half the tidal radius. We relax this constraint since there is absolutely no evidence for it and show that the dark matter paradigm, with either cusped or cored dark matter profiles, has no trouble sustaining the orbits of the two least massive globular clusters for a Hubble time almost regardless of their initial distance from Fornax. The three most massive globulars can remain in orbit as long as their starting distances are marginally outside the tidal radius. The outlook for modified Newtonian dynamics is also not nearly as bleak as previously reported. Although dynamical friction inside the tidal radius is far stronger in MOND, outside dynamical friction is negligible due to the absence of stars. This allows highly radial orbits to survive, but more importantly circular orbits at distances more than 85% of Fornax's tidal radius to survive indefinitely. The probability of the globular clusters being on circular orbits at this distance compared with their current projected distances is discussed and shown to be plausible. Finally, if we ignore the presence of the most massive globular (giving it a large line of sight distance) we demonstrate that the remaining four globulars can survive within the tidal radius for the Hubble time with perfectly sensible orbits.

Traditional estimators of the mass of galaxy clusters assume that the cluster components (galaxies, intracluster medium, and dark matter) are in dynamical equilibrium. Two additional estimators, that do not require this assumption, were proposed in the 1990s: gravitational lensing and the caustic technique. With these methods, we can measure the cluster mass within radii much larger than the virial radius. In the caustic technique, the mass measurement is only based on the celestial coordinates and redshifts of the galaxies in the cluster field of view; therefore, unlike lensing, it can be, in principle, applied to clusters at any redshift. Here, we review the origin, the basics and the performance of the caustic method.

The study of the broad-band emission of GHz-Peaked-Spectrum (GPS) radio galaxies is a powerful tool to investigate the physical processes taking place in the central, kpc-sized region of their active hosts, where the jets propagate and the lobes expand, interacting with the surrounding interstellar medium (ISM). We recently developed a new dynamical-radiative model to describe the evolution of the GPS phenomenon (Stawarz et al. 2008): as the relativistic jets propagate through the ISM, gradually engulfing narrow-line emitting gas clouds along their way, the electron population of the expanding lobes evolves, emitting synchrotron light, as well as inverse-Compton radiation via up-scattering of the photon fields from the host galaxy and its active nucleus. The model, which successfully reproduces the key features of the GPS radio sources as a class, provides a description of the evolution of their spectral energy distribution (SED) with the lobes' expansion, predicting significant and complex X-ray to gamma-ray emission. We apply here the model to the broad-band SED's of a sample of known, X-ray emitting GPS galaxies, and show that: (i) the free-free absorption mechanism enables us to reproduce the radio continuum at frequencies below the turnover; (ii) the lobes' non-thermal, inverse-Compton emission can account for the observed X-ray spectra, providing a viable alternative to the thermal, accretion-dominated scenario. We also show that, in our sample, the relationship between the X-ray and radio hydrogen column densitities, N_H and N_HI, is suggestive of a positive correlation, which, if confirmed, would support the scenario of high-energy emitting lobes.

We run adiabatic N-body/hydrodynamical simulations of isolated self-gravitating gas clouds to test whether conformal gravity, an alternative theory to General Relativity, is able to explain the properties of X-ray galaxy clusters without resorting to dark matter. We show that the gas clouds rapidly reach equilibrium with a density profile which is well fit by a beta-model whose normalization and slope are in approximate agreement with observations. However, conformal gravity fails to yield the observed thermal properties of the gas cloud: (i) the mean temperature is at least an order of magnitude larger than observed; (ii) the temperature profiles increase with the square of the distance from the cluster center, in clear disagreement with real X-ray clusters. These results depend on a gravitational potential whose parameters reproduce the velocity rotation curves of spiral galaxies. However, this parametrization stands on an arbitrarily chosen conformal factor. It remains to be seen whether a different conformal factor, specified by a spontaneous breaking of the conformal symmetry, can reconcile this theory with observations.

We use scalar-field Lagrangians with a non-canonical kinetic term to obtain unified dark matter models where both the dark matter and the dark energy, the latter mimicking a cosmological constant, are described by the scalar field itself. In this framework, we propose a technique to reconstruct models where the effective speed of sound is small enough that the scalar field can cluster. These models avoid the strong time evolution of the gravitational potential and the large Integrated Sachs-Wolfe effect which have been a serious drawback of previously considered models. Moreover, these unified dark matter scalar field models can be easily generalized to behave as dark matter plus a dark energy component behaving like any type of quintessence fluid.

The recently revised cosmological constraints from the Five-Year WMAP data ameliorate previous tension between cosmological constraints from the microwave background and from cluster abundances. We demonstrate that the revised estimates of cosmological parameters are in excellent agreement with the mass function of X-ray clusters in the Sloan Digital Sky Survey. Velocity segregation between galaxies and the underlying dark matter could cause virial mass estimates to be biased, causing the mass scale of the mass function to be offset from the true value. Modest velocity segregation ($\sigma_{gxy}/\sigma_{DM}$=1.13$^{+0.06}_{-0.05}$) is sufficient to match the mass function to the Five-Year WMAP results. When the new WMAP results are combined with constraints from supernovae and baryon acoustic oscillations, there is no need for velocity segregation ($\sigma_{gxy}/\sigma_{DM}$=1.05$\pm$0.05). This result agrees with expectations for velocity segregation from state-of-the-art numerical simulations of clusters. Together with the improved agreement between the new WMAP results and recent cosmic shear measurements, this result demonstrates that the amplitude of large-scale structure in the nearby universe matches that predicted from the structure seen in the microwave background. The new constraint we place on velocity segregation in clusters indicates that virial mass estimates for clusters are reasonably accurate. This result suggests that future cluster surveys will be able to probe both cosmological parameters and fundamental cluster physics.

Under the assumption that General Relativity (GR) correctly describes the phenomenology of our Universe, astronomical observations provide compelling evidence that (1) the dynamics of cosmic structure is dominated by dark matter (DM), an exotic matter mostly made of hypothetical elementary particles, and (2) the expansion of the Universe is currently accelerating because of the presence of a positive cosmological constant Lambda. The DM particles have not yet been detected and there is no theoretical justification for the tiny positive Lambda implied by observations. Therefore, over the last decade, the search for extended or alternative theories of gravity has flourished.

We present the work of an international team at the International Space Science Institute (ISSI) in Bern that worked together to review the current observational and theoretical status of the non-virialised X-ray emission components in clusters of galaxies. The subject is important for the study of large-scale hierarchical structure formation and to shed light on the "missing baryon" problem. The topics of the team work include thermal emission and absorption from the warm-hot intergalactic medium, non-thermal X-ray emission in clusters of galaxies, physical processes and chemical enrichment of this medium and clusters of galaxies, and the relationship between all these processes. One of the main goals of the team is to write and discuss a series of review papers on this subject. These reviews are intended as introductory text and reference for scientists wishing to work actively in this field. The team consists of sixteen experts in observations, theory and numerical simulations.

The study of the metal enrichment of the intra-cluster and inter-galactic media (ICM and IGM) represents a direct means to reconstruct the past history of star formation, the role of feedback processes and the gas-dynamical processes which determine the evolution of the cosmic baryons. In this paper we review the approaches that have been followed so far to model the enrichment of the ICM in a cosmological context. While our presentation will be focused on the role played by hydrodynamical simulations, we will also discuss other approaches based on semi-analytical models of galaxy formation, also critically discussing pros and cons of the different methods. We will first review the concept of the model of chemical evolution to be implemented in any chemo-dynamical description. We will emphasise how the predictions of this model critically depend on the choice of the stellar initial mass function, on the stellar life-times and on the stellar yields. We will then overview the comparisons presented so far between X-ray observations of the ICM enrichment and model predictions. We will show how the most recent chemo-dynamical models are able to capture the basic features of the observed metal content of the ICM and its evolution. We will conclude by highlighting the open questions in this study and the direction of improvements for cosmological chemo-dynamical models of the next generation.

There are many processes that can transport gas from the galaxies to their environment and enrich the environment in this way with metals. These metal enrichment processes have a large influence on the evolution of both the galaxies and their environment. Various processes can contribute to the gas transfer: ram-pressure stripping, galactic winds, AGN outflows, galaxy-galaxy interactions and others. We review their observational evidence, corresponding simulations, their efficiencies, and their time scales as far as they are known to date. It seems that all processes can contribute to the enrichment. There is not a single process that always dominates the enrichment, because the efficiencies of the processes vary strongly with galaxy and environmental properties.

Non-thermal components are key ingredients for understanding clusters of galaxies. In the hierarchical model of structure formation, shocks and large-scale turbulence are unavoidable in the cluster formation processes. Understanding the amplification and evolution of the magnetic field in galaxy clusters is necessary for modelling both the heat transport and the dissipative processes in the hot intra-cluster plasma. The acceleration, transport and interactions of non-thermal energetic particles are essential for modelling the observed emissions. Therefore, the inclusion of the non-thermal components will be mandatory for simulating accurately the global dynamical processes in clusters. In this review, we summarise the results obtained with the simulations of the formation of galaxy clusters which address the issues of shocks, magnetic field, cosmic ray particles and turbulence.

Modern hydrodynamical simulations offer nowadays a powerful means to trace the evolution of the X-ray properties of the intra-cluster medium (ICM) during the cosmological history of the hierarchical build up of galaxy clusters. In this paper we review the current status of these simulations and how their predictions fare in reproducing the most recent X-ray observations of clusters. After briefly discussing the shortcomings of the self-similar model, based on assuming that gravity only drives the evolution of the ICM, we discuss how the processes of gas cooling and non-gravitational heating are expected to bring model predictions into better agreement with observational data. We then present results from the hydrodynamical simulations, performed by different groups, and how they compare with observational data. As terms of comparison, we use X-ray scaling relations between mass, luminosity, temperature and pressure, as well as the profiles of temperature and entropy. The results of this comparison can be summarised as follows: (a) simulations, which include gas cooling, star formation and supernova feedback, are generally successful in reproducing the X-ray properties of the ICM outside the core regions; (b) simulations generally fail in reproducing the observed ``cool core'' structure, in that they have serious difficulties in regulating overcooling, thereby producing steep negative central temperature profiles. This discrepancy calls for the need of introducing other physical processes, such as energy feedback from active galactic nuclei, which should compensate the radiative losses of the gas with high density, low entropy and short cooling time, which is observed to reside in the innermost regions of galaxy clusters.

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.

Clusters of galaxies are self-gravitating systems of mass ~10^14-10^15 Msun. They consist of dark matter (~80 %), hot diffuse intracluster plasma (< 20 %) and a small fraction of stars, dust, and cold gas, mostly locked in galaxies. In most clusters, scaling relations between their properties testify that the cluster components are in approximate dynamical equilibrium within the cluster gravitational potential well. However, spatially inhomogeneous thermal and non-thermal emission of the intracluster medium (ICM), observed in some clusters in the X-ray and radio bands, and the kinematic and morphological segregation of galaxies are a signature of non-gravitational processes, ongoing cluster merging and interactions. In the current bottom-up scenario for the formation of cosmic structure, clusters are the most massive nodes of the filamentary large-scale structure of the cosmic web and form by anisotropic and episodic accretion of mass. In this model of the universe dominated by cold dark matter, at the present time most baryons are expected to be in a diffuse component rather than in stars and galaxies; moreover, ~50 % of this diffuse component has temperature ~0.01-1 keV and permeates the filamentary distribution of the dark matter. The temperature of this Warm-Hot Intergalactic Medium (WHIM) increases with the local density and its search in the outer regions of clusters and lower density regions has been the quest of much recent observational effort. Over the last thirty years, an impressive coherent picture of the formation and evolution of cosmic structures has emerged from the intense interplay between observations, theory and numerical experiments. Future efforts will continue to test whether this picture keeps being valid, needs corrections or suffers dramatic failures in its predictive power.

We present a new determination of the cluster mass function and velocity dispersion function in a volume $\sim10^7 h^3$Mpc$^{-3}$ using the Fourth Data Release of the Sloan Digital Sky Survey (SDSS). We use the caustic technique to remove foreground and background galaxies. The cluster virial mass function agrees well with recent estimates from both X-ray observations and cluster richnesses. The mass function lies between those predicted by the First-Year and Three-Year WMAP data. We constrain the cosmological parameters $\Omega_m$ and $\sigma_8$ and find good agreement with WMAP and constraints from other techniques. With the CIRS mass function alone, we estimate $\Omega_m=0.24^{+0.14}_{-0.09}$ and $\sigma_8=0.92^{+0.24}_{-0.19}$, or $\sigma_8=0.84\pm$0.03 when holding $\Omega_m=0.3$ fixed. We also use the WMAP parameters as priors and constrain velocity segregation in clusters. Using the First and Third-Year results, we infer velocity segregation of $\sigma_{gxy}/\sigma_{DM}\approx0.94\pm$0.05 or 1.28$\pm$0.06 respectively. We compare the velocity dispersion function of clusters to that of early-type galaxies and conclude that clusters comprise the high-velocity end of the velocity dispersion function of dark matter haloes. The evolution of cluster abundances provides constraints on dark energy models; the mass function presented here offers an important low redshift calibration benchmark.

We compare the properties of galaxy groups extracted from the Updated Zwicky Catalogue (UZC) with those of groups extracted from N-body simulations of the local Universe, in a LambdaCDM and a tauCDM cosmology. In the simulations, the initial conditions of the dark matter density field are set to reproduce the present time distribution of the galaxies within 80 Mpc/h from the Milky Way. These initial conditions minimize the uncertainty originated by cosmic variance, which has affected previous analyses of this small volume of the Universe. The simulations also model the evolution of the photometric properties of the galaxy population with semi-analytic prescriptions. The models yield a galaxy luminosity function sensibly different from that of the UZC and are unable to reproduce the distribution of groups and their luminosity content. The discrepancy between the model and the UZC reduces substantially, if we redistribute the luminosity among the galaxies in the simulation according to the UZC luminosity function while preserving the galaxy luminosity rank. The modified LambdaCDM model provides the best match to the UZC: the abundances of groups by harmonic radius, velocity dispersion, mass and luminosity are consistent with observations. We find that this model also reproduces the halo occupation number of groups and clusters. However, the large-scale distribution of groups is marginally consistent with the UZC and the redshift-space correlation function of galaxies on scales larger than 6 Mpc/h is still more than 3-sigma smaller than observed. We conclude that reproducing the properties of the observed groups certainly requires a more sophisticated treatment of galaxy formation, and possibly an improvement of the dark matter model.

We use the Fourth Data Release of the Sloan Digital Sky Survey to test the ubiquity of infall patterns around galaxy clusters and measure cluster mass profiles to large radii. We match X-ray cluster catalogs with SDSS, search for infall patterns, and compute mass profiles for a complete sample of X-ray selected clusters. Very clean infall patterns are apparent in most of the clusters, with the fraction decreasing with increasing redshift due to shallower sampling. All 72 clusters in a well-defined sample limited by redshift (ensuring good sampling) and X-ray flux (excluding superpositions) show infall patterns sufficient to apply the caustic technique. This sample is by far the largest sample of cluster mass profiles extending to large radii to date. Similar to CAIRNS, cluster infall patterns are better defined in observations than in simulations. Further work is needed to determine the source of this difference. We use the infall patterns to compute mass profiles for 72 clusters and compare them to model profiles. Cluster scaling relations using caustic masses agree well with those using X-ray or virial mass estimates, confirming the reliability of the caustic technique. We confirm the conclusion of CAIRNS that cluster infall regions are well fit by NFW and Hernquist profiles and poorly fit by singular isothermal spheres. This much larger sample enables new comparisons of cluster properties with those in simulations. The shapes (specifically, NFW concentrations) of the mass profiles agree well with the predictions of simulations. The mass inside the turnaround radius is on average 2.19$\pm$0.18 times that within the virial radius. This ratio agrees well with recent predictions from simulations of the final masses of dark matter haloes.

We study the efficiency and reliability of cluster mass estimators that are based on the projected phase-space distribution of galaxies in a cluster region. To this aim, we analyse a data-set of 62 clusters extracted from a concordance LCDM cosmological hydrodynamical simulation. Galaxies (or Dark Matter particles) are first selected in cylinders of given radius (from 0.5 to 1.5 Mpc/h) and ~200 Mpc/h length. Cluster members are then identified by applying a suitable interloper removal algorithm. Two cluster mass estimators are considered: the virial mass estimator (Mvir), and a mass estimator (Msigma) based entirely on the cluster velocity dispersion estimate. Mvir overestimates the true mass by ~10%, and Msigma underestimates the true mass by ~15%, on average, for sample sizes of > 60 cluster members. For smaller sample sizes, the bias of the virial mass estimator substantially increases, while the Msigma estimator becomes essentially unbiased. The dispersion of both mass estimates increases by a factor ~2 as the number of cluster members decreases from ~400 to ~20. The bias in the Mvir estimates is reduced in clusters without significant evidence for subclustering, and when only early-type galaxies are selected. Radially-dependent incompleteness can drastically affect Mvir estimates, but leaves the Msigma estimates almost unaffected. Other observational effects, like centering and velocity errors, and different observational apertures, have little effect on the mass estimates. (Abridged)

The angular-diameter distance D_A of a galaxy cluster can be measuread by combining its X-ray emission with the cosmic microwave background fluctution due to the Sunyaev-Zeldovich effect. The application of this distance indicator usually assumes that the cluster is spherically symmetric, the gas is distributed according to the isothermal beta-model, and the X-ray temperature is an unbiased measure of the electron temperature. We test these assumptions with galaxy clusters extracted from an extended set of cosmological N-body/hydrodynamical simulations of a LCDM concordance cosmology, which include the effect of radiative cooling, star formation and energy feedback from supernovae. We find that, due to the steep temperature gradients which are present in the central regions of simulated clusters, the assumption of isothermal gas leads to a significant underestimate of D_A. This bias is efficiently corrected by using the polytropic version of the beta-model to account for the presence of temperature gradients. In this case, once irregular clusters are removed, the correct value of D_A is recovered with a ~ 5 per cent accuracy on average, with a ~ 20 per cent intrinsic scatter due to cluster asphericity. This result is valid when using either the electron temperature or a spectroscopic-like temperature. When using instead the emission-weighted definition for the temperature of the simulated clusters, D_A is biased low by \~ 20 per cent. We discuss the implications of our results for an accurate determination of the Hubble constant H_0 and of the density parameter Omega_m. We find that H_0 can be potentially recovered with exquisite precision, while the resulting estimate of Omega_m, which is unbiased, has typical errors Delta(Omega_m) ~ 0.05.

We examine the systematics affecting the X-ray mass estimators applied to a set of five simulated galaxy clusters. They have been processed through the X-ray Map Simulator, X-MAS, to provide Chandra-like long exposures that are analyzed to reconstruct the gas temperature, density, and mass profiles used as input. We find that at R_2500 the mass profile obtained via a direct application of the hydrostatic equilibrium equation is consistent within 1 sigma with the actual mass; although we notice this estimator shows high statistical errors due to high level of Chandra background. Instead, the poorness of the beta-model in describing the gas density profile makes the evaluated masses to be underestimated by \sim 40 per cent with respect to the true mass, both with an isothermal and a polytropic temperature profile. We also test ways to recover the mass by adopting an analytic mass model, such as those proposed by Navarro et al. (1997) and Rasia et al. (2004), and fitting the temperature profile expected from the hydrostatic equilibrium equation to the observed one. These methods and the one of the hydrostatic equilibrium equation provide a more robust mass estimation than the ones based on the beta-model. In the present work the main limitation for a precise mass reconstruction is to ascribe to the relatively high level of the background chosen to reproduce the Chandra one. After artificially reducing it by a factor of 100, we find that the estimated mass significantly underestimates the true mass profiles. This is manly due (i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to the bias towards lower values of the X-ray temperature measurements because of the complex thermal structure of the emitting plasma.

We study the properties of the diffuse X-ray background by using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model. The simulation follows gravitational and gas dynamics and includes a treatment of physical processes like radiative cooling, star formation and supernova feedback. From the simulation outputs, we produce a set of two-dimensional maps of the intergalactic medium X-ray emission integrated over redshift. We find that the signal in the soft (0.5-2 keV) band is lognormally distributed with a mean intensity of about 4 10^-12 erg s^-1 cm^-2 deg^-2; approximately 40 per cent of the emission originates from warm-hot gas (defined as baryons with 10^5<T<10^7 K), and 90 per cent comes from structures at z<0.9. Since the spectrum is soft, being mostly provided by the intergalactic medium at low temperature, the total mean intensity in the hard (2-10 keV) X-ray band is smaller by a factor of about 4. In order to constrain the physical processes included in our simulation, we compare our results with the observed upper limit (1.2 +/- 0.3) 10^-12 erg s^-1 cm^-2 deg^-2 of the soft X-ray emission due to diffuse gas. To this purpose, we remove the contributions of observable extended objects (groups and clusters of galaxies) from the simulated maps by adopting different detectability criteria which are calibrated on the properties of systems at intermediate redshifts observed by Chandra. We show that the simulated diffuse soft X-ray emission is consistent with the present observed upper limit. However, if future measurements will decrease the level of the unresolved X-ray background by a factor of two, a more efficient feedback mechanism should be required to suppress the soft emission of the gas residing within filaments and group-size haloes.

We discuss an extended set of Tree+SPH simulations of galaxy clusters, with the goal of investigating the interplay between numerical resolution effects and star-formation/feedback processes. The simulated clusters span the mass range (0.1-2.3) 10^{15}Msun/h, with mass resolution varying by several decades. At the highest achieved resolution, we resolve the virial region of a Virgo-like cluster with more than 2 million gas particles and with at least as many dark-matter (DM) particles. Our resolution study confirms that, in the absence of an efficient feedback mechanism, runaway cooling leads to about 35 per cent of baryons in clusters to be locked up in long lived stars at our highest resolution, with no evidence of convergence. However, including feedback causes the fraction of cooled baryons to converge at about 15 per cent already at modest resolution. Feedback also stabilizes other gas-related quantities, such as radial profiles of entropy, gas density and temperature, against variations due to changes in resolution. We also investigate the influence of the gravitational force softening length, and that of numerical heating of the gas induced by two-body encounters between DM and lighter gas particles. We show that simulations where more DM than gas particles are used, show a significantly enhanced efficiency of star formation at z>3. Our results are important for establishing and delineating the regime of numerical reliability of the present generation of hydrodynamical simulations of galaxy clusters.

We investigate the environmental dependence of star formation in cluster virial regions and infall regions as part of CAIRNS (Cluster And Infall Region Nearby Survey), a large spectroscopic survey of the infall regions surrounding nine nearby rich clusters of galaxies. Our long-slit spectroscopy yields estimates of star formation rates in environments from cluster cores to the general large-scale structure. The fraction of galaxies with current star formation in their inner disks as traced by H-alpha emission increases with distance from the cluster and converges to the ``field'' value only at 2-3 virial radii, in agreement with other investigations. However, among galaxies with significant current star formation (EW[Ha]geq2\AA), there is no difference in the distribution of EW[Ha] inside and outside the virial radius. This surprising result, first seen by Carter et al., suggests that (1) star formation is truncated on either very short timescales or only at moderate and high redshifts or (2) that projection effects contaminate the measurement. The number density profiles of star-forming and non-star-forming galaxies indicate that, among galaxies projected inside the virial radius, at least half of the former and 20% of the latter are ``infall interlopers,'' galaxies in the infall region but outside the virial region. The kinematics of star-forming galaxies in the infall region closely match those of absorption-dominated galaxies. This result shows that the star forming galaxies in the infall regions are not interlopers from the field and excludes one model of the backsplash scenario of galaxy transformation. Finally, we quantify systematic uncertainties in estimating the global star formation in galaxies from their inner disks.

There are only two methods for estimating the mass distribution in the outer regions of galaxy clusters, where virial equilibrium does not hold: weak gravitational lensing and identification of caustics in redshift space. For the first time, we apply both methods to three clusters: A2390, MS1358 and Cl 0024. The two measures are in remarkably good agreement out to ~2 Mpc/h from the cluster centers. This result demonstrates that the caustic technique is a valuable complement to weak lensing. With a few tens of redshifts per (Mpc/h)^2 within the cluster, the caustic method is applicable for any z<~0.5.

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $\Lambda$ Cold Dark Matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of -0.8 or less. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement. Voids very clearly correspond to minima in the smoothed initial density field. We find a universal void mass profile of the form $\rho(<r)/\rho(r_{\rm eff}) \propto \exp[(r/r_{\rm eff})^\alpha]$ where $r_{\rm eff}$ is the effective radius of a void and $\alpha\sim 2$. The mass function of haloes in voids is steeper than that of haloes that populate denser regions. The abundances of void haloes seem to evolve somewhat more strongly between redshifts $\sim 1$ and 0 than the global abundances of haloes. (abridged)

The fluctuations in the Cosmic Microwave Background (CMB) intensity due to the Sunyaev-Zeldovich (SZ) effect are the sum of a thermal and a kinetic contribution. Separating the two components to measure the peculiar velocity of galaxy clusters requires radio and microwave observations at three or more frequencies, and knowledge of the temperature T_e of the intracluster medium weighted by the electron number density. To quantify the systematics of this procedure, we extract a sample of 117 massive clusters at redshift z=0 from an N-body hydrodynamical simulation, with 2x480^3 particles, of a cosmological volume 192 Mpc/h on a side of a flat Cold Dark Matter model with Omega_0=0.3 and Lambda=0.7. Our simulation includes radiative cooling, star formation and the effect of feedback and galactic winds from supernovae. We find that (1) our simulated clusters reproduce the observed scaling relations between X-ray and SZ properties; (2) bulk flows internal to the intracluster medium affect the velocity estimate by less than 200 km/s in 93 per cent of the cases; (3) using the X-ray emission weighted temperature, as an estimate of T_e, can overestimate the peculiar velocity by 20-50 per cent, if the microwave observations do not spatially resolve the cluster. For spatially resolved clusters, the assumptions on the spatial distribution of the ICM, required to separate the two SZ components, still produce a velocity overestimate of 10-20 per cent, even with an unbiased measure of T_e. Thanks to the large size of our cluster samples, these results set a robust lower limit of 200 km/s to the systematic errors that will affect upcoming measures of cluster peculiar velocities with the SZ effect.

The detection of excess of soft X-ray or Extreme Ultraviolet (EUV) radiation, above the thermal contribution from the hot intracluster medium (ICM), has been a controversial subject ever since the initial discovery of this phenomenon. We use a large--scale hydrodynamical simulation of a concordance $\Lambda$CDM model, to investigate the possible thermal origin for such an excess in a set of 20 simulated clusters having temperatures in the range 1--7 keV. Simulated clusters are analysed by mimicking the observational procedure applied to ROSAT--PSPC data, which for the first time showed evidences for the soft X-ray excess. For cluster--centric distances $0.4< R/R_{\rm vir}< 0.7$ we detect a significant excess in most of the simulated clusters, whose relative amount changes from cluster to cluster and, for the same cluster, by changing the projection direction. In about 30 per cent of the cases, the soft X-ray flux is measured to be at least 50 per cent larger than predicted by the one--temperature plasma model. We find that this excess is generated in most cases within the cluster virialized regions. It is mainly contributed by low--entropy and high--density gas associated with merging sub--halos, rather than to diffuse warm gas. Only in a few cases the excess arises from fore/background groups observed in projection, while no evidence is found for a significant contribution from gas lying within large--scale filaments. We compute the distribution of the relative soft excess, as a function of the cluster--centric distance, and compare it with the observational result by Bonamente et al. (2003) for the Coma cluster. Similar to observations, we find that the relative excess increases with the distance from the cluster center, with no significant excess detected for $R<0.4R_{\rm vir}$. (abridged)

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is thermally complex, T_ew fails at reproducing T_spec, and propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T_sl relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and \beta--model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X-ray temperature function could bias low \sigma_8 by 10-20 per cent. This may alleviate the tension between the value of \sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

(ABRIDGED) We analyze the X-ray properties of a sample of local and high redshift galaxy clusters extracted from a large cosmological hydrodynamical simulation. This simulation has been realized using the Tree+SPH code GADGET-2 for a LambdaCDM model. In our analysis, we consider only objects with T_ew >2 keV and adopt an approach that mimics observations, associating with each measurement an error comparable with recent observations and providing best-fit results via robust techniques. Within the clusters, baryons are distributed among (i) a cold neutral phase, with a relative contribution that increases from less than 1 to 3 per cent at higher redshift, (ii) stars which contribute with about 20 per cent and (iii) the X-ray emitting plasma that contributes by 80 (76) per cent at z=0 (1) to the total baryonic budget. A depletion of the cosmic baryon fraction of ~7 (at z=0) and 5 (at z=1) per cent is measured at the virial radius, R_vir, in good agreement with adiabatic hydrodynamical simulations. We confirm that, also at redshift >0.5, power-law relations hold between gas temperature, T, bolometric luminosity, L, central entropy, S, gas mass, M_gas, and total gravitating mass, M_tot and that these relations are steeper than predicted by simple gravitational collapse. A significant, negative evolution in the L-T and L-M_tot relations and positive evolution in the S-T relation are detected at 0.5 < z < 1 in this set of simulated galaxy clusters. This is partially consistent with recent analyses of the observed properties of z>0.5 X-ray galaxy clusters. By fixing the slope to the values predicted by simple gravitational collapse, we measure at high redshift normalizations lower by 10-40 per cent in the L-T, M_tot-T, M_gas-T, f_gas-T and L-M_tot relations than the observed estimates.

In this paper we present the analysis of Beppo-SAX observations of the cluster pairs A2061-A2067 and A2122-A2124, located in the Corona Borealis supercluster, which have been selected as candidate merging clusters. The aim of this work is to study the physics of the intracluster medium and to look for the possible presence of merging signatures. We derived the global temperatures and abundances and the temperature profiles and maps for these clusters. We do not find any significant evidence of interaction between the clusters forming the pairs, but we detect a candidate shock inside A2061. On the basis of the X-ray and optical properties of this cluster we propose a scenario in which a group is falling inside A2061. This interaction is in the phase in which the cores have not encountered yet and in which the formation of a shock is expected.

CAIRNS (Cluster And Infall Region Nearby Survey) is a spectroscopic survey of the infall regions surrounding nine nearby rich clusters of galaxies. In Paper I, we used redshifts within $\sim 10\Mpc$ of the centers of the clusters to determine the mass profiles of the clusters based on the phase space distribution of the galaxies. Here, we use 2MASS photometry and an additional \ncziinew redshifts to investigate the environmental dependence of near-infrared mass-to-light ratios. In the virial regions, the halo occupation function is non-linear; the number of bright galaxies per halo increases more slowly than the mass of the halo. On larger scales, the light contained in galaxies is less clustered than the mass in rich clusters. Specifically, the mass-to-light ratio inside the virial radius is a factor of $1.8\pm0.3$ larger than that outside the virial radius. This difference could result from changing fractions of baryonic to total matter or from variations in the efficiency of galaxy formation or disruption with environment. The average mass-to-light ratio $M/L_K = 53\pm 5 h$ implies $\Omega_m = 0.18\pm 0.03$ (statistical) using the luminosity density based on 2dFGRS data. These results are difficult to reconcile with independent methods which suggest higher $\Omega_m$. Reconciling these values by invoking bias requires that the typical value of $M/L_K$ changes significantly at densities of $\lesssim3\rho_c$.

We study the properties of the diffuse light in galaxy clusters forming in a large hydrodynamical cosmological simulation of the Lambda-CDM cosmology. The simulation includes a model for radiative cooling, star formation in dense cold gas, and feedback by SN-II explosions. We select clusters having mass M>10^(14) h^(-1) Msun and study the spatial distribution of their star particles. While most stellar light is concentrated in gravitationally bound galaxies orbiting in the cluster potential, we find evidence for a substantial diffuse component, which may account for the extended halos of light observed around central cD galaxies. We find that more massive simulated clusters have a larger fraction of stars in the diffuse light than the less massive ones. The intracluster light is more centrally concentrated than the galaxy light, and the stars in the diffuse component are on average older than the stars in cluster galaxies, supporting the view that the diffuse light is not a random sampling of the stellar population in the cluster galaxies. We thus expect that at least ~10% of the stars in a cluster may be distributed as intracluster light, largely hidden thus far due to its very low surface brightness.

We present results on the X-ray properties of clusters and groups of galaxies, extracted from a large hydrodynamical simulation. We used the GADGET code to simulate a LambdaCDM model within a box of 192 Mpc/h on a side, with 480^3 dark matter particles and as many gas particles. The simulation includes radiative cooling, star formation and supernova feedback. The simulated M-T relation is consistent with observations once we mimic the procedure for mass estimates applied to real clusters. Also, with the adopted choices of Omega_m=0.3 and sigma_8=0.8 the resulting XTF agrees with observational determinations. The L-T relation also agrees with observations for clusters with T>2 keV, with no change of slope at the scale of groups. The entropy in central cluster regions is higher than predicted by gravitational heating alone, the excess being almost the same for clusters and groups. The simulated clusters appear to have suffered some overcooling, with f*~0.2, thus about twice as large as the value observed. Interestingly, temperature profiles are found to steadily increase toward cluster centers. They decrease in the outer regions, much like observational data do at r>0.2r_vir, while not showing an isothermal regime followed by a smooth temperature decline in the innermost regions.

The CAIRNS (Cluster And Infall Region Nearby Survey) project is a spectroscopic survey of the infall regions surrounding eight nearby, rich, X-ray luminous clusters of galaxies. We collect 15665 redshifts (3471 new or remeasured) within \sim 5-10 Mpc of the centers of the clusters, making it the largest study of the infall regions of clusters. We determine cluster membership and the mass profiles of the clusters based on the phase space distribution of the galaxies. All of the clusters display decreasing velocity dispersion profiles. The mass profiles are fit well by functional forms based on numerical simulations but exclude an isothermal sphere. Specifically, NFW and Hernquist models provide good descriptions of cluster mass profiles to their turnaround radii. Our sample shows that the predicted infall pattern is ubiquitous in rich, X-ray luminous clusters over a large mass range. The caustic mass estimates are in excellent agreement with independent X-ray estimates at small radii and with virial estimates at intermediate radii. The mean ratio of the caustic mass to the X-ray mass is 1.03\pm0.11 and the mean ratio of the caustic mass to the virial mass (when corrected for the surface pressure term) is 0.93\pm0.07. We further demonstrate that the caustic technique provides reasonable mass estimates even in merging clusters.

We use a simple analytic model to compute the angular correlation function of clusters identified in upcoming thermal SZ effect surveys. We then compute the expected fraction of close pairs of clusters on the sky that are also close along the line of sight. We show how the expected number of cluster pairs as a function of redshift is sensitive to the assumed biasing relation between the cluster and the mass distribution. We find that, in a LambdaCDM model, the fraction of physically associated pairs is 70% for angular separations smaller than 20 arcmin and clusters with specific flux difference larger than 200 mJy at 143 GHz. The agreement of our analytic results with the Hubble volume N-body simulations is satisfactory. These results quantify the feasibility of using SZ surveys to compile catalogues of superclusters at any redshifts.

Using a redshift survey of 1323 galaxies (1092 new or remeasured) in a region of 95 square degrees centered on the nearby galaxy cluster Abell 2199, we analyze the supercluster containing A2199, A2197, and an X-ray group. The caustic technique accurately reproduces the true mass profiles of simulated simple superclusters (i.e., superclusters where the virial mass of one cluster is 2-10 times the virial mass of all other clusters in the supercluster). We calculate the masses of the two main components of A2197 (A2197W and A2197E) using archival X-ray observations and demonstrate that the A2199 supercluster is simple and thus that the caustic technique should yield an accurate mass profile. The mass profile is uncertain by ~30% within 3 Mpc/h and by a factor of two within 8 Mpc/h and is one of only a few for a supercluster on such large scales. Independent X-ray mass estimates agree with our results at all radii where they overlap. The mass profile strongly disagrees with an isothermal sphere profile but agrees with profiles suggested by simulations. We discuss the interplay of the supercluster dynamics and the dynamics of the bound subclusters. The agreement between the infall mass profile and other techniques shows that the caustic technique is surprisingly robust for simple superclusters (abridged).

Using a redshift survey of 1779 galaxies and photometry from the 2-Micron All-Sky Survey (2MASS) covering 200 square degrees, we calculate independent mass and light profiles for the infall region of the Coma cluster of galaxies. The redshift survey is complete to $K_s=12.2$ (622 galaxies), 1.2 magnitudes fainter than $M^*_{K_s}$ at the distance of Coma. We confirm the mass profile obtained by Geller, Diaferio, & Kurtz. The enclosed mass-to-light ratio measured in the $K_s$ band is approximately constant to a radius of $10 \Mpc$, where $M/L_{K_s}= 75\pm 23\mlsun$, in agreement with weak lensing results on similar scales. Within $2.5\Mpc$, X-ray estimates yield similar mass-to-light ratios (67$\pm32h$). The constant enclosed mass-to-light ratio with radius suggests that K-band light from bright galaxies in clusters traces the total mass on scales $\lesssim10 \Mpc$. Uncertainties in the mass profile imply that the mass-to-light ratio inside $r_{200}$ may be as much as a factor of 2.5 larger than that outside $r_{200}$. These data demonstrate that K-band light is not positively biased with respect to the mass; we cannot rule out antibias. These results imply $\Omega_m = 0.17 \pm 0.05$. Estimates of possible variations in $M/L_{K_s}$ with radius suggest that the density parameter is no smaller than $\Omega_m \approx 0.08$.

We use numerical simulations of a (480 Mpc/h)^3 volume to show that the distribution of peak heights in maps of the temperature fluctuations from the kinematic and thermal Sunyaev-Zeldovich effects will be highly non-Gaussian, and very different from the peak height distribution of a Gaussian random field. We then show that it is a good approximation to assume that each peak in either SZ effect is associated with one and only one dark matter halo. This allows us to use our knowledge of the properties of haloes to estimate the peak height distributions. At fixed optical depth, the distribution of peak heights due to the kinematic effect is Gaussian, with a width which is approximately proportional to optical depth; the non-Gaussianity comes from summing over a range of optical depths. The optical depth is an increasing function of halo mass, and the distribution of halo speeds is Gaussian, with a dispersion which is approximately independent of halo mass. This means that observations of the kinematic effect can be used to put constraints on how the abundance of massive clusters evolves, and on the evolution of cluster velocities. The non-Gaussianity of the thermal effect, on the other hand, comes primarily from the fact that, on average, the effect is larger in more massive haloes, and the distribution of halo masses is highly non-Gaussian. We also show that because haloes of the same mass may have a range of density and velocity dispersion profiles, the relation between halo mass and the amplitude of the thermal effect is not deterministic, but has some scatter.

We discuss H-band (1.65 micrometer) near-infrared photometry of the central 3x3h^{-2} Mpc^2 of Abell 1644 to a limiting M_H~M^*_H+3 (throughout this paper H_0=100 h km s^{-1} Mpc^{-1}). There are 861 galaxies in the photometric survey region. We also measured radial velocities of 155 galaxies; 141 of these are cluster members within 2.44 h^{-1} Mpc of the cluster center. The cluster velocity dispersion of sigma~1000 km s^{-1} remains constant out to the limiting radius. We find no evidence for substructure in the cluster. The cluster mass within R=2.4 h^{-1} Mpc is 7.6+-1.3x10^{14} h^{-1} M_\odot. We compute the cluster luminosity function; the Schechter parameters \alpha=-1.14\pm0.08 and M^*_H=-24.3\pm0.2 (with h=0.5) agree well with other H-band luminosity functions. From the virial theorem and the caustic method we compute one of the first mass-to-light ratios at H; the result is M/L_H=82-127 h M_\odot/L_\odot within 1.5 h^{-1} Mpc. This ratio corresponds to 374-579 h M_\odot/L_\odot at R. The agreement of our IR measurement with previous $M/L$ determinations indicates that at low redshift dust and young stellar populations may produce only negligible systematic errors in optical mass-to-light ratios.

Using a large redshift survey covering 95 square degrees, we demonstrate that the infall region of Abell 2199 contains Abell 2197, one or two X-ray emitting groups, and up to five additional groups identified in redshift surveys. Our survey shows that the X-ray emitting systems, located at projected radii of $1.^\circ4, 1.^\circ9$, and $5.^\circ1$ (2.2, 3.1, and $8.0 h^{-1}$Mpc), are connected kinematically to A2199. A2197 is itself an optically rich cluster; its weak X-ray emission suggests that it is much less massive than A2199. The absence of a sharp peak in the infall pattern at the position of A2197 supports this hypothesis. The outermost group is well outside the virial region of A2199 and it distorts the infall pattern in redshift space. The two X-ray emitting groups are roughly colinear, suggesting the existence of an extended ($8.0 h^{-1}$Mpc) filament. The identification of these infalling groups provides direct support of hierarchical structure formation; studies of these systems will provide insights into structure evolution. Groups in the infall regions of nearby clusters may offer a unique probe of the physics of the warm/hot ionized medium (WHIM) which is difficult to observe directly with current instruments.

A simple model of how objects of different masses stream towards each other as they cluster gravitationally is described. The model shows how the mean streaming velocity of dark matter particles is related to the motions of the parent dark matter haloes. It also provides a reasonably accurate description of how the pairwise velocity dispersion of dark matter particles differs from that of the parent haloes. The analysis is then extended to describe the streaming motions of galaxies. This shows explicitly that the streaming motions measured in a given galaxy sample depend on how the sample was selected, and shows how to account for this dependence on sample selection. In addition,we show that the pairwise dispersion should also depend on sample type. Our model predicts that, on small scales, redshift space distortions should affect red galaxies more strongly than blue.

We write the correlation function of dark matter particles, xi(r), as the sum of two terms - one which accounts for nonlinear evolution, and dominates on small scales, and another which is essentially the term from linear theory, and dominates on large scales. We use models of the number and spatial distribution of haloes and halo density profiles to describe the nonlinear term and its evolution. The result provides a good description of the evolution of xi(r) in simulations. We then use this decomposition to provide simple and accurate models of how the single particle velocity dispersion evolves with time, and how the first and second moments of the pairwise velocity distribution depend on scale. The key idea is to use the simple physics of linear theory on large scales, the simple physics of the virial theorem on small scales, and our model for the correlation function to tell us how to weight the two types of contributions (linear and nonlinear) to the pairwise velocity statistics. When incorporated into the streaming model, our results will allow a simple accurate description of redshift-space distortions over the entire range of linear to highly nonlinear regimes.

We present a simple model for the shape of the distribution function of galaxy peculiar velocities. We show how both nonlinear and linear theory terms combine to produce a distribution which has an approximately Gaussian core with exponential wings. The model is easily extended to study how the statistic depends on the type of particle used to trace the velocity field (dark matter particles, dark matter haloes, galaxies), and on the density of the environment in which the test particles are. Comparisons with simulations suggest that our model is accurate. We also show that the evolution of the peculiar velocities depends on the local, rather than the global density. Since clusters populate denser regions on average, using cluster velocities with the linear theory scaling may lead to an overestimate of the global value of Omega. Conversely, using linear theory with the global value of Omega to scale cluster velocities from the initial to the present time results in an underestimate of their true velocities. In general, however, the directions of motions of haloes are rather well described by linear theory. Our results help to simplify models of redshift-space distortions considerably.

We describe observations of the nearby cluster of galaxies A576 beyond the virial radius and into the infall region. Using 1057 redshifts, we use the infall pattern in redshift space to determine the mass profile of A576 to a radius of ~4 Mpc/h. This mass estimation technique makes no assumptions about the equilibrium state of the cluster. Within 1 Mpc/h, the mass profile we derive exceeds that determined from X-ray observations by a factor of 2.5. At \~2.5 Mpc/h, however, the mass profile agrees with virial mass estimates. Our mass profile is consistent with a NFW or Hernquist profile, but it is inconsistent with an isothermal sphere. R-band images of a $3^\circ x 3^\circ$ region centered on the cluster allow an independent determination of the cluster light profile. We calculate the integrated mass-to-light ratio as a function of cluster radius; it decreases smoothly from the core to ~4 Mpc/h. The differential dM/dL_R profile decreases more steeply; we find M/L_R ~100 h at ~4 Mpc/h, in good agreement with the mass-to-light ratios of individual galaxies. This value implies $\Omega_m \lesssim 0.4$ at 95% confidence. For a Hernquist model, the best-fit mass profiles differ from the observed surface number density of galaxies; the galaxies have a larger scale radius than the mass. This result is consistent with the centrally peaked $M/L_R$ profile. Similarly, the scale radius of the light profile is larger than that of the mass profile. We discuss some potential systematic effects; none can easily reconcile our results with a constant mass-to-light ratio. (abstract edited)

We combine dissipationless N-body simulations and semi-analytic models of galaxy formation to study the spatial and kinematic distributions of cluster galaxies in a LCDM cosmology. We investigate how the star formation rates, colours and morphologies of galaxies vary as a function of distance from the cluster centre and compare our results with the CNOC1 survey of galaxies from 15 X-ray luminous clusters in the redshift range 0.18 to 0.55. In our model, gas no longer cools onto galaxies after they fall into the cluster and their star formation rates decline on timescales of 1-2 Gyr. Galaxies in cluster cores have lower star formation rates and redder colours than galaxies in the outer regions because they were accreted earlier. Our colour and star formation gradients agree with those those derived from the data. The difference in velocity dispersions between red and blue galaxies observed in the CNOC1 clusters is also well reproduced by the model. We assume that the morphologies of cluster galaxies are determined solely by their merging histories. Morphology gradients in clusters arise naturally, with the fraction of bulge- dominated galaxies highest in cluster cores. We compare these gradients with the CNOC1 data and find excellent agreement for bulge-dominated galaxies. The simulated clusters contain too few galaxies of intermediate bulge-to-disk ratio, suggesting that additional processes may influence the morphological evolution of disk-dominated galaxies in clusters. Although the properties of the cluster galaxies in our model agree extremely well with the data, the same is not true of field galaxies. Both the star formation rates and the colours of bright field galaxies appear to evolve much more strongly from redshift 0.2 to 0.4 in the CNOC1 field sample than in our simulations.

We identify high density regions of supercluster size in high resolution N-body simulations of a representative volume of three Cold Dark Matter Universes. By assuming that (1) the density and peculiar velocities of baryons trace those of the dark matter, and (2) the temperature of plasma is proportional to the velocity dispersion of the dark matter particles in regions where the crossing times is smaller than the supercluster free-fall time, we investigate how thermal motions of electrons in the intra-cluster medium and peculiar velocity of clusters can affect the secondary anisotropies in the cosmic microwave background (CMB). We show that the thermal effect dominates the kinematic effect and that the largest thermal decrements are associated with the most massive clusters in superclusters. Thus, searching for the presence of two or more close large CMB decrements represents a viable strategy for identifying superclusters at cosmological distances. Moreover, maps of the kinematic effect in superclusters are characterized by neighboring large peaks of opposite signs. These peaks can be as high as ~ 10 microK at the arcminute angular resolution. Simultaneous pointed observations of superclusters in the millimeter and submillimeter bands with upcoming sensitive CMB experiments can separate between the thermal and kinematic effect contributions and constrain the evolution of the velocity field in large overdense regions.

We analyse the redshift space topology and geometry of the nearby Universe by computing the Minkowski functionals of the Updated Zwicky Catalogue (UZC). The UZC contains the redshifts of almost 20,000 galaxies, is 96% complete to the limiting magnitude m_Zw=15.5 and includes the Center for Astrophysics (CfA) Redshift Survey (CfA2). From the UZC we can extract volume limited samples reaching a depth of 70 hMpc before sparse sampling dominates. We quantify the shape of the large-scale galaxy distribution by deriving measures of planarity and filamentarity from the Minkowski functionals. The nearby Universe shows a large degree of planarity and a small degree of filamentarity. This quantifies the sheet-like structure of the Great Wall which dominates the northern region (CfA2N) of the UZC. We compare these results with redshift space mock catalogues constructed from high resolution N-body simulations of two Cold Dark Matter models with either a decaying massive neutrino (tauCDM) or a non-zero cosmological constant (LambdaCDM). We use semi-analytic modelling to form and evolve galaxies in these dark matter-only simulations. We are thus able, for the first time, to compile redshift space mock catalogues which contain galaxies, along with their observable properties, rather than dark matter particles alone. In both models the large scale galaxy distribution is less coherent than the observed distribution, especially with regard to the large degree of planarity of the real survey. However, given the small volume of the region studied, this disagreement can still be a result of cosmic variance.

For the first time, the combination of semi-analytic modelling of galaxy formation and N-body simulations of cosmic structure formation enables us to model, at the same time, both the photometric and the clustering properties of galaxies. Two Cold Dark Matter Universes provide a reasonable fit to the observed properties of galaxies, groups and clusters, including luminosities, colours, density and velocity biases. We show how the properties of galaxies and groups on small scales are inextricably connected with the global properties of the Universe.

We present a technique for estimating the mass in the outskirts of galaxy clusters where the usual assumption of dynamical equilibrium is not valid. The method assumes that clusters form through hierarchical clustering and requires only galaxy redshifts and positions on the sky. We apply the method to dissipationless cosmological N-body simulations where galaxies form and evolve according to semi-analytic modelling. The method recovers the actual cluster mass profile within a factor of two to several megaparsecs from the cluster centre. This error originates from projection effects, sparse sampling, and contamination by foreground and background galaxies. In the absence of velocity biases, this method can provide an estimate of the mass-to-light ratio on scales ~1-10 Mpc/h where this quantity is still poorly known.

This is the third paper in a series which combines N-body simulations and semi-analytic modelling to provide a fully spatially resolved simulation of the galaxy formation and clustering processes. Here we extract mock redshift surveys from our simulations: a Cold Dark Matter model with either Omega_0=1 (tauCDM) or Omega_0=0.3 and Lambda=0.7 (LambdaCDM). We compare the mock catalogues with the northern region (CfA2N) of the Center for Astrophysics (CfA) Redshift Surveys. We study the properties of galaxy groups and clusters identified using standard observational techniques and we study the relation of these groups to real virialised systems. Most features of CfA2N groups are reproduced quite well by both models with no obvious dependence on Omega_0. Redshift space correlations and pairwise velocities are also similar in the two cosmologies. The luminosity functions predicted by our galaxy formation models depend sensitively on the treatment of star formation and feedback. For the particular choices of Paper I they agree poorly with the CfA survey. To isolate the effect of this discrepancy on our mock redshift surveys, we modify galaxy luminosities in our simulations to reproduce the CfA luminosity function exactly. This adjustment improves agreement with the observed abundance of groups, which depends primarily on the galaxy luminosity density, but other statistics, connected more closely with the underlying mass distribution, remain unaffected. Regardless of the luminosity function adopted, modest differences with observation remain. These can be attributed to the presence of the ``Great Wall'' in the CfA2N. It is unclear whether the greater coherence of the real structure is a result of cosmic variance, given the relatively small region studied, or reflects a physical deficiency of the models.

We use a new redshift survey complete to m_R=15.4 within 4.25 deg from the center of the Coma cluster to measure the mass profile of the cluster to r=5.5 Mpc/h. We extend the profile to r=10 Mpc/h with a further sample complete to m_R=15.4 in 42% of the area within a 10 deg radius and to m_{Zw}=15.5 in the remaining area. Galaxies within this region are falling onto the cluster on moderately radial orbits and thus do not satisfy virial equilibrium. Nonetheless, identification of the caustics in redshift space provides an estimate of the gravitational potential at radius r and hence of the system mass, M(<10 Mpc/h)=(1.65+-0.41)10^{15} M_\odot/h (1-sigma error). Previous mass estimates derived from optical and X-ray observations are limited to r<2.5 Mpc/h. Our mass profile is consistent with these estimates but extends to distances four times as large. Over the entire range, the mass increases with r at the rate expected for a Navarro, Frenk & White (1997) density profile.

In hierarchical cosmologies the evolution of galaxy clustering depends both on cosmological quantities such as Omega and Lambda, which determine how dark matter halos form and evolve, and on the physical processes - cooling, star formation and feedback - which drive the formation of galaxies within these merging halos. In this paper, we combine dissipationless cosmological N-body simulations and semi-analytic models of galaxy formation in order to study how these two aspects interact. We focus on the differences in clustering predicted for galaxies of differing luminosity, colour, morphology and star formation rate and on what these differences can teach us about the galaxy formation process. We show that a "dip" in the amplitude of galaxy correlations between z=0 and z=1 can be an important diagnostic. Such a dip occurs in low-density CDM models because structure forms early and dark matter halos of 10**12 solar masses, containing galaxies with luminosities around L*, are unbiased tracers of the dark matter over this redshift range; their clustering amplitude then evolves similarly to that of the dark matter. At higher redshifts bright galaxies become strongly biased and the clustering amplitude increases again. In high density models, structure forms late and bias evolves more rapidly. As a result, the clustering amplitude of L* galaxies remains constant from z=0 to 1. The strength of these effects is sensitive to sample selection. The dip is weaker for galaxies with lower star formation rates, redder colours, higher luminosities and earlier morphological types. We explain why this is the case and how it is related to the variation with redshift of the abundance and and environment of the observed galaxies. Studies of clustering evolution as a function of galaxy properties should place strong constraints on models.

We introduce a new technique for following the formation and evolution of galaxies in cosmological N-body simulations. Dissipationless simulations are used to track the formation and merging of dark matter halos as a function of redshift. Simple prescriptions, taken directly from semi-analytic models of galaxy formation, are adopted for cooling, star formation, supernova feedback and the merging of galaxies within the halos. This scheme enables us to study the clustering properties of galaxies and to investigate how selection by type, colour or luminosity influences the results. In this paper, we study properties of the galaxy distribution at z=0. These include luminosity functions, colours, correlation functions, pairwise peculiar velocities, cluster M/L ratios and star formation rates. We focus on two variants of a CDM cosmology: a high- density model with Gamma=0.21 (TCDM) and a low-density model with Omega=0.3 and Lambda=0.7 (LCDM). Both are normalized to reproduce the I-band Tully-Fisher relation near a circular velocity of 220 km/s. Our results depend strongly both on this normalization and on the prescriptions for star formation and feedback. Very different assumptions are required to obtain an acceptable model in the two cases. For TCDM, efficient feedback is required to suppress the growth of galaxies low-mass field halos. Without it, there are too many galaxies and the correlation function turns over below 1 Mpc. For LCDM, feedback must be weak, otherwise too few L* galaxies are produced and the correlation function is too steep. Given the uncertainties in modelling some of the key physical processes, we conclude that it is not yet possible to draw conclusions about the values of cosmological parameters from studies of this kind. Further work on global star formation and feedback effects is required to narrow the range of possibilities

Using high resolution cosmological N-body simulations, we investigate the survival of dark matter satellites falling into larger haloes. Satellites preserve their identity for some time after merging. We compute their loss of mass, energy and angular momentum as dynamical friction, tidal forces and collisions with other satellites dissolve them. We also analyse the evolution of their internal structure. Satellites with less than a few per cent the mass of the main halo may survive for several billion years, whereas larger satellites rapidly sink into the center of the main halo potential well and lose their identity. Penetrating encounters between satellites are frequent and may lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10 per cent on average) is bound to substructure at most redshifts of interest. We discuss the application of these results to the survival and extent of dark matter haloes associated with cluster galaxies, and to interactions between galaxies in clusters. We find that 35-40 per cent of galaxy dark matter haloes are disrupted by the present time. The fraction of satellites undergoing close encounters is similar to the fraction of interacting or merging galaxies in clusters at moderate redshift.

In hierarchical clustering, galaxy clusters accrete mass through the aggregation of smaller systems. Thus, the velocity field of the infall regions of clusters contains significant random motion superimposed on radial infall. Because the purely spherical infall model does not predict the amplitude of the velocity field correctly, methods estimating the cosmological density parameter Omega_0 based on this model yield unreliable biased results. In fact, the amplitude of the velocity field depends on local dynamics and only very weakly on the global properties of the universe. We use N-body simulations of flat and open universes to show that the amplitude of the velocity field of the infall regions of dark matter halos is a direct measure of the escape velocity within these regions. We can use this amplitude to estimate the mass of dark matter halos within a few megaparsecs from the halo center. In this region dynamical equilibrium assumptions do not hold. The method yields a mass estimate with better than 30% accuracy. If galaxies trace the velocity field of the infall regions of clusters reliably, this method provides a straightforward way to estimate the amount of mass surrounding rich galaxy clusters from redshift data alone.

The redshift-space correlation function $\xi_s$ for projected galaxy separations $\simless 1h^{-1}$ Mpc can be expressed as the convolution of the real-space correlation function with the galaxy pairwise velocity distribution function (PVDF). An exponential PVDF yields the best fit to the $\xi_s$ measured from galaxy samples of different redshift surveys. We show that this exponential PVDF is not merely a fitting function but arises from well defined gravitational processes. Two ingredients conspire to yield a PVDF with a nearly exponential shape: (i) the number density $n(\sigma)$ of systems with velocity dispersion $\sigma$; (ii) the unrelaxed dynamical state of most galaxy systems. The former ingredient determines the exponential tail and the latter the central peak of the PVDF. We compare the $\xi_s$ measured on non-linear scales from galaxy samples of the CfA redshift surveys with different models of the PVDF convolved with the measured real-space correlation function. This preliminary comparison indicates that the agreement between model and observations depends strongly on both the underlying cosmological model and the internal dynamics of galaxy systems. Neither parameter dominates. Moreover, the agreement depends sensitively on the accuracy of the galaxy position and velocity measurements.