Search for: All records

While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time$t_{\mathrm{damp}}\sim \rho v^2/{\dot{E}}_{\mathrm{CR}}\propto E_{\mathrm{CR}}^{-1}$becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating at${\dot{E}}_{\mathrm{CR}}\sim \widetilde{ϵ}$, where$\widetilde{ϵ}$is the turbulent energy input rate. In that case, almost all of the energy in large-scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticE_CR∼E_gscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e.,$\widetilde{ϵ}\gg \rho v^3/L$.

 

We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As mass dropout from thermal instability causes the background density and hence plasma β ≡ Pg/PB to fall, the CGM becomes globally unstable. At the cool disc-to-hot−halo interface, a sharp drop in density boosts Alfven speeds and CR gradients, driving a transition from diffusive to streaming transport. CR forces and heating strengthen, while countervailing gravitational forces and radiative cooling weaken, resulting in a loss of both hydrostatic and thermal equilibrium. In lower β haloes, CR heating drives a hot, single-phase diffuse wind with velocities v ∝ (theat/tff)−1, which exceeds the escape velocity when theat/tff ≲ 0.4. In higher β haloes, where the Alfven Mach number is higher, CR forces drive multi-phase winds with cool, dense fountain flows and significant turbulence. These flows are CR dominated due to ‘trapping’ of CRs by weak transverse B-fields, and have the highest mass loading factors. Thus, local thermal instability can result in winds or fountain flows where either the heat or momentum input of CRs dominates.

 

Spurred by rich, multiwavelength observations and enabled by new simulations, ranging from cosmological to subparsec scales, the past decade has seen major theoretical progress in our understanding of the circumgalactic medium (CGM). We review key physical processes in the CGM. Our conclusions include the following: ▪ The properties of the CGM depend on a competition between gravity-driven infall and gas cooling. When cooling is slow relative to free fall, the gas is hot (roughly virial temperature), whereas the gas is cold ( T ∼ 10⁴K) when cooling is rapid. ▪ Gas inflows and outflows play crucial roles, as does the cosmological environment. Large-scale structure collimates cold streams and provides angular momentum. Satellite galaxies contribute to the CGM through winds and gas stripping. ▪ In multiphase gas, the hot and cold phases continuously exchange mass, energy, and momentum. The interaction between turbulent mixing and radiative cooling is critical. A broad spectrum of cold gas structures, going down to subparsec scales, arises from fragmentation, coagulation, and condensation onto gas clouds. ▪ Magnetic fields, thermal conduction, and cosmic rays can substantially modify how the cold and hot phases interact, although microphysical uncertainties are presently large. Key open questions for future work include the mutual interplay between small-scale structure and large-scale dynamics, and how the CGM affects the evolution of galaxies.

 

Astrophysical gases such as the interstellar-, circumgalactic-, or intracluster-medium are commonly multiphase, which poses the question of the structure of these systems. While there are many known processes leading to fragmentation of cold gas embedded in a (turbulent) hot medium, in this work, we focus on the reverse process: coagulation. This is often seen in wind-tunnel and shearing layer simulations, where cold gas fragments spontaneously coalesce. Using 2D and 3D hydrodynamical simulations, we find that sufficiently large (≫cstcool), perturbed cold gas clouds develop pulsations which ensure cold gas mass growth over an extended period of time (≫r/cs). This mass growth efficiently accelerates hot gas which in turn can entrain cold droplets, leading to coagulation. The attractive inverse square force between cold gas droplets has interesting parallels with gravity; the ‘monopole’ is surface area rather than mass. We develop a simple analytic model which reproduces our numerical findings.

 

Understanding the survival, growth, and dynamics of cold gas is fundamental to galaxy formation. While there has been a plethora of work on ‘wind tunnel’ simulations that study such cold gas in winds, the infall of this gas under gravity is at least equally important, and fundamentally different since cold gas can never entrain. Instead, velocity shear increases and remains unrelenting. If these clouds are growing, they can experience a drag force due to the accretion of low-momentum gas, which dominates over ram pressure drag. This leads to subvirial terminal velocities, in line with observations. We develop simple analytic theory and predictions based on turbulent radiative mixing layers. We test these scalings in 3D hydrodynamic simulations, both for an artificial constant background and a more realistic stratified background. We find that the survival criterion for infalling gas is more stringent than in a wind, requiring that clouds grow faster than they are destroyed ($t_{\rm grow} \lt 4\, t_{\rm cc}$). This can be translated to a critical pressure, which for Milky Way-like conditions is $P \sim 3000 \, {k}_\mathrm{ B} \, {\rm K}\, {\rm cm}^{-3}$. Cold gas that forms via linear thermal instability (tcool/tff < 1) in planar geometry meets the survival threshold. In stratified environments, larger clouds need only survive infall until cooling becomes effective. We discuss applications to high-velocity clouds and filaments in galaxy clusters.

 

Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasmaβ. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-βenvironments like the interstellar medium but remain unchanged in high-βenvironments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss ratev_A· ∇P_CR/P_CR∼v_A/L₀, due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength overβ∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent spectrum.

 

Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔPc is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures, on ∼kpc scales.

 

