Linking Electron Density with Elevated Star Formation Activity from z = 0 to z = 10

We present a [C ii] 158 μm map of the entire M51 (including M51b) grand design spiral galaxy observed with the Far Infrared Field-Imaging Line Spectrometer (FIFI-LS) instrument on board the Stratospheric Observatory For Infrared Astronomy (SOFIA). We compare the [C ii] emission with the total far-infrared (TIR) intensity and star formation rate (SFR) surface density maps (derived using Hα and 24 μm emission) to study the relationship between [C ii] and the star formation activity in a variety of environments within M51 on scales of 16″ corresponding to ∼660 pc. We find that [C ii] and the SFR surface density are well correlated in the central, spiral arm, and inter-arm regions. The correlation is in good agreement with that found for a larger sample of nearby galaxies at kpc scales. We find that the SFR, and [C ii] and TIR luminosities in M51, are dominated by the extended emission in M51's disk. The companion galaxy M51b, however, shows a deficit of [C ii] emission compared with the TIR emission and SFR surface density, with [C ii] emission detected only in the SW part of this galaxy. The [C ii] deficit is associated with an enhanced dust temperature in this galaxy. We interpret the faint [C ii] emission in M51b to be a result of suppressed star formation in this galaxy, while the bright mid- and far-infrared emission, which drive the TIR and SFR values, are powered by other mechanisms. A similar but less-pronounced effect is seen at the location of the black hole in M51's center. The observed [C ii] deficit in M51b suggests that this galaxy is a valuable laboratory to study the origin of the apparent [C ii] deficit observed in ultra-luminous galaxies.

We have derived the column densities of heavy elements in three gamma-ray burst (GRB) optical transients, associated with the circumburst or interstellar medium (ISM) of the host galaxy. In comparison with the same elements observed in damped Lyα (DLA) systems along QSO sight lines, we find evidence for much higher column densities of Zn II. The gap between the QSO-DLA and GRB-DLA distributions is smoothly bridged by observations of the interstellar absorption in the Milky Way and the Magellanic Clouds. Very small [Fe/Zn], [Si/Zn], and [Cr/Zn] values in GRB-DLAs indicate large dust depletions. Once the dust-to-metals ratios are determined, we find an optical extinction AV ≈ 1 mag, to be compared with typical AV 0.1 in most QSO-DLAs. Our inference of high dust content appears to be in contradiction with the typical low reddening previously found in GRBs. One possible way to reconcile is a dust grain size distribution biased toward big grains, which would give a gray extinction. Possibly, the small dust grains have been destroyed by the GRBs themselves. Our findings support the idea that primarily optically selected QSOs probe mainly low-gas/dust regions of high-redshift galaxies, while the more powerful GRBs can be detected through denser regions of their ISM (molecular clouds and star-forming regions). Therefore, GRB-DLAs and QSO-DLAs together provide a more complete picture of the global properties of the interstellar medium in high-redshift galaxies.

We investigate the nature of the star formation law at low gas surface densities using a sample of 19 low surface brightness (LSB) galaxies with existing HI maps in the literature, UV imaging from the Galaxy Evolution Explorer satellite, and optical images from the Sloan Digital Sky Survey. All of the LSB galaxies have (NUV-r) colors similar to those for higher surface brightness star-forming galaxies of similar luminosity indicating that their average star formation histories are not very different. Based upon four LSB galaxies with both UV and FIR data, we find FIR/UV ratios significantly less than one, implying low amounts of internal UV extinction in LSB galaxies. We use the UV images and HI maps to measure the star formation rate and hydrogen gas surface densities within the same region for all of the galaxies. The LSB galaxy star formation rate surface densities lie below the extrapolation of the power law fit to the star formation rate surface density as a function of the total gas density for higher surface brightness galaxies. Although there is more scatter, the LSB galaxies also lie below a second version of the star formation law in which the star formation rate surface density is correlated with the gas density divided by the orbital time in the disk. The downturn seen in both star formation laws is consistent with theoretical models that predict lower star formation efficiencies in LSB galaxies due to the declining molecular fraction with decreasing density.

