Axisymmetric disks of eccentric Kepler orbits are vulnerable to an instability which causes orbits to exponentially grow in inclination, decrease in eccentricity, and cluster in their angle of pericenter. Geometrically, the disk expands to a cone shape which is asymmetric about the mid-plane. In this paper, we describe how secular gravitational torques between individual orbits drive this "inclination instability". We derive growth timescales for a simple two-orbit model using a Gauss $N$-ring code, and generalize our result to larger $N$ systems with $N$-body simulations. We find that two-body relaxation slows the growth of the instability at low $N$ and that angular phase coverage of orbits in the disk is important at higher $N$. As $N \to \infty$, the e-folding timescale converges to that expected from secular theory.

Cold ($T\sim 10^{4} \ \mathrm{K}$) gas is very commonly found in both galactic and cluster halos. There is no clear consensus on its origin. Such gas could be uplifted from the central galaxy by galactic or AGN winds. Alternatively, it could form in situ by thermal instability. Fragmentation into a multi-phase medium has previously been shown in hydrodynamic simulations to take place once $t_\mathrm{cool}/t_\mathrm{ff}$, the ratio of the cooling time to the free-fall time, falls below a threshold value. Here, we use 3D plane-parallel MHD simulations to investigate the influence of magnetic fields. We find that because magnetic tension suppresses buoyant oscillations of condensing gas, it destabilizes all scales below $l_\mathrm{A}^\mathrm{cool} \sim v_\mathrm{A} t_\mathrm{cool}$, enhancing thermal instability. This effect is surprisingly independent of magnetic field orientation or cooling curve shape, and sets in even at very low magnetic field strengths. Magnetic fields critically modify both the amplitude and morphology of thermal instability, with $\delta \rho/\rho \propto \beta^{-1/2}$, where $\beta$ is the ratio of thermal to magnetic pressure. In galactic halos, magnetic fields can render gas throughout the entire halo thermally unstable, and may be an attractive explanation for the ubiquity of cold gas, even in the halos of passive, quenched galaxies.

In some galaxies, the stars orbiting the supermassive black hole take the form of an eccentric nuclear disk, in which every star is on a coherent, apsidally-aligned orbit. The most famous example of an eccentric nuclear disk is the double nucleus of Andromeda, and there is strong evidence for many more in the local universe. Despite their apparent ubiquity however, a dynamical explanation for their longevity has remained a mystery: differential precession should wipe out large-scale apsidal-alignment on a short timescale. Here we identify a new dynamical mechanism which stabilizes eccentric nuclear disks, and explain for first time the negative eccentricity gradient seen in the Andromeda nucleus. The stabilizing mechanism drives oscillations of the eccentricity vectors of individual orbits, both in direction (about the mean body of the disk) and in magnitude. Combined with the negative eccentricity gradient, the eccentricity oscillations push some stars near the inner edge of the disk extremely close to the black hole, potentially leading to tidal disruption events. Order of magnitude calculations predict extremely high rates in recently-formed eccentric nuclear disks ($\sim0.1 - 1$ ${\rm yr}^{-1} {\rm gal}^{-1}$). Unless the stellar disks are replenished, these rates should decrease with time as the disk depletes in mass. If eccentric nuclear disks form during gas-rich major mergers, this may explain the preferential occurrence of tidal disruption events in recently-merged and post-merger (E+A/K+A) galaxies.

Motivated by observational and theoretical work which both suggest very small scale ($\lesssim 1\,$pc) structure in the circum-galactic medium of galaxies and in other environments, we study Lyman-$\alpha$ (Ly$\alpha$) radiative transfer in an extremely clumpy medium with many "clouds" of neutral gas along the line of sight. While previous studies have typically considered radiative transfer through sightlines intercepting $\lesssim 10$ clumps, we explore the limit of a very large number of clumps per sightline (up to $f_{\mathrm{c}} \sim 1000$). Our main finding is that, for covering factors greater than some critical threshold, a multiphase medium behaves similar to a homogeneous medium in terms of the emergent Ly$\alpha$ spectrum. The value of this threshold depends on both the clump column density and on the movement of the clumps. We estimate this threshold analytically and compare our findings to radiative transfer simulations with a range of covering factors, clump column densities, radii, and motions. Our results suggest that (i) the success in fitting observed Ly$\alpha$ spectra using homogeneous "shell models" (and the corresponding failure of multiphase models) hints towards the presence of very small-scale structure in neutral gas, in agreement within a number of other observations; and (ii) the recurrent problems of reproducing realistic line profiles from hydrodynamical simulations may be due to their inability to resolve small-scale structure, which causes simulations to underestimate the effective covering factor of neutral gas clouds.

