Source author record

Ian J. Parrish

Ian J. Parrish appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Resonance Broadening and Heating of Charged Particles in Magnetohydrodynamic Turbulence

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Thermal Instability in Gravitationally-Stratified Plasmas: Implications for Multi-Phase Structure in Clusters and Galaxy 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.

preprint2011arXiv

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

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.

preprint2011arXiv

Turbulent Pressure Support in the Outer Parts of Galaxy Clusters

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.

preprint2010arXiv

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

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.

preprint2010arXiv

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

Observations of the cores of nearby galaxy clusters show H$α$ 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.

preprint2010arXiv

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

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.

preprint2009arXiv

The Heating of Test Particles in Numerical Simulations of Alfvenic Turbulence

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.