Neutral hydrogen (HI) velocity dispersions are believed to be set by turbulence in the interstellar medium (ISM). Although turbulence is widely believed to be driven by star formation (SF), recent studies have shown that this driving mechanism may not be dominant in regions of low SF rate surface density (SFRSD), such as found in dwarf galaxies or the outer regions of spirals. We have generated average HI line profiles in a number of nearby dwarfs and low-mass spirals by co-adding HI spectra in regions with either a common radius or SFRSD. We find that the spatially-resolved superprofiles are composed of a central narrow peak (5-15 km/s) with higher velocity wings to either side. With the assumption that the central peak reflects the turbulent velocity dispersion, we compare HI kinematics to local ISM properties, including surface mass densities and measures of SF. The HI velocity dispersion is correlated most strongly with surface mass density, which points at a gravitational origin for turbulence, but it is unclear which instabilities can operate efficiently in these systems. SF energy is produced at a level sufficient to drive HI turbulent motions where SFRSD > 10^-4 Msun yr^-1 kpc^-2. At low SF intensities, SF does not supply enough energy for turbulence, nor does it uniquely determine the velocity dispersion. Nevertheless, SF appears to provide a lower threshold for HI velocity dispersions. We find that coupling efficiency decreases with increasing SFRSD, consistent with a picture where SF couples to the ISM with constant efficiency, but that less of that energy is found in HI at higher SFRSD. We examine a number of potential drivers of HI turbulence, including SF, gravitational instabilities, the magnetorotational instability, and accretion, and find that no single mechanism can drive the observed levels of turbulence at low SFRSD. We discuss possible solutions to this conundrum.

Star formation cannot truly be understood from observational data alone; only with simulations is it possible to assemble the complete picture. Observations guide the physics we build into our simulations, yet the impact of different star formation and feedback models can only be investigated with simulations. Synthetic observations allow us to make a realistic comparison to true observations as well as teach us about the emission tracers we depend upon. Through coupling the stellar population synthesis code SLUG2 to galaxy simulations, we can generate synthetic star formation rate tracer maps. These maps assume different stellar metallicities, star formation rate surface densities, and suffer from varied amounts of extinction. This allows us to explore and constrain the environmental effects on the characteristic emission lifetimes — the duration for which a tracer is visible. With these emission lifetimes and in conjunction with a new statistical method, the ‘uncertainty principle for star formation’, constraints can be placed upon the durations of different evolutionary phases of the star formation process, allowing us to better understand the physics of star formation and feedback on sub-galactic scales. Studying the interstellar medium can also reveal information about stellar feedback: the gas density structure is altered as a result of the injected energy, momentum, and matter. Surveys of the CO emission in galaxies can tell us how the properties of this medium have evolved over cosmic time. Using DESPOTIC to model CO line emission of gas found within the IllustrisTNG50 cosmological simulation, we produce an equivalent synthetic survey. This synthetic survey can be used as a basis for comparison and predictor of observational trends.

Building on our previous research involving multiwavelength data from UV to IR, we employ spectroscopic observations of ionized gas, as well as neutral hydrogen gas obtained from the Five-hundred Meter Aperture Spherical Telescope (FAST), to explore the intrinsic processes of star formation and chemical enrichment within NGC 628. Our analysis focuses on several key properties, including gas-phase extinction, star formation rate (SFR) surface density, oxygen abundance, and H i mass surface density. The azimuthal distributions of these parameters in relation to the morphological and kinematic features revealed by FAST H i reveal that NGC 628 is an isolated galaxy that has not undergone recent interactions. We observe a mild radial extinction gradient accompanied by a notable dispersion. The SFR surface density also shows a gentle radial gradient, characteristic of typical spiral galaxies. Additionally, we find a negative radial metallicity gradient of −0.44 dex R 25 − 1 , supporting the “inside-out” scenario of galaxy formation. We investigate the resolved mass–metallicity relation (MZR) and the resolved star formation main sequence (SFMS) alongside their dependences on the physical properties of both ionized and neutral hydrogen gas. Our findings indicate no secondary dependence of the resolved MZR on SFR surface density or H i mass surface density. Furthermore, we observe that gas-phase extinction and the equivalent width of Hα both increase with SFR surface density in the resolved SFMS.