We use idealized 3D hydrodynamic simulations to study the dynamics and thermal structure of the circumgalactic medium (CGM). Our simulations quantify the role of cooling, stellar feedback driven galactic winds and cosmological gas accretion in setting the properties of the CGM in dark matter haloes ranging from $10^{11}$ to $10^{12}$ M$_\odot$. Our simulations support a conceptual picture in which the properties of the CGM, and the key physics governing it, change markedly near a critical halo mass of M$_{\rm crit} \approx 10^{11.5}$ M$_\odot$. As in calculations without stellar feedback, above M$_{\rm crit}$ halo gas is supported by thermal pressure created in the virial shock. The thermal properties at small radii are regulated by feedback triggered when $t_{\rm cool}/t_{\rm ff}\lesssim10$ in the hot gas. Below M$_{\rm crit}$, however, there is no thermally supported halo and self-regulation at $t_{\rm cool}/t_{\rm ff}\sim10$ does not apply. Instead, the gas is out of hydrostatic equilibrium and largely supported against gravity by bulk flows (turbulence and coherent inflow/outflow) arising from the interaction between cosmological gas inflow and outflowing galactic winds. In these lower mass haloes, the phase structure depends sensitively on the outflows' energy per unit mass and mass-loading, which may allow measurements of the CGM thermal state to constrain the nature of galactic winds. Our simulations account for some of the properties of the multiphase halo gas inferred from quasar absorption line observations, including the presence of significant mass at a wide range of temperatures, and the characteristic OVI and CIV column densities and kinematics. However, we underpredict the neutral hydrogen content of the $z\sim0$ CGM.

We perform the first magnetohydrodynamical simulations of tidal disruptions of stars by supermassive black holes. We consider stars with both tangled and ordered magnetic fields, for both grazing and deeply disruptive encounters. When the star survives disruption, we find its magnetic field amplifies by a factor of up to twenty, but see no evidence for the a self-sustaining dynamo that would yield arbitrary field growth. For stars that do not survive, and within the tidal debris streams produced in partial disruptions, we find that the component of the magnetic field parallel to the direction of stretching along the debris stream only decreases slightly with time, eventually resulting in a stream where the magnetic pressure is in equipartition with the gas. Our results suggest that the returning gas in most (if not all) stellar tidal disruptions is already highly magnetized by the time it returns to the black hole.

Lyman-Alpha (Ly$\alpha$) is the strongest emission line in the Universe and is frequently used to detect and study the most distant galaxies. Because Lya is a resonant line, photons typically scatter prior to escaping; this scattering process complicates the interpretation of Ly$\alpha$ spectra, but also encodes a wealth of information about the structure and kinematics of neutral gas in the galaxy. Modeling the Ly$\alpha$ line therefore allows us to study tiny-scale features of the gas, even in the most distant galaxies. Curiously, observed Ly$\alpha$ spectra can be modeled successfully with very simple, homogeneous geometries (such as an expanding, spherical shell), whereas more realistic, multiphase geometries often fail to reproduce the observed spectra. This seems paradoxical since the gas in galaxies is known to be multiphase. In this Letter, we show that spectra emerging from extremely clumpy geometries with many clouds along the line of sight converge to the predictions from simplified, homogeneous models. We suggest that this resolves the apparent discrepancy, and may provide a way to study the gas structure in galaxies on scales far smaller than can be probed in either cosmological simulations or direct (i.e., spatially-resolved) observations.

