No Arabic abstract
We present an analysis of the global and spatially-resolved Kennicutt-Schmidt (KS) star formation relation in the FIRE (Feedback In Realistic Environments) suite of cosmological simulations, including halos with $z = 0$ masses ranging from $10^{10}$ -- $10^{13}$ M$_{odot}$. We show that the KS relation emerges and is robustly maintained due to the effects of feedback on local scales regulating star-forming gas, independent of the particular small-scale star formation prescriptions employed. We demonstrate that the time-averaged KS relation is relatively independent of redshift and spatial averaging scale, and that the star formation rate surface density is weakly dependent on metallicity and inversely dependent on orbital dynamical time. At constant star formation rate surface density, the `Cold & Dense gas surface density (gas with $T < 300$~K and $n > 10$~cm$^{-3}$, used as a proxy for the molecular gas surface density) of the simulated galaxies is $sim$0.5~dex less than observed at $sim$kpc scales. This discrepancy may arise from underestimates of the local column density at the particle-scale for the purposes of shielding in the simulations. Finally, we show that on scales larger than individual giant molecular clouds, the primary condition that determines whether star formation occurs is whether a patch of the galactic disk is thermally Toomre-unstable (not whether it is self-shielding): once a patch can no longer be thermally stabilized against fragmentation, it collapses, becomes self-shielding, cools, and forms stars, regardless of epoch or environment.
The surface densities of molecular gas, $Sigma_{rm H_2}$, and the star formation rate (SFR), $dotSigma_star$, correlate almost linearly on kiloparsec scales in observed star-forming (non-starburst) galaxies. We explore the origin of the linear slope of this correlation using a suite of isolated $L_star$ galaxy simulations. We show that in simulations with efficient feedback, the slope of the $dotSigma_star$-$Sigma_{rm H_2}$ relation on kiloparsec scales is insensitive to the slope of the $dotrho_star$-$rho$ relation assumed at the resolution scale. We also find that the slope on kiloparsec scales depends on the criteria used to identify star-forming gas, with a linear slope arising in simulations that identify star-forming gas using a virial parameter threshold. This behavior can be understood using a simple theoretical model based on conservation of interstellar gas mass as the gas cycles between atomic, molecular, and star-forming states under the influence of feedback and dynamical processes. In particular, we show that the linear slope emerges when feedback efficiently regulates and stirs the evolution of dense, molecular gas. We show that the model also provides insights into the likely origin of the relation between the SFR and molecular gas in real galaxies on different scales.
We use a new method to trace backwards the star formation history of the Milky Way disk, using a sample of M dwarfs in the solar neighbourhood which is representative for the entire solar circle. M stars are used because they show H_alpha emission until a particular age which is a well calibrated function of their absolute magnitudes. This allows us to reconstruct the rate at which disk stars have been born over about half the disks lifetime. Our star formation rate agrees well with those obtained by using other, independent, methods and seems to rule out a constant star formation rate. The principal result of this study is to show that a relation of the Schmidt-Kennicut type (which relates the star formation rate to the interstellar gas content of galaxy disks) has pertained in the Milky Way disk during the last 5 Gyr. The star formation rate we derive from the M dwarfs and the interstellar gas content of the disk can be inferred as a function of time from a model of the chemical enrichment of the disk, which is well constrained by the observations indicating that the metallicity of the Galactic disk has remained nearly constant over the timescales involved. We demonstrate that the star formation rate and gas surface densities over the last 5 Gyrs can be accurately described by a Schmidt-Kennicutt law with an index of Gamma = 1.45 (+0.22,-0.09). This is, within statistical uncertainties, the same value found for other galaxies.
We have mapped the northern area (30 times 20) of a local group spiral galaxy M33 in 12CO(J=1-0) line with the 45-m telescope at the Nobeyama Radio Observatory. Along with Halpha and Spitzer 24-micron data, we have investigated the relationship between the surface density of molecular gas mass and that of star formation rate (SFR) in an external galaxy (Kennicutt-Schmidt law) with the highest spatial resolution (~80pc) to date, which is comparable to scales of giant molecular clouds (GMCs). At positions where CO is significantly detected, the SFR surface density exhibits a wide range of over four orders of magnitude, from Sigma(SFR)<10^{-10} to ~10^{-6}M_solar yr^{-1} pc^{-2}, whereas the Sigma(H2) values are mostly within 10 to 40 M_solar pc^{-2}. The surface density of gas and that of SFR correlate well at a 1-kpc resolution, but the correlation becomes looser with higher resolution and breaks down at GMC scales. The scatter of the Sigma(SFR)-Sigma(H2) relationship in the 80-pc resolution results from the variety of star forming activity among GMCs, which is attributed to the various evolutionary stages of GMCs and to the drift of young clusters from their parent GMCs. This result shows that the Kennicutt-Schmidt law is valid only in scales larger than that of GMCs, when we average the spatial offset between GMCs and star forming regions, and their various evolutionary stages.
We compile observations of molecular gas contents and infrared-based star formation rates (SFRs) for 112 circumnuclear star forming regions, in order to re-investigate the form of the disk-averaged Schmidt surface density star formation law in starbursts. We then combine these results with total gas and SFR surface densities for 153 nearby non-starbursting disk galaxies from de los Reyes & Kennicutt (2019), to investigate the properties of the combined star formation law, following Kennicutt (1998; K98). We confirm that the combined Schmidt law can be fitted with a single power law with slope $n = 1.5pm0.05$ (including fitting method uncertainties), somewhat steeper than the value $n = 1.4pm0.15$ found by K98. Fitting separate power laws to the non-starbursting and starburst galaxies, however, produces very different slopes ($n = 1.34pm0.07$ and $0.98pm0.07$, respectively), with a pronounced offset in the zeropoint ($sim$0.6,dex) of the starburst relation to higher SFR surface densities. This offset is seen even when a common conversion factor between CO intensity and molecular hydrogen surface density is applied, and is confirmed when disk surface densities of interstellar dust are used as proxies for gas measurements. Tests for possible systematic biases in the starburst data fail to uncover any spurious sources for such a large offset. We tentatively conclude that the global Schmidt law in galaxies, at least as it is conventionally measured, is bimodal or possibly multi-modal. Possible causes may include changes in the small-scale structure of the molecular ISM or the stellar initial mass function. A single $n sim 1.5$ power law still remains as a credible approximation or recipe for analytical or numerical models of galaxy formation and evolution.
We address a simple model where the Kennicutt-Schmidt (KS) relation between the macroscopic densities of star-formation rate (SFR, $rho_{rm sfr}$) and gas ($n$) in galactic discs emerges from self-regulation of the SFR via supernova feedback. It arises from the physics of supernova bubbles, insensitive to the microscopic SFR recipe and not explicitly dependent on gravity. The key is that the filling factor of SFR-suppressed supernova bubbles self-regulates to a constant, $fsim 0.5$. Expressing the bubble fading radius and time in terms of $n$, the filling factor is $f propto S,n^{-s}$ with $ssim 1.5$, where $S$ is the supernova rate density. A constant $f$ thus refers to $rho_{rm sfr} propto n^{1.5}$, with a density-independent SFR efficiency per free-fall time $sim 0.01$. The self-regulation to $f sim 0.5$ and the convergence to a KS relation independent of the local SFR recipe are demonstrated in cosmological and isolated-galaxy simulations using different codes and recipes. In parallel, the spherical analysis of bubble evolution is generalized to clustered supernovae, analytically and via simulations, yielding $s simeq 1.5 pm 0.5$. An analysis of photo-ionized bubbles about pre-supernova stars yields a range of KS slopes but the KS relation is dominated by the supernova bubbles. Superbubble blowouts may lead to an alternative self-regulation by outflows and recycling. While the model is over-simplified, its simplicity and validity in the simulations may argue that it captures the origin of the KS relation.