Search for: All records

Arkenstone is a new scheme that allows multiphase, stellar feedback-driven winds to be included in coarse resolution cosmological simulations. The evolution of galactic winds and their subsequent impact on the circumgalactic medium are altered by exchanges of mass, energy, momentum, and metals between their component phases. These exchanges are governed by complex, small-scale physical processes that cannot be resolved in cosmological simulations. In this second presentation paper, we describe Arkenstone’s novel cloud particle approach for modelling unresolvable cool clouds entrained in hot, fast winds. This general framework allows models of the cloud–wind interaction, derived from state-of-the-art high-resolution simulations, to be applied in a large-scale context. In this work, we adopt a cloud evolution model that captures simultaneous cloud mass loss to and gain from the ambient hot phase via turbulent mixing and radiative cooling, respectively. We demonstrate the scheme using non-cosmological idealized simulations of a galaxy with a realistic circumgalactic medium component, using the arepo code. We show that the ability of a high-specific energy wind component to perform preventative feedback may be limited by heavy loading of cool clouds coupled into it. We demonstrate that the diverging evolution of clouds of initially differing masses leads to a complex velocity field for the cool phase and a cloud mass function that varies both spatially and temporally in a non-trivial manner. These latter two phenomena can manifest in the simulation because of our choice of a Lagrangian discretization of the cloud population, in contrast to other proposed schemes.

 

The end of supernova remnant (SNR) evolution is characterized by a so-called “radiative” stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to gamma rays, including radio and gamma-ray brightening by nearly 2 orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array).

 

We present a suite of six high-resolution chemodynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback among the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation Ω =v_circ/Rvia the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with Ω, so that the bulge-dominated galaxies have higher midplane densities, lower velocity dispersions, and higher molecular gas fractions than their main-sequence counterparts. The star formation rate in five out of six galaxies is independent of Ω and is consistent with regulation by the midplane gas pressure alone. In the sixth galaxy, which has the most centrally concentrated bulge and thus the highest Ω, we reproduce dynamical suppression of the star formation efficiency in agreement with observations. This produces a transition away from pressure-regulated star formation.

 

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.

 

