Search

Refracted Gravity (RG) is a a classical theory of gravity where a gravitational permittivity $ a monotonically-increasing function of the local density rho , is introduced in the Poisson equation to mimic the effect of dark matter at astrophysical scales. We use high precision spectroscopic data of two massive galaxy clusters, MACS J1206.2-0847 at redshift z=0.44, and Abell S1063 (RXC J2248.7-4431) at z=0.35, to determine the total gravitational potential in the context of RG and to constrain the three, supposedly universal, free parameters of this model. Using an upgraded version of the MG-MAMPOSSt algorithm, we perform a kinematic analysis which combines the velocity distribution of the cluster galaxies and the velocity dispersion profile of the stars within the Brightest Cluster Galaxy (BCG). The unprecedented dataset used has been obtained by an extensive spectroscopic campaign carried out with the VIMOS and MUSE spectrographs at the ESO VLT. We found that RG describes the kinematics of these two clusters as well as Newtonian gravity, although the latter is slightly preferred. However, (i) each cluster requires a different set of the three free RG parameters, and (ii) the two sets are inconsistent with other results in the literature at different scales. We discuss the limitation of the method used to constrain the RG parameters as well as possible systematic effects which can give rise to the observed tension, notably deviations from the spherical symmetry and from the dynamical equilibrium of the clusters., Comment: 10 pages, 6 figures. Submitted to A&A

Ram pressure stripping is perhaps the most efficient mechanism for removing gas and quenching galaxies in dense environments as they move through the intergalactic medium. Extreme examples of on-going ram pressure stripping are known as jellyfish galaxies, characterized by a tail of stripped material that can be directly observed in multiple wavelengths. Using the largest homogeneous broad-band optical jellyfish candidate sample in local clusters known to date, we measure the angle between the direction of the tails visible in the galaxies, and the direction towards the host cluster center. We find that $33\%$ of the galaxy tails point away from the cluster center, $18\%$ point towards the cluster center, and $49\%$ point elsewhere. Moreover, we find stronger signatures of ram pressure stripping happening on galaxies with a tail pointing away and towards the cluster center, and larger velocity dispersion profiles for galaxies with tails pointing away. These results are consistent with a scenario where ram pressure stripping has a stronger effect for galaxies following radial orbits on first infall. The results also suggest that in many cases, radially infalling galaxies are able to retain their tails after pericenter and continue to experience significant on-going ram pressure stripping. We further constrain the lifespan of the optical tails from the moment they first appear to the moment they disappear, by comparing the observed tail directions with matched N-body simulations through Bayesian parameter estimation. We obtain that galaxy tails appear for the first time at $\sim 1.16$ R$_{200}$ and disappear $\sim660$ Myr after pericenter., Comment: 19 pages, 17 figures

Observational studies have widely demonstrated that galaxy physical properties are strongly affected by the surrounding environment. On one side, gas inflows provide galaxies with new fuel for star formation. On the other side, the high temperatures and densities of the medium are expected to induce quenching in the star formation. Observations of large structures, in particular filaments at the cluster outskirts (r>2r$_{200}$), are currently limited to the low redshift Universe. We present a multi-band dataset for the cluster MACS J0416.1-2403 (z=0.397), observed in the context of the Galaxy Assembly as a function of Mass and Environment with VST (VST-GAME) survey. The project aims at gathering deep ($r$<24.4) and wide (20x20Mpc$^2$) observations at optical wavelengths for six massive galaxy clusters at 0.2

The study of cluster substructures is important for the determination of the cluster dynamical status, assembly history, and the evolution of cluster galaxies, and it allows to set of constraints on the nature of dark matter and cosmological parameters. We present and test DS+, a new method for the identification and characterization of group-sized substructures in clusters. Our new method is based on the projected positions and line-of-sight velocities of cluster galaxies, and it is an improvement and extension of the traditional method of Dressler & Shectman (1988). We test it on cluster-size cosmological halos extracted from the IllustrisTNG simulations, with virial masses $\rm{14 \lesssim \log (M_{200}/M_{\odot}) \lesssim 14.6}$, that contain on average $\sim 190$ galaxies. We also present an application of our method on a real data set, the Bullet cluster. DS+ is able to identify $\sim 80\%$ of real group galaxies as members of substructures, and at least 60\% of the galaxies assigned to substructures belong to real groups. The physical properties of the real groups are significantly correlated with those of the corresponding detected substructures, albeit with significant scatter, and overestimated on average. Application of the DS+ method to the Bullet cluster confirms the presence and main properties of the high-speed collision and identifies other substructures along the main cluster axis. DS+ proves to be a reliable method for the identification of substructures in clusters. The method is made freely available to the community as a Python code., Comment: 13 pages, 8 figures, accepted for publication in A&A

Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ~ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest-mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses (Mstar <~ 1E10 Msun), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help to shed light on the low-mass tension., Comment: 20 pages, 10 figures, resubmitted to MNRAS following minor revisions

