The disruption of substructure in galaxy clusters likely plays an important role in shaping the cluster population as a significant fraction of cluster galaxies today have spent time in a previous host system, and thus may have been pre-processed. Once inside the cluster, group galaxies face the combined environmental effects from group and cluster - so called 'post-processing'. We investigate these concepts, by tracking the evolution of satellites and their hosts after entering the cluster and find that tidal forces during their first pericentric passage are very efficient at breaking up groups, preferentially removing satellites at larger distances from their hosts. 92.2% of satellites whose host has passed pericentre will leave their host by $z=0$, typically no later than half a Gyr after pericentric passage. We find satellites leave with high velocities, and quickly separate to large distances from their hosts, making their identification within the cluster population challenging. Those few satellites ($\sim$7.8%) that remain bound to their hosts after a pericentric passage are typically found close to their host centres. This implies that substructure seen in clusters today is very likely on first infall into the cluster, and yet to pass pericentre. This is even more likely if the substructure is extended, with satellites beyond R$_{200}$ of their host. We find the cluster dominates the tidal mass loss and destruction of satellites, and is responsible for rapidly halting the accretion of new satellites onto hosts once they reach 0.6-0.8 R$_{200}$ radii from the cluster.

Galaxies within galaxy clusters are known to be subject to a wide variety of environmental effects, both gravitational and hydrodynamical. In this study, we examine the purely gravitational interaction of idealised galaxy models falling into a galaxy cluster in the context of Modified Newtonian Dynamics (MOND). This modification of gravity gives rise to an external field effect (EFE), where the internal dynamics of a system are affected by the presence of external gravitational fields. We examine the consequences of the EFE on low and high mass disk galaxies in time-evolving analytic background cluster potentials, considering orbits with weak and strong tidal fields. By varying the orbital plane of the galaxies we also test the effect of having the tidal interaction orthogonal or parallel to the disk. Furthermore, we consider as a control sample models where the EFE has been removed and they are only affected by tides. Our results suggest that MOND cluster galaxies should exhibit clear asymmetries in their isophotes, suffer increased mass loss and a reduction in their rotation curves due to the combined effect of cluster tides and the external field. In particular, low mass galaxies are hit hard by the EFE, becoming dominated by dispersion rather than rotation even in the absence of tides.

Studies of cluster galaxies are increasingly finding galaxies with spectacular one-sided tails of gas and young stars, suggestive of intense ram-pressure stripping. These so-called "jellyfish" galaxies typically have late-type morphology. In this paper, we present MUSE observations of an elliptical galaxy in Abell 2670 with long tails of material visible in the optical spectra, as well as blobs with tadpole-like morphology. The spectra in the central part of the galaxy reveals a stellar component as well as ionized gas. The stellar component does not have significant rotation, while the ionized gas defines a clear star-forming gas disk. We argue, based on deep optical images of the galaxy, that the gas was most likely acquired during a past wet merger. It is possible that the star-forming blobs are also remnants of the merger. In addition, the direction and kinematics of the one-sided ionized tails, combined with the tadpole morphology of the star-forming blobs, strongly suggests that the system is undergoing ram pressure from the intracluster medium. In summary, this paper presents the discovery of a post-merger elliptical galaxy undergoing ram pressure 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.

The recently developed code for N-body/hydrodynamics simulations in Modified Newtonian Dynamics (MOND), known as RAyMOND, is used to investigate the consequences of MOND on structure formation in a cosmological context, with a particular focus on the velocity field. This preliminary study investigates the results obtained with the two formulations of MOND implemented in RAyMOND, as well as considering the effects of changing the choice of MOND interpolation function, and the cosmological evolution of the MOND acceleration scale. The simulations are contrived such that structure forms in a background cosmology that is similar to $\Lambda$CDM, but with a significantly lower matter content. Given this, and the fact that a fully consistent MOND cosmology is still lacking, we compare our results with a standard $\Lambda$CDM simulation, rather than observations. As well as demonstrating the effectiveness of using RAyMOND for cosmological simulations, it is shown that a significant enhancement of the velocity field is likely an unavoidable consequence of the gravitational modification implemented in MOND, and may represent a clear observational signature of such a modification. It is further suggested that such a signal may be clearest in intermediate density regions such as cluster outskirts and filaments.

