We investigate the effect of ram-pressure from the intracluster medium on the stripping of HI gas in galaxies in a massive, relaxed, X-ray bright, galaxy cluster at z = 0.2 from the Blind Ultra Deep HI Environmental Survey (BUDHIES). We use cosmological simulations, and velocity vs. position phase-space diagrams to infer the orbital histories of the cluster galaxies. In particular, we embed a simple analytical description of ram-pressure stripping in the simulations to identify the regions in phase-space where galaxies are more likely to have been sufficiently stripped of their HI gas to fall below the detection limit of our survey. We find a striking agreement between the model predictions and the observed location of HI-detected and non-detected blue (late-type) galaxies in phase-space, strongly implying that ram-pressure plays a key role in the gas removal from galaxies, and that this can happen during their first infall into the cluster. However, we also find a significant number of gas-poor, red (early-type) galaxies in the infall region of the cluster that cannot easily be explained with our model of ram-pressure stripping alone. We discuss different possible additional mechanisms that could be at play, including the pre-processing of galaxies in their previous environment. Our results are strengthened by the distribution of galaxy colours (optical and UV) in phase-space, that suggests that after a (gas-rich) field galaxy falls into the cluster, it will lose its gas via ram-pressure stripping, and as it settles into the cluster, its star formation will decay until it is completely quenched. Finally, this work demonstrates the utility of phase-space diagrams to analyze the physical processes driving the evolution of cluster galaxies, in particular HI gas stripping.

Using the latest cosmological hydrodynamic N-body simulations of groups and clusters, we study how location in phase-space coordinates at z=0 can provide information on environmental effects acting in clusters. We confirm the results of previous authors showing that galaxies tend to follow a typical path in phase-space as they settle into the cluster potential. As such, different regions of phase-space can be associated with different times since first infalling into the cluster. However, in addition, we see a clear trend between total mass loss due to cluster tides, and time since infall. Thus we find location in phase-space provides information on both infall time, and tidal mass loss. We find the predictive power of phase-space diagrams remains even when projected quantities are used (i.e. line-of-sight velocities, and projected distances from the cluster). We provide figures that can be directly compared with observed samples of cluster galaxies and we also provide the data used to make them as supplementary data, in order to encourage the use of phase-space diagrams as a tool to understand cluster environmental effects. We find that our results depend very weakly on galaxy mass or host mass, so the predictions in our phase-space diagrams can be applied to groups or clusters alike, or to galaxy populations from dwarfs up to giants.

We conduct a comprehensive numerical study of the orbital dependence of harassment on early-type dwarfs consisting of 168 different orbits within a realistic, Virgo-like cluster, varying in eccentricity and pericentre distance. We find harassment is only effective at stripping stars or truncating their stellar disks for orbits that enter deep into the cluster core. Comparing to the orbital distribution in cosmological simulations, we find that the majority of the orbits (more than three quarters) result in no stellar mass loss. We also study the effects on the radial profiles of the globular cluster systems of early-type dwarfs. We find these are significantly altered only if harassment is very strong. This suggests that perhaps most early-type dwarfs in clusters such as Virgo have not suffered any tidal stripping of stars or globular clusters due to harassment, as these components are safely embedded deep within their dark matter halo. We demonstrate that this result is actually consistent with an earlier study of harassment of dwarf galaxies, despite the apparent contradiction. Those few dwarf models that do suffer stellar stripping are found out to the virial radius of the cluster at redshift=0, which mixes them in with less strongly harassed galaxies. However when placed on phase-space diagrams, strongly harassed galaxies are found offset to lower velocities compared to weakly harassed galaxies. This remains true in a cosmological simulation, even when halos have a wide range of masses and concentrations. Thus phase-space diagrams may be a useful tool for determining the relative likelihood that galaxies have been strongly or weakly harassed.

The dynamical stability of disk galaxies is sensitive to whether their anomalous rotation curves are caused by dark matter halos or Milgromian dynamics (MOND). We investigate this by setting up a MOND model of M33. We first simulate it in isolation for 6 Gyr, starting from an initial good match to the rotation curve (RC). Too large a bar and bulge form when the gas is too hot, but this is avoided by reducing the gas temperature. A strong bar still forms in 1 Gyr, but rapidly weakens and becomes consistent with the observed weak bar. Previous work showed this to be challenging in Newtonian models with a live dark matter halo, which developed strong bars. The bar pattern speed implies a realistic corotation radius of 3 kpc. However, the RC still rises too steeply, and the central line-of-sight velocity dispersion (LOSVD) is too high. We then add a constant external acceleration field of 8.4 × 10−12 m s−2 at 30° to the disk as a first-order estimate for the gravity exerted by M31. This suppresses buildup of material at the center, causing the RC to rise more slowly and reducing the central LOSVD. Overall, this simulation bears good resemblance to several global properties of M33, and highlights the importance of including even a weak external field on the stability and evolution of disk galaxies. Further simulations with a time-varying external field, modeling the full orbit of M33, will be needed to confirm its resemblance to observations.

We present accurate models of the gravitational potential produced by a radially exponential disk mass distribution. The models are produced by combining three separate Miyamoto-Nagai disks. Such models have been used previously to model the disk of the Milky Way, but here we extend this framework to allow its application to disks of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disk treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disk by <0.4 % out to 4 disk scalelengths, and <1.9 % out to 10 disk scalelengths. We tabulate fitting parameters which facilitate construction of exponential disks for any scalelength, and a wide range of disk thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disk galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disk, and a disky dwarf galaxy.