We present deep JVLA observations at 1.4 GHz and 2.7 GHz (full polarization), as well as optical OmegaWINGS/WINGS and X-ray observations of two extended radio galaxies in the IIZW108 galaxy cluster at z = 0.04889. They show a bent tail morphology in agreement with a radio lobed galaxy falling into the cluster potential. Both galaxies are found to possess properties comparable with {narrow-angle} tail galaxies in the literature even though they are part of a low mass cluster. We find a spectral index steepening and an increase in fractional polarization through the galaxy jets and an ordered magnetic field component mostly aligned with the jet direction. This is likely caused by either shear due to the velocity difference of the intracluster medium and the jet fluid and/or magnetic draping of the intracluster medium across the galaxy jets. We find clear evidence that one source is showing two active galactic nuclei (AGN) outbursts from which we expect the AGN has never turned off completely. We show that pure standard electron cooling cannot explain the jet length. We demonstrate therefore that these galaxies can be used as a laboratory to study gentle re-acceleration of relativistic electrons in galaxy jets via transition from laminar to turbulent motion., Comment: accepted for publication in MNRAS

We measure the rate of environmentally-driven star formation quenching in galaxies at $z\sim 1$, using eleven massive ($M\approx 2\times10^{14}\,\mathrm{M}_\odot$) galaxy clusters spanning a redshift range $1.010.5$) we do not find any significant excess of transition galaxies in clusters, relative to a comparison field sample at the same redshift. It is likely that such galaxies were quenched prior to their accretion in the cluster, in group, filament or protocluster environments. For lower stellar mass galaxies ($9.5<\log{M/\mathrm{M}_\odot}<10.5$) there is a small but significant excess of transition galaxies in clusters, accounting for an additional $\sim 5-10$ per cent of the population compared with the field. We show that our data are consistent with a scenario in which 20--30 per cent of low-mass, star-forming galaxies in clusters are environmentally quenched every Gyr, and that this rate slowly declines from $z=1$ to $z=0$. While environmental quenching of these galaxies may include a long delay time during which star formation declines slowly, in most cases this must end with a rapid ($\tau<1$ Gyr) decline in star formation rate., Comment: Accepted for publication in MNRAS Sept 6, 2021

We aim to study how the orbits of galaxies in clusters depend on the prominence of the corresponding central galaxies. We divided our data set of $\sim$ 100 clusters and groups into four samples based on their magnitude gap between the two brightest members, $\Delta m_{12}$. We then stacked all the systems in each sample, in order to create four stacked clusters, and derive the mass and velocity anisotropy profiles for the four groups of clusters using the MAMPOSSt procedure. Once the mass profile is known, we also obtain the (non parametric) velocity anisotropy profile via the inversion of the Jeans equation. In systems with the largest $\Delta m_{12}$, galaxy orbits are prevalently radial, except near the centre, where orbits are isotropic (or tangential when also the central galaxies are considered in the analysis). In the other three samples with smaller $\Delta m_{12}$, galaxy orbits are isotropic or only mildly radial. Our study supports the results of numerical simulations that identify radial orbits of galaxies as the cause of an increasing $\Delta m_{12}$ in groups., Comment: 12 pages, 5 figures, and 6 tables. Accepted for publication in Astronomy and Astrophysics

We present results on the measured shapes of 832 galaxies in 11 galaxy clusters at 1.0 < z <1.4 from the GOGREEN survey. We measure the axis ratio ($q$), the ratio of the minor to the major axis, of the cluster galaxies from near-infrared Hubble Space Telescope imaging using S\'ersic profile fitting and compare them with a field sample. We find that the median $q$ of both star-forming and quiescent galaxies in clusters increases with stellar mass, similar to the field. Comparing the axis ratio distributions between clusters and the field in four mass bins, the distributions for star-forming galaxies in clusters are consistent with those in the field. Conversely, the distributions for quiescent galaxies in the two environments are distinct, most remarkably in $10.1\leq\log(M/{\rm M}_{\odot})<10.5$ where clusters show a flatter distribution, with an excess at low $q$. Modelling the distribution with oblate and triaxial components, we find that the cluster and field sample difference is consistent with an excess of flattened oblate quiescent galaxies in clusters. The oblate population contribution drops at high masses, resulting in a narrower $q$ distribution in the massive population than at lower masses. Using a simple accretion model, we show that the observed $q$ distributions and quenched fractions are consistent with a scenario where no morphological transformation occurs for the environmentally quenched population in the two intermediate mass bins. Our results suggest that environmental quenching mechanism(s) likely produce a population that has a different morphological mix than those resulting from the dominant quenching mechanism in the field., Comment: Accepted for publication in ApJ. 25 pages, 15 figures

We use photometric redshifts and statistical background subtraction to measure stellar mass functions in galaxy group-mass ($4.5-8\times10^{13}~\mathrm{M}_\odot$) haloes at $11$., Comment: Accepted July 6, 2021, MNRAS

