No Arabic abstract
The cloud-scale physics of star formation and feedback represent the main uncertainty in galaxy formation studies. Progress is hampered by the limited empirical constraints outside the restricted environment of the Local Group. In particular, the poorly-quantified time evolution of the molecular cloud lifecycle, star formation, and feedback obstructs robust predictions on the scales smaller than the disc scale height that are resolved in modern galaxy formation simulations. We present a new statistical method to derive the evolutionary timeline of molecular clouds and star-forming regions. By quantifying the excess or deficit of the gas-to-stellar flux ratio around peaks of gas or star formation tracer emission, we directly measure the relative rarity of these peaks, which allows us to derive their lifetimes. We present a step-by-step, quantitative description of the method and demonstrate its practical application. The methods accuracy is tested in nearly 300 experiments using simulated galaxy maps, showing that it is capable of constraining the molecular cloud lifetime and feedback time-scale to $<0.1$ dex precision. Access to the evolutionary timeline provides a variety of additional physical quantities, such as the cloud-scale star formation efficiency, the feedback outflow velocity, the mass loading factor, and the feedback energy or momentum coupling efficiencies to the ambient medium. We show that the results are robust for a wide variety of gas and star formation tracers, spatial resolutions, galaxy inclinations, and galaxy sizes. Finally, we demonstrate that our method can be applied out to high redshift ($zlesssim4$) with a feasible time investment on current large-scale observatories. This is a major shift from previous studies that constrained the physics of star formation and feedback in the immediate vicinity of the Sun.
We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale by reconstructing the underlying evolutionary timeline. However, by itself this new method only recovers the relative durations of different evolutionary phases. To enable observational applications, it therefore requires knowledge of an absolute reference time-scale to convert relative time-scales into absolute values. The logical choice for this reference time-scale is the duration over which the star formation rate (SFR) tracer is visible because it can be characterised using stellar population synthesis (SPS) models. In this paper, we calibrate this reference time-scale using synthetic emission maps of several SFR tracers, generated by combining the output from a hydrodynamical disc galaxy simulation with the SPS model SLUG2. We apply our statistical method to obtain self-consistent measurements of each tracers reference time-scale. These include H${alpha}$ and 12 ultraviolet (UV) filters (from GALEX, Swift, and HST), which cover a wavelength range 150-350 nm. At solar metallicity, the measured reference time-scales of H${alpha}$ are ${4.32^{+0.09}_{-0.23}}$ Myr with continuum subtraction, and 6-16 Myr without, where the time-scale increases with filter width. For the UV filters we find 17-33 Myr, nearly monotonically increasing with wavelength. The characteristic time-scale decreases towards higher metallicities, as well as to lower star formation rate surface densities, owing to stellar initial mass function sampling effects. We provide fitting functions for the reference time-scale as a function of metallicity, filter width, or wavelength, to enable observational applications of our statistical method across a wide variety of galaxies.
Recent observational studies aiming to quantify the molecular cloud lifecycle require the use of known reference time-scales to turn the relative durations of different phases of the star formation process into absolute time-scales. We previously constrained the characteristic emission time-scales of different star formation rate (SFR) tracers, as a function of the SFR surface density and metallicity. However, we omitted the effects of dust extinction. Here, we extend our suite of SFR tracer emission time-scales by accounting for extinction, using synthetic emission maps of a high-resolution hydrodynamical simulation of an isolated, Milky-Way-like disc galaxy. The stellar feedback included in the simulation is inefficient compared to observations, implying that it represents a limiting case in which the duration of embedded star formation (and the corresponding effect of extinction) is overestimated. Across our experiments, we find that extinction mostly decreases the SFR tracer emission time-scale, changing the time-scales by factors of 0.04-1.74, depending on the gas column density. UV filters are more strongly affected than H$alpha$ filters. We provide the limiting correction factors as a function of the gas column density and flux sensitivity limit for a wide variety of SFR tracers. Applying these factors to observational characterisations of the molecular cloud lifecycle produces changes that broadly fall within the quoted uncertainties, except at high kpc-scale gas surface densities ($Sigma_{rm g}gtrsim20~{mathrm{M_{odot},pc^{-2}}}$). Under those conditions, correcting for extinction may decrease the measured molecular cloud lifetimes and feedback time-scales, which further strengthens previous conclusions that molecular clouds live for a dynamical time and are dispersed by early, pre-supernova feedback.
We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.
Diffuse emission is observed in galaxies in many tracers across the electromagnetic spectrum, including tracers of star formation, such as H$alpha$ and ultraviolet (UV), and tracers of gas mass, such as carbon monoxide (CO) transition lines and the 21-cm line of atomic hydrogen (HI). Its treatment is key to extracting meaningful information from observations such as cloud-scale star formation rates. Finally, studying diffuse emission can reveal information about the physical processes taking place in the ISM, such as chemical transitions and the nature of stellar feedback (through the photon escape fraction). We present a physically-motivated method for decomposing astronomical images containing both diffuse emission and compact regions of interest, such as HII regions or molecular clouds, into diffuse and compact component images through filtering in Fourier space. We have previously presented a statistical method for constraining the evolutionary timeline of star formation and mean separation length between compact star forming regions with galaxy-scale observations. We demonstrate how these measurements are biased by the presence of diffuse emission in tracer maps and that by using the mean separation length as a critical length scale to separate diffuse emission from compact emission, we are able to filter out this diffuse emission, thus removing its biasing effect. Furthermore, this method provides, without the need for interferometry or ancillary spectral data, a measurement of the diffuse emission fraction in input tracer maps and decomposed diffuse and compact emission maps for further analysis.
We investigate the balance of power between stars and AGN across cosmic history, based on the comparison between the infrared (IR) galaxy luminosity function (LF) and the IR AGN LF. The former corresponds to emission from dust heated by stars and AGN, whereas the latter includes emission from AGN-heated dust only. We find that at all redshifts (at least up to z~2.5), the high luminosity tails of the two LFs converge, indicating that the most infrared-luminous galaxies are AGN-powered. Our results shed light to the decades-old conundrum regarding the flatter high-luminosity slope seen in the IR galaxy LF compared to that in the UV and optical. We attribute this difference to the increasing fraction of AGN-dominated galaxies with increasing total infrared luminosity (L_IR). We partition the L_IR-z parameter space into a star-formation and an AGN-dominated region, finding that the most luminous galaxies at all epochs lie in the AGN-dominated region. This sets a potential `limit to attainable star formation rates, casting doubt on the abundance of `extreme starbursts: if AGN did not exist, L_IR>10^13 Lsun galaxies would be significantly rarer than they currently are in our observable Universe. We also find that AGN affect the average dust temperatures (T_dust) of galaxies and hence the shape of the well-known L_IR-T_dust relation. We propose that the reason why local ULIRGs are hotter than their high redshift counterparts is because of a higher fraction of AGN-dominated galaxies amongst the former group.