We present ALMA CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in the 2A 0335+096 galaxy cluster (z = 0.0346). The total molecular gas mass of (1.13+/-0.15) x 10^9 M_sun is divided into two components: a nuclear region and a 7 kpc long dusty filament. The central molecular gas component accounts for (3.2+/-0.4) x 10^8 M_sun of the total supply of cold gas. Instead of forming a rotationally-supported ring or disk, it is composed of two distinct, blueshifted clumps south of the nucleus and a series of low-significance redshifted clumps extending toward a nearby companion galaxy. The velocity of the redshifted clouds increases with radius to a value consistent with the companion galaxy, suggesting that an interaction between these galaxies <20 Myr ago disrupted a pre-existing molecular gas reservoir within the BCG. Most of the molecular gas, (7.8+/-0.9) x 10^8 M_sun, is located in the filament. The CO emission is co-spatial with a 10^4 K emission-line nebula and soft X-rays from 0.5 keV gas, indicating that the molecular gas has cooled out of the intracluster medium over a period of 25-100 Myr. The filament trails an X-ray cavity, suggesting that the gas has cooled from low entropy gas that has been lifted out of the cluster core and become thermally unstable. We are unable to distinguish between inflow and outflow along the filament with the present data. Cloud velocities along the filament are consistent with gravitational free-fall near the plane of the sky, although their increasing blueshifts with radius are consistent with outflow.

We present ALMA observations of the CO(1-0) and CO(3-2) line emission tracing filaments of cold molecular gas in the central galaxy of the cluster PKS0745-191. The total molecular gas mass of 4.6 +/- 0.3 x 10^9 solar masses, assuming a Galactic X_{CO} factor, is divided roughly equally between three filaments each extending radially 3-5 kpc from the galaxy centre. The emission peak is located in the SE filament roughly 1 arcsec (2 kpc) from the nucleus. The velocities of the molecular clouds in the filaments are low, lying within +/-100 km/s of the galaxy's systemic velocity. Their FWHMs are less than 150 km/s, which is significantly below the stellar velocity dispersion. Although the molecular mass of each filament is comparable to a rich spiral galaxy, such low velocities show that the filaments are transient and the clouds would disperse on <10^7 yr timescales unless supported, likely by the indirect effect of magnetic fields. The velocity structure is inconsistent with a merger origin or gravitational free-fall of cooling gas in this massive central galaxy. If the molecular clouds originated in gas cooling even a few kpc from their current locations their velocities would exceed those observed. Instead, the projection of the N and SE filaments underneath X-ray cavities suggests they formed in the updraft behind bubbles buoyantly rising through the cluster atmosphere. Direct uplift of the dense gas by the radio bubbles appears to require an implausibly high coupling efficiency. The filaments are coincident with low temperature X-ray gas, bright optical line emission and dust lanes indicating that the molecular gas could have formed from lifted warmer gas that cooled in situ.

We use analytic estimates and numerical simulations of test particles interacting with magnetohydrodynamic (MHD) turbulence to show that subsonic MHD turbulence produces efficient second-order Fermi acceleration of relativistic particles. This acceleration is not well-described by standard quasi-linear theory but is a consequence of resonance broadening of wave-particle interactions in MHD turbulence. We provide momentum diffusion coefficients that can be used for astrophysical and heliospheric applications and discuss the implications of our results for accretion flows onto black holes. In particular, we show that particle acceleration by subsonic turbulence in radiatively inefficient accretion flows can produce a non-thermal tail in the electron distribution function that is likely important for modeling and interpreting the emission from low luminosity systems such as Sgr A* and M87.

We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed clusters there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are dominated by buoyancy, down to small, turbulent scales: $(\delta\rho_k/\rho)^2 = \eta_1^2 (V_{1,k}/c_s)^2$, where $\delta\rho_k/\rho$ is the spectral amplitude of the density perturbations at wave number $k$, $V_{1,k}^2=V_k^2/3$ is the mean square component of the velocity field, $c_s$ is the sound speed, and $\eta_1$ is a dimensionless constant of order unity. Using cosmological simulations of relaxed galaxy clusters, we calibrate this relation and find $\eta_1\approx 1 \pm 0.3$. We argue that this value is set at large scales by buoyancy physics, while at small scales the density and velocity power spectra are proportional because the former are a passive scalar advected by the latter. This opens an interesting possibility to use gas density power spectra as a proxy for the velocity power spectra in relaxed clusters, across a wide range of scales.