Bubbles driven by energy and mass injection from small scales are ubiquitous in astrophysical fluid systems and essential to feedback across multiple scales. In particular, O stars in young clusters produce high-velocity winds that create hot bubbles in the surrounding gas. We demonstrate that the dynamical evolution of these bubbles is critically dependent upon the geometry of their interfaces with their surroundings and the nature of heat transport across these interfaces. These factors together determine the amount of energy that can be lost from the interior through cooling at the interface, which in turn determines the ability of the bubble to do work on its surroundings. We further demonstrate that the scales relevant to physical dissipation across this interface are extremely difficult to resolve in global numerical simulations of bubbles for parameter values of interest. This means the dissipation driving evolution of these bubbles in numerical simulations is often of a numerical nature. We describe the physical and numerical principles that determine the level of dissipation in these simulations; we use this, along with a fractal model for the geometry of the interfaces, to explain differences in convergence behavior between hydrodynamical and magnetohydrodynamical simulations presented here. We additionally derive an expression for momentum as a function of bubble radius expected when the relevant dissipative scales are resolved and show that it still results in efficiently cooled solutions, as postulated in previous work.

 

Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of V. Springel & L. Hernquist (hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially resolved simulations of the turbulent, multiphase, star-forming interstellar medium (ISM) have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local-Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the pressure-regulated feedback-modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the “native” simulation outputs with postprocessed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step toward implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations.

 

Measuring the properties of the cold neutral medium (CNM) in low-metallicity galaxies provides insight into heating and cooling mechanisms in early Universe-like environments. We report detections of two localized atomic neutral hydrogen (Hi) absorption features in NGC 6822, a low-metallicity (0.2 Z⊙) dwarf galaxy in the Local Group. These are the first unambiguous CNM detections in a low-metallicity dwarf galaxy outside the Magellanic Clouds. The Local Group L-Band Survey (LGLBS) enabled these detections due to its high spatial (15 pc for Hi emission) and spectral (0.4 km s−1) resolution. We introduce LGLBS and describe a custom pipeline to search for Hi absorption at high angular resolution and extract associated Hi emission. A detailed Gaussian decomposition and radiative transfer analysis of the NGC 6822 detections reveals five CNM components, with key properties: a mean spin temperature of 32±6 K, a mean CNM column density of 3.1×1020 cm−2, and CNM mass fractions of 0.33 and 0.12 for the two sightlines. Stacking non-detections does not reveal low-level signals below our median optical depth sensitivity of 0.05. One detection intercepts a star-forming region, with the Hi absorption profile encompassing the CO (2−1) emission, indicating coincident molecular gas and a depression in high-resolution Hi emission. We also analyze a nearby sightline with deep, narrow Hi self-absorption dips, where the background warm neutral medium is attenuated by intervening CNM. The association of CNM, CO, and Hα emissions suggests a close link between the colder, denser Hi phase and star formation in NGC 6822. 

Measuring the properties of the cold neutral medium (CNM) in low-metallicity galaxies provides insights into heating and cooling mechanisms in early Universe-like environments. We report detections of two localized atomic neutral hydrogen (Hi) absorption features in NGC 6822, a low-metallicity (0.2Z_⊙) dwarf galaxy in the Local Group. These are the first unambiguous CNM detections in a low-metallicity dwarf galaxy outside the Magellanic Clouds. The Local GroupL-band Survey (LGLBS) enabled these detections, due to its high spatial (15 pc for Hiemission) and spectral (0.4 km s⁻¹) resolution. We introduce LGLBS and describe a custom pipeline for searching for Hiabsorption at high angular resolution and extracting associated Hiemission. A detailed Gaussian decomposition and radiative transfer analysis of the NGC 6822 detections reveals five CNM components, with key properties: a mean spin temperature of 32 ± 6 K, a mean CNM column density of 3.1 × 10²⁰cm⁻², and CNM mass fractions of 0.33 and 0.12 for the two sightlines. Stacking nondetections does not reveal low-level signals below our median optical depth sensitivity of 0.05. One detection intercepts a star-forming region, with the Hiabsorption profile encompassing the CO (2−1) emission, indicating coincident molecular gas and a depression in high-resolution Hiemission. We also analyze a nearby sightline with deep, narrow Hiself-absorption dips, where the background warm neutral medium is attenuated by intervening CNM. The association of CNM, CO, and Hαemissions suggests a close link between the colder, denser Hiphase and star formation in NGC 6822.

 

Abstract We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and they span more than six orders of magnitude in radius, reaching scales similar to that of the black hole horizon. The key physical ingredients are radiative cooling and a phenomenological heating model. We find that the accretion flow takes the form of multiphase gas at radii less than about a kpc. The cold gas accretion includes two dynamically distinct stages: the typical disk stage in which the cold gas resides in a rotationally supported disk, and relatively rare chaotic stages (≲10% of the time) in which the cold gas inflows via chaotic streams. Though cold gas accretion dominates the time-averaged accretion rate at intermediate radii, accretion at the smallest radii is dominated by hot virialized gas at most times. The accretion rate scales with radius as M ̇ ∝ r 1 / 2 when hot gas dominates, and we obtain M ̇ ≃ 10 − 4 – 10 − 3 M ⊙ yr − 1 near the event horizon, similar to what is inferred from EHT observations. The orientation of the cold gas disk can differ significantly on different spatial scales. We propose a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by ∼ ( r g / r Bondi ) 1 / 2 . Our results can also provide more realistic initial conditions for simulations of black hole accretion at the event horizon scale. 

Abstract In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, ultra-diffuse galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 H i -selected UDGs and 21 “normal” dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint spectral energy distribution fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime—although the UDGs in our sample show unusually low star formation rate surface densities given their H i content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the “normal” dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the H i -rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H 2 poor compared even to the overall dwarf population.