Pressure Regulated Formation of Molecular Clouds and Stars: The case of the Milky Way

We present a steady-state analytical model for pressure-regulated formation of molecular clouds (MC) and stars (SF) in gaseous galactic disks and apply it to the Milky Way (MW). MC formation depends on midplane interstellar pressure $P_{\text{ISM}}$ and metallicity $Z$, and for galactocentric distances $R\gtrsim 5$ kpc, $P_{\text{ISM}}(R)$ scales approximately linearly with molecular gas surface density ${\mathrm{\Sigma }}_{\mathrm{m}\mathrm{o}\mathrm{l}}(R)$. The molecularization of the cold neutral medium (CNM) is due to the opacity of small dust grains that protect the center of the cloud from dissociating radiation when the column density is ${\mathrm{\Sigma }}_d\ge 5(Z_{\odot }\mathrm{/}Z)M_{\odot }{\text{ pc}}^{-2}$. The H${}_2$ formation rate per hydrogen atom is $F\sim {10}^{-15}(P_{\text{ISM}}\mathrm{/}{P}_{\odot })T_{100}^{-1\mathrm{/}2}{\text{s}}^{-1}$, and the corresponding formation rate per unit area is ${\dot{\mathrm{\Sigma }}}_{\mathrm{m}\mathrm{o}\mathrm{l}}^{+}\sim 5\times {10}^{-2}\left(P_{\text{ISM}}\mathrm{/}{P}_{\odot }\right)T_{100}^{-1\mathrm{/}2}{M}_{\odot }{\text{kpc}}^{-2}{\text{yr}}^{-1}$, where $P_{\odot }$ is the pressure at the solar circle and $T_{100}=T\mathrm{/}100\text{ K}$ is the temperature of the cloud. In equilibrium, this equals the molecular gas destruction rate ${\dot{\mathrm{\Sigma }}}_{\mathrm{m}\mathrm{o}\mathrm{l}}^{-}$ due to SF. Self-gravity sets in when the column density of a cloud reaches ${\mathrm{\Sigma }}_{\mathrm{s}\mathrm{g}}={\mathrm{\Sigma }}_{\mathrm{s}\mathrm{g},\odot }(P_{\text{ISM}}\mathrm{/}{P}_{\odot }{)}^{1\mathrm{/}2}$, with ${\mathrm{\Sigma }}_{\mathrm{s}\mathrm{g},\odot }\sim 30M_{\odot }{\text{pc}}^{-2}$. Given the distribution of $P_{\text{ISM}}(R)$ and $Z(R)$ in the MW, the SF process at $5\lesssim R\lesssim 11$ kpc follows a two-step track: first, MCs form from CNM gas and then they form stars when self-gravity sets in. The resulting SFR surface density is ${\mathrm{\Sigma }}_{\text{SFR}}(R)\approx (1.6-4)\times {10}^{-3}\left(P_{\text{ISM}}\mathrm{/}{P}_{\odot }\right){\text{M}}_{\odot }{\text{kpc}}^{-2}{\text{yr}}^{-1}$ with an average final SF efficiency of ${\epsilon }_{\mathrm{s}\mathrm{f}}\sim (3-8)\times {10}^{-2}$.