We investigate the connection between galactic outflows and star formation using two independent data sets covering a sample of 22 galaxies between 1 ≲ z ≲ 1.5. The Hubble Space Telescope WFC3/G141 grism provides low spectral resolution, high spatial resolution spectroscopy yielding Hα emission-line maps from which we measure the spatial extent and strength of star formation. In the rest-frame near-UV, Keck/DEIMOS observes Fe ii and Mg ii interstellar absorption lines, which provide constraints on the intensity and velocity of the outflows. We compare outflow properties from individual and composite spectra with the star formation rate (SFR) and SFR surface density (ΣSFR), as well as the stellar mass and specific SFR (sSFR). The Fe ii and Mg ii equivalent widths (EWs) increase with both SFR and ΣSFR at ≳3σ significance, while the composite spectra show larger Fe ii EWs and outflow velocities in galaxies with higher SFR, ΣSFR, and sSFR. Absorption-line profiles of the composite spectra further indicate that the differences between subsamples are driven by outflows rather than the interstellar medium. While these results are consistent with those of previous studies, the use of Hα images makes them the most direct test of the relationship between star formation and outflows at z > 1 to date. Future facilities such as the James Webb Space Telescope and the upcoming Extremely Large Telescopes will extend these direct, Hα-based studies to lower masses and SFRs, probing galactic feedback across orders of magnitude in galaxy properties and augmenting the correlations we find here.

Context. Long gamma-ray bursts (GRBs) offer a promising tool for tracing the cosmic history of star formation, especially at high redshift, where conventional methods are known to suffer from intrinsic biases. Previous studies of GRB host galaxies at low redshift showed that high surface density of stellar mass and high surface density of star formation rate (SFR) can potentially enhance the GRB production. Evaluating the effect of such stellar densities at high redshift is therefore crucial to fully control the ability of long GRBs for probing the activity of star formation in the distant Universe. Aims. We assess how the size, stellar mass, and star formation rate surface densities of distant galaxies affect the probability of their hosting a long GRB, using a sample of GRB hosts at z > 1 and a control sample of star-forming sources from the field. Methods. We gathered a sample of 45 GRB host galaxies at 1 < z < 3.1 observed with the Hubble Space Telescope WFC3 camera in the near-infrared. Our subsample at 1 < z < 2 has cumulative distributions of redshift and stellar mass consistent with the host galaxies of known unbiased GRB samples, while our GRB host selection at 2 < z < 3.1 has lower statistics and is probably biased toward the high end of the stellar mass function. Using the GALFIT parametric approach, we modeled the GRB host light profile with a Sérsic component and derived the half-light radius for 35 GRB hosts, which we used to estimate the star formation rate and stellar mass surface densities of each object. We compared the distribution of these physical quantities to the SFR-weighted properties of a complete sample of star-forming galaxies from the 3D-HST deep survey at a comparable redshift and stellar mass. Results. We show that similarly to z < 1, GRB hosts are smaller in size and they have higher stellar mass and star formation rate surface densities than field galaxies at 1 < z < 2. Interestingly, this result is robust even when separately considering the hosts of GRBs with optically bright afterglows and the hosts of dark GRBs, as the two subsamples share similar size distributions. At z > 2, however, GRB hosts appear to have sizes and stellar mass surface densities more consistent with those characterizing the field galaxies. This may reveal an evolution with redshift of the bias between GRB hosts and the overall population of star-forming sources, although we cannot exclude that our result at z > 2 is also affected by the prevalence of dark GRBs in our selection. Conclusions. In addition to a possible trend toward a low-metallicity environment, other environmental properties such as stellar density appear to play a role in the formation of long GRBs, at least up to z ∼ 2. This might suggest that GRBs require special environments to enhance their production.

