Turbulence In Intracluster Medium Research Articles

Multiscale analysis of turbulence evolution in the density-stratified intracluster medium

The diffuse hot medium inside clusters of galaxies typically exhibits turbulent motions whose amplitude increases with radius, as revealed by cosmological hydrodynamical simulations. However, its physical origin remains unclear. It could either be due to an excess injection of turbulence at large radii, or faster turbulence dissipation at small radii. We investigate this by studying the time evolution of turbulence in the intracluster medium (ICM) after major mergers, using the Omega500 non-radiative hydrodynamical cosmological simulations. By applying a novel wavelet analysis to study the radial dependence of the ICM turbulence spectrum, we discover that faster turbulence dissipation in the inner high density regions leads to the increasing turbulence amplitude with radius. We also find that the ICM turbulence at all radii decays in two phases after a major merger: an early fast decay phase followed by a slow secular decay phase. The buoyancy effects resulting from the ICM density stratification becomes increasingly important during turbulence decay, as revealed by a decreasing turbulence Froude number $Fr \sim \mathcal{O}(1)$. Our results indicate that the stronger density stratification and smaller eddy turn-over time are the likely causes of the faster turbulence dissipation rate in the inner regions of the cluster.

Three-dimensional Magnetohydrodynamical Simulations of the Morphology of Head–Tail Radio Galaxies Based on the Magnetic Tower Jet Model

Abstract The distinctive morphology of head–tail radio galaxies reveals strong interactions between the radio jets and their intra-cluster environment, the general consensus on the morphology origin of head–tail sources is that radio jets are bent by violent intra-cluster weather. We demonstrate in this paper that such strong interactions provide a great opportunity to study the jet properties and also the dynamics of the intra-cluster medium (ICM). By three-dimensional magnetohydrodynamical simulations, we analyze the detailed bending process of a magnetically dominated jet, based on the magnetic tower jet model. We use stratified atmospheres modulated by wind/shock to mimic the violent intra-cluster weather. Core sloshing is found to be inevitable during the wind-cluster core interaction, which induces significant shear motion and could finally drive ICM turbulence around the jet, making it difficult for the jet to survive. We perform a detailed comparison between the behavior of pure hydrodynamical jets and the magnetic tower jet and find that the jet-lobe morphology could not survive against the violent disruption in all of our pure hydrodynamical jet models. On the other hand, the head–tail morphology is well reproduced by using a magnetic tower jet model bent by wind, in which hydrodynamical instabilities are naturally suppressed and the jet could always keep its integrity under the protection of its internal magnetic fields. Finally, we also check the possibility for jet bending by shock only. We find that shock could not bend the jet significantly, and thus could not be expected to explain the observed long tails in head–tail radio galaxies.

Features of collisionless turbulence in the intracluster medium from simulated Faraday rotation maps – II. The effects of instabilities feedback

Statistical analysis of Faraday Rotation Measure (RM) maps of the intracluster medium (ICM) of galaxy clusters provides a unique tool to evaluate some spatial features of the magnetic fields there. Its combination with numerical simulations of magnetohydrodynamic (MHD) turbulence allows the diagnosis of the ICM turbulence. Being the ICM plasma weakly collisional, the thermal velocity distribution of the particles naturally develops anisotropies as a consequence of the large scale motions and the conservation of the magnetic moment of the charged particles. A previous study (Paper I) analyzed the impact of large scale thermal anisotropy on the statistics of RM maps synthesized from simulations of turbulence; these simulations employed a collisionless MHD model which considered a tensor pressure with uniform anisotropy. In the present work, we extend that analysis to a collisionless MHD model in which the thermal anisotropy develops according to the conservation of the magnetic moment of the thermal particles. We also consider the effect of anisotropy relaxation caused by the micro-scale mirror and firehose instabilities. We show that if the relaxation rate is fast enough to keep the anisotropy limited by the threshold values of the instabilities, the dispersion and power spectrum of the RM maps are indistinguishable from those obtained from collisional MHD. Otherwise, there is a reduction in the dispersion and steepening of the power spectrum of the RM maps (compared to the collisional case). Considering the first scenario, the use of collisional MHD simulations for modeling the RM statistics in the ICM becomes better justified.

Thermal SZ fluctuations in the ICM: probing turbulence and thermodynamics in Coma cluster withPlanck

