Normalized to: Parrish, I.

**30** article(s) in total. **82** co-authors, from **1** to **18** common article(s). Median position in authors list is **2,5**.

[1] oai:arXiv.org:1610.00716 [pdf] - 1531849

Molecular Gas Along a Bright H-alpha Filament in 2A 0335+096 Revealed by
ALMA

Comments: 21 pages, 14 figures, accepted to ApJ

Submitted: **2016-10-03**, last modified: **2016-10-07**

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.

[2] oai:arXiv.org:1602.05962 [pdf] - 1370417

ALMA observations of cold molecular gas filaments trailing rising radio
bubbles in PKS0745-191

Comments: 16 pages, 10 figures, accepted by MNRAS

Submitted: **2016-02-18**

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.

[3] oai:arXiv.org:1403.3123 [pdf] - 857783

Acceleration of Relativistic Electrons by MHD Turbulence: Implications
for Non-thermal Emission from Black Hole Accretion Disks

Comments: 7 pages, 4 figures; accepted by the Astrophysical Journal

Submitted: **2014-03-12**, last modified: **2014-07-27**

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.

[4] oai:arXiv.org:1404.5306 [pdf] - 1771440

The relation between gas density and velocity power spectra in galaxy
clusters: qualitative treatment and cosmological simulations

Comments: 6 pages, 3 figures, submitted to ApJ Letters

Submitted: **2014-04-21**

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.

[5] oai:arXiv.org:1307.0021 [pdf] - 796476

The Efficiency of Second-Order Fermi Acceleration by Weakly Compressible
MHD Turbulence

Comments: 12 pages, 9 figures; accepted for publication in Astrophysical
Journal

Submitted: **2013-06-28**, last modified: **2013-08-23**

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.

[6] oai:arXiv.org:1212.0543 [pdf] - 1158235

What Sets Temperature Gradients in Galaxy Clusters? Implications for
non-thermal pressure support and mass-observable scaling relations

Comments: accepted for publication in MNRAS

Submitted: **2012-12-03**, last modified: **2013-03-15**

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.

[7] oai:arXiv.org:1105.2563 [pdf] - 598094

Thermal Instability in Gravitationally-Stratified Plasmas: Implications
for Multi-Phase Structure in Clusters and Galaxy Halos

Comments: This is an electronic version of an article published in Monthly
Notices of the Royal Astronomical Society, Volume 419, Issue 4, pp. 3319-3337

Submitted: **2011-05-12**, last modified: **2012-12-04**

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.

[8] oai:arXiv.org:1206.4314 [pdf] - 1124241

On the Structure of Hot Gas in Halos: Implications for the Lx-Tx
Relation & Missing Baryons

Comments: 13 pages, 4 figures; accepted in MNRAS; a new figure and changes
based on referee's report; critical density using H0=70 km/s/Mpc instead of
the earlier value based on H0=100 km/s/Mpc

Submitted: **2012-06-19**, last modified: **2012-09-16**

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.

[9] oai:arXiv.org:1204.0155 [pdf] - 710350

Resonance Broadening and Heating of Charged Particles in
Magnetohydrodynamic Turbulence

Comments: 20 pages, 15 figures; accepted by The Astrophysical Journal; added
clarifying appendices, but no major changes to results

Submitted: **2012-03-31**, last modified: **2012-09-06**

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.

[10] oai:arXiv.org:1110.5875 [pdf] - 1085160

X-Ray surface brightness and gas density fluctuations in the Coma
cluster

Comments: published in MNRAS

Submitted: **2011-10-26**, last modified: **2012-04-07**

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.

[11] oai:arXiv.org:1202.2117 [pdf] - 529349

The Dynamics of Rayleigh-Taylor Stable and Unstable Contact
Discontinuities with Anisotropic Thermal Conduction

Comments: 18 pages, 15 figures, submitted to MNRAS

Submitted: **2012-02-09**

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.

[12] oai:arXiv.org:1201.0754 [pdf] - 483873

The Effects of Anisotropic Viscosity on Turbulence and Heat Transport in
the Intracluster Medium

Comments: Submitted to MNRAS. Comments are welcome

Submitted: **2012-01-03**

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.

[13] oai:arXiv.org:1109.1285 [pdf] - 1083881

Turbulent Pressure Support in the Outer Parts of Galaxy Clusters

Comments: Accepted to MNRAS Letters. In Press

Submitted: **2011-09-06**, last modified: **2011-10-07**

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.

[14] oai:arXiv.org:1105.3015 [pdf] - 447588

An unusual stellar death on Christmas Day

Comments: 41 pages, 21 figures, accepted for publication in the original
journal!

Submitted: **2011-05-16**, last modified: **2011-10-03**

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).

[15] oai:arXiv.org:1106.4816 [pdf] - 1077505

Thermal Instability & the Feedback Regulation of Hot Halos in Clusters,
Groups, and Galaxies

Comments: 26 pages, 18 figures; abstract abridged; to be submitted to MNRAS