It is widely known that the gas in galaxy discs is highly turbulent, but there is much debate on which mechanism can energetically maintain this turbulence. Among the possible candidates, supernova (SN) explosions are likely the primary drivers but doubts remain on whether they can be sufficient in regions of moderate star formation activity, in particular in the outer parts of discs. Thus, a number of alternative mechanisms have been proposed. In this paper, we measure the SN efficiencyη, namely the fraction of the total SN energy needed to sustain turbulence in galaxies, and verify that SNe can indeed be the sole driving mechanism. The key novelty of our approach is that we take into account the increased turbulence dissipation timescale associated with the flaring in outer regions of gaseous discs. We analyse the distribution and kinematics of HI and CO in ten nearby star-forming galaxies to obtain the radial profiles of the kinetic energy per unit area for both the atomic gas and the molecular gas. We use a theoretical model to reproduce the observed energy with the sum of turbulent energy from SNe, as inferred from the observed star formation rate (SFR) surface density, and the gas thermal energy. For the atomic gas, we explore two extreme cases in which the atomic gas is made either of cold neutral medium or warm neutral medium, and the more realistic scenario with a mixture of the two phases. We find that the observed kinetic energy is remarkably well reproduced by our model across the whole extent of the galactic discs, assumingηconstant with the galactocentric radius. Taking into account the uncertainties on the SFR surface density and on the atomic gas phase, we obtain that the median SN efficiencies for our sample of galaxies are ⟨ηatom⟩ = 0.015−0.008+0.018for the atomic gas and ⟨ηmol⟩ = 0.003−0.002+0.006for the molecular gas. We conclude that SNe alone can sustain gas turbulence in nearby galaxies with only few percent of their energy and that there is essentially no need for any further source of energy.

In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, ultra-diffuse galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 H i-selected UDGs and 21 “normal” dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint spectral energy distribution fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime—although the UDGs in our sample show unusually low star formation rate surface densities given their H i content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the “normal” dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the H i-rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H2 poor compared even to the overall dwarf population.

In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, Ultra-Diffuse Galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 HI-selected UDGs and 21 "normal" dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint SED fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime - although the UDGs in our sample show unusually low star formation rate surface densities given their HI content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the "normal" dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the HI-rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H$_2$-poor compared even to the overall dwarf population.

We present a sub-grid model for star formation in galaxy simulations, incorporating molecular hydrogen (H 2 ) production via dust grain condensation and its destruction through star formation and photodissociation. Implemented within the magnetohydrodynamical code AREPO , our model tracks the non-equilibrium mass fractions of molecular, ionised, and atomic hydrogen, as well as a stellar component, by solving a system of differential equations governing mass exchange between these phases. Star formation is treated with a variable rate dependent on the local H 2 abundance, which itself varies in a complex way with key quantities such as gas density and metallicity. Testing the model in a cosmological simulation of a Milky Way-mass galaxy, we obtain a well-defined spiral structure at z =0, including a gas disc twice the size of the stellar one, alongside a realistic star formation history. Our results show a broad range of star formation efficiencies per free-fall time, from as low as 0.001% at high redshift to values between 0.1% and 10% for ages ≳ 3−4 Gyr. These findings align well with observational estimates and simulations of a turbulent interstellar medium. Notably, our model reproduces a star formation rate versus molecular hydrogen surface densities relation akin to the molecular Kennicutt-Schmidt law. Furthermore, we find that the star formation efficiency varies with density and metallicity, providing an alternative to fixed-efficiency assumptions and enabling comparisons with more detailed star formation models. Comparing different star formation prescriptions, we find that in models that link star formation to H 2 , star formation onset is ∼ 500 Myr later than those relying solely on total or cold gas density.