We attempt to find a progenitor for the ultra-faint object Segue 1 under the assumption that it formed as a dark matter free star cluster in the past. We look for orbits, using the elongation of Segue 1 on the sky as a tracer of its path. Those orbits are followed backwards in time to find the starting points of our N-body simulations. The successful orbit, with which we can reproduce Segue 1 has a proper motion of mu_alpha = -0.19 mas/yr and mu_delta = -1.9 mas/yr, placing Segue 1 near its apo-galacticon today. Our best fitting model has an initial mass of 6224 Msun and an initial scale-length of 5.75 pc.

Using three-dimensional stellar kinematic data from simulated galaxies, we examine the efficacy of a Jeans equation analysis in reconstructing the total disk surface density, including the dark matter, at the "Solar" radius. Our simulation dataset includes galaxies formed in a cosmological context using state-of-the-art high resolution cosmological zoom simulations, and other idealised models. The cosmologically formed galaxies have been demonstrated to lie on many of the observed scaling relations for late-type spirals, and thus offer an interesting surrogate for real galaxies with the obvious advantage that all the kinematical data are known perfectly. We show that the vertical velocity dispersion is typically the dominant kinematic quantity in the analysis, and that the traditional method of using only the vertical force is reasonably effective at low heights above the disk plane. At higher heights the inclusion of the radial force becomes increasingly important. We also show that the method is sensitive to uncertainties in the measured disk parameters, particularly the scale lengths of the assumed double exponential density distribution, and the scale length of the radial velocity dispersion. In addition, we show that disk structure and low number statistics can lead to significant errors in the calculated surface densities. Finally we examine the implications of our results for previous studies of this sort, suggesting that more accurate measurements of the scale lengths may help reconcile conflicting estimates of the local dark matter density in the literature.

We present MUSE observations of AM1354-250, confirming its status as a collisional ring galaxy which has recently undergone an interaction, creating its distinctive shape. We analyse the stellar and gaseous emission throughout the galaxy finding direct evidence that the gaseous ring is expanding with a velocity of $\sim$70km.s$^{-1}$ and that star formation is occurring primarily in HII regions associated with the ring. This star formation activity is likely triggered by the interaction. We find evidence for several excitation mechanisms in the gas, including emission consistent with shocked gas in the expanding ring and a region of LINER-like emission in the central core of the galaxy. Evidence of kinematic disturbance in both the stars and gas, possibly also triggered by the interaction, can be seen in all of the velocity maps. The ring galaxy retains weak spiral structure, strongly suggesting the progenitor galaxy was a massive spiral prior to the collision with its companion an estimated $140 \pm 12$ Myr ago.

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.

We examine the effects of gas-expulsion on initially substructured distributions of stars. We perform N-body simulations of the evolution of these distributions in a static background potential to mimic the gas. We remove the static potential instantaneously to model gas-expulsion. We find that the exact dynamical state of the cluster plays a very strong role in affecting a cluster's survival, especially at early times: they may be entirely destroyed or only weakly affected. We show that knowing both detailed dynamics and relative star-gas distributions can provide a good estimate of the post-gas expulsion state of the cluster, but even knowing these is not an absolute way of determining the survival or otherwise of the cluster.

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.

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.

The LCDM concordance cosmological model is supported by a wealth of observational evidence, particularly on large scales. At galactic scales, however, the model is poorly constrained and recent observations suggest a more complex behaviour in the dark sector than may be accommodated by a single cold dark matter component. Furthermore, a modification of the gravitational force in the very weak field regime may account for at least some of the phenomenology of dark matter. A well-known example of such an approach is MOdified Newtonian Dynamics (MOND). While this idea has proven remarkably successful in the context of stellar dynamics in individual galaxies, the effects of such a modification of gravity on galaxy interactions and environmental processes deserves further study. To explore this arena we modify the parallel adaptive mesh refinement code RAMSES to use two formulations of MOND. We implement both the fully non-linear aquadratic Lagrangian (AQUAL) formulation as well as the simpler quasi-linear formulation (QUMOND). The relevant modifications necessary for the Poisson solver in RAMSES are discussed in detail. Using idealised tests, in both serial and parallel runs, we demonstrate the effectiveness of the code.