Astrophysical gases are commonly multiphase and highly turbulent. In this work, we investigate the survival and growth of cold gas in such a turbulent, multiphase medium using three-dimensional hydrodynamical simulations. Similar to previous work simulating coherent flow (winds), we find that cold gas survives if the cooling time of the mixed gas is shorter than the Kelvin–Helmholtz time of the cold gas clump (with some weak additional Mach number dependence). However, there are important differences. Near the survival threshold, the long-term evolution is highly stochastic, and subject to the existence of sufficiently large clumps. In a turbulent flow, the cold gas continuously fragments, enhancing its surface area. This leads to exponential mass growth, with a growth time given by the geometric mean of the cooling and the mixing time. The fragmentation process leads to a large number of small droplets which follow a scale-free dN/dm ∝ m−2 mass distribution, and dominate the area covering fraction. Thus, whilst survival depends on the presence of large ‘clouds’, these in turn produce a ‘fog’ of smaller droplets tightly coupled to the hot phase which are probed by absorption line spectroscopy. We show with the aid of Monte Carlo simulations that the simulated mass distribution emerges naturally due to the proportional mass growth and the coagulation of droplets. We discuss the implications of our results for convergence criteria of larger scale simulations and observations of the circumgalactic medium.

 

ABSTRACT Cosmic ray (CR)-modified shocks are a demanding test of numerical codes. We use them to test and validate the two-moment method for CR hydrodynamics, as well as characterize the realism of CR shock acceleration in two-fluid simulations which inevitably arises. Previously, numerical codes were unable to incorporate streaming in this demanding regime, and have never been compared against analytic solutions. First, we find a new analytic solution highly discrepant in acceleration efficiency from the standard solution. It arises from bi-directional streaming of CRs away from the subshock, similar to a Zeldovich spike in radiative shocks. Since fewer CRs diffuse back upstream, this favours a much lower acceleration efficiency, typically ${\lesssim}10{{\ \rm per\ cent}}$ (even for Mach number > 10) as opposed to ${\gtrsim}50{{\ \rm per\ cent}}$ found in previous analytic work. At Mach number ≳10, the new solution bifurcates into three branches, with efficient, intermediate, and inefficient CR acceleration. Our two-moment code accurately recovers these solutions across the entire parameter space probed, with no ad hoc closure relations. For generic initial conditions, the inefficient branch is robustly chosen by the code; the intermediate branch is unstable. The preferred branch is very weakly modified by CRs. At high Mach numbers (≳10), the gas jump conditions approach that of a purely hydrodynamic shock, and a sub-grid prescription for thermal injection is required for reasonable acceleration efficiencies ${\sim}10{{\ \rm per\ cent}}$. CR-modified shocks have very long equilibration times (∼1000 diffusion time) required to develop the precursor, which must be resolved by ≳10 cells for convergence. Non-equilibrium effects, poor resolution, and obliquity of the magnetic field all reduce CR acceleration efficiency. Shocks in galaxy-scale simulations will generally contribute little to CR acceleration without sub-grid modification. 

ABSTRACT Active galactic nuclei (AGNs) feedback is responsible for maintaining plasma in global thermal balance in extended haloes of elliptical galaxies and galaxy clusters. Local thermal instability in the hot gas leads to the formation of precipitating cold gas clouds that feed the central supermassive black holes, thus heating the hot gas and maintaining global thermal equilibrium. We perform 3D magnetohydrodynamical (MHD) simulations of self-regulated AGNs feedback in a Perseus-like galaxy cluster with the aim of understanding the impact of the feedback physics on the turbulence properties of the hot and cold phases of the intracluster medium (ICM). We find that, in general, the cold phase velocity structure function (VSF) is steeper than the prediction from Kolmogorov’s theory. We attribute the physical origin of the steeper slope of the cold phase VSF to the driving of turbulent motions primarily by the gravitational acceleration acting on the ballistic clouds. We demonstrate that, in the pure hydrodynamical case, the precipitating cold filaments may be the dominant agent driving turbulence in the hot ICM. The arguments in favour of this hypothesis are that: (i) the cold phase mass dominates over hot gas mass in the inner cool core; (ii) hot and cold gas velocities are spatially correlated; (iii) both the cold and hot phase velocity distributions are radially biased. We show that, in the MHD case, the turbulence in the ambient hot medium (excluding the jet cone regions) can also be driven by the AGN jets. The driving is then facilitated by enhanced coupling due to magnetic fields of the ambient gas and the AGN jets. In the MHD case, turbulence may thus be driven by a combination of AGN jet stirring and filament motions. We conclude that future observations, including those from high spatial and spectral resolution X-ray missions, may help to constrain self-regulated AGN feedback by quantifying the multitemperature VSF in the ICM.