We find that clouds of optically-thin, pressure-confined gas are prone to fragmentation as they cool below $\sim10^6$ K. This fragmentation follows the lengthscale $\sim{c}_{\text{s}}\,t_{\text{cool}}$, ultimately reaching very small scales ($\sim{0.1} \text{pc}/n$) as they reach the temperature $\sim10^4$ K at which hydrogen recombines. While this lengthscale depends on the ambient pressure confining the clouds, we find that the column density through an individual fragment $N_{\text{cloudlet}}\sim10^{17} \text{cm}^{-3}$ is essentially independent of environment; this column density represents a characteristic scale for atomic gas at $10^4$ K. We therefore suggest that "clouds" of cold, atomic gas may in fact have the structure of a mist or a fog, composed of tiny fragments dispersed throughout the ambient medium. We show that this scale emerges in hydrodynamic simulations, and that the corresponding increase in the surface area may imply rapid entrainment of cold gas. We also apply it to a number of observational puzzles, including the large covering fraction of diffuse gas in galaxy halos, the broad line widths seen in quasar and AGN spectra, and the entrainment of cold gas in galactic winds. While our simulations make a number of assumptions and thus have associated uncertainties, we show that this characteristic scale is consistent with a number of observations, across a wide range of astrophysical environments. We discuss future steps for testing, improving, and extending our model.

We study the dynamical evolution of the putative gas clouds G1 and G2 recently discovered in the Galactic center. Following earlier studies suggesting that these two clouds are part of a larger gas streamer, we combine their orbits into a single trajectory. Since the gas clouds experience a drag force from background gas, this trajectory is not exactly Keplerian. By assuming the G1 and G2 clouds trace this trajectory, we fit for the drag force they experience and thus extract information about the accretion flow at a distance of thousands of Schwarzschild radii from the black hole. This range of radii is important for theories of black hole accretion, but is currently unconstrained by observations. In this paper we extend our previous work by accounting for radial forces due to possible inflow or outflow of the background gas. Such radial forces drive precession in the orbital plane, allowing a slightly better fit to the G1 and G2 data. This precession delays the pericenter passage of G2 by 4-5 months relative to estimates derived from a Keplerian orbital fit; if it proves possible to identify the pericenter time observationally, this enables an immediate test of whether G1 and G2 are gas clouds part of a larger gas streamer. If G2 is indeed a gas cloud, its closest approach likely occurred in late summer 2014, after many of the observing campaigns monitoring G2's anticipated pericenter passage ended. We discuss how this affects interpretation of the G2 observations.

Disks of bodies orbiting a much more massive central object are extremely common in astrophysics. When the orbits comprising such disks are eccentric, we show they are susceptible to a new dynamical instability. Gravitational forces between bodies in the disk drive exponential growth of their orbital inclinations and clustering in their angles of pericenter, expanding an initially thin disk into a conical shape by giving each orbit an identical 'tilt' with respect to the disk plane. This new instability dynamically produces the unusual distribution of orbits observed for minor planets beyond Neptune, suggesting that the instability has shaped the outer Solar System. It also implies a large initial disk mass (1-10 Earth masses) of scattered bodies at hundreds of AU; we predict increasing numbers of detections of minor planets clustered in their angles of pericenter with high inclinations.

The kinetic energy of a star in orbit about a supermassive black hole is a significant fraction of its rest mass energy when its periapse is comparable to its tidal radius. Upon its destruction, a fraction of this energy is extracted and injected into the stellar debris, half of which becomes unbound from the black hole, with the fastest material moving at $\sim 0.03 c$. In this paper, we present a formalism for determining the fate of these unbound debris streams (UDSs) as they depart from the black hole and interact with the surrounding gas. As the density and velocity varies along the length of a UDS, we find that hydrodynamical drag quickly shapes UDSs into loop-like structures, with the densest portions of the streams leading portions of lower density. As UDSs travel outwards, their drag against the ISM increases quadratically with distance, which causes UDSs to deposit their momentum and energy into the ambient medium before the surrounding shocked ISM has a chance to cool. This sudden injection of $\sim 10^{50}$ erg into the ambient medium generates a Sedov-like unbound debris remnant (UDR) that mimics supernova remnants (SNRs) in energetics and appearance, accelerates particles which will produce cosmic rays and synchrotron emission, and provides momentum feedback into the molecular clouds surrounding a black hole. We estimate that a few of these UDRs might be present within a couple degrees of the Galactic Center masquerading as SNRs, and that the UDR scenario is a plausible explanation for Sgr A East.

