Search for: All records

Massive elliptical galaxies harbor large amounts of hot gas (T≳ 10⁶K) in their interstellar medium (ISM) but are typically quiescent in star formation. The jets of active galactic nuclei (AGNs) and Type Ia supernovae (SNe Ia) inject energy into the ISM, which offsets its radiative losses and keeps it hot. SNe Ia deposit their energy locally within the galaxy compared to the larger few ×10 kiloparsec-scale AGN jets. In this study, we perform high-resolution (512³) hydrodynamic simulations of a local (1 kpc³) density-stratified patch of the ISM of massive galaxies. We include radiative cooling and shell-averaged volume heating, as well as randomly exploding SN Ia. We study the effect of different fractions of supernova (SN) heating (with respect to the net cooling rate), different initial ISM density/entropy (which controls the growth timet_tiof the thermal instability), and different degrees of stratification (which affect the freefall timet_ff). We find that SNe Ia drive predominantly compressive turbulence in the ISM with a velocity dispersion ofσ_vup to 40 km s⁻¹and logarithmic density dispersion ofσ_s∼ 0.2–0.4. These fluctuations trigger multiphase condensation in regions of the ISM, where$\min (t_{\mathrm{ti}})/t_{\mathrm{ff}}\lesssim 0.6\exp (6{\sigma }_s)$, in agreement with theoretical expectations that large density fluctuations efficiently trigger multiphase gas formation. Since the SN Ia rate is not self-adjusting, when the net cooling drops below the net heating rate, SNe Ia drive a hot wind which sweeps out most of the mass in our local model. Global simulations are required to assess the ultimate fate of this gas.

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

Observed accretion rates onto the Milky Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is a rotating inflow in the volume-filling hot phase ($\sim 10^6\, {\rm K}$) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic, and quasi-spherical. Within the radius of angular momentum support ($\sim 15\, {\rm kpc}$ for the Milky Way) the hot flow flattens into a disc geometry and then cools from $\sim 10^6$ to $\sim 10^4\, {\rm K}$ at the disc–halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with ‘precipitation’ models for galaxy accretion. Prior to cooling and accretion the inflow completes ≈tcool/tff radians of rotation, where tcool/tff is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio tcool/tff may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We show that if rotating hot inflows are common in Milky-Way-size disc galaxies, as predicted, then signatures of the expected hot gas rotation profile should be observable with X-ray telescopes and fast radio burst surveys.

Cosmic rays (CRs) with energies ≪ TeV comprise a significant component of the interstellar medium (ISM). Major uncertainties in CR behaviour on observable scales (much larger than CR gyroradii) stem from how magnetic fluctuations scatter CRs in pitch angle. Traditional first-principles models, which assume these magnetic fluctuations are weak and uniformly scatter CRs in a homogeneous ISM, struggle to reproduce basic observables such as the dependence of CR residence times and scattering rates on rigidity. We therefore explore a new category of ‘patchy’ CR scattering models, wherein CRs are pre-dominantly scattered by intermittent strong scattering structures with small volume-filling factors. These models produce the observed rigidity dependence with a simple size distribution constraint, such that larger scattering structures are rarer but can scatter a wider range of CR energies. To reproduce the empirically inferred CR scattering rates, the mean free path between scattering structures must be $\ell _{\rm mfp}\sim 10\, {\rm pc}$ at GeV energies. We derive constraints on the sizes, internal properties, mass/volume-filling factors, and the number density any such structures would need to be both physically and observationally consistent. We consider a range of candidate structures, both large scale (e.g. H ii regions) and small scale (e.g. intermittent turbulent structures, perhaps even associated with radio plasma scattering) and show that while many macroscopic candidates can be immediately ruled out as the primary CR scattering sites, many smaller structures remain viable and merit further theoretical study. We discuss future observational constraints that could test these models.

Low-collisionality plasma in a magnetic field generically develops anisotropy in its distribution function with respect to the magnetic field direction. Motivated by the application to radiation from accretion flows and jets, we explore the effect of temperature anisotropy on synchrotron emission. We derive analytically and provide numerical fits for the polarized synchrotron emission and absorption coefficients for a relativistic bi-Maxwellian plasma (we do not consider Faraday conversion/rotation). Temperature anisotropy can significantly change how the synchrotron emission and absorption coefficients depend on observing angle with respect to the magnetic field. The emitted linear polarization fraction does not depend strongly on anisotropy, while the emitted circular polarization does. We apply our results to black hole imaging of Sgr A* and M87* by ray tracing a GRMHD simulation and assuming that the plasma temperature anisotropy is set by the thresholds of kinetic-scale anisotropy-driven instabilities. We find that the azimuthal asymmetry of the 230 GHz images can change by up to a factor of 3, accentuating (T_⊥>T_∥) or counteracting (T_⊥<T_∥) the image asymmetry produced by Doppler beaming. This can change the physical inferences from observations relative to models with an isotropic distribution function, e.g., by allowing for larger inclination between the line of sight and spin direction in Sgr A*. The observed image diameter and the size of the black hole shadow can also vary significantly due to plasma temperature anisotropy. We describe how the anisotropy of the plasma can affect future multifrequency and photon ring observations. We also calculate kinetic anisotropy-driven instabilities (mirror, whistler, and firehose) for relativistically hot plasmas.