Studies of cluster mass and velocity anisotropy profiles are useful tests of dark matter models, and of the assembly history of clusters of galaxies. These studies might be affected by unknown systematics caused by projection effects. We aim at testing observational methods for the determination of mass and velocity anisotropy profiles of clusters of galaxies. Particularly, we focus on the MAMPOSSt technique (Mamon et al. 2013). We use results from two semi-analytic models of galaxy formation coupled with high-resolution N-body cosmological simulations, the catalog of De Lucia & Blaizot (2007) and the FIRE catalog based on the new GAlaxy Evolution and Assembly model. We test the reliability of the Jeans equation in recovering the true mass profile when full projected phase-space information is available. We examine the reliability of the MAMPOSSt method in estimating the true mass and velocity anisotropy profiles of the simulated halos when only projected phase-space information is available, as in observations. The spherical Jeans equation provides a reliable tool for the determination of cluster mass profiles, also for subsamples of tracers separated by galaxy color. Results are equally good for prolate and oblate clusters. Using only projected phase-space information, MAMPOSSt provides estimates of the mass profile with a standard deviation of 35-69 %, and a negative bias of 7-17 %, nearly independent of radius, and that we attribute to the presence of interlopers in the projected samples. The bias changes sign, that is, the mass is over-estimated, for prolate clusters with their major axis aligned along the line-of-sight. MAMPOSSt measures the velocity anisotropy profiles accurately in the inner cluster regions, with a slight overestimate in the outer regions, both for the whole sample of observationally-identified cluster members and separately for red and blue galaxies., Comment: 14 pages, 17 figures

We present the first public data release of the GOGREEN and GCLASS surveys of galaxies in dense environments, spanning a redshift range $0.8

We study the star formation histories (SFHs) and mass-weighted ages of 331 UVJ-selected quiescent galaxies in 11 galaxy clusters and in the field at 11 has been driven by different physical processes than those at play at z=0., Comment: accepted Sept 7 2020, MNRAS

We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0

We present results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs), we find that cluster galaxies have suppressed SFRs at fixed stellar mass in comparison to their field counterparts by a factor of 1.4 $\pm$ 0.1 ($\sim3.3\sigma$) across the stellar mass range: $9.0 < \log(M_{*} /M_{\odot}) < 11.2$. We also find that this modest suppression in the cluster galaxy star-forming main sequence is mass and redshift dependent: the difference between cluster and field increases towards lower stellar masses and lower redshift. When comparing the distribution of cluster and field galaxy SFRs to the star-forming main sequence, we find an overall shift towards lower SFRs in the cluster population, and note the absence of a tail of high SFR galaxies as seen in the field. Given this observed suppression in the cluster galaxy star-forming main sequence, we explore the implications for several scenarios such as formation time differences between cluster and field galaxies, and environmentally-induced star formation quenching and associated timescales., Comment: Accepted for publication in MNRAS on February 25 2020, appended erratum accepted in MNRAS on September 17 2020

With MUSE, Chandra, VLA, ALMA and UVIT data from the GASP programme we study the multiphase baryonic components in a jellyfish galaxy (JW100) with a stellar mass 3.2 X 10^{11} M_sun hosting an AGN. We present its spectacular extraplanar tails of ionized and molecular gas, UV stellar light, X-ray and radio continuum emission. This galaxy represents an excellent laboratory to study the interplay between different gas phases and star formation, and the influence of gas stripping, gas heating, and AGN. We analyze the physical origin of the emission at different wavelengths in the tail, in particular in-situ star formation (related to Halpha, CO and UV emission), synchrotron emission from relativistic electrons (producing the radio continuum) and heating of the stripped interstellar medium (ISM) (responsible for the X-ray emission). We show the similarities and differences of the spatial distributions of ionized gas, molecular gas and UV light, and argue that the mismatch on small scales (1kpc) is due to different stages of the star formation process. We present the relation Halpha--X-ray surface brightness, which is steeper for star-forming regions than for diffuse ionised gas regions with high [OI]/Halpha ratio. We propose that ISM heating due to interaction with the intracluster medium (either for mixing, thermal conduction or shocks) is responsible for the X-ray tail, the observed [OI]-excess and the lack of star formation in the northern part of the tail. We also report the tentative discovery in the tail of the most distant (and among the brightest) currently known ULX, a point-like ultraluminous X-ray source commonly originating in a binary stellar system powered either by an intermediate-mass black hole or a magnetized neutron star., Comment: accepted for publication in ApJ

Stacked analyses of galaxy clusters at low-to-intermediate redshift show signatures attributable to dust, but the origin of this dust is uncertain. We test the hypothesis that the bulk of cluster dust derives from galaxy ejecta. To do so, we employ dust abundances obtained from detailed chemical evolution models of galaxies. We integrate the dust abundances over cluster luminosity functions (one-slope and two-slope Schechter functions). We consider both a hierarchical scenario of galaxy formation and an independent evolution of the three main galactic morphologies: elliptical/S0, spiral and irregular. We separate the dust residing within galaxies from the dust ejected in the intracluster medium. To the latter, we apply thermal sputtering. The model results are compared to low-to-intermediate redshift observations of dust masses. We find that in any of the considered scenarios, elliptical/S0 galaxies contribute negligibly to the present-time intracluster dust, despite producing the majority of gas-phase metals in galaxy clusters. Spiral galaxies, instead, provide both the bulk of the spatially-unresolved dust and of the dust ejected into the intracluster medium. The total dust-to-gas mass ratio in galaxy clusters amounts to $10^{-4}$, while the intracluster medium dust-to-gas mass ratio amounts to $10^{-6}$ at most. These dust abundances are consistent with the estimates of cluster observations at $0.2 < z <1$. We propose that galactic sources, spiral galaxies in particular, are the major contributors to the cluster dust budget., Comment: MNRAS. Accepted on Feb 11th, 2020. 13 pages, 4 figures