Hercules is a dwarf spheroidal satellite of the Milky Way, found at a distance of about 138 kpc, and showing evidence of tidal disruption. It is very elongated and exhibits a velocity gradient of 16 +/- 3 km/s/kpc. Using this data a possible orbit of Hercules has previously been deduced in the literature. In this study we make use of a novel approach to find a best fit model that follows the published orbit. Instead of using trial and error, we use a systematic approach in order to find a model that fits multiple observables simultaneously. As such, we investigate a much wider parameter range of initial conditions and ensure we have found the best match possible. Using a dark matter free progenitor that undergoes tidal disruption, our best-fit model can simultaneously match the observed luminosity, central surface brightness, effective radius, velocity dispersion, and velocity gradient of Hercules. However, we find it is impossible to reproduce the observed elongation and the position angle of Hercules at the same time in our models. This failure persists even when we vary the duration of the simulation significantly, and consider a more cuspy density distribution for the progenitor. We discuss how this suggests that the published orbit of Hercules is very likely to be incorrect.

As star clusters are expected to form with low star formation efficiencies, the gas in the cluster is expelled quickly and early in their development: the star cluster "pops." This leads to an unbound stellar system, evolving in the Galactic potential. Previous N-body simulations have demonstrated the existence of a stepped number density distribution of cluster stars after popping, both in vertical position and vertical velocity, with a passing resemblance to a Christmas tree. Using numerical and analytical methods, we investigate the source of this structure, which arises due to the phase mixing of the out-of-equilibrium stellar system, determined entirely by the background analytic potential. Considering only the vertical motions, we construct a theoretical model to describe the time evolution of the phase space distribution of stars in a Miyamoto-Nagai disk potential and a full Milky-Way type potential comprising bulge, halo and disk components, which is then compared with N-body simulations. Using our theoretical model, we investigate the possible observational signatures and the feasibility of detection.

The formation mechanism of tidal dwarf galaxies means they are expected to contain little or no dark matter. As such, they might be expected to be very sensitive to their environment. We investigate the impact of ram pressure on tidal dwarf galaxies in a parameter study, varying dwarf galaxy properties and ram pressures. We submit model tidal dwarf galaxies to wind-tunnel style tests using a toy ram pressure model. The effects of ram pressure are found to be substantial. If tidal dwarf galaxies have their gas stripped, they may be completely destroyed. Ram pressure drag causes acceleration of our dwarf galaxy models, and this further enhances stellar losses. The dragging can also cause stars to lie in a low surface brightness stellar stream that points in the opposite direction to the stripped gas, in a manner distinctive from tidal streams. We investigate the effects of ram pressure on surface density profiles, the dynamics of the stars, and discuss the consequences for dynamical mass measurements.

Recent deep photometry of the dwarf spheroidal Ursa Major II's morphology, and spectroscopy of individual stars, have provided a number of new constraints on its properties. With a velocity dispersion $\sim$6 km s$^{-1}$, and under the assumption that the galaxy is virialised, the mass-to-light ratio is found to be approaching $\sim$2000 - apparently heavily dark matter dominated. Using N-Body simulations, we demonstrate that the observed luminosity, ellipticity, irregular morphology, velocity gradient, and the velocity dispersion can be well reproduced through processes associated with tidal mass loss, and in the absence of dark matter. These results highlight the considerable uncertainty that exists in measurements of the dark matter content of Ursa Major II. The dynamics of the inner tidal tails, and tidal stream, causes the observed velocity dispersion of stars to be boosted to values of $>$5 km s$^{-1}$ ($>$20 km s$^{-1}$ at times). This effect is responsible for raising the velocity dispersion of our model to the observed values in UMaII. We test an iterative rejection technique for removing unbound stars from samples of UMaII stars whose positions on the sky, and line-of-sight velocities, are provided. We find this technique is very effective at providing an accurate bound mass from this information, and only fails when the galaxy has a bound mass less than 10$%$ of its initial mass. However when $<2%$ mass remains bound, mass overestimation by $>$3 orders of magnitude are seen. Additionally we find that mass measurements are sensitive to measurement uncertainty in line-of-sight velocities. Measurement uncertainties of 1-4 km s$^{-1}$ result in mass overestimates by a factor of $\sim$1.3-5.7.