Metallicity Dependence of Pressure-regulated Feedback-modulated Star Formation in the TIGRESS-NCR Simulation Suite

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 ′ ≡ Z / Z ⊙ ∼ 0.1 - 3 , gas surface density Σ_gas ∼ 5-150 M _⊙ pc⁻², and stellar surface density Σ_star ∼ 1-50 M _⊙ pc⁻². The range of emergent star formation rate surface density is Σ_SFR ∼ 10⁻⁴-0.5 M _⊙ kpc⁻² yr⁻¹, and ISM total midplane pressure is P _tot/k _B = 10³-10⁶ cm⁻³ K, with P _tot equal to the ISM weight W . For given Σ_gas and Σ_star, we find Σ SFR ∝ Z ′ 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 ϒ th ∝ W − 0.46 Z ′ − 0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒ turb + mag ∝ W − 0.22 Z ′ − 0.18 . The reduction in Σ_SFR at 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.