We investigate the effects of pitch-angle scattering on the efficiency of particle heating and acceleration by MHD turbulence using phenomenological estimates and simulations of non-relativistic test particles interacting with strong, subsonic MHD turbulence. We include an imposed pitch-angle scattering rate, which is meant to approximate the effects of high frequency plasma waves and/or velocity space instabilities. We focus on plasma parameters similar to those found in the near-Earth solar wind, though most of our results are more broadly applicable. An important control parameter is the size of the particle mean free path lambda_{mfp} relative to the scale of the turbulent fluctuations L. For small scattering rates, particles interact quasi-resonantly with turbulent fluctuations in magnetic field strength. Scattering increases the long-term efficiency of this resonant heating by factors of a few-10, but the distribution function does not develop a significant non-thermal power-law tail. For higher scattering rates, the interaction between particles and turbulent fluctuations becomes non-resonant, governed by particles heating and cooling adiabatically as they encounter turbulent density fluctuations. Rapid pitch-angle scattering can produce a power-law tail in the proton distribution function but this requires fine-tuning of parameters. Moreover, in the near-Earth solar wind, a significant power-law tail cannot develop by this mechanism because the particle acceleration timescales are longer than the adiabatic cooling timescale set by the expansion of the solar wind. Our results thus imply that MHD-scale turbulent fluctuations are unlikely to be the origin of the v^{-5} tail in the proton distribution function observed in the solar wind.

We present a spherically symmetric model for the origin and evolution of the temperature profiles in the hot plasma filling galaxy groups and clusters. We find that the gas in clusters is generically not isothermal, and that the temperature declines with radius at large distances from the cluster center (outside the core- and scale radii). This temperature profile is determined by the accretion history of the halo, and is not quantitatively well-described by a polytropic model. We explain quantitatively how the large-scale temperature gradient persists in spite of thermal conduction and convection. These results are a consequence of the cosmological assembly of clusters and cannot be reproduced with non-cosmological simulations of isolated halos. We show that the variation in halo assembly histories produces a ~10% scatter in temperature at fixed mass. On top of this scatter, conduction decreases the temperature of the gas near the scale radius in massive clusters, which may bias hydrostatic mass estimates inferred from x-ray and SZ observations. As an example application of our model profiles, we use mixing-length theory to estimate the turbulent pressure support created by the magnetothermal instability (MTI): in agreement with our earlier MHD simulations, we find that the convection produced by the MTI can provide ~5% non-thermal pressure support near r_500. The magnitude of this turbulent pressure support is likely to be non-monotonic in halo mass, peaking in ~10^14.5 M_sun halos.

We study the interplay among cooling, heating, conduction, and magnetic fields in gravitationally stratified plasmas using simplified, plane-parallel numerical simulations. Since the physical heating mechanism remains uncertain in massive halos such as groups or clusters, we adopt a simple, observationally-motivated prescription which enforces global thermal equilibrium when averaged over large scales. The plasma remains susceptible to local thermal instability, however, and cooling drives an inward flow of material. In contrast to previous results, we argue that the thermal stability of the plasma is independent of its convective stability. We find that the ratio of the cooling timescale to the dynamical timescale t_cool/t_ff controls the saturation of the thermal instability: when t_cool/t_ff < 1, the plasma develops extended multi-phase structure, whereas when t_cool / t_ff > 1 it does not. (In a companion paper, we show that the criterion for thermal instability in a spherical potential is somewhat less stringent, t_cool / t_ff < 10.) When thermal conduction is anisotropic with respect to the magnetic field, the criterion for multi-phase structure is essentially independent of the thermal conductivity of the plasma. Our criterion for local thermal instability to produce multi-phase structure is an extension of the cold vs. hot accretion modes in galaxy formation that applies at all radii in hot halos, not just to the virial shock. We show that this criterion is consistent with data on multi-phase gas in the ACCEPT sample of clusters; in addition, when t_cool / t_ff > 1, the net cooling rate to low temperatures and the mass flux to small radii are suppressed enough relative to models without heating to be qualitatively consistent with star formation rates and x-ray line emission in groups and clusters.

We present one-dimensional models of the hot gas in dark-matter halos, which both predict the existence of cool cores and explain their structure. Our models are directly applicable to semi-analytic models (SAMs) of galaxy formation. We have previously argued that filaments of cold (~10^4 K) gas condense out of the intracluster medium (ICM) in hydrostatic and thermal equilibrium when the ratio of the thermal instability timescale to the free-fall time $t_{TI}/t_{ff}$ falls below 5-10. This criterion corresponds to an upper limit on the density of the ICM and motivates a model in which a density core forms wherever $t_{TI}/t_{ff} \lesssim 10$. Consistent with observations and numerical simulations, this model predicts larger and more tenuous cores for lower-mass halos---while the core density in a cluster may be as large as ~ 0.1 cm^{-3}, the core density in the Galactic halo should not exceed ~ 10^{-4} cm^{-3}. We can also explain the large densities in smaller mass halos (galactic `coronae') if we include the contribution of the central galaxy to the gravitational potential. Our models produce a favorable match to the observational X-ray luminosity-temperature (Lx-Tx) relation. For halo masses $\lesssim 10^{13}$ Msun the core size approaches the virial radius. Thus, most of the baryons in such halos cannot be in the hot ICM, but either in the form of stars or in the form of hot gas beyond the virial radius. Because of the smaller mass in the ICM and much larger mass available for star formation, the majority of the baryons in low mass halos ($\lesssim 10^{13}$ Msun) can be expelled beyond the virial radius due to supernova feedback. This can account for the baryons `missing' from low mass halos, such as the Galactic halo.

The heating, acceleration, and pitch-angle scattering of charged particles by MHD turbulence are important in a wide range of astrophysical environments, including the solar wind, accreting black holes, and galaxy clusters. We simulate the interaction of high-gyrofrequency test particles with fully dynamical simulations of subsonic MHD turbulence, focusing on the parameter regime with beta ~ 1, where beta is the ratio of gas to magnetic pressure. We use the simulation results to calibrate analytical expressions for test particle velocity-space diffusion coefficients and provide simple fits that can be used in other work. The test particle velocity diffusion in our simulations is due to a combination of two processes: interactions between particles and magnetic compressions in the turbulence (as in linear transit-time damping; TTD) and what we refer to as Fermi Type-B (FTB) interactions, in which charged particles moving on field lines may be thought of as beads spiralling around moving wires. We show that test particle heating rates are consistent with a TTD resonance which is broadened according to a decorrelation prescription that is Gaussian in time. TTD dominates the heating for v_s >> v_A (e.g. electrons), where v_s is the thermal speed of species s and v_A is the Alfven speed, while FTB dominates for v_s << v_A (e.g. minor ions). Proton heating rates for beta ~ 1 are comparable to the turbulent cascade rate. Finally, we show that velocity diffusion of collisionless, large gyrofrequency particles due to large-scale MHD turbulence does not produce a power-law distribution function.

X-ray surface brightness fluctuations in the core ($650 \times 650$ kpc) region of the Coma cluster observed with XMM-Newton and Chandra are analyzed using a 2D power spectrum approach. The resulting 2D spectra are converted to 3D power spectra of gas density fluctuations. Our independent analyses of the XMM-Newton and Chandra observations are in excellent agreement and provide the most sensitive measurements of surface brightness and density fluctuations for a hot cluster. We find that the characteristic amplitude of the volume filling density fluctuations relative to the smooth underlying density distribution varies from 7-10% on scales of $\sim$500 kpc down to $\sim$5% at scales $\sim$ 30 kpc. On smaller spatial scales, projection effects smear the density fluctuations by a large factor, precluding strong limits on the fluctuations in 3D. On the largest scales probed (hundreds of kpc), the dominant contributions to the observed fluctuations most likely arise from perturbations of the gravitational potential by the two most massive galaxies in Coma, NGC4874 and NGC4889, and the low entropy gas brought to the cluster by an infalling group. Other plausible sources of X-ray surface brightness fluctuations are discussed, including turbulence, metal abundance variations, and unresolved sources. Despite a variety of possible origins for density fluctuations, the gas in the Coma cluster core is remarkably homogeneous on scales from $\sim$ 500 to $\sim$30 kpc.

We study the effects of anisotropic thermal conduction along magnetic field lines on an accelerated contact discontinuity in a weakly collisional plasma. We first perform a linear stability analysis similar to that used to derive the Rayleigh-Taylor instability (RTI) dispersion relation. We find that anisotropic conduction is only important for compressible modes, as incompressible modes are isothermal. Modes grow faster in the presence of anisotropic conduction, but growth rates do not change by more than a factor of order unity. We next run fully non-linear numerical simulations of a contact discontinuity with anisotropic conduction. The non-linear evolution can be thought of as a superposition of three physical effects: temperature diffusion due to vertical conduction, the RTI, and the heat flux driven buoyancy instability (HBI). In simulations with RTI-stable contact discontinuities, the temperature discontinuity spreads due to vertical heat conduction. This occurs even for initially horizontal magnetic fields due to the initial vertical velocity perturbation and numerical mixing across the interface. The HBI slows this temperature diffusion by reorienting initially vertical magnetic field lines to a more horizontal geometry. In simulations with RTI-unstable contact discontinuities, the dynamics are initially governed by temperature diffusion, but the RTI becomes increasingly important at late times. We discuss the possible application of these results to supernova remnants, solar prominences, and cold fronts in galaxy clusters.

In the intracluster medium (ICM) of galaxy clusters, heat and momentum are transported almost entirely along (but not across) magnetic field lines. We perform the first fully self-consistent Braginskii-MHD simulations of galaxy clusters including both of these effects. Specifically, we perform local and global simulations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability (HBI) and assess the effects of viscosity on their saturation and astrophysical implications. We find that viscosity has only a modest effect on the saturation of the MTI. As in previous calculations, we find that the MTI can generate nearly sonic turbulent velocities in the outer parts of galaxy clusters, although viscosity somewhat suppresses the magnetic field amplification. At smaller radii in cool-core clusters, viscosity can decrease the linear growth rates of the HBI. However, it has less of an effect on the HBI's nonlinear saturation, in part because three-dimensional interchange motions (magnetic flux tubes slipping past each other) are not damped by anisotropic viscosity. In global simulations of cool core clusters, we show that the HBI robustly inhibits radial thermal conduction and thus precipitates a cooling catastrophe. The effects of viscosity are, however, more important for higher entropy clusters. We argue that viscosity can contribute to the global transition of cluster cores from cool-core to non cool-core states: additional sources of intracluster turbulence, such as can be produced by AGN feedback or galactic wakes, suppress the HBI, heating the cluster core by thermal conduction; this makes the ICM more viscous, which slows the growth of the HBI, allowing further conductive heating of the cluster core and a transition to a non cool-core state.

We use three-dimensional MHD simulations with anisotropic thermal conduction to study turbulence due to the magnetothermal instability (MTI) in the intracluster medium (ICM) of galaxy clusters. The MTI grows on timescales of ~1 Gyr and is capable of driving vigorous, sustained turbulence in the outer parts of galaxy clusters if the temperature gradient is maintained in spite of the rapid thermal conduction. If this is the case, turbulence due to the MTI can provide up to 5-30% of the pressure support beyond r_500 in galaxy clusters, an effect that is strongest for hot, massive clusters. The turbulence driven by the MTI is generally additive to other sources of turbulence in the ICM, such as that produced by structure formation. This new source of non-thermal pressure support reduces the observed Sunyaev-Zel'dovich (SZ) signal and X-ray pressure gradient for a given cluster mass and introduces a cluster mass and temperature gradient-dependent bias in SZ and X-ray mass estimates of clusters. This additional physics may also need to be taken into account when estimating the matter power spectrum normalization, sigma-8, through simulation templates from the observed amplitude of the SZ power spectrum.

Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- sociated with supernovae (Woosley et al. 2006). They release ultra-relativistic jets producing non-thermal emission through synchrotron radiation as they interact with the surrounding medium (Zhang et al. 2004). Here we report observations of the peculiar GRB 101225A (the "Christmas burst"). Its gamma-ray emission was exceptionally long and followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which gets thermalized by interacting with the dense, previously ejected material and thus creating the observed black-body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy (Campana et al. 2011).

Observations of clusters and groups imply that such halos are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally-motivated heating prescriptions, we show that local thermal instability (TI) can produce a multi-phase medium---with ~ 10000 K cold filaments condensing out of the hot ICM---only when the ratio of the TI timescale in the hot plasma (t_{TI}) to the free-fall timescale (t_{ff}) satisfies t_{TI}/t_{ff} <~ 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling, and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of ~ 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling-flow is ~ 0.001 for clusters and decreases for lower mass halos; supernova heating may be energetically sufficient to balance cooling in galactic halos. We further argue that the ICM self-adjusts so that t_{TI}/t_{ff} >~ 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives t_{TI}/t_{ff} >~ 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminosity-temperature (L_X-T_X) relation. The high-velocity clouds observed in the Galactic halo can also be due to local TI producing multi-phase gas close to the virial radius.

(abridged) The ICM has been suggested to be buoyantly unstable in the presence of magnetic field and anisotropic thermal conduction. We perform first cosmological simulations of galaxy cluster formation that simultaneously include magnetic fields, radiative cooling and anisotropic thermal conduction. In isolated and idealized cluster models, the magnetothermal instability (MTI) tends to reorient the magnetic fields radially. Using cosmological simulations of the Santa Barbara cluster we detect radial bias in the velocity and magnetic fields. Such radial bias is consistent with either the inhomogeneous radial gas flows due to substructures or residual MTI-driven field rearangements that are expected even in the presence of turbulence. Although disentangling the two scenarios is challenging, we do not detect excess bias in the runs that include anisotropic thermal conduction. The anisotropy effect is potentially detectable via radio polarization measurements with LOFAR and SKA and future X-ray spectroscopic studies with the IXO. We demonstrate that radiative cooling boosts the amplification of the magnetic field by about two orders of magnitude beyond what is expected in the non-radiative cases. At z=0 the field is amplified by a factor of about 10^6 compared to the uniform magnetic field evolved due to the universal expansion alone. Interestingly, the runs that include both radiative cooling and anisotropic thermal conduction exhibit stronger magnetic field amplification than purely radiative runs at the off-center locations. In these runs, shallow temperature gradients away from the cluster center make the ICM neutrally buoyant. The ICM is more easily mixed in these regions and the winding up of the frozen-in magnetic field is more efficient resulting in stronger magnetic field amplification.

We study the effects of anisotropic thermal conduction on low-collisionality, astrophysical plasmas using two and three-dimensional magnetohydrodynamic simulations. For weak magnetic fields, dilute plasmas are buoyantly unstable for either sign of the temperature gradient: the heat-flux-driven buoyancy instability (HBI) operates when the temperature increases with radius while the magnetothermal instability (MTI) operates in the opposite limit. In contrast to previous results, we show that, in the presence of a sustained temperature gradient, the MTI drives strong turbulence and operates as an efficient magnetic dynamo (akin to standard, adiabatic convection). Together, the turbulent and magnetic energies contribute up to ~10% of the pressure support in the plasma. In addition, the MTI drives a large convective heat flux, ~1.5% of rho c_s^3. These findings are robust even in the presence of an external source of strong turbulence. Our results on the nonlinear saturation of the HBI are consistent with previous studies but we explain physically why the HBI saturates quiescently by re-orienting the magnetic field (suppressing the conductive heat flux through the plasma), while the MTI saturates by generating sustained turbulence. We also systematically study how an external source of turbulence affects the saturation of the HBI: such turbulence can disrupt the HBI only on scales where the shearing rate of the turbulence is faster than the growth rate of the HBI. In particular, our results provide a simple mapping between the level of turbulence in a plasma and the effective isotropic thermal conductivity. We discuss the astrophysical implications of these findings, with a particular focus on the intracluster medium of galaxy clusters.

We study the effects of externally imposed turbulence on the thermal properties of galaxy cluster cores, using three-dimensional numerical simulations including magnetic fields, anisotropic thermal conduction, and radiative cooling. The imposed "stirring" crudely approximates the effects of galactic wakes, waves generated by galaxies moving through the intracluster medium (ICM), and/or turbulence produced by a central active galactic nucleus. The simulated clusters exhibit a strong bimodality. Modest levels of turbulence, ~100 km/s (~10% of the sound speed), suppress the heat-flux-driven buoyancy instability (HBI), resulting in an isotropically tangled magnetic field and a quasi-stable, high entropy, thermal equilibrium with no cooling catastrophe. Thermal conduction dominates the heating of the cluster core, but turbulent mixing is critical because it suppresses the HBI and (to a lesser extent) the thermal instability. Lower levels of turbulent mixing (approximately less than 100 km/s) are insufficient to suppress the HBI, rapidly leading to a thermal runaway and a cool-core cluster. Remarkably, then, small fluctuations in the level of turbulence in galaxy cluster cores can initiate transitions between cool-core (low entropy) and non cool-core (high entropy) states.

Observations of the cores of nearby galaxy clusters show H$\alpha$ and molecular emission line filaments. We argue that these are the result of {\em local} thermal instability in a {\em globally} stable galaxy cluster core. We present local, high resolution, two-dimensional magnetohydrodynamic simulations of thermal instability for conditions appropriate to the intracluster medium (ICM); the simulations include thermal conduction along magnetic field lines and adiabatic cosmic rays. Thermal conduction suppresses thermal instability along magnetic field lines on scales smaller than the Field length ($\gtrsim$10 kpc for the hot, diffuse ICM). We show that the Field length in the cold medium must be resolved both along and perpendicular to the magnetic field in order to obtain numerically converged results. Because of negligible conduction perpendicular to the magnetic field, thermal instability leads to fine scale structure in the perpendicular direction. Filaments of cold gas along magnetic field lines are thus a natural consequence of thermal instability with anisotropic thermal conduction. Nonlinearly, filaments of cold ($\sim 10^4$ K) gas should have lengths (along the magnetic field) comparable to the Field length in the cold medium $\sim 10^{-4}$ pc! Observations show, however, that the atomic filaments in clusters are far more extended, $\sim 10$ kpc. Cosmic ray pressure support (or a small scale turbulent magnetic pressure) may resolve this discrepancy: even a small cosmic ray pressure in the diffuse ICM, $\sim 10^{-4}$ of the thermal pressure, can be adiabatically compressed to provide significant pressure support in cold filaments. This is qualitatively consistent with the large population of cosmic rays invoked to explain the atomic and molecular line ratios observed in filaments.

The intracluster medium (ICM) is stably stratified in the hydrodynamic sense with the entropy $s$ increasing outwards. However, thermal conduction along magnetic field lines fundamentally changes the stability of the ICM, leading to the "heat-flux buoyancy instability" when $dT/dr>0$ and the "magnetothermal instability" when $dT/dr<0$. The ICM is thus buoyantly unstable regardless of the signs of $dT/dr$ and $ds/dr$. On the other hand, these temperature-gradient-driven instabilities saturate by reorienting the magnetic field (perpendicular to $\hat{\bf r}$ when $dT/dr>0$ and parallel to $\hat{\bf r}$ when $dT/dr<0$), without generating sustained convection. We show that after an anisotropically conducting plasma reaches this nonlinearly stable magnetic configuration, it experiences a buoyant restoring force that resists further distortions of the magnetic field. This restoring force is analogous to the buoyant restoring force experienced by a stably stratified adiabatic plasma. We argue that in order for a driving mechanism (e.g, galaxy motions or cosmic-ray buoyancy) to overcome this restoring force and generate turbulence in the ICM, the strength of the driving must exceed a threshold, corresponding to turbulent velocities $\gtrsim 10 -100 {km/s}$. For weaker driving, the ICM remains in its nonlinearly stable magnetic configuration, and turbulent mixing is effectively absent. We discuss the implications of these findings for the turbulent diffusion of metals and heat in the ICM.

We study the heating of charged test particles in three-dimensional numerical simulations of weakly compressible magnetohydrodynamic (MHD) turbulence (``Alfvenic turbulence''); these results are relevant to particle heating and acceleration in the solar wind, solar flares, accretion disks onto black holes, and other astrophysics and heliospheric environments. The physics of particle heating depends on whether the gyrofrequency of a particle is comparable to the frequency of a turbulent fluctuation that is resolved on the computational domain. Particles with these frequencies nearly equal undergo strong perpendicular heating (relative to the local magnetic field) and pitch angle scattering. By contrast, particles with large gyrofrequency undergo strong parallel heating. Simulations with a finite resistivity produce additional parallel heating due to parallel electric fields in small-scale current sheets. Many of our results are consistent with linear theory predictions for the particle heating produced by the Alfven and slow magnetosonic waves that make up Alfvenic turbulence. However, in contrast to linear theory predictions, energy exchange is not dominated by discrete resonances between particles and waves; instead, the resonances are substantially ``broadened.'' We discuss the implications of our results for solar and astrophysics problems, in particular the thermodynamics of the near-Earth solar wind. We conclude that Alfvenic turbulence produces significant parallel heating via the interaction between particles and magnetic field compressions (``slow waves''). However, on scales above the proton Larmor radius, Alfvenic turbulence does not produce significant perpendicular heating of protons or minor ions.

We perform a suite of simulations of cooling cores in clusters of galaxies in order to investigate the effect of the recently discovered heat flux buoyancy instability (HBI) on the evolution of cores. Our models follow the 3-dimensional magnetohydrodynamics (MHD) of cooling cluster cores and capture the effects of anisotropic heat conduction along the lines of magnetic field, but do not account for the cosmological setting of clusters or the presence of AGN. Our model clusters can be divided into three groups according to their final thermodynamical state: catastrophically collapsing cores, isothermal cores, and an intermediate group whose final state is determined by the initial configuration of magnetic field. Modeled cores that are reminiscent of real cluster cores show evolution towards thermal collapse on a time scale which is prolonged by a factor of ~2-10 compared with the zero-conduction cases. The principal effect of the HBI is to re-orient field lines to be perpendicular to the temperature gradient. Once the field has been wrapped up onto spherical surfaces surrounding the core, the core is insulated from further conductive heating (with the effective thermal conduction suppressed to less than 1/100th of the Spitzer value) and proceeds to collapse. We speculate that, in real clusters, the central AGN and possibly mergers play the role of "stirrers," periodically disrupting the azimuthal field structure and allowing thermal conduction to sporadically heat the core.

We examine the long-standing cooling flow problem in galaxy clusters with 3D MHD simulations of isolated clusters including radiative cooling and anisotropic thermal conduction along magnetic field lines. The central regions of the intracluster medium (ICM) can have cooling timescales of ~200 Myr or shorter--in order to prevent a cooling catastrophe the ICM must be heated by some mechanism such as AGN feedback or thermal conduction from the thermal reservoir at large radii. The cores of galaxy clusters are linearly unstable to the heat-flux-driven buoyancy instability (HBI), which significantly changes the thermodynamics of the cluster core. The HBI is a convective, buoyancy-driven instability that rearranges the magnetic field to be preferentially perpendicular to the temperature gradient. For a wide range of parameters, our simulations demonstrate that in the presence of the HBI, the effective radial thermal conductivity is reduced to less than 10% of the full Spitzer conductivity. With this suppression of conductive heating, the cooling catastrophe occurs on a timescale comparable to the central cooling time of the cluster. Thermal conduction alone is thus unlikely to stabilize clusters with low central entropies and short central cooling timescales. High central entropy clusters have sufficiently long cooling times that conduction can help stave off the cooling catastrophe for cosmologically interesting timescales.

Using a linear stability analysis and two and three-dimensional nonlinear simulations, we study the physics of buoyancy instabilities in a combined thermal and relativistic (cosmic ray) plasma, motivated by the application to clusters of galaxies. We argue that cosmic ray diffusion is likely to be slow compared to the buoyancy time on large length scales, so that cosmic rays are effectively adiabatic. If the cosmic ray pressure $p_{cr}$ is $\gtrsim 25 %$ of the thermal pressure, and the cosmic ray entropy ($p_{\rm cr}/\rho^{4/3}$; $\rho$ is the thermal plasma density) decreases outwards, cosmic rays drive an adiabatic convective instability analogous to Schwarzschild convection in stars. Global simulations of galaxy cluster cores show that this instability saturates by reducing the cosmic ray entropy gradient and driving efficient convection and turbulent mixing. At larger radii in cluster cores, the thermal plasma is unstable to the heat flux-driven buoyancy instability (HBI), a convective instability generated by anisotropic thermal conduction and a background conductive heat flux. Cosmic-ray driven convection and the HBI may contribute to redistributing metals produced by Type 1a supernovae in clusters. Our calculations demonstrate that adiabatic simulations of galaxy clusters can artificially suppress the mixing of thermal and relativistic plasma; anisotropic thermal conduction allows more efficient mixing, which may contribute to cosmic rays being distributed throughout the cluster volume.

The electron mean free path in the intracluster medium (ICM) of galaxy clusters is much larger than the gyroradius; thus, heat is transported anisotropically along magnetic field lines. We show that the intracluster medium is unstable to the magnetothermal instability (MTI) using MHD simulations with anisotropic thermal conduction. As a result of the MTI, we find that the temperature profile of the ICM can be substantially modified on timescales of several billion years while the magnetic field is amplified by dynamo action up to more than fifty times the original energy. We also show that the instability drives field lines to become preferentially radial leading to conduction that is a highly efficient fraction of the Spitzer conductivity. As such, we present the first self-consistent calculation of the effective thermal conductivity in the ICM.

In low collisionality plasmas heat flows almost exclusively along magnetic field lines, and the condition for stability to convection is modified from the standard Schwarzschild criterion. We present local two and three-dimensional simulations of a new heat flux driven buoyancy instability (the HBI) that occurs when the temperature in a plasma decreases in the direction of gravity. We find that the HBI drives a convective dynamo that amplifies an initially weak magnetic field by a factor of ~20. In simulations that begin with the magnetic field aligned with the temperature gradient, the HBI saturates by rearranging the magnetic field lines to be almost purely perpendicular to the initial temperature gradient. This magnetic field reorientation results in a net heat flux through the plasma that is less than 1% of the field-free (Spitzer) value. We show that the HBI is likely to be present in the cool cores of clusters of galaxies between ~0.1-100 kpc, where the temperature increases outwards. The saturated state of the HBI suggests that inward thermal conduction from large radii in clusters is unlikely to solve the cooling flow problem. Finally, we also suggest that the HBI may contribute to suppressing conduction across cold fronts in galaxy clusters.

We perform a linear magnetohydrodynamic perturbation analysis for a stratified magnetized envelope where the diffusion of heat is mediated by charged particles that are confined to flow along magnetic field lines. We identify an instability, the ``coulomb bubble instability,'' which may be thought of as standard magnetosonic fast and slow waves, driven by the rapid diffusion of heat along the direction of the magnetic field. We calculate the growth rate and stability criteria for the coulomb bubble instability for various choices of equilibrium conditions. The coulomb bubble instability is most strongly driven for weakly magnetized atmospheres that are strongly convectively stable. We briefly discuss a possible application of astrophysical interest: diffusion of interstellar cosmic rays in the hot T ~ 10^6 K Galactic corona. We show that for commonly accepted values of the cosmic ray and gas pressure as well as its overall characteristic dimensions, the Galactic corona is in a marginal state of stability with respect to a cosmic ray coulomb bubble instability. The implication being that a cosmic ray coulomb bubble instability plays a role regulating both the pressure and transport properties of interstellar cosmic rays, while serving as a source of acoustic power above the galactic disk. (abridged)

In dilute astrophysical plasmas, thermal conduction is primarily along magnetic field lines, and therefore highly anisotropic. As a result, the usual convective stability criterion is modified from a condition on entropy to a condition on temperature. For small magnetic fields or small wavenumbers, instability occurs in any atmosphere where the temperature and pressure gradients point in the same direction. We refer to the resulting convective instability as the magnetothermal instability (MTI). We present fully three-dimensional simulations of the MTI and show that saturation results in an atmosphere with different vertical structure, dependent upon the boundary conditions. When the temperature at the boundary of the unstable layer is allowed to vary, the temperature gradient relaxes until the unstable region is almost isothermal. When the temperature at the boundary of the unstable region is fixed, the magnetic field is reoriented to an almost vertical geometry as a result of buoyant motions. This case exhibits more vigorous turbulence. In both cases the resulting saturated heat flux is almost one-half of the value expected if the conduction were purely isotropic. The action of the MTI results in dynamical processes that lead to significant transport perpendicular to the initial direction of the magnetic field. The resulting magnetoconvection in both cases amplifies the magnetic field until it is almost in equipartition with sustained subsonic turbulence. These results are relevant to understanding measurements of the temperature profiles of the intracluster medium of clusters of galaxies as well as the structure of radiatively inefficient accretion flows.

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows.

In weakly magnetized, dilute plasmas in which thermal conduction along magnetic field lines is important, the usual convective stability criterion is modified. Instead of depending on entropy gradients, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) > 0, which we refer to as the Balbus criterion. We refer to the convective instability that results in this regime as the magnetothermal instability (MTI). We use numerical MHD simulations which include anisotropic electron heat conduction to follow the growth and saturation of the MTI in two-dimensional, plane parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured in the simulations agree with the weak field dispersion relation. We investigate the effect of strong fields and isotropic conduction on the linear properties and nonlinear regime of the MTI. In the nonlinear regime, the instability saturates and convection decays away, when the atmosphere becomes isothermal. Sustained convective turbulence can be driven if there is a fixed temperature difference between the top and bottom edges of the simulation domain, and if isotropic conduction is used to create convectively stable layers that prevent the formation of unresolved, thermal boundary layers. The largest component of the time-averaged heat flux is due to advective motions. These results have implications for a variety of astrophysical systems, such as the temperature profile of hot gas in galaxy clusters, and the structure of radiatively inefficient accretion flows.