Search for: All records

Abstract We present new simulations of local star-forming disks that self-consistently evolve cosmic rays (CRs) and multiphase gas using Tigress++. To isolate the role of CRs, we conduct paired simulations under solar-neighborhood conditions: a magnetohydrodynamics (MHD) model following the standardTigress-classicframework with far-UV heating and supernova (SN) feedback from star clusters formed via gravitational collapse, and a CR + MHD (CRMHD) model in which an additional 10% of each SN’s energy is injected as CRs. These CRs are transported anisotropically along magnetic field lines via a two-moment solver, with the CR scattering rate set by balancing Alfvén-wave growth and damping based on the self-confinement paradigm. The CRMHD model develops a characteristic two-zone vertical CR profile: uniform pressure in the diffusion-dominated, high-density midplane gas, and an exponential atmosphere shaped primarily by advection and streaming in low-density extraplanar gas. The CR pressure is comparable to the total thermal gas pressure in the midplane, but is too uniform to affect gas dynamics, leaving the star formation rates unchanged. In contrast, the vertical CR pressure gradient at ∣z∣ ≳ 1 kpc accelerates warm outflowing gas, resulting in an approximately 4 times higher mass loading factor than in the MHD model. CR–gas interactions increase the CR energy near the midplane through compressional work, while CR streaming heats low-density warm–hot gas. Tigress++opens a path toward investigating CR transport and CR-regulated interstellar medium and outflow dynamics at high resolution across diverse galactic environments. 

Abstract We use models of spectrally resolved cosmic-ray (CR) transport in TIGRESS MHD simulations of the local interstellar medium to produce synthetic synchrotron emission and to test, on scales from a few kiloparsecs down to ∼10 pc, the traditional estimate of magnetic field strength based on the assumption of equipartition between the magnetic and total CR energy densities. Our analysis shows that the traditional equipartition estimate works well at the kiloparsec scale of the simulation box, but breaks down at smaller scales. We find that the predicted magnetic field strength can be improved at small scales by assuming a constant CR energy density across each mock radio observation. The large-scale mean CR energy density can be estimated by assuming equipartition with the large-scale mean magnetic energy density, or as a function of additional observable quantities such as the star formation rate surface density or gas weight. In addition to estimating the magnetic field strength, we use synthetic polarized emission to create maps of the magnetic field direction. We find that the true magnetic field direction can be recovered well from the mock observations. 

Abstract We investigate the transport of spectrally resolved cosmic-ray (CR) protons with kinetic energies between 1 and 100 GeV within the multiphase, dynamic interstellar medium (ISM), using a two-moment CR fluid solver applied to a TIGRESS MHD simulation with conditions similar to the solar neighborhood. Our CR transport prescription incorporates space- and momentum-dependent CR scattering coefficientsσ = κ⁻¹, computed from the local balance between streaming-driven Alfvèn wave growth and damping processes. We find that advection combines with momentum-dependent diffusion to produce a CR distribution functionf(p) ∝ p^−γwithγ ≈ 4.6 that agrees with observations, steepened from an injected power-law slopeγ_inj = 4.3. The CR pressure is uniform in the highly diffusive, mostly neutral midplane region, but decreases exponentially in the ionized extraplanar region where scattering is efficient. To interpret these numerical results, we develop a two-zone analytic model that captures and links the two (physically and spatially) distinct regimes of CR transport in the multiphase, dynamic ISM. At low momenta, CR transport is dominated by gas advection, while at high momenta, both advection and diffusion contribute. At high momentum, the analytic prediction for the spectral slope approachesγ = (4/3)γ_inj − 1, and the predicted scaling of grammage with momentum is $X\propto p^{1-{\gamma }_{\mathrm{inj}}/3}$ , consistent with the simulations. These results support a physical picture in which CRs are confined within the neutral midplane by the surrounding ionized gas, with their escape regulated by both the CR scattering rate in the ionized extraplanar gas and the velocity and Alfvén speed of that gas, at effective speed $v_{\mathrm{c},\mathrm{eff}}\approx (1/2){\left[{\kappa }_{\parallel }d\right(v+v_{\mathrm{A},\mathrm{i}})/dz]}^{1/2}$ . 

Abstract We model the transport and spectral evolution of 1–100 GeV cosmic-ray electrons (CREs) in TIGRESS magnetohydrodynamic simulations of the magnetized, multiphase interstellar medium. We postprocess a kpc-sized galactic disk patch representative of the solar neighborhood using a two-moment method for cosmic ray (CR) transport that includes advection, streaming, and diffusion. The diffusion coefficient is set by balancing wave growth via the CR streaming instability against wave damping (nonlinear Landau and ion–neutral collisions), depending on local gas and CR properties. Implemented energy loss mechanisms include synchrotron, inverse Compton, ionization, and bremsstrahlung. We evaluate CRE losses by different mechanisms as a function of energy and distance from the midplane, and compare loss timescales to transport and diffusion timescales. This comparison shows that CRE spectral steepening abovep= 1 GeV c⁻¹is due to a combination of energy-dependent transport and losses. Our evolved CRE spectra are consistent with direct observations in the solar neighborhood, with a spectral index that steepens from an injected value of −2.3 to an energy-dependent value between −2.7 and −3.3. We also show that the steepening is independent of the injection spectrum. Finally, we present potential applications of our models, including to the production of synthetic synchrotron emission. Our simulations demonstrate that the CRE spectral slope can be accurately recovered from pairs of radio observations in the range 1.5–45 GHz. 

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

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

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

Abstract We study the propagation of mildly relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta et al. to postprocess three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the postprocessed simulations to make predictions regarding the potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extraplanar region (∣ z ∣ > 0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR. 

ABSTRACT The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (∼0.01−0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momentum of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ∼0.01 pc scales with <0.05 km s−1 velocity resolution using the dense gas tracer N2H+. Here, we present the first data sets from DiSCo towards the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core’s internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first data sets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample.