Balmer Decrement and IRX Break in Tracing Dust Attenuation at Scales of Individual Star-forming Regions in NGC 628

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.

We analyze the optical continuum of star-forming galaxies in the Sloan Digital Sky Survey by fitting stacked spectra with stellar population synthesis models to investigate the relation between stellar mass, stellar metallicity, dust attenuation, and star formation rate. We fit models calculated with star formation and chemical evolution histories that are derived empirically from multi-epoch observations of the stellar mass–star formation rate and the stellar mass–gas-phase metallicity relations, respectively. We also fit linear combinations of single-burst models with a range of metallicities and ages. Star formation and chemical evolution histories are unconstrained for these models. The stellar mass–stellar metallicity relations obtained from the two methods agree with the relation measured from individual supergiant stars in nearby galaxies. These relations are also consistent with the relation obtained from emission-line analysis of gas-phase metallicity after accounting for systematic offsets in the gas-phase metallicity. We measure dust attenuation of the stellar continuum and show that its dependence on stellar mass and star formation rate is consistent with previously reported results derived from nebular emission lines. However, stellar continuum attenuation is smaller than nebular emission line attenuation. The continuum-to-nebular attenuation ratio depends on stellar mass and is smaller in more massive galaxies. Our consistent analysis of stellar continuum and nebular emission lines paves the way for a comprehensive investigation of stellar metallicities of star-forming and quiescent galaxies.

Context. Understanding the wavelength dependence of dust attenuation is vital for inferring the properties of galaxies from their spectral energy distribution (SED) fitting. The dust attenuation curves in star-forming galaxies depend on the complex interplay between the intrinsic physical dust properties and dust-to-star geometry. Due to the lack of observational constraints at high redshift, dust attenuation and extinction laws measured in the local Universe (e.g., the Calzetti attenuation law and the Small Magellanic Cloud and Milky Way extinction laws) have been employed to describe the dust attenuation at early epochs. Aims. We exploit the high sensitivity and spectral resolution of the James Webb Space Telescope (JWST) to constrain dust attenuation laws in z ∼ 7–8 galaxies. Our goals are to: i) check whether dust attenuation curves at high-z differ from the ones measured in the local Universe and ii) quantify the dependence of the inferred galaxy properties on the assumed dust attenuation law. Methods. We developed a modified version of the SED fitting code BAGPIPES by including a detailed dust attenuation curve parameterization. We applied our method to the JWST Near Infrared Spectrograph (NIRSpec) spectra in the ∼0.6–5.3 µm range to probe the nebular line (Hα, Hβ, Hγ, [O II] λ3727, [O III] λλ4959, 5007, [Ne III] λ3869) and continuum emissions of three star-forming galaxies at z = 7–8. Dust attenuation parameters and global galaxy properties are derived from the fit to the data. Results. We find that the attenuation curves of the analyzed high-z galaxies differ from local templates. One out of the three galaxies shows a characteristic 2175Å bump, typically associated with the presence of small carbonaceous dust grains such as polycyclic aromatic hydrocarbons (PAHs). This is among the first pieces of evidence suggesting the presence of PAHs in early galaxies. Galaxy properties such as the stellar mass (M*) and star formation rate (SFR) inferred from the SED fitting are affected by the assumed attenuation curve (with deviations of up to ∼0.35 dex), however, the adopted star formation history plays the dominant role (up to ∼0.4 dex for the same galaxy properties). Conclusions. Our results highlight the importance of accounting for the potential diversity among dust attenuation laws when analyzing the spectra of high-z galaxies, whose dust properties and dust-to-star geometry are still poorly understood. The application of our method to a larger sample of galaxies observed with JWST can provide important insights into the properties of dust and galaxies in the early Universe.

We compare the infrared excess (IRX) and Balmer decrement ( ) as dust attenuation indicators in relation to other galaxy parameters using a sample of ∼32,000 local star-forming galaxies (SFGs) carefully selected from the Sloan Digital Sky Survey, the Galaxy Evolution Explorer, and the Wide-field Infrared Survey Explore. While at fixed , IRX turns out to be independent on galaxy stellar mass, the Balmer decrement does show a strong mass dependence at fixed IRX. We find the discrepancy, parameterized by the color excess ratio , is not dependent on the gas-phase metallicity and axial ratio, but on the specific star formation rate (SSFR) and galaxy size (R e) following . This finding reveals that the nebular attenuation as probed by the Balmer decrement becomes increasingly larger than the global (stellar) attenuation of SFGs with decreasing SSFR surface density. This can be understood in the context of an enhanced fraction of intermediate-age stellar populations that are less attenuated by dust than the H ii region-traced young population, in conjunction with a decreasing dust opacity of the diffuse interstellar matter when spreading over a larger spatial extent. Once the SSFR surface density of an SFG is known, the conversion between attenuation of nebular and stellar emission can be well estimated using our scaling relation.

We exploit the unprecedented depth of integral field data from the KMOS Ultra-deep Rotational Velocity Survey (KURVS) to analyse the strong (Hα) and forbidden ([N ii], [S ii]) emission line ratios in 22 main-sequence galaxies at $z\, \approx \, 1.5$. Using the [N ii]/Hα emission-line ratio, we confirm the presence of the stellar mass – gas-phase metallicity relation at this epoch, with galaxies exhibiting on average 0.13 ± 0.04 dex lower gas-phase metallicity (12 + log(O/H)M13 = 8.40 ± 0.03) for a given stellar mass (log10(M*[M⊙] = 10.1 ± 0.1) .than local main-sequence galaxies. We determine the galaxy-integrated [S ii] doublet ratio, with a median value of [S ii]λ6716/λ6731 = 1.26 ± 0.14 equivalent to an electron density of log10(ne[cm−3]) = 1.95 ± 0.12. Utilising CANDELS HST multi-band imaging we define the pixel surface-mass and star-formation rate density in each galaxy and spatially resolve the fundamental metallicity relation at $z\, \approx \, 1.5$, finding an evolution of 0.05 ± 0.01 dex compared to the local relation. We quantify the intrinsic gas-phase metallicity gradient within the galaxies using the [N ii]/Hα calibration, finding a median annuli-based gradient of ΔZ/ΔR = −0.015 ± 0.005 dex kpc−1. Finally, we examine the azimuthal variations in gas-phase metallicity, which show a negative correlation with the galaxy integrated star-formation rate surface density ($r_{\rm s}\,$ = −0.40, ps = 0.07) but no connection to the galaxies kinematic or morphological properties nor radial variations in stellar mass surface density or star formation rate surface density. This suggests both the radial and azimuthal variations in interstellar medium properties are connected to the galaxy integrated density of recent star formation.

Star formation rates (SFRs), gas-phase metallicities, and stellar masses are crucial for studying galaxy evolution. The different relations resulting from these properties give insights into the complex interplay of gas inside galaxies and their evolutionary trajectory and current characteristics. We aim to characterize these relations at $z 0.3$, corresponding to a 3-4 Gyr lookback time, to gather insight into the galaxies' redshift evolution. We utilized optical integral field spectroscopy data from 65 emission-line galaxies from the MUSE large program MAGPI at a redshift of $0.28<z<0.35$ (average redshift of $z 0.3$) and spanning a total stellar mass range of $8.2< /M_ odot ) < 11.4$. We measured emission line fluxes and stellar masses, allowing us to determine spatially resolved SFRs, gas-phase metallicities, and stellar mass surface densities. We derived the resolved star formation main sequence (rSFMS), resolved mass metallicity relation (rMZR), and resolved fundamental metallicity relation (rFMR) at $z 0.3$, and compared them to results for the local Universe. We find a relatively shallow rSFMS slope of $ 0.014$ compared to the expected slope at this redshift for an ordinary least square (OLS) fitting routine. For an orthogonal distance regression (ODR) routine, a much steeper slope of $ 0.022$ is measured. We confirm the existence of an rMZR at $z 0.3$ with an average metallicity located $ 0.03$ dex above the local Universe's metallicity. Via partial correlation coefficients, evidence is found that the local metallicity is predominantly determined by the stellar mass surface density and has a weak secondary (inverse) dependence on the SFR surface density $ SFR $. Additionally, a significant dependence of the local metallicity on the total stellar mass $M_ $ is found. Furthermore, we find that the stellar mass surface density $ $ and $M_ $ have a significant influence in determining the strength with which $ SFR $ correlates with the local metallicity. We observe that at lower stellar masses, there is a tighter correlation between $ SFR $ and the gas-phase metallicity, resulting in a more pronounced rFMR.

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at 1.3 ≤ z ≤ 2.6 from the MOSFIRE Deep Evolution Field survey. We divide galaxies with a detected Hα emission line and coverage of Hβ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), and we then stack their spectra. From each stack, we measure the Balmer decrement and gas-phase metallicity, and then we compute the median A V and UV continuum spectral slope (β). First, we find that none of the dust properties (Balmer decrement, A V, or β) varies with the axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model in which attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases the SFR and lowers the metallicity, but leaves the dust column density mostly unchanged. We quantify this idea using the Kennicutt–Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation.