We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 \mu$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we derive a shallow faint end slope of $ \alpha \sim -0.35^{+0.15}_{-0.15} $ and $ M_{H}^{*} \sim -23.52^{+0.15}_{-0.17} $, suggesting that there is a deficit of faint red sequence galaxies in clusters at high redshift. By comparing the stacked red sequence LF of our sample with a sample of clusters at z~0.6, we find an evolution in the faint end of the red sequence over the ~2.6 Gyr between the two samples, with the mean faint end red sequence luminosity growing by more than a factor of two. The faint-to-luminous ratio of our sample ($0.78^{+0.19}_{-0.15}$) is consistent with the trend of decreasing ratio with increasing redshift as proposed in previous studies. A comparison with the field shows that the faint-to-luminous ratios in clusters are consistent with the field at z~1.15 and exhibit a stronger redshift dependence. Our results support the picture that the build up of the faint red sequence galaxies occurs gradually over time and suggest that faint cluster galaxies, similar to bright cluster galaxies, experience the quenching effect induced by environment already at z~1.15., Comment: Accepted for publication in ApJ. 24 pages, 13 figures

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter., Comment: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

The gravitational slip parameter is an important discriminator between large classes of gravity theories at cosmological and astrophysical scales. In this work we use a combination of simulated information of galaxy cluster mass profiles, inferred by Strong+Weak lensing analyses and by the study of the dynamics of the cluster member galaxies, to reconstruct the gravitational slip parameter $\eta$ and predict the accuracy with which it can be constrained with current and future galaxy cluster surveys. Performing a full-likelihood statistical analysis, we show that galaxy cluster observations can constrain $\eta$ down to the percent level already with a few tens of clusters. We discuss the significance of possible systematics, and show that the cluster masses and numbers of galaxy members used to reconstruct the dynamics mass profile have a mild effect on the predicted constraints., Comment: 18 pages plus references, 8 Figures

We study the dependency of the concentration on mass and redshift using three large N-body cosmological hydrodynamic simulations carried out by the Magneticum project. We constrain the slope of the mass-concentration relation with an unprecedented mass range for hydrodynamic simulations and find a negative trend on the mass-concentration plane and a slightly negative redshift dependency, in agreement with observations and other numerical works. We also show how the concentration correlates with the fossil parameter, defined as the stellar mass ratio between the central galaxy and the most massive satellite, in agreement with observations. We find that haloes with high fossil parameter have systematically higher concentration and investigate the cause in two different ways. First we study the evolution of haloes that lives unperturbed for a long period of time, where we find that the internal region keeps accreting satellites as the fossil parameter increases and the scale radius decreases (which increases the concentration). We also study the dependency of the concentration on the virial ratio and the energy term from the surface pressure $E_s$. We conclude that fossil objects have higher concentration because they are dynamically relaxed, with no in-fall/out-fall material and had time to accrete their satellites., Comment: 13 pages, 10 figures

The effect of galactic orbits on a galaxy's internal evolution within a galaxy cluster environment has been the focus of heated debate in recent years. To understand this connection, we use both the $(0.5 \,$Gpc)$^3$ and the Gpc$^3$ boxes from the cosmological hydrodynamical simulation set Magneticum Pathfinder. We investigate the velocity-anisotropy, phase space, and the orbital evolution of up to $\sim 5 \cdot 10^{5}$ resolved satellite galaxies within our sample of 6776 clusters with $M_{\mathrm{vir}} > 10^{14} \, \mathrm{M_{\odot}}$ at low redshift, which we also trace back in time. In agreement with observations, we find that star-forming satellite galaxies inside galaxy clusters are characterised by more radially dominated orbits, independent of cluster mass. Furthermore, the vast majority of star-forming satellite galaxies stop forming stars during their first passage. We find a strong dichotomy both in line-of-sight and radial phase space between star-forming and quiescent galaxies, in line with observations. The tracking of individual orbits shows that the star-formation of almost all satellite galaxies drops to zero within $1 \, \mathrm{Gyr}$ after in-fall. Satellite galaxies that are able to remain star-forming longer are characterised by tangential orbits and high stellar mass. All this indicates that in galaxy clusters the dominant quenching mechanism is ram-pressure stripping., Comment: 22 pages, 16 figures, accepted by MNRAS