We report the detection of thermal Sunyaev-Zeldovich effect (SZ) fluctuations in the intracluster medium (ICM) of Coma cluster observed with ${\it Planck}$. The SZ data links the maximum observable X-ray scale to the large Mpc scale, extending our knowledge of the power spectrum of ICM fluctuations. Deprojecting the 2-d SZ perturbations into 3-d pressure fluctuations, we find an amplitude spectrum which peaks at $\delta P/P = 33\pm 12\%$ and $74\pm19\%$ in the $15'$ and $40'$ radius region, respectively. We perform tests to ensure fluctuations are intrinsic to the cluster and not due to noise contamination. By using high-resolution hydrodynamical models, we improve the ICM turbulence constraints in Coma, finding 3-d Mach number ${\rm Ma_{3d}}= 0.8\pm0.3$ (15' region), increasing to supersonic values at larger radii (40'), and an injection scale $L_{\rm inj}\approx 500$ kpc. Such properties are consistent with driving due to mergers, in particular tied to internal galaxy groups. The large pressure fluctuations show that Coma is in adiabatic mode (mediated by sound waves), rather than isobaric mode (mediated by buoyancy waves). As predicted by turbulence models, the distribution of SZ fluctuations is log-normal with mild non-Gaussianities (heavy tails). The substantial non-thermal pressure support implies hydrostatic mass bias $b_M=-15\%$ to $-45\%$ from the core to the outskirt region, respectively. While total SZ power probes the thermal energy content, the SZ fluctuations constrain the non-thermal deviations important for precision cosmology. The proposed, novel approach can be exploited by multifrequency observations using ground based interferometers and future space CMB missions.

Stochastic reacceleration of relativistic electrons by turbulent reconnection: a mechanism for cluster-scale radio emission?

In this paper we investigate a situation where relativistic particles are reaccelerated diffusing across regions of reconnection and magnetic dynamo in super-Alfvenic, incompressible large-scale turbulence. We present an exploratory study of this mechanism in the intra-cluster-medium (ICM). In view of large-scale turbulence in the ICM we adopt a reconnection scheme that is based on turbulent reconnection and MHD turbulence. In this case particles are accelerated and decelerated in a systematic way in reconnecting and magnetic-dynamo regions, respectively, and on longer time-scales undergo a stochastic process diffusing across these sites (similar to second-order Fermi). Our study extends on larger scales numerical studies that focused on the acceleration in and around turbulent reconnecting regions. We suggest that this mechanism may play a role in the reacceleration of relativistic electrons in galaxy clusters providing a new physical scenario to explain the origin of cluster-scale diffuse radio emission. Indeed differently from current turbulent reacceleration models proposed for example for radio halos this mechanism is based on the effect of large-scale incompressible and super-Alfvenic turbulence. In this new model turbulence governs the interaction between relativistic particles and magnetic field lines that diffuse, reconnect and are stretched in the turbulent ICM.

Viscous Kelvin–Helmholtz instabilities in highly ionized plasmas

Transport coefficients in highly ionized plasmas like the intracluster medium (ICM) are still ill-constrained. They influence various processes, among them the mixing at shear flow interfaces due to the Kelvin–Helmholtz instability (KHI). The observed structure of potential mixing layers can be used to infer the transport coefficients, but the data interpretation requires a detailed knowledge of the long-term evolution of the KHI under different conditions. Here we present the first systematic numerical study of the effect of constant and temperature-dependent isotropic viscosity over the full range of possible values. We show that moderate viscosities slow down the growth of the KHI and reduce the height of the KHI rolls and their rolling-up. Viscosities above a critical value suppress the KHI. The effect can be quantified in terms of the Reynolds number Re=Uλν, where U is the shear velocity, λ the perturbation length and ν the kinematic viscosity. We derive the critical Re for constant and temperature-dependent Spitzer-like viscosities, an empirical relation for the viscous KHI growth time as a function of Re and density contrast, and describe special behaviours for Spitzer-like viscosities and high density contrasts. Finally, we briefly discuss several astrophysical situations where the viscous KHI could play a role, i.e. sloshing cold fronts, gas stripping from galaxies, buoyant cavities, ICM turbulence and high-velocity clouds.

Detection of a radio bridge in Abell 3667

We have detected a radio bridge of unpolarized synchrotron emission connecting the NW relic of the galaxy cluster Abell 3667 to its central regions. We used data at 2.3 GHz from the S-band Polarization All Sky Survey (S-PASS) and at 3.3 GHz from a follow up observation, both conducted with the Parkes Radio Telescope. This emission is further aligned with a diffuse X-ray tail, and represents the most compelling evidence for an association between intracluster medium turbulence and diffuse synchrotron emission. This is the first clear detection of a bridge associated both with an outlying cluster relic and X-ray diffuse emission. All the indicators point toward the synchrotron bridge being related to the post-shock turbulent wake trailing the shock front generated by a major merger in a massive cluster. Although predicted by simulations, this is the first time such emission is detected with high significance and clearly associated with the path of a confirmed shock. Although the origin of the relativistic electrons is still unknown, the turbulent re-acceleration model provides a natural explanation for the large-scale emission. The equipartition magnetic field intensity of the bridge is B_eq = 2.2 +/- 0.3 \mu G. We further detect diffuse emission coincident with the central regions of the cluster for the first time.

COMPARISONS OF COSMOLOGICAL MAGNETOHYDRODYNAMIC GALAXY CLUSTER SIMULATIONS TO RADIO OBSERVATIONS