The cosmic time evolution of the radial structure is one of the key topics in the investigation of disc galaxies. In the build-up of galactic discs, gas infall is an important ingredient and it produces radial gas inflows as a physical consequence of angular momentum conservation since the infalling gas onto the disc at a specific radius has lower angular momentum than the circular motions of the gas at the point of impact. NGC,300 is a well-studied isolated, bulgeless, and low-mass disc galaxy ideally suited for an investigation of galaxy evolution with radial gas inflows. Our aim is to investigate the effects of radial gas inflows on the physical properties of NGC,300, for example the radial profiles of H i gas mass and star formation rate (SFR) surface densities, specific star formation rate (sSFR), and metallicity, and to study how the metallicity gradient evolves with cosmic time. A chemical evolution model for NGC,300 was constructed by assuming its disc builds up progressively by the infalling of metal-poor gas and the outflowing of metal-enriched gas. Radial gas inflows were also considered in the model. We used the model to build a bridge between the available data (e.g. gas content, SFR, and chemical abundances) observed today and the galactic key physical processes. Our model including the radial gas inflows and an inside-out disc formation scenario can simultaneously reproduce the present-day observed radial profiles of H i gas mass surface density, SFR surface density, sSFR, gas-phase, and stellar metallicity. We find that, although the value of radial gas inflow velocity is as low as -0.1 the radial gas inflows steepen the present-day radial profiles of H i gas mass surface density, SFR surface density, and metallicity, but flatten the radial sSFR profile. Incorporating radial gas inflows significantly improves the agreement between our model predicted present-day sSFR profile and the observations of NGC,300. Our model predictions are also in good agreement with the star-forming galaxy main sequence and the mass-metallicity relation of star-forming galaxies. It predicts a significant flattening of the metallicity gradient with cosmic time. We also find that the model predicted star formation has been more active recently, indicating that the radial gas inflows may help to sustain star formation in local spirals, at least in NGC,300.

Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the ISM. We present a multiband RC survey of the Local Group galaxy M 33 down to ≃30 pc linear resolution observed with the Karl G. Jansky Very Large Array (VLA). We calibrate the star formation rate surface density and investigate the impact of diffuse emission on this calibration using a structural decomposition. Separating the thermal and non-thermal emission components, the correlation between different phases of the ISM, and the impact of massive star formation are being investigated. Radio sources with sizes ≲200 pc constitute about 36 per cent (46 per cent) of the total RC emission at 1.5 GHz (6.3 GHz) in the inner 18 × 18 arcmin2 (or 4 kpc × 4 kpc) disc of M 33. The non-thermal spectral index becomes flatter with increasing star formation rate surface density, indicating the escape of cosmic ray electrons from their birth places. The magnetic field strength also increases with star formation rate following a bi-modal relation, indicating that the small-scale turbulent dynamo acts more efficiently at higher luminosities and star formation rates. Although the correlations are tighter in star-forming regions, the non-thermal emission is also correlated with the more quiescent molecular gas in the ISM. An almost linear molecular star formation law exists in M 33 when excluding diffuse structures. Massive star formation amplifies the magnetic field and increases the number of high-energy cosmic ray electrons, which can help the onset of winds and outflows.

The initial mass function (IMF) is a construct that describes the distribution of stellar masses for a newly formed population of stars. It is a fundamental element underlying all of star and galaxy formation and has been the subject of extensive investigation for more than 60 yr. In the past few decades, there has been a growing, and now substantial, body of evidence supporting the need for a variable IMF. In this light, it is crucial to investigate the IMF’s characteristics across different spatial scales and to understand the factors driving its variability. We make use of spatially resolved spectroscopy to examine the high-mass IMF slope of star-forming galaxies within the SAMI survey. By applying the Kennicutt method and stellar population synthesis models, we estimated both the spaxel-resolved ( $\alpha_{res}$ ) and galaxy-integrated ( $\alpha_{int}$ ) high-mass IMF slopes of these galaxies. Our findings indicate that the resolved and integrated IMF slopes exhibit a near 1:1 relationship for $\alpha_{int}\gtrsim -2.7$ . We observe a wide range of $\alpha_{res}$ distributions within galaxies. To explore the sources of this variability, we analyse the relationships between the resolved and integrated IMF slopes and both the star formation rate (SFR) and SFR surface density ( $\Sigma_{\textrm{SFR}}$ ). Our results reveal a strong correlation where flatter/steeper slopes are associated with higher/lower SFR and $\Sigma_{\textrm{SFR}}$ . This trend is qualitatively similar for resolved and global scales. Additionally, we identify a mass dependency in the relationship with SFR, though none was found in the relation between the resolved slope and $\Sigma_{\textrm{SFR}}$ . These findings suggest an scenario where the formation of high-mass stars is favoured in regions with more concentrated star formation. This may be a consequence of the reduced fragmentation of molecular clouds, which nonetheless accrete more material.