It is well known that galaxies falling into clusters can experience gas stripping due to ram-pressure by the intra-cluster medium (ICM). The most spectacular examples are galaxies with extended tails of optically-bright stripped material known as "jellyfish". We use the first large homogeneous compilation of jellyfish galaxies in clusters from the WINGS and OmegaWINGS surveys, and follow-up MUSE observations from the GASP MUSE programme to investigate the orbital histories of jellyfish galaxies in clusters and reconstruct their stripping history through position vs. velocity phase- space diagrams. We construct analytic models to define the regions in phase-space where ram-pressure stripping is at play. We then study the distribution of cluster galaxies in phase-space and find that jellyfish galaxies have on average higher peculiar velocities (and higher cluster velocity dispersion) than the overall population of cluster galaxies at all clustercentric radii, which is indicative of recent infall into the cluster and radial orbits. In particular, the jellyfish galaxies with the longest gas tails reside very near the cluster cores (in projection) and are moving at very high speeds, which coincides with the conditions of the most intense ram-pressure. We conclude that many of the jellyfish galaxies seen in clusters likely formed via fast (~1- 2 Gyr), incremental, outside-in ram-pressure stripping during first infall into the cluster in highly radial orbits., Comment: 14 pages, 9 figures, MNRAS in press

We present the first results of a new spectroscopic survey of the cluster Abell 85 targeting 1466 candidate cluster members within the central $\sim$1 deg$^2$ of the cluster and having magnitudes $m_r < 20.5$ using VIMOS/VLT and HYDRA/WIYN. A total of 520 galaxies are confirmed as either relaxed cluster members or part of an infalling population. A significant fraction are low mass; the median stellar mass of the sample is $10^{9.6} M_{\odot} $, and 25% have stellar masses below $10^9 M_{\odot}$ (i.e. 133 dwarf galaxies). We also identify seven active galactic nuclei (AGN), four of which reside in dwarf host galaxies. We probe the evolution of star formation rates, based on H$\alpha$ emission and continuum modeling, as a function of both mass and environment. We find that more star forming galaxies are observed at larger clustercentric distances, while infalling galaxies show evidence for recently enhanced star forming activity. Main sequence galaxies, defined by their continuum star formation rates, show different evolutionary behavior based on their mass. At the low mass end, the galaxies have had their star formation recently quenched, while more massive galaxies show no significant change. The timescales probed here favor fast quenching mechanisms, such as ram-pressure stripping. Galaxies within the green valley, defined similarly, do not show evidence of quenching. Instead, the low mass galaxies maintain their levels of star forming activity, while the more massive galaxies have experienced a recent burst., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

In the context of the GAs Stripping Phenomena in galaxies with Muse (GASP) survey, we present the characterization of JO204, a jellyfish galaxy in A957, a relatively low-mass cluster with $M=4.4 \times10^{14}M_\odot$. This galaxy shows a tail of ionized gas that extends up to 30 kpc from the main body in the opposite direction of the cluster center. No gas emission is detected in the galaxy outer disk, suggesting that gas stripping is proceeding outside-in. The stellar component is distributed as a regular disk galaxy; the stellar kinematics shows a symmetric rotation curve with a maximum radial velocity of 200km/s out to 20 kpc from the galaxy center. The radial velocity of the gas component in the central part of the disk follows the distribution of the stellar component; the gas kinematics in the tail retains the rotation of the galaxy disk, indicating that JO204 is moving at high speed in the intracluster medium. Both the emission and radial velocity maps of the gas and stellar components indicate ram-pressure as the most likely primary mechanism for gas stripping, as expected given that JO204 is close to the cluster center and it is likely at the first infall in the cluster. The spatially resolved star formation history of JO204 provides evidence that the onset of ram-pressure stripping occurred in the last 500 Myr, quenching the star formation activity in the outer disk, where the gas has been already completely stripped. Our conclusions are supported by a set of hydrodynamic simulations., Comment: accepted for publication in ApJ

[Abridged] We use the WINGS database to select a sample of 67 nearby galaxy clusters with at least 30 spectroscopic members each. 53 of these clusters do not show evidence of substructures in phase-space, while 14 do. We estimate the virial radii and circular velocities of the 67 clusters by a variety of proxies (velocity dispersion, X-ray temperature, and richness) and use these estimates to build stack samples from these 53 and 14 clusters ('Reg' and 'Irr' stacks, respectively). We determine the number-density and velocity-dispersion profiles (VDPs) of E, S0, and Sp+Irr (S) galaxies in the Reg and Irr samples, separately, and fit models to these profiles. The number density profiles of E, S0, and S galaxies are adequately described by either a NFW or a cored King model, both for the Reg and Irr samples, with a slight preference for the NFW model. The spatial distribution concentration increases from the S to the S0 and to the E populations, both in the Reg and the Irr stacks, reflecting the well-known morphology-radius relation. Reg clusters have a more concentrated spatial distribution of E and S0 galaxies than Irr clusters, while the spatial distributions of S galaxies in Reg and Irr clusters are similar. We propose a new phenomenological model that provides acceptable fits to the VDP of all our galaxy samples. The VDPs become steeper and with a higher normalization from E to S0 to S galaxies. The S0 VDP is close to that of E galaxies in Reg clusters, and intermediate between those of E and S galaxies in Irr clusters. Our results suggest that S galaxies are a recently accreted cluster population, that take less than 3 Gyr to evolve into S0 galaxies after accretion, and in doing so modify their phase-space distribution, approaching that of cluster ellipticals. While in Reg clusters this evolutionary process is mostly completed, it is still ongoing in Irr clusters., Comment: A&A, in press - 11 pages, 9 figures, 4 tables

We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $110$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius., Comment: Accepted for publication in MNRAS May 26, 2017

GASP (GAs Stripping Phenomena in galaxies with MUSE) is a new integral-field spectroscopic survey with MUSE at the VLT aiming at studying gas removal processes in galaxies. We present an overview of the survey and show a first example of a galaxy undergoing strong gas stripping. GASP is obtaining deep MUSE data for 114 galaxies at z=0.04-0.07 with stellar masses in the range 10^9.2-10^11.5 M_sun in different environments (galaxy clusters and groups, over more than four orders of magnitude in halo mass). GASP targets galaxies with optical signatures of unilateral debris or tails reminiscent of gas stripping processes ("jellyfish galaxies"), as well as a control sample of disk galaxies with no morphological anomalies. GASP is the only existing Integral Field Unit (IFU) survey covering both the main galaxy body and the outskirts and surroundings, where the IFU data can reveal the presence and the origin of the outer gas. To demonstrate GASP's ability to probe the physics of gas and stars, we show the complete analysis of a textbook case of a "jellyfish" galaxy, JO206. This is a massive galaxy (9 x 10^10 M_sun in a low-mass cluster (sigma ~500 km/s), at a small projected clustercentric radius and a high relative velocity, with >=90kpc-long tentacles of ionized gas stripped away by ram pressure. We present the spatially resolved kinematics and physical properties of gas and stars, and depict the evolutionary history of this galaxy., Comment: accepted for publication in ApJ

The so-called jellyfish galaxies are objects exhibiting disturbed morphology, mostly in the form of tails of gas stripped from the main body of the galaxy. Several works have strongly suggested ram pressure stripping to be the mechanism driving this phenomenon. Here, we focus on one of these objects, drawn from a sample of optically selected jellyfish galaxies, and use it to validate SINOPSIS, the spectral fitting code that will be used for the analysis of the GASP (GAs Stripping Phenomena in galaxies with MUSE) survey, and study the spatial distribution and physical properties of gas and stellar populations in this galaxy. We compare the model spectra to those obtained with GANDALF, a code with similar features widely used to interpret the kinematic of stars and gas in galaxies from IFU data. We find that SINOPSIS can reproduce the pixel-by-pixel spectra of this galaxy at least as good as GANDALF does, providing reliable estimates of the underlying stellar absorption to properly correct the nebular gas emission. Using these results, we find strong evidences of a double effect of ram pressure exerted by the intracluster medium onto the gas of the galaxy. A moderate burst of star formation, dating between 20 and 500 Myr ago and involving the outer parts of the galaxy more strongly than the inner regions, was likely induced by a first interaction of the galaxy with the intracluster medium. Stripping by ram pressure, plus probable gas depletion due to star formation, contributed to create a truncated ionized gas disk. The presence of an extended stellar tail on only one side of the disk, points instead to another kind of process, likely a gravitational interaction by a fly-by or a close encounter with another galaxy in the cluster., Comment: ApJ in press, 26 pages, 18 figures

