No Arabic abstract
Deep observations of nearby galaxy clusters with Chandra have revealed concave bay structures in a number of systems (Perseus, Centaurus and Abell 1795), which have similar X-ray and radio properties. These bays have all the properties of cold fronts, where the temperature rises and density falls sharply, but are concave rather than convex. By comparing to simulations of gas sloshing, we find that the bay in the Perseus cluster bears a striking resemblance in its size, location and thermal structure, to a giant ($approx$50 kpc) roll resulting from Kelvin-Helmholtz instabilities. If true, the morphology of this structure can be compared to simulations to put constraints on the initial average ratio of the thermal and magnetic pressure, $beta= p_{rm th} / p_{rm B}$, throughout the overall cluster before the sloshing occurs, for which we find $beta=200$ to best match the observations. Simulations with a stronger magnetic field ($beta=100$) are disfavoured, as in these the large Kelvin-Helmholtz rolls do not form, while in simulations with a lower magnetic field ($beta=500$) the level of instabilities is much larger than is observed. We find that the bay structures in Centaurus and Abell 1795 may also be explained by such features of gas sloshing.
Sloshing cold fronts (CFs) arise from minor merger triggered gas sloshing. Their detailed structure depends on the properties of the intra-cluster medium (ICM): hydrodynamical simulations predict the CFs to be distorted by Kelvin-Helmholtz instabilities (KHIs), but aligned magnetic fields, viscosity, or thermal conduction can suppress the KHIs. Thus, observing the detailed structure of sloshing CFs can be used to constrain these ICM properties. Both smooth and distorted sloshing CFs have been observed, indicating that the KHI is suppressed in some clusters, but not in all. Consequently, we need to address at least some sloshing clusters individually before drawing general conclusions about the ICM properties. We present the first detailed attempt to constrain the ICM properties in a specific cluster from the structure of its sloshing CF. Proximity and brightness make the Virgo cluster an ideal target. We combine observations and Virgo-specific hydrodynamical sloshing simulations. Here we focus on a Spitzer-like temperature dependent viscosity as a mechanism to suppress the KHI, but discuss the alternative mechanisms in detail. We identify the CF at 90 kpc north and north-east of the Virgo center as the best location in the cluster to observe a possible KHI suppression. For viscosities $gtrsim$ 10% of the Spitzer value KHIs at this CF are suppressed. We describe in detail the observable signatures at low and high viscosities, i.e. in the presence or absence of KHIs. We find indications for a low ICM viscosity in archival XMM-Newton data and demonstrate the detectability of the predicted features in deep Chandra observations.
We investigate the origin and nature of the multiple sloshing cold fronts in the core of Abell 496 by direct comparison between observations and dedicated hydrodynamical simulations. Our simulations model a minor merger with a 4{times}10^13M{circ} subcluster crossing A496 from the south-west to the north-north-east, passing the cluster core in the south-east at a pericentre distance 100 to a few 100 kpc about 0.6 to 0.8 Gyr ago. The gas sloshing triggered by the merger can reproduce almost all observed features, e.g. the characteristic spiral-like brightness residual distribution in the cluster centre and its asymmetry out to 500 kpc, also the positions of and contrasts across the cold fronts. If the subcluster passes close (100 kpc) to the cluster core, the resulting shear flows are strong enough to trigger Kelvin-Helmholtz instabilities that in projection resemble the peculiar kinks in the cold fronts of Abell 496. Finally, we show that sloshing does not lead to a significant modification of the global ICM profiles but a mild oscillation around the initial profiles.
Sloshing cold fronts in clusters, produced as the dense cluster core moves around in the cluster potential in response to in-falling subgroups, provide a powerful probe of the physics of the intracluster medium (ICM), and the magnetic fields permeating it. These sharp discontinuities in density and temperature rise gradually outwards with age in a characteristic spiral pattern, embedding into the intracluster medium a record of the minor merging activity of clusters: the further from the cluster centre a cold front is, the older it is. Recently it has been discovered that these cold fronts can survive out to extremely large radii in the Perseus cluster. Here we report on high spatial resolution Chandra observations of the large scale cold front in Perseus. We find that rather than broadening through diffusion, the cold front remains extremely sharp (consistent with abrupt jumps in density) but instead is split into two sharp edges. These results show that magnetic draping can suppress diffusion for vast periods of time, around ~5 Gyr, even as the cold front expands out to nearly half the cluster virial radius.
There has been interest in recent years to assess the ability of astrophysical hydrodynamics codes to correctly model the Kelvin-Helmholtz instability. Smoothed particle hydrodynamics (SPH), in particular, has received significant attention, though there has yet to be a clear demonstration that SPH yields converged solutions that are in agreement with other methods. We have performed SPH simulations of the Kelvin-Helmholtz instability using the test problem put forward by Lecoanet et al (2016). We demonstrate that the SPH solutions converge to the reference solution in both the linear and non-linear regimes. Quantitative convergence in the strongly non-linear regime is achieved by using a physical Navier-Stokes viscosity and thermal conductivity. We conclude that standard SPH with an artificial viscosity can correctly capture the Kelvin-Helmholtz instability.
We perform simulations of the Kelvin-Helmholtz instability using smoothed particle hydrodynamics (SPH). The instability is studied both in the linear and strongly non-linear regimes. The smooth, well-posed initial conditions of Lecoanet et al. (2016) are used, along with an explicit Navier-Stokes viscosity and thermal conductivity to enforce the evolution in the non-linear regime. We demonstrate convergence to the reference solution using SPH. The evolution of the vortex structures and the degree of mixing, as measured by a passive scalar `colour field, match the reference solution. Tests with an initial density contrast produce the correct qualitative behaviour. The L2 error of the SPH calculations decreases as the resolution is increased. The primary source of error is numerical dissipation arising from artificial viscosity, and tests with reduced artificial viscosity have reduced L2 error. A high-order smoothing kernel is needed in order to resolve the initial velocity amplitude of the seeded mode and inhibit excitation of spurious modes. We find that standard SPH with an artificial viscosity has no difficulty in correctly modelling the Kelvin-Helmholtz instability and yields convergent solutions.