We investigate the energy sources of random turbulent motions of ionised gas from H$\alpha$ emission in eight local star-forming galaxies from the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. These galaxies satisfy strict pure star-forming selection criteria to avoid contamination from active galactic nuclei (AGN) or strong shocks/outflows. Using the relatively high spatial and spectral resolution of SAMI, we find that -- on sub-kpc scales our galaxies display a flat distribution of ionised gas velocity dispersion as a function of star formation rate (SFR) surface density. A major fraction of our SAMI galaxies shows higher velocity dispersion than predictions by feedback-driven models, especially at the low SFR surface density end. Our results suggest that additional sources beyond star formation feedback contribute to driving random motions of the interstellar medium (ISM) in star-forming galaxies. We speculate that gravity, galactic shear, and/or magnetorotational instability (MRI) may be additional driving sources of turbulence in these galaxies.

Context. Feedback from massive stars plays a crucial role in regulating the growth of young star-forming galaxies (SFGs) and in shaping their interstellar medium (ISM). This feedback contributes to the removal and mixing of metals via galactic outflows and to the clearance of neutral gas, which facilitates the escape of ionizing photons. Aims. Our goal is to study the impact of stellar feedback on the chemical abundances of the ISM in a sample of SFGs with strong emission lines at z ∼ 3. Methods. We selected 35 low-mass SFGs (7.9 < log(M⋆/M⊙) < 10.3) from deep spectroscopic surveys based on their CIII]λ1908 emission. We used new follow-up near-infrared (NIR) observations to examine their rest-optical emission lines and to identify ionized outflow signatures through broad emission line wings detected after Gaussian modeling of [OIII]λλ4959,5007 profiles. We characterized the gas-phase metallicity and carbon-to-oxygen (C/O) abundance of the galaxies using a Te-based method via the OIII]λ1666/[OIII]λ5007 ratio and photoionization models. Results. We find line ratios and rest-frame equivalent widths (EWs) characteristic of high-ionization conditions powered by massive stars. Our sample displays a mean rest-frame EW([OIII]λ5007) of ∼560 Å, while about 15% of the SFGs show EW([OIII]λλ4959,5007) > 1000 Å and EW(CIII]) > 5 Å, closely resembling those now seen in epoch of reionization (EoR) galaxies with the James Webb Space Telescope. We find high Te values, which imply low gas-phase metallicities 12+log(O/H) ∼ 7.5–8.5 (mean of 17% solar) and C/O abundances from 23% to 128% solar, with no apparent increasing trend with metallicity. Our sample follows the mass-metallicity relation at z ∼ 3, with some galaxies showing lower gas-phase metallicities. This results in significant deviations from the fundamental metallicity relation. From our [OIII]λλ4959,5007 line profile modeling, we find that 65% of our sample shows an outflow component, which is found both blue- or redshifted relative to the ionized gas systemic velocity, and the mean maximum velocities are vmax ∼ 280 km s−1. We find a weak correlation between vmax and the star formation rate surface density (ΣSFR) of vmax = (2.41 ± 0.03) × ΣSFR(0.06 ± 0.03). Moreover, we find that the mass-loading factor μ of our galaxy sample is typically lower than in more massive galaxies from the literature, but it is higher than in typical local dwarf galaxies. In the stellar mass range covered by our sample, we find that μ increases with ΣSFR. This suggests that for a given stellar mass, denser starbursts in low-mass galaxies produce stronger outflows. Our results complement the picture drawn by similar studies at lower redshift, suggesting that the removal of ionized gas in low-mass SFGs driven by stellar feedback is regulated by their stellar mass and by the strength and concentration of their star formation, that is, ΣSFR.