The Ophiuchus cluster, at a redshift z=0.0296, is known from X-rays to be one of the most massive nearby clusters, but due to its very low Galactic latitude its optical properties have not been investigated in detail. We discuss the optical properties of the galaxies in the Ophiuchus cluster, in particular with the aim of understanding better its dynamical properties. We have obtained deep optical imaging in several bands with various telescopes, and applied a sophisticated method to model and subtract the contributions of stars in order to measure galaxy magnitudes as accurately as possible. The colour-magnitude relations obtained show that there are hardly any blue galaxies in Ophiuchus (at least brighter than r'<=19.5), and this is confirmed by the fact that we only detect two galaxies in Halpha. We also obtained a number of spectra with ESO-FORS2, that we combined with previously available redshifts. Altogether, we have 152 galaxies with spectroscopic redshifts in the 0.02<=z<=0.04 range, and 89 galaxies with both a redshift within the cluster redshift range and a measured r' band magnitude (limited to the Megacam 1x1 deg^2 field). A complete dynamical analysis based on the galaxy redshifts available shows that the overall cluster is relaxed and has a mass of 1.1x10^15 solar masses. The Sernal-Gerbal method detects a main structure and a much smaller substructure that are not separated in projection. From its dynamical properties derived from optical data, the Ophiuchus cluster seems to be overall a relaxed structure, or at most a minor merger, though in X-rays the central region (radius ~ 150 kpc) may show evidence for merging effects., Comment: Accepted for publication in Astronomy & Astrophysics

The next generation of large sky photometric surveys will finally be able to use arc statistics as a cosmological probe. Here we present the first of a series of papers on this topic. In particular, we study how arc counts are sensitive to the variation of two cosmological parameters: the (total) matter density parameter, $\Omega_m$, and the normalisation of the primordial power spectrum, expressed in terms of $\sigma_8$. Both these parameters influence the abundances of collapsed structures and their internal structure. We compute the expected number of gravitational arcs with various length-to-width ratios in mock light cones, by varying these cosmological parameters in the ranges $0.1\leq\Omega_m\leq0.5$ and $0.6\leq\sigma_8\leq 1$. We find that the arc counts dependence on $\Omega_m$ and $\sigma_8$ is similar, but not identical, to that of the halo counts. We investigate how the precision of the constraints on the cosmological parameters based on arc counts depends on the survey area. We find that the constraining power of arc statistics degrades critically only for surveys covering an area smaller than $10\%$ of the whole sky. Finally, we consider the case in which the search for arcs is done only in frames where galaxy clusters have been previously identified. Adopting the selection function for galaxy clusters expected to be detected from photometric data in future wide surveys, we find that less than $10\%$ of the arcs will be missed, with only a small degradation of the corresponding cosmological constraints., Comment: 12 pages, 12 figures, accepted by MNRAS

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., Comment: 23 pages, 13 figures, accepted for publication in A&A. Astro-ph abstract abridged

We compute optical galaxy luminosity functions (GLFs) in the B, V, R, and I rest-frame bands for one of the largest medium-to-high-redshift (0.4 < z < 0.9) cluster samples to date in order to probe the abundance of faint galaxies in clusters. We also study how the GLFs depend on cluster redshift, mass, and substructure, and compare the GLFs of clusters with those of the field. We separately investigate the GLFs of blue and red-sequence (RS) galaxies to understand the evolution of different cluster populations. We find that the shapes of our GLFs are similar for the B, V, R, and I bands with a drop at the red GLF faint end that is more pronounced at high-redshift: alpha(red) ~ -0.5 at 0.40 < z < 0.65 and alpha(red) > 0.1 at 0.65 < z < 0.90. The blue GLFs have a steeper faint end (alpha(blue) ~ -1.6) than the red GLFs, that appears to be independent of redshift. For the full cluster sample, blue and red GLFs intersect at M(V) = -20, M(R) = -20.5, and M(I) = -20.3. A study of how galaxy types evolve with redshift shows that late type galaxies appear to become early types between z ~ 0.9 and today. Finally, the faint ends of the red GLFs of more massive clusters appear to be richer than less massive clusters, which is more typical of the lower redshift behaviour. Our results indicate that our clusters form at redshifts higher than z = 0.9 from galaxy structures that already have an established red sequence. Late type galaxies then appear to evolve into early types, enriching the red-sequence between this redshift and today. This effect is consistent with the evolution of the faint end slope of the red-sequence and the galaxy type evolution that we find. Finally, faint galaxies accreted from the field environment at all redshifts might have replaced the blue late type galaxies that converted into early types, explaining the lack of evolution in the faint end slopes of the blue GLFs., Comment: accepted for publication in A&A

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

We present a joint shear-and-magnification weak-lensing analysis of a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters at 0.19

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2\rho)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program., Comment: 5 pages, 1 figures, accepted for publication in the Astrophysical Journal Letters

Aim: We aim to compute the mass and velocity anisotropy profiles of Abell 2142 and, from there, the pseudo phase--space density profile $Q(r)$ and the density slope - velocity anisotropy $\beta - \gamma$ relation, and then to compare them with theoretical expectations. Methods: The mass profiles were obtained by using three techniques based on member galaxy kinematics, namely the caustic method, the method of Dispersion - Kurtosis, and MAMPOSSt. Through the inversion of the Jeans equation, it was possible to compute the velocity anisotropy profiles. Results: The mass profiles, as well as the virial values of mass and radius, computed with the different techniques agree with one another and with the estimates coming from X-ray and weak lensing studies. A concordance mass profile is obtained by averaging the lensing, X-ray, and kinematics determinations. The cluster mass profile is well fitted by an NFW profile with $c=4.0 \pm 0.5$. The population of red and blue galaxies appear to have a different velocity anisotropy configuration, since red galaxies are almost isotropic, while blue galaxies are radially anisotropic, with a weak dependence on radius. The $Q(r)$ profile for the red galaxy population agrees with the theoretical results found in cosmological simulations, suggesting that any bias, relative to the dark matter particles, in velocity dispersion of the red component is independent of radius. The $\beta - \gamma$ relation for red galaxies matches the theoretical relation only in the inner region. The deviations might be due to the use of galaxies as tracers of the gravitational potential, unlike the non--collisional tracer used in the theoretical relation., Comment: 14 pages, 14 figures. Consolidated version including the Corrigendum published on A&A

We provide a new observational test for a key prediction of the \Lambda CDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $\sigma$ agreement with \Lambda CDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated., Comment: 25 pages, 16 figures, 6 tables, 2 machine readable tables, accepted for publication in ApJ, updated acknowledgements and data table format modifications made