Radio observations of galaxy clusters show that there are $\mu$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($\sigma_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $\sigma_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

The role of pressure anisotropy in the turbulent intracluster medium

In low-density plasma environments, such as the intracluster medium (ICM), the Larmour frequency is much larger than the ion-ion collision frequency. In such a case, the thermal pressure becomes anisotropic with respect to the magnetic field orientation and the evolution of the turbulent gas is more correctly described by a kinetic approach. A possible description of these collisionless scenarios is given by the so-called kinetic magnetohydrodynamic (KMHD) formalism, in which particles freely stream along the field lines, while moving with the field lines in the perpendicular direction. In this way a fluid-like behavior in the perpendicular plane is restored. In this work, we study fast growing magnetic fluctuations in the smallest scales which operate in the collisionless plasma that fills the ICM. In particular, we focus on the impact of a particular evolution of the pressure anisotropy and its implications for the turbulent dynamics of observables under the conditions prevailing in the ICM. We present results from numerical simulations and compare the results which those obtained using an MHD formalism.

EVOLUTION AND DISTRIBUTION OF MAGNETIC FIELDS FROM ACTIVE GALACTIC NUCLEI IN GALAXY CLUSTERS. I. THE EFFECT OF INJECTION ENERGY AND REDSHIFT

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.

TURBULENCE AND DYNAMO IN GALAXY CLUSTER MEDIUM: IMPLICATIONS ON THE ORIGIN OF CLUSTER MAGNETIC FIELDS

We present self-consistent cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and the magnetic field ejection by an active galactic nucleus (AGN). We find that the magnetic fields ejected by the AGNs, though initially distributed in relatively small volumes, can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process. The ICM turbulence is shown to be generated and sustained by the frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process is shown to exist in the ICM and amplify the magnetic field energy and flux. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs while micro Gauss ($\mu$G) fields can distribute over $\sim$ Mpc scales throughout the whole cluster. This finding shows that magnetic fields from AGNs, being further amplified by the ICM turbulence through small-scale dynamo processes, can be the origin of cluster-wide magnetic fields.

Three‐dimensional Magnetohydrodynamic Simulations of Cold Fronts in Magnetically Turbulent ICM

Steep gradients of temperature and density, called cold fronts, are observed by Chandra in a leading edge of subclusters moving through the intracluster medium (ICM). The presence of cold fronts indicates that thermal conduction across the front is suppressed by magnetic fields. We carried out three-dimensional magnetohydrodynamic (MHD) simulations including anisotropic thermal conduction of a subcluster moving through a magnetically turbulent ICM. We found that turbulent magnetic fields are stretched and amplified by shear flows along the interface between the subcluster and the ambient ICM. Since magnetic fields reduce the efficiency of thermal conduction across the front, the cold front survives at least 1 Gyr. We also found that a moving subcluster works as an amplifier of magnetic fields. Numerical results indicate that stretched turbulent magnetic fields accumulate behind the subcluster and are further amplified by vortex motions. The moving subcluster creates a long tail of ordered magnetic fields, in which the magnetic field strength attains plasma beta < 10.

AGN Feedback and Cooling Flows: Problems with Simple Hydrodynamic Models

In recent years it has become increasingly clear that active galactic nuclei, and radio galaxies in particular, have an impact on large-scale structure and galaxy formation. In principle, radio galaxies are energetic enough to halt the cooling of the virialized intracluster medium (ICM) in the inner regions of galaxy clusters, solving the cooling flow problem and explaining the high-mass truncation of the galaxy luminosity function. We explore this process through a series of high-resolution, three-dimensional hydrodynamic simulations of jetted active galaxies that act in response to cooling-mediated accretion of an ICM atmosphere. We find that our models are incapable of producing a long-term balance of heating and cooling; catastrophic cooling can be delayed by the jet action but inevitably takes hold. At the heart of the failure of these models is the formation of a low-density channel through which the jet can freely flow, carrying its energy out of the cooling core. It is possible that this failure is due to an oversimplified treatment of the fast jet (which may underestimate the dentist drill effect). However, it seems likely that additional complexity (large-angle jet precession or ICM turbulence) or additional physics (magnetohydrodynamic effects and plasma transport processes) is required to produce a spatial distribution of jet heating that can prevent catastrophic cooling. This work also underscores the importance of including jet dynamics in any feedback model, as opposed to the isotropically inflated bubble approach taken in some previous works.

Turbulence generation by substructure motion in clusters of galaxies

Abstract Clusters of galaxies form through major merger and/or absorption of smaller groups. In fact, some characteristic structures such as cold fronts, which are likely relevant to moving substructures, are found by Chandra. It is expected that moving substructures generate turbulence in the intracluster medium (ICM). Such turbulence probably plays a crucial role in mixture and transport of gas energy and heavy elements, and particle acceleration. The Astro-E2 satellite, which is planned to be launched in 2005, will detect broadened lines due to turbulent motion. In order to explore the above-mentioned issues, it is important to investigate the generation processes and structure of ICM turbulence. We investigate the ICM dynamical evolution in and around a moving substructure with three-dimensional hydrodynamical simulations. Eddy-like structures develop near the boundary between the substructure and the ambient ICM through Kelvin–Helmholtz instabilities. Because of these structures, characteristic patterns appear in the line-of-sight velocity distribution of the ICM.