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

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 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi>Z</mml:mi> <mml:mo accent="false">′</mml:mo> <mml:mo>≡</mml:mo> <mml:mi>Z</mml:mi> <mml:mo stretchy="true">/</mml:mo> <mml:msub> <mml:mi>Z</mml:mi> <mml:mo>⊙</mml:mo> </mml:msub> <mml:mo>∼</mml:mo> <mml:mn>0.1</mml:mn> <mml:mo>-</mml:mo> <mml:mn>3</mml:mn> </mml:math> , gas surface density Σ gas ∼ 5–150 M ⊙ pc −2 , and stellar surface density Σ star ∼ 1–50 M ⊙ pc −2 . The range of emergent star formation rate surface density is Σ SFR ∼ 10 −4 –0.5 M ⊙ kpc −2 yr −1 , and ISM total midplane pressure is P tot / k B = 10 3 –10 6 cm −3 K, with P tot equal to the ISM weight <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi class="MJX-tex-calligraphic" mathvariant="script">W</mml:mi> </mml:math> . For given Σ gas and Σ star , we find <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal">Σ</mml:mi> </mml:mrow> <mml:mrow> <mml:mi>SFR</mml:mi> </mml:mrow> </mml:msub> <mml:mo>∝</mml:mo> <mml:mi>Z</mml:mi> <mml:msup> <mml:mrow> <mml:mo accent="true">′</mml:mo> </mml:mrow> <mml:mrow> <mml:mn>0.3</mml:mn> </mml:mrow> </mml:msup> </mml:math> . 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 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal">ϒ</mml:mi> </mml:mrow> <mml:mrow> <mml:mi>th</mml:mi> </mml:mrow> </mml:msub> <mml:mo>∝</mml:mo> <mml:msup> <mml:mrow> <mml:mi class="MJX-tex-calligraphic" mathvariant="script">W</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.46</mml:mn> </mml:mrow> </mml:msup> <mml:mi>Z</mml:mi> <mml:msup> <mml:mrow> <mml:mo accent="true">′</mml:mo> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.53</mml:mn> </mml:mrow> </mml:msup> </mml:math> , while the combined turbulent and magnetic feedback yield shows weaker dependence <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal">ϒ</mml:mi> </mml:mrow> <mml:mrow> <mml:mi>turb</mml:mi> <mml:mo>+</mml:mo> <mml:mi>mag</mml:mi> </mml:mrow> </mml:msub> <mml:mo>∝</mml:mo> <mml:msup> <mml:mrow> <mml:mi class="MJX-tex-calligraphic" mathvariant="script">W</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.22</mml:mn> </mml:mrow> </mml:msup> <mml:mi>Z</mml:mi> <mml:msup> <mml:mrow> <mml:mo accent="true">′</mml:mo> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.18</mml:mn> </mml:mrow> </mml:msup> </mml:math> . 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.