Submitted: **2011-06-23**

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.

[16] oai:arXiv.org:1010.2277 [pdf] - 1041232

Cosmological MHD simulations of cluster formation with anisotropic
thermal conduction

Comments: submitted to ApJ, higher resolution figures available at:
http://www.astro.lsa.umich.edu/~mateuszr/

Submitted: **2010-10-11**

(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.

[17] oai:arXiv.org:1009.2498 [pdf] - 598093

Can conduction induce convection? The non-linear saturation of buoyancy
instabilities in dilute plasmas

Comments: 18 pages, 14 figures. Submitted to MNRAS

Submitted: **2010-09-13**

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.

[18] oai:arXiv.org:0912.1851 [pdf] - 902689

Turbulence in Galaxy Cluster Cores: a Key to Cluster Bimodality?

Comments: Published in ApJ Letters

Submitted: **2009-12-09**, last modified: **2010-03-26**

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.

[19] oai:arXiv.org:1003.5546 [pdf] - 1025985

Thermal Instability with Anisotropic Thermal Conduction and Adiabatic
Cosmic Rays: Implications for Cold Filaments in Galaxy Clusters

Comments: submitted to ApJ; 13 figs. 31 pages; abstract shortened; figures
reduced in size; see http://astro.berkeley.edu/~psharma/TI-v6.pdf for a copy
with high resolution figures

Submitted: **2010-03-25**

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.

[20] oai:arXiv.org:0909.0270 [pdf] - 1017333

Turbulence and Mixing in the Intracluster Medium

Comments: 8 pages, 2 figs., submitted to the conference proceedings of "The
Monster's Fiery Breath;" a follow up of arXiv:0901.4786 focusing on the
general mixing properties of the ICM

Submitted: **2009-09-01**

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.

[21] oai:arXiv.org:0908.4078 [pdf] - 1017267

The Heating of Test Particles in Numerical Simulations of Alfvenic
Turbulence

Comments: Submitted to ApJ

Submitted: **2009-08-28**

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.

[22] oai:arXiv.org:0905.4508 [pdf] - 24700

Simulations of MHD Instabilities in Intracluster Medium Including
Anisotropic Thermal Conduction

Comments: 16 pages, 3 tables, 17 figures, accepted to ApJ with minor revisions,
to appear in Volume 704, Oct 20, 2009 issue

Submitted: **2009-05-27**, last modified: **2009-08-26**

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.

[23] oai:arXiv.org:0905.4500 [pdf] - 24694

Anisotropic Thermal Conduction and the Cooling Flow Problem in Galaxy
Clusters

Comments: Submitted to ApJ, 14 pages, 14 figures

Submitted: **2009-05-27**

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.

[24] oai:arXiv.org:0901.4786 [pdf] - 900592

Buoyancy Instabilities in Galaxy Clusters: Convection Due to Adiabatic
Cosmic Rays and Anisotropic Thermal Conduction

Comments: submitted to ApJ; 15 pages and 12 figures; abstract shortened to < 24
lines; for high resolution movies see
http://astro.berkeley.edu/~psharma/clustermovie.html

Submitted: **2009-01-29**

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.

[25] oai:arXiv.org:0806.2148 [pdf] - 13511

The Magnetothermal Instability in the Intracluster Medium

Comments: Submitted to ApJ

Submitted: **2008-06-12**

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.

[26] oai:arXiv.org:0712.3048 [pdf] - 8287

Nonlinear Simulations of the Heat Flux Driven Buoyancy Instability and
its Implications for Galaxy Clusters

Comments: 4 pages, 4 figures, submitted to ApJ Letters

Submitted: **2007-12-19**

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.

[27] oai:arXiv.org:astro-ph/0703097 [pdf] - 89906

Coulomb Bubbles: Over-stable Driving of Magnetoacoustic Waves Due to the
Rapid and Anisotropic Diffusion of Energy

Comments: 16 pages. Submitted to ApJ

Submitted: **2007-03-06**

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)

[28] oai:arXiv.org:astro-ph/0612195 [pdf] - 87533

Saturation of the Magnetothermal Instability in Three Dimensions

Comments: 15 pages, 15 figures, submitted to ApJ. High resolution figures and
movies are available at http://www.astro.princeton.edu/~iparrish/MTI3D/

Submitted: **2006-12-07**

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.

[29] oai:arXiv.org:astro-ph/0609732 [pdf] - 85324

Simulation of the Magnetothermal Instability

Comments: Accepted for publication in Astrophysics and Space Science as
proceedings of the 6th High Energy Density Laboratory Astrophysics (HEDLA)
Conference

Submitted: **2006-09-26**

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.

[30] oai:arXiv.org:astro-ph/0507212 [pdf] - 74381

Nonlinear Evolution of the Magnetothermal Instability in Two Dimensions

Comments: 37 pages, 17 figures, accepted for publication in ApJ Corrected sign
typo in instability criterion

Submitted: **2005-07-08**, last modified: **2005-07-13**

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.