We study the scaling relations between gas-phase metallicity, stellar mass surface density (Σ*), star formation rate surface density (ΣSFR), and molecular gas surface density ($\Sigma _{{\rm H}_2}$) in local star-forming galaxies on scales of a kpc. We employ optical integral field spectroscopy from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, and ALMA data for a subset of MaNGA galaxies. We use partial correlation coefficients and Random Forest regression to determine the relative importance of local and global galactic properties in setting the gas-phase metallicity. We find that the local metallicity depends primarily on Σ* (the resolved mass–metallicity relation, rMZR), and has a secondary anticorrelation with ΣSFR (i.e. a spatially resolved version of the ‘Fundamental Metallicity Relation’, rFMR). We find that $\Sigma _{{\rm H}_2}$ is less important than ΣSFR in determining the local metallicity. This result indicates that gas accretion, resulting in local metallicity dilution and local boosting of star formation, is unlikely to be the primary origin of the rFMR. The local metallicity depends also on the global properties of galaxies. We find a strong dependence on the total stellar mass (M*) and a weaker (inverse) dependence on the total SFR. The global metallicity scaling relations, therefore, do not simply stem out of their resolved counterparts; global properties and processes, such as the global gravitational potential well, galaxy-scale winds and global redistribution/mixing of metals, likely contribute to the local metallicity, in addition to local production and retention.

We present an analysis of the chemical abundance properties of ≈650 star-forming galaxies at z ≈ 0.6 – 1.8. Using integral-field observations from the K - band Multi-Object Spectrograph (KMOS), we quantify the [N ii]/Hα emission-line ratio, a proxy for the gas-phase Oxygen abundance within the interstellar medium. We define the stellar mass – metallicity relation at z ≈ 0.6 – 1.0 and z ≈ 1.2 – 1.8 and analyse the correlation between the scatter in the relation and fundamental galaxy properties (e.g. Hα star-formation rate, Hα specific star-formation rate, rotation dominance, stellar continuum half-light radius and Hubble-type morphology). We find that for a given stellar mass, more highly star-forming, larger and irregular galaxies have lower gas-phase metallicities, which may be attributable to their lower surface mass densities and the higher gas fractions of irregular systems. We measure the radial dependence of gas-phase metallicity in the galaxies, establishing a median, beam smearing-corrected, metallicity gradient of ΔZ/ΔR= 0.002 ± 0.004 dex kpc−1, indicating on average there is no significant dependence on radius. The metallicity gradient of a galaxy is independent of its rest-frame optical morphology, whilst correlating with its stellar mass and specific star-formation rate, in agreement with an inside-out model of galaxy evolution, as well as its rotation dominance. We quantify the evolution of metallicity gradients, comparing the distribution of ΔZ/ΔR in our sample with numerical simulations and observations at z ≈ 0 – 3. Galaxies in our sample exhibit flatter metallicity gradients than local star-forming galaxies, in agreement with numerical models in which stellar feedback plays a crucial role redistributing metals.

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.

