Search for: All records

Abstract Turbulent radiative mixing layers play an important role in many astrophysical contexts where cool (≲10⁴K) clouds interact with hot flows (e.g., galactic winds, high-velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; however, turbulence in these multiphase flows remains poorly understood. We have investigated the emergent turbulence arising in the interaction between clouds and supersonic winds in hydrodynamicenzo-esimulations. In order to obtain robust results, we employed multiple metrics to characterize the turbulent velocity,v_turb. We find four primary results when cooling is sufficient for cloud survival. First,v_turbmanifests clear temperature dependence. Initially,v_turbroughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time untilv_turbis constant. Second, the relative velocity between the cloud and wind initially drives rapid growth ofv_turb. As it drops (from entrainment),v_turbstarts to decay before it stabilizes at roughly half its maximum. At late times, cooling flows appear to support turbulence. Third, the magnitude ofv_turbscales with the ratio between the hot phase sound-crossing time and the minimum cooling time. Finally, we find tentative evidence for a length scale associated with resolving turbulence. Underresolving this scale may cause violent shattering and affect the cloud’s large-scale morphological properties. 

Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over $Z\prime \equiv Z/Z_{\odot }\sim 0.1-3$, gas surface density Σ_gas∼ 5–150M_⊙pc⁻², and stellar surface density Σ_star∼ 1–50M_⊙pc⁻². The range of emergent star formation rate surface density is Σ_SFR∼ 10⁻⁴–0.5M_⊙kpc⁻²yr⁻¹, and ISM total midplane pressure isP_tot/k_B= 10³–10⁶cm⁻³K, withP_totequal to the ISM weight $\mathcal{W}$. For given Σ_gasand Σ_star, we find ${\mathrm{\Sigma }}_{\mathrm{SFR}}\propto Z{\prime }^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to Σ_SFR. The thermal feedback yield varies sensitively with both weight and metallicity as ${\mathrm{\Upsilon }}_{\mathrm{th}}\propto {\mathcal{W}}^{-0.46}Z{\prime }^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence ${\mathrm{\Upsilon }}_{\mathrm{turb}+\mathrm{mag}}\propto {\mathcal{W}}^{-0.22}Z{\prime }^{-0.18}$. The reduction in Σ_SFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 

Abstract We make an in-depth analysis of different active galactic nuclei (AGN) jet models’ signatures, inducing quiescence in galaxies with a halo mass of 10¹²M_⊙. Three jet models, including cosmic-ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. Each of our simulations is idealized isolated galaxy simulations with AGN jet powers that are constant in time and generated using GIZMO and with FIRE stellar feedback. We examine the distribution of Mgii, Ovi, and Oviiiions, alongside gas temperature and density profiles. Low-energy ions, like Mgii, concentrate in the interstellar medium (ISM), while higher energy ions, e.g., Oviii, prevail at the AGN jet cocoon’s edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic-ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic-ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in Oviand Oviiiis demonstrated in the ISM-to-circumgalactic medium (CGM) transition (between 10 and 30 kpc), growing smoothly toward the CGM at greater distances. 

ABSTRACT In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. Active galactic nucleus jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work concludes from simulations of a halo with $$10^{14} \,\mathrm{ M}_\odot$$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $$10^{12}-10^{15}\, {\rm M}_{\odot }$$ haloes using high-resolution, non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without violating observational constraints. On the contrary, if we scale the energy flux based on the total cooling rate within the cooling radius, strong interstellar medium cooling dominates this scaling, resulting in a jet flux exceeding what is needed. Among the three jet types, the CR-dominant jet is most effective in suppressing cooling flows across all surveyed halo masses due to enhanced CR pressure support. We confirm that the criteria for a successful jet model work across a wider range, encompassing halo masses of $$10^{12}-10^{15} {\rm M_\odot }$$. 

ABSTRACT The origin of the cold phase in the circumgalactic medium (CGM) is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the CGM of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressure stripping and induced cooling in the mixing layer of the stripped cold gas. These two mechanisms contribute a comparable amount of cold gas to the host CGM. We find that the less massive satellites (≤109M⊙) not only lose all of their cold gas in a short period (∼ 0.5–1 Gyr), but their stripped cold clouds also mix with the hot CGM gas and get heated up quickly. However, stellar feedback from these less massive satellites can hugely alter the fate of their stripped gas. Feedback speeds up the destruction of the stripped cold clouds from these satellites by making them more diffuse with more surface area. On the other hand, the more massive satellites (LMC or SMC-like ∼1010M⊙) can add cold gas to the total gas budget of the host CGM for several Gyr. 

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

Abstract The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develop across a broad range of scales in the turbulent interstellar medium (ISM) remain unclear. Using unprecedentedly high-resolution 3D magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated reconnection form regularly throughout the volume. The plasmoids form in three distinct environments: (i) within cold clumps, (ii) at the asymmetric interface of the cold and warm phases, and (iii) within the warm, volume-filling phase. We then show that the complex magnetothermal phase structure is characterized by a predominantly highly magnetized cold phase, but that regions of high magnetic curvature, which are the sites of reconnection, span a broad range in temperature. Furthermore, we show that thermal instabilities change the scale-dependent anisotropy of the turbulent magnetic field, reducing the increase in eddy elongation at smaller scales. Finally, we show that most of the mass is contained in one contiguous cold structure surrounded by smaller clumps that follow a scale-free mass distribution. These clumps tend to be highly elongated and exhibit a size versus internal velocity relation consistent with supersonic turbulence and a relative clump distance–velocity scaling consistent with subsonic motion. We discuss the striking similarity of cold plasmoids to observed tiny-scale atomic and ionized structures and H i fibers and consider how the presence of plasmoids will modify the motion of charged particles, thereby impacting cosmic-ray transport and thermal conduction in the ISM and other similar systems. 

Abstract The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM coevolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (“thermalization”) from a cool, turbulence-supported phase to a virial-temperature, thermally supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high redshift in ultrafaint progenitors and at low redshift in classicalM_vir∼ 10¹¹M_⊙dwarfs, while Milky Way–mass halos undergo the transition atz≈ 0.5. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multiphase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions. 

Abstract Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy’s star-forming efficiency. Mass, energy, and metal outflows (η_M,η_E, andη_Z, here normalized by the star formation rate, the SNe energy, and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas, which enters the interstellar medium (ISM), although whether such high mass loadings agree with observations is still unclear. To better understand how the relative importance of ejective (i.e., high mass loading) versus preventative (i.e., high energy loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from 10¹⁰to 10¹²M_⊙, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low-mass-loaded (η_M∼ 0.1–10) but high-energy-loaded (η_E∼ 0.1–1) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass loading of outflows but are quite sensitive to our choice of the energy loading, preferringη_E∼ 1 for the lowest-mass halos and ∼0.1 for Milky Way–like halos.