In the local Universe, galaxy properties show a strong dependence on environment. In cluster cores, early type galaxies dominate, whereas star-forming galaxies are more and more common in the outskirts. At higher redshifts and in somewhat less dense environments (e.g. galaxy groups), the situation is less clear. One open issue is that of whether and how the star formation rate (SFR) of galaxies in groups depends on the distance from the centre of mass. To shed light on this topic, we have built a sample of X-ray selected galaxy groups at 010^10.3 M_sun in order to have a high spectroscopic completeness. As we use only spectroscopic redshifts, our results are not affected by uncertainties due to projection effects. We use several SFR indicators to link the star formation (SF) activity to the galaxy environment. Taking advantage of the extremely deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations, we have an accurate, broad-band measure of the SFR for the bulk of the star-forming galaxies. We use multi-wavelength SED fitting techniques to estimate the stellar masses of all objects and the SFR of the MIPS and PACS undetected galaxies. We analyse the dependence of the SF activity, stellar mass and specific SFR on the group-centric distance, up to z~1.6, for the first time. We do not find any correlation between the mean SFR and group-centric distance at any redshift. We do not observe any strong mass segregation either, in agreement with predictions from simulations. Our results suggest that either groups have a much smaller spread in accretion times with respect to the clusters and that the relaxation time is longer than the group crossing time., Comment: Accepted for publication in MNRAS

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5\sigma_\kappa. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2{\sigma} departure above the maximum simulated halo mass at this redshift., Comment: 24 pages, 16 pages, 6 tables; matches version accepted for publication in ApJ

[Abridged] We present an analysis of the relation between the masses of cluster- and group-sized halos, extracted from $\Lambda$CDM cosmological N-body and hydrodynamic simulations, and their velocity dispersions, at different redshifts from $z=2$ to $z=0$. The main aim of this analysis is to understand how the implementation of baryonic physics in simulations affects such relation, i.e. to what extent the use of the velocity dispersion as a proxy for cluster mass determination is hampered by the imperfect knowledge of the baryonic physics. In our analysis we use several sets of simulations with different physics implemented. Velocity dispersions are determined using three different tracers, DM particles, subhalos, and galaxies. We confirm that DM particles trace a relation that is fully consistent with the theoretical expectations based on the virial theorem and with previous results presented in the literature. On the other hand, subhalos and galaxies trace steeper relations, and with larger values of the normalization. Such relations imply that galaxies and subhalos have a $\sim10$ per cent velocity bias relative to the DM particles, which can be either positive or negative, depending on halo mass, redshift and physics implemented in the simulation. We explain these differences as due to dynamical processes, namely dynamical friction and tidal disruption, acting on substructures and galaxies, but not on DM particles. These processes appear to be more or less effective, depending on the halo masses and the importance of baryon cooling, and may create a non-trivial dependence of the velocity bias and the $\soneD$--$\Mtwo$ relation on the tracer, the halo mass and its redshift. These results are relevant in view of the application of velocity dispersion as a proxy for cluster masses in ongoing and future large redshift surveys., Comment: 13 pages, 16 figures. Minor modifications to match the version in press on MNRAS

Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and thoroughly tested. MAMPOSSt assumes a shape for the gravitational potential, but instead of postulating a shape for the distribution function in terms of energy and angular momentum, or supposing Gaussian line-of-sight velocity distributions, MAMPOSSt assumes a VA profile and a shape for the 3D velocity distribution, here Gaussian. MAMPOSSt requires no binning, differentiation, nor extrapolation of the observables. Tests on cluster-mass haloes from LambdaCDM cosmological simulations show that, with 500 tracers, MAMPOSSt is able to jointly recover the virial radius, tracer scale radius, dark matter scale radius and outer or constant VA with small bias (<10% on scale radii and <2% on the two other quantities) and inefficiencies of 10%, 27%, 48% and 20%, respectively. MAMPOSSt does not perform better when some parameters are frozen, and even worse when the virial radius is set to its true value, which appears to be the consequence of halo triaxiality. The accuracy of MAMPOSSt depends weakly on the adopted interloper removal scheme, including an efficient iterative Bayesian scheme that we introduce here, which can directly obtain the virial radius with as good precision as MAMPOSSt. Our tests show that MAMPOSSt with Gaussian 3D velocities is very competitive with, and up to 1000x faster than other methods. Hence, MAMPOSSt is a very powerful and rapid tool for the mass and anisotropy modeling of systems such as clusters and groups of galaxies, elliptical and dwarf spheroidal galaxies., Comment: MNRAS in press (fixes to eqs [20] and [30] plus minor touch-ups). 21 pages, 10 figures, 6 Tables

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters., Comment: Accepted by ApJ (30 pages, 17 figures), one new figure (Figure 10) added, minor text changes; a version with high resolution figures available at http://www.asiaa.sinica.edu.tw/~keiichi/upfiles/MACS1206/ms_highreso.pdf