We use local stratified shearing-box simulations to elucidate the impact of two-temperature thermodynamics on the thermal structure of coronae in radiatively efficient accretion flows. Rather than treating the coronal plasma as an isothermal fluid, we use a simple, parametrized cooling function that models the collisional transfer of energy from the ions to the rapidly cooling leptons. Two-temperature models naturally form temperature inversions, with a hot, magnetically dominated corona surrounding a cold disc. Simulations with net vertical flux (NF) magnetic fields launch powerful magnetocentrifugal winds that would enhance accretion in a global system. The outflow rates are much better converged with increasing box height than analogous isothermal simulations, suggesting that the winds into two-temperature coronae may be sufficiently strong to evaporate a thin disc and form a radiatively inefficient accretion flow under some conditions. We find evidence for multiphase structure in the corona, with broad density and temperature distributions, and we propose criteria for the formation of a multiphase corona. The fraction of cooling in the surface layers of the disc is substantially larger for NF fields compared to zero net-flux configurations, with moderate NF simulations radiating ≳30 per cent of the flow’s total luminosity above two mid-plane scale heights. Our work shows that NF fields may efficiently power the coronae of luminous Seyfert galaxies and quasars, providing compelling motivation for future studies of the heating mechanisms available to NF fields and the interplay of radiation with two-temperature thermodynamics.

Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological fire-2 CR-magnetohydrodynamic simulations of L* galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics-driven hysteresis. SC models show a higher tendency to undergo ‘ejective’ feedback events due to a runaway buildup of CR pressure in dense gas due to the behaviour of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, interstellar medium gas, and B properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L* galaxies, and indicating the prospect for non-thermal radio continuum observations to constrain CR transport physics.

In 1977, Blandford and Znajek showed that the electromagnetic field surrounding a rotating black hole can harvest its spin energy and use it to power a collimated astrophysical jet, such as the one launched from the center of the elliptical galaxy M87. Today, interferometric observations with the Event Horizon Telescope (EHT) are delivering high-resolution, event-horizon-scale, polarimetric images of the supermassive black hole M87* at the jet launching point. These polarimetric images offer an unprecedented window into the electromagnetic field structure around a black hole. In this paper, we show that a simple polarimetric observable—the phase ∠β₂of the second azimuthal Fourier mode of the linear polarization in a near-horizon image—depends on the sign of the electromagnetic energy flux and therefore provides a direct probe of black hole energy extraction. In Boyer–Lindquist coordinates, the Poynting flux for axisymmetric electromagnetic fields is proportional to the productB^ϕB^r. The phase ∠β₂likewise depends on the ratioB^ϕ/B^r, thereby enabling an observer to determine the direction of electromagnetic energy flow in the near-horizon environment experimentally. Data from the 2017 EHT observations of M87* are consistent with electromagnetic energy outflow. Currently envisioned multifrequency observations of M87* will achieve higher dynamic range and angular resolution, and hence deliver measurements of ∠β₂closer to the event horizon as well as better constraints on Faraday rotation. Such observations will enable a definitive test for energy extraction from the black hole M87*.

When collapse of the iron core in a massive red or yellow supergiant does not lead to an energetic supernova, a significant fraction of the convective hydrogen envelope will fall in towards the black hole formed from the collapsing core. The random velocity field in the convective envelope results in finite specific angular momentum in each infalling shell. Using 3D hydrodynamical simulations, we follow the infall of this material to small radii, resolving the circularization radii of the flow. We show that infall of the convective envelope leads to nearly complete envelope ejection in a ≳1048 erg explosion with outflow speeds of ≳200 km s−1. The light curve of such an explosion would show a characteristic, red plateau as the ejecta cools and a hydrogen recombination front recedes through the expanding ejecta. Adopting supernova IIp scalings, the event would have a plateau luminosity of ≳1040 erg s−1 and a duration of several hundreds of days. These events would appear quite similar to luminous red novae with red or yellow supergiant progenitors; some luminous red novae may, in fact, be signposts of black hole formation. The mechanism studied here produces more energetic explosions than the weak shock generated from radiation of neutrino energy during the protoneutron star phase. Because we cannot simulate all the way to the horizon, our results are likely lower limits on the energy and luminosity of transients produced during the collapse of a red or yellow supergiant to form a black hole.

Gas in the central regions of cool-core clusters and other massive haloes has a short cooling time (≲1 Gyr). Theoretical models predict that this gas is susceptible to multiphase condensation, in which cold gas is expected to condense out of the hot phase if the ratio of the thermal instability growth time-scale (tti) to the free-fall time (tff) is tti/tff ≲ 10. The turbulent mixing time tmix is another important time-scale: if tmix is short enough, the fluctuations are mixed before they can cool. In this study, we perform high-resolution (5122 × 768–10242 × 1536 resolution elements) hydrodynamic simulations of turbulence in a stratified medium, including radiative cooling of the gas. We explore the parameter space of tti/tff and tti/tmix relevant to galaxy and cluster haloes. We also study the effect of the steepness of the entropy profile, the strength of turbulent forcing and the nature of turbulent forcing (natural mixture versus compressive modes) on multiphase gas condensation. We find that larger values of tti/tff or tti/tmix generally imply stability against multiphase gas condensation, whereas larger density fluctuations (e.g. due to compressible turbulence) promote multiphase gas condensation. We propose a new criterion min (tti/min (tmix, tff)) ≲ c2 × exp (c1σs) for when the halo becomes multiphase, where σs denotes the amplitude of logarithmic density fluctuations and c1 ≃ 6, c2 ≃ 1.8 from an empirical fit to our results.