No Arabic abstract
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 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.
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.
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.
Does star formation proceed in the same way in large spirals such as the Milky Way and in smaller chemically younger galaxies? Earlier work suggests a more rapid transformation of H$_2$ into stars in these objects but (1) a doubt remains about the validity of the H$_2$ mass estimates and (2) there is currently no explanation for why star formation should be more efficient. M~33, a local group spiral with a mass $sim 10$% and a metallicity half that of the Galaxy, represents a first step towards the metal poor Dwarf Galaxies. We have searched for molecular clouds in the outer disk of M~33 and present here a set of detections of both $^{12}$CO and $^{13}$CO, including the only detections (for both lines) beyond the R$_{25}$ radius in a subsolar metallicity galaxy. The spatial resolution enables mass estimates for the clouds and thus a measure of the $N({rm H}_2) / I_{rm CO}$ ratio, which in turn enables a more reliable calculation of the H$_2$ mass. Our estimate for the outer disk of M~33 is $N({rm H}_2) / I_{rm CO(1-0)} sim 5 times 10^{20} ,{rm cm^{-2}/(K{rm km s^{-1}})}$ with an estimated uncertainty of a factor $le 2$. While the $^{12/13}$CO line ratios do not provide a reliable measure of $N({rm H}_2) / I_{rm CO}$, the values we find are slightly greater than Galactic and corroborate a somewhat higher $N({rm H}_2) / I_{rm CO}$ value. Comparing the CO observations with other tracers of the interstellar medium, no reliable means of predicting where CO would be detected was identified. In particular, CO detections were often not directly on local HI or FIR or H$alpha$ peaks, although generally in regions with FIR emission and high HI column density. The results presented here provide support for the quicker transformation of H$_2$ into stars in M~33 than in large local universe spirals.
There is a remarkably tight relation between the observationally inferred dust masses and star-formation rates (SFRs) of SDSS galaxies, Mdust $propto$ SFR$^{1.11}$ (Da Cunha et al. 2010). Here we extend the Mdust-SFR relation to the high end and show that it bends over at very large SFRs (i.e., dust masses are lower than predicted for a given SFR). We identify several distinct evolutionary processes in the diagram: (1) A star-bursting phase in which dust builds up rapidly at early times. The maximum attainable dust mass in this process is the cause of the bend-over of the relation. A high dust-formation efficiency, a bottom-light initial mass function, and negligible supernova shock dust destruction are required to produce sufficiently high dust masses. (2) A quiescent star-forming phase in which the subsequent parallel decline in dust mass and SFR gives rise to the Mdust-SFR relation, through astration and dust destruction. The dust-to-gas ratio is approximately constant along the relation. We show that the power-law slope of the Mdust-SFR relation is inversely proportional to the global Schmidt-Kennicutt law exponent (i.e., $sim 0.9$) in simple chemical evolution models. (3) A quenching phase which causes star formation to drop while the dust mass stays roughly constant or drops proportionally. Combined with merging, these processes, as well as the range in total baryonic mass, give rise to a complex population of the diagram which adds significant scatter to the original Mdust-SFR relation. (4) At very high redshifts, a population of galaxies located significantly below the local relation is predicted.