The question of what regulates star formation is a long standing issue. To investigate this issue, we run simulations of a kiloparsec cube section of a galaxy with three kinds of stellar feedback: the formation of HII regions, the explosion of supernovae, and the UV heating. We show that stellar feedback is sufficient to reduce the averaged star formation rate (SFR) to the level of the Schmidt- Kennicutt law in Milky-Way like galaxies but not in high-redshift gas rich galaxies suggesting that another type of support should be added. We investigate whether an external driving of the turbulence such as the one created by the large galactic scales could diminish the SFR at the observed level. Assuming that the Toomre parameter is close to 1 as suggested by the observations, we infer a typical turbulent forcing that we argue should be applied parallel to the plane of the galactic disc. When this forcing is applied in our simulations, the SFR within our simulations closely follows the Schmidt- Kennicutt relation. We found that the velocity dispersion is strongly anisotropic with the velocity dispersion alongside the galactic plane being up to 10 times larger than the perpendicular velocity.

The interstellar medium (ISM) is a complex non-linear system governed by gravity and magneto-hydrodynamics, as well as radiative, thermodynamical, and chemical processes. Our understanding of it mostly progresses through observations and numerical simulations, and a quantitative comparison between these two approaches requires a generic and comprehensive statistical description. The goal of this paper is to build such a description, with the purpose to permit an efficient comparison independent of any specific prior or model. We start from the Wavelet Scattering Transform (WST), a low-variance statistical description of non-Gaussian processes, developed in data science, that encodes long-range interactions through a hierarchical multiscale approach based on the Wavelet transform. We perform a reduction of the WST through a fit of its angular dependencies, allowing to gather most of the information it contains into a few components whose physical meanings are identified, and that describe, e.g., isotropic and anisotropic behaviours. The result of this paper is the Reduced Wavelet Scattering Transform (RWST), a statistical description with a small number of coefficients that characterizes complex structures arising from non-linear phenomena, free from any specific prior. The RWST coefficients encode moments of order up to four, have reduced variances, and quantify the couplings between scales. To show the efficiency and generality of this description, we apply it successfully to three kinds of processes: fractional Brownian motions, MHD simulations, and Herschel observations in a molecular cloud. With fewer than 100 coefficients when probing 6 scales and 8 angles on 256*256 maps, we were able with the RWST to perform quantitative comparisons, to infer relevant physical properties, and to produce realistic synthetic fields.

A numerical shearing box is used to perform three-dimensional simulations of a 1 kpc stratified cubic box of turbulent and self-gravitating interstellar medium (in a rotating frame) with supernovae and HII feedback. We vary the value of the velocity gradient induced by the shear and the initial value of the galactic magnetic field. Finally the different star formation rates and the properties of the structures associated with this set of simulations are computed. We first confirm that the feedback has a strong limiting effect on star formation. The galactic shear has also a great influence: the higher the shear, the lower the SFR. Taking the value of the velocity gradient in the solar neighbourhood, the SFR is too high compared to the observed Kennicutt law, by a factor approximately three to six. This discrepancy can be solved by arguing that the relevant value of the shear is not the one in the solar neighbourhood, and that in reality the star formation efficiency within clusters is not 100%. Taking into account the fact that star-forming clouds generally lie in spiral arms where the shear can be substantially higher (as probed by galaxy-scale simulations), the SFR is now close to the observed one. Different numerical recipes have been tested for the sink particles, giving a numerical incertitude of a factor of about two on the SFR. Finally we have also estimated the velocity dispersions in our dense clouds and found that they lie below the observed Larson law by a factor of about two. Conclusions. In our simulations, magnetic field, shear, HII regions, and supernovae all contribute significantly to reduce the SFR. In this numerical setup with feedback from supernovae and HII regions and a relevant value of galactic shear, the SFRs are compatible with those observed, with a numerical incertitude factor of about two.