We examine the role of thermal conduction and magnetic fields in cores of galaxy clusters through global simulations of the intracluster medium (ICM). In particular, we study the influence of thermal conduction, both isotropic and anisotropic, on the condensation of multiphase gas in cluster cores. Previous hydrodynamic simulations have shown that cold gas condenses out of the hot ICM in thermal balance only when the ratio of the cooling time ($t_{\rm cool}$) and the free-fall time ($t_{\rm ff}$) is less than $\approx 10$. Since thermal conduction is significant in the ICM and it suppresses local cooling at small scales, it is imperative to include thermal conduction in such studies. We find that anisotropic (along local magnetic field lines) thermal conduction does not influence the condensation criterion for a general magnetic geometry, even if thermal conductivity is large. However, with isotropic thermal conduction cold gas condenses only if conduction is suppressed (by a factor $\lesssim 0.3$) with respect to the Spitzer value.

The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.

The massive black hole in our galactic center, Sgr A*, accretes only a small fraction of the gas available at its Bondi radius. The physical processes determining this accretion rate remain unknown, partly due to a lack of observational constraints on the gas at distances between ~10 and ~10$^5$ Schwarzschild radii (Rs) from the black hole. Recent infrared observations identify low-mass gas clouds, G1 and G2, moving on highly eccentric, nearly co-planar orbits through the accretion flow around Sgr A*. Although it is not yet clear whether these objects contain embedded stars, their extended gaseous envelopes evolve independently as gas clouds. In this paper we attempt to use these gas clouds to constrain the properties of the accretion flow at ~10$^3$ Rs. Assuming that G1 and G2 follow the same trajectory, we model the small differences in their orbital parameters as evolution resulting from interaction with the background flow. We find evolution consistent with the G-clouds originating in the clockwise disk. Our analysis enables the first unique determination of the rotation axis of the accretion flow: we localize the rotation axis to within 20 degrees, finding an orientation consistent with the parsec-scale jet identified in x-ray observations and with the circumnuclear disk, a massive torus of molecular gas ~1.5 pc from Sgr A*. This suggests that the gas in the accretion flow comes predominantly from the circumnuclear disk, rather than the winds of stars in the young clockwise disk. This result will be tested by the Event Horizon Telescope within the next year. Our model also makes testable predictions for the orbital evolution of G1 and G2, falsifiable on a 5-10 year timescale.

We present three-dimensional magnetohydrodynamic simulations of magnetized gas clouds accelerated by hot winds. We initialize gas clouds with tangled internal magnetic fields and show that this field suppresses the disruption of the cloud: rather than mixing into the hot wind as found in hydrodynamic simulations, cloud fragments end up co-moving and in pressure equilibrium with their surroundings. We also show that a magnetic field in the hot wind enhances the drag force on the cloud by a factor ~(1+v_A^2/v_wind^2)$, where v_A is the Alfven speed in the wind and v_wind measures the relative speed between the cloud and the wind. We apply this result to gas clouds in several astrophysical contexts, including galaxy clusters, galactic winds, the Galactic center, and the outskirts of the Galactic halo. Our results can explain the prevalence of cool gas in galactic winds and galactic halos and how such cool gas survives in spite of its interaction with hot wind/halo gas. We also predict that drag forces can lead to a deviation from Keplerian orbits for the G2 cloud in the galactic center.

We report the discovery of bent double jets in the elliptical galaxy NGC 1272, a member of the Perseus cluster. The jets have a radius of curvature of ~2 kpc, much smaller than the galaxy half-light radius of ~11 kpc. This bending is likely a result of ram pressure, and indicates that the intracluster gas enters deep within the galaxy and that the interstellar medium of NGC 1272 has been significantly removed. X-ray observations of the surrounding intracluster medium allow us to constrain the pressure within the jet. We find that the standard assumptions of equipartition often used in interpreting other jets underestimate the pressure in the jets by a factor ~30.

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.

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.

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.

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.

The galaxy cluster RX J1347-1145 is one of the most X-ray luminous and most massive clusters known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and cosmology. In this paper we present new high-resolution HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing analysis. Combining this information with shape measurements of weak lensing sources in the central regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map provides the best available quantification of the total mass of the central part of the cluster to date. We compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude that both mass estimates agree extremely well in the observed region, namely within 400 / h_70 kpc of the cluster center. In addition we study the major baryonic components (gas and stars) and hence derive the dark matter distribution in the center of the cluster. We find that the dark matter and baryons are both centered on the BCG within the uncertainties (alignment is better than <10 kpc). We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with both NFW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint the concentration parameter, and hence the inner slope of the dark matter profile.