Search for: All records

The property of star formation rate (SFR) is tightly connected to the amount of dense gas in molecular clouds. However, it is not fully understood how the relationship between dense molecular gas and star formation varies within galaxies and in different morphological environments. Most previous studies have typically been limited to kiloparsec-scale resolution such that different environments could not be resolved. In this work, we present new ALMA observations of HCN(1−0) at 260 pc scale to test how the amount of dense gas and its ability to form stars varies with environmental properties. Combined with existing CO(2−1) observations from ALMA and Hαfrom MUSE, we measured the HCN/CO line ratio, a proxy for the dense gas fraction, and SFR/HCN, a proxy for the star formation efficiency of the dense gas. We find a systematic > 1 dex increase (decreases) of HCN/CO (SFR/HCN) towards the centre of the galaxy, and roughly flat trends of these ratios (average variations < 0.3 dex) throughout the disc. While spiral arms, interarm regions, and bar ends show similar HCN/CO and SFR/HCN, on the bar, there is a significantly lower SFR/HCN at a similar HCN/CO. The strong environmental influence on dense gas and star formation in the centre of NGC 4321, suggests either that clouds couple strongly to the surrounding pressure or that HCN emission traces more of the bulk molecular gas that is less efficiently converted into stars. Across the disc, where the ISM pressure is typically low, SFR/HCN is more constant, indicating a decoupling of the clouds from their surrounding environment. The low SFR/HCN on the bar suggests that gas dynamics (e.g. shear and streaming motions) can have a large effect on the efficiency with which dense gas is converted into stars. In addition, we show that HCN/CO is a good predictor of the mean molecular gas surface density at 260 pc scales across environments and physical conditions. 

ABSTRACT The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at z ≲ 2.5 leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z > 3) JWST observations have challenged the fundamental (i.e. redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and Evolution and Assembly of GaLaxies and their Environment (EAGLE) all predict MZRs that exhibit scatter with a secondary dependence on SFR up to z = 8. We introduce the concept of a ‘strong’ FMR, where the strength of correlated scatter does not evolve with time, and a ‘weak’ FMR, where there is some time evolution. We find that each simulation analysed has a statistically significant weak FMR – there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ∼10–40 per cent at z ≳ 3 when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR. 

ABSTRACT The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (MZ*R) in the hydrodynamic simulations Illustris, TNG, and EAGLE (Evolution and Assembly of GaLaxies and their Environment) to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulations. Interestingly, we observe significant variations in the overall normalization and redshift evolution of the MZ*R across the three simulations. However, all simulations consistently demonstrate a tertiary dependence on the specific star formation rate (sSFR) of galaxies. This finding parallels the relationship seen in both simulations and observations between stellar mass, gas-phase metallicity, and some proxy of galaxy gas content (e.g. SFR, gas fraction, and atomic gas mass). Since we find this correlation exists in all three simulations, each employing a subgrid treatment of the dense, star-forming interstellar medium (ISM) to simulate smooth stellar feedback, we interpret this result as a fairly general feature of simulations of this kind. Furthermore, with a toy analytic model, we propose that the tertiary correlation in the stellar component is sensitive to the extent of the ‘burstiness’ of feedback within galaxies. 

Abstract We present the largest catalog to date of star clusters and compact associations in nearby galaxies. We have performed aV-band-selected census of clusters across the 38 spiral galaxies of the PHANGS–Hubble Space Telescope (HST) Treasury Survey, and measured integrated, aperture-corrected near-ultraviolet-U-B-V-Iphotometry. This work has resulted in uniform catalogs that contain ∼20,000 clusters and compact associations, which have passed human inspection and morphological classification, and a larger sample of ∼100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color–color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS–Atacama Large Millimeter/submillimeter Array. We explore new ways of understanding this large data set in a multiscale context by bringing together once-separate techniques for the characterization of clusters (color–color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution fitting techniques. 

Abstract We use 0.1″ observations from the Atacama Large Millimeter Array (ALMA), Hubble Space Telescope (HST), and JWST to study young massive clusters (YMCs) in their embedded “infant” phase across the central starburst ring in NGC 3351. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (∼10⁵M_⊙), small radii (≲ 5 pc), large escape velocities (6–10 km s⁻¹), and short freefall times (0.5–1 Myr). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed Hiiregion–cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline, ranging from ∼1–2 Myr before cluster formation as starless clumps, to ∼4–6 Myr after as exposed Hiiregion–cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit a nonuniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows. 

ABSTRACT We present radial gas-phase metallicity profiles, gradients, and break radii at redshift z = 0–3 from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, Rbreak, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of ∼2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that Rbreak occurs where the gas mixing time-scale is ∼10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic discs can be thought of in two-parts: a star-forming inner disc with a steep gradient and a mixing-dominated outer disc with a flat gradient, with the break radius marking the region of transition between them. 

Abstract The CO-to-H 2 conversion factor ( α CO ) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (3), (6), and (7) observations toward the inner ∼2 kpc of NGC 3627 and NGC 4321 tracing 12 CO, 13 CO, and C 18 O lines on ∼100 pc scales. Our multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

Abstract Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions. 

ABSTRACT Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar dust lanes, which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate on to the nuclear ring of the barred galaxy NGC 1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position–position–velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $$\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot \, yr^{-1}$$ averaged over a time span of 40 Myr, which varies by a factor of a few over time-scales of ∼10 Myr. Most of the inflow appears to be consumed by star formation in the ring, which is currently occurring at a star formation rate (SFR) of $$\simeq\!1.8\!-\!2 \, \rm M_\odot \, yr^{-1}$$, suggesting that the inflow is causally controlling the SFR in the ring as a function of time. 

Abstract Using recently acquired Hubble Space Telescope NIR observations ( J , Pa β , and H bands) of the nearby galaxy NGC 1313, we investigate the timescales required by a young star cluster to emerge from its natal cloud. We search for extincted star clusters, potentially embedded in their natal cloud as either (1) compact sources in regions with high H α /Pa β extinctions or (2) compact H ii regions that appear as point-like sources in the Pa β emission map. The NUV–optical–NIR photometry of the candidate clusters is used to derive their ages, masses, and extinctions via a least- χ 2 spectral energy distribution broad- and narrowband fitting process. The 100 clusters in the final samples have masses in the range and moderate extinctions, E ( B − V ) ≲ 1.0 mag. Focusing on the young clusters (0–6 Myr), we derive a weak correlation between extinction and age of the clusters. Almost half of the clusters have low extinctions, E ( B − V ) < 0.25 mag, already at very young ages (≤3 Myr), suggesting that dust is quickly removed from clusters. A stronger correlation is found between the morphology of the nebular emission (compact, partial or absent, both in H α and Pa β ) and cluster age. Relative fractions of clusters associated with a specific nebular morphology are used to estimate the typical timescales for clearing the natal gas cloud, resulting in between 3 and 5 Myr, ∼1 Myr older than what was estimated from NUV–optical-based cluster studies. This difference hints at a bias for optical-only-based studies, which James Webb Space Telescope will address in the coming years.