Cosmic rays (CRs) are a plausible mechanism for launching winds of cool material from the discs of star-forming galaxies. However, there is no consensus on what types of galaxies likely host CR-driven winds, or what role these winds might play in regulating galaxies’ star formation rates. Using a detailed treatment of the transport and losses of hadronic CRs developed in the previous paper in this series, here we develop a semi-analytical model that allows us to assess the viability of using CRs to launch cool winds from galactic discs. In particular, we determine the critical CR fluxes – and corresponding star formation rate surface densities – above which hydrostatic equilibrium within a given galaxy is precluded because CRs drive the gas off in a wind or otherwise render it unstable. Our model demonstrates that catastrophic, CR-driven wind loss is a possibility at galactic mean surface densities below ${\lesssim}10^2 \ \mathrm{ M}_{\odot }$ pc−2. In this regime – encompassing the Galaxy and local dwarfs – the locus of the CR-stability curve patrols the high side of the observed distribution of galaxies in the Kennicutt–Schmidt parameter space of star formation rate versus gas surface density. However, hadronic losses render CRs unable to drive global winds in galaxies with surface densities above the ∼102−103 M⊙ pc−2 transition region. Our results show that quiescent, low surface density galaxies like the Milky Way are poised on the cusp of instability, such that small changes to interstellar mass (ISM) parameters can lead to the launching of CR-driven outflows, and we suggest that, as a result, CR feedback sets an ultimate limit to the star formation efficiency of most modern galaxies.

The dust attenuation for a sample of ∼10,000 local (z ≲ 0.1) star-forming galaxies is constrained as a function of their physical properties. We utilize aperture-matched multiwavelength data available from the Galaxy Evolution Explorer and the Sloan Digital Sky Survey to ensure that regions of comparable size in each galaxy are being analyzed. We follow the method of Calzetti et al. and characterize the dust attenuation through the UV power-law index, β, and the dust optical depth, which is quantified using the difference in Balmer emission line optical depth, . The observed linear relationship between β and is similar to the local starburst relation, but the large scatter (σ int = 0.44) suggests that there is significant variation in the local universe. We derive a selective attenuation curve over the range 1250 Å < λ < 8320 Å and find that a single attenuation curve is effective for characterizing the majority of galaxies in our sample. This curve has a slightly lower selective attenuation in the UV compared to previously determined curves. We do not see evidence to suggest that a 2175 Å feature is significant in the average attenuation curve. Significant positive correlations are seen between the amount of UV and optical reddening and galaxy metallicity, mass, star formation rate (SFR), and SFR surface density. This provides a potential tool for gauging attenuation where the stellar population is unresolved, such as at high z.

The carbon monoxide (CO) spectral line energy distributions (SLEDs) of galaxies contain a wealth of information about conditions in their cold interstellar gas. Here, we use galaxy-scale observations of the three lowest energy CO lines to determine SLEDs and line ratios in a sample of 47 nearby, predominantly star-forming galaxies. We find a systematic trend of higher gas excitation with increasing star formation rate (SFR) and SFR surface density (Σ SFR ), with the range of variations being even larger than predicted by simulations. Power-law fits of the CO line ratios as a function of SFR and Σ SFR provide a good description of the trends seen in our sample and also accurately predict values for a wide range of galaxy types compiled from the literature. Based on these fits, we provide prescriptions for estimating CO(1–0) luminosities and molecular gas masses using CO(3–2) or CO(2–1) in cases where CO(1–0) is not observed directly. We compare our observed SLEDs with molecular cloud models in order to examine how the physical properties of cold gas vary across the galaxy population. We find that gas conditions in star-forming and starburst galaxies lie on a continuum with increasing gas density in more actively star-forming systems.