In the outskirts of the intracluster medium (ICM) in galaxy clusters, the temperature decreases with radius. Due to the weakly collisional nature of the plasma, these regions are susceptible to the magneto-thermal instability (MTI), which can sustain turbulence and provide turbulent pressure support in the ICM. This instability arises due to heat conduction directed along the magnetic field, with a heat conductivity which is normally assumed to be given by the Spitzer value. Recent numerical studies of the ion mirror and the electron whistler instability using particle-in-cell codes have shown that microscale instabilities can lead to a reduced value for the heat conductivity in the ICM. This could in turn influence the efficiency with which the MTI drives turbulence. In this paper we investigate the influence of reduced heat transport on the nonlinear evolution of the MTI by employing a subgrid model that mimics the influence of the mirror instability on the heat conductivity. We use this subgrid model to assess the effect of microscales on the large scale dynamics of the ICM. We find that the nonlinear saturation of the MTI is surprisingly robust. Over a factor of $\sim 10^3$ the thermal-to-magnetic pressure ratio and collisionality we find at most modest changes to the saturation of the MTI with respect to reference simulations where heat transport is unsuppressed.

The assumption of a spatially uniform helium distribution in the intracluster medium can lead to biases in the estimates of key cluster parameters if composition gradients are present. The helium concentration profile in galaxy clusters is unfortunately not directly observable. Current models addressing the putative sedimentation are one-dimensional and parametrize the presence of magnetic fields in a crude way, ignoring the weakly-collisional, magnetized nature of the medium. When these effects are considered, a wide variety of instabilities can play an important role in the plasma dynamics. In a series of recent papers, we have developed the local, linear theory of these instabilities and addressed their non-linear development with a modified version of Athena. Here, we extend our study by developing a quasi-global approach that we use to simulate the mixing of helium as induced by generalizations of the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI), which feed off thermal and composition gradients. In the inner region of the ICM, mixing can occur on few Gyrs, after which the average magnetic field inclination angle is $\sim 45^{\circ}$ resulting in an averaged Spitzer parameter higher by about 20 % than the value obtained in homogeneous simulations. In the cluster outskirts the instabilities are rather inefficient, due to the shallow gradients. This suggests that compositions gradients in cluster cores might be shallower than one-dimensional models predict. More quantitative statements demand more refined models that can incorporate the physics driving the sedimentation process and simultaneously account for the weakly-collisional nature of the plasma.

The distribution of Helium in the intracluster medium (ICM) permeating galaxy clusters is not well constrained due to the very high plasma temperature. Therefore, the plasma is often assumed to be homogeneous. A non-uniform Helium distribution can however lead to biases when measuring key cluster parameters. This has motivated one-dimensional models that evolve the ICM composition assuming that the effects of magnetic fields can be parameterized or ignored. Such models for non-isothermal clusters show that Helium can sediment in the cluster core leading to a peak in concentration offset from the cluster center. The resulting profiles have recently been shown to be linearly unstable when the weakly-collisional character of the magnetized plasma is considered. In this paper, we present a modified version of the MHD code Athena, which makes it possible to evolve a weakly-collisional plasma subject to a gravitational field and stratified in both temperature and composition. We thoroughly test our implementation and confirm excellent agreement against several analytical results. In order to isolate the effects of composition, in this initial study we focus our attention on isothermal plasmas. We show that plasma instabilities, feeding off gradients in composition, can induce turbulent mixing and saturate by re-arranging magnetic field lines and alleviating the composition gradient. Composition profiles that increase with radius lead to instabilities that saturate by driving the average magnetic field inclination to roughly $45^{\circ}$. We speculate that this effect may alleviate the core insulation observed in homogeneous settings, with potential consequences for the associated cooling flow problem.

High redshift, massive halos are observed to have sustained, high star formation rates, which require that the amount of cold gas in the halo is continuously replenished. The cooling time scale for the hot virialized halo gas is too long to provide the source of cold gas. Supersonic, cold streams have been invoked as a mechanism for feeding massive halos at high redshift and deliver the cold gas required for continued star formation at the rates observed. This mechanism for replenishing the cold gas reservoir is motivated by some cosmological simulations. However, the cold streams are likely to be subject to the supersonic version of the Kelvin-Helmholtz instability (KHI), which eventually leads to stream disruption. Cosmological simulations have yet to obtain the spatial resolution required for understanding the detailed stability properties of cold streams. In this paper, we consider instead an idealized model of magnetized cold streams that we spatially resolve. Using linear theory we show how magnetic fields with dynamically important field strengths do not inhibit the KHI but rather enhance its growth rate. We perform nonlinear simulations of magnetized stream disruption and find that magnetic fields can nevertheless increase stream survival times by suppressing the mixing rate of cold gas with the circumgalactic medium. We find that magnetic fields can allow streams to survive $\sim 2-8$ times longer and, consequently, that streams $\sim 2-8$ times thinner can reach the central galaxy if the magnetic field strength is $\sim 0.3-0.8 \mu$G.

We present a method for efficiently modelling Braginskii viscosity on an unstructured, moving mesh. Braginskii viscosity, i.e., anisotropic transport of momentum with respect to the direction of the magnetic field, is thought to be of prime importance for studies of the weakly collisional plasma that comprises the intracluster medium (ICM) of galaxy clusters. Here anisotropic transport of heat and momentum has been shown to have profound consequences for the stability properties of the ICM. Our new method for modelling Braginskii viscosity has been implemented in the moving mesh code Arepo. We present a number of examples that serve to test the implementation and illustrate the modified dynamics found when including Braginskii viscosity in simulations. These include (but are not limited to) damping of fast magneto-sonic waves, interruption of linearly polarized Alfv\'en waves by the firehose instability and the inhibition of the Kelvin-Helmholtz instability by Braginskii viscosity. An explicit update of Braginskii viscosity is associated with a severe time step constraint that scales with $(\Delta x)^2$ where $\Delta x$ is the grid size. In our implementation, this restrictive time step constraint is alleviated by employing 2nd order accurate Runge-Kutta-Legendre super-time-stepping. We envision including Braginskii viscosity in future large-scale simulations of Kelvin-Helmholtz unstable cold fronts in cluster mergers and AGN-generated bubbles in central cluster regions.

The Kelvin-Helmholtz instability (KHI) is a standard test of hydrodynamic and magnetohydrodynamic (MHD) simulation codes and finds many applications in astrophysics. The classic linear theory considers a discontinuity in density and velocity at the interface of two fluids. However, for numerical simulations of the KHI such initial conditions do not yield converged results even at the linear stage of the instability. Instead, smooth profiles of velocity and density are required for convergence. This renders analytical theory to be only approximately valid and hinders quantitative comparisons between the classical theory and simulations. In this paper we derive a linear theory for the KHI with smooth profiles and illustrate code testing with the MHD code Athena. We provide the linear solution for the KHI with smooth initial conditions in three different limits: inviscid hydrodynamics, ideal MHD and Braginskii-MHD. These linear solutions are obtained numerically with the framework Psecas (Pseudo-Spectral Eigenvalue Calculator with an Automated Solver), which generates and solves numerical eigenvalue problems using an equation-parser and pseudo-spectral methods. The Athena simulations are carried out on a periodic, Cartesian domain which is useful for code testing purposes. Using Psecas and analytic theory, we outline the differences between this artificial numerical setup and the KHI on an infinite Cartesian domain and the KHI in cylindrical geometry. We discuss several astrophysical applications, such as cold flows in galaxy formation and cold fronts in galaxy cluster mergers. Psecas, and the linear solutions used for code testing, are publicly available and can be downloaded from the web.

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable, to different kinds of instabilities depending on the magnetic field orientation, at all radii. The fastest growing modes are usually related to generalizations of the Magnetothermal Instability (MTI) and the Heat-flux-driven Buoyancy Instability (HBI) which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow fast compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.