Search for: All records

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 Current methods of identifying the ionizing source of nebular emission in galaxies are well defined for the era of single-fiber spectroscopy, but still struggle to differentiate the complex and overlapping ionization sources in some galaxies. With the advent of integral field spectroscopy, the limits of these previous classification schemes are more apparent. We propose a new method for distinguishing the ionizing source in resolved galaxy spectra by use of a multidimensional diagnostic diagram that compares emission-line ratios with velocity dispersion on a spaxel-by-spaxel basis within a galaxy. This new method is tested using the Sydney-Australian-Astronomical-Observatory Multi-object Integral-Field Spectrograph Galaxy Survey (SAMI) Data Release 3 (DR3), which contains 3068 galaxies atz< 0.12. Our results are released as ionization maps available alongside the SAMI DR3 public data. Our method accounts for a more diverse range of ionization sources than the standard suite of emission-line diagnostics; we find 1433 galaxies with a significant contribution from non-star-forming ionization using our improved method as compared to 316 galaxies identified using only emission-line ratio diagnostics. Within these galaxies, we further identify 886 galaxies hosting unique signatures inconsistent with standard ionization by Hiiregions, active galactic nuclei, or shocks. These galaxies span a wide range of masses and morphological types and comprise a sizable portion of the galaxies used in our sample. With our revised method, we show that emission-line diagnostics alone do not adequately differentiate the multiple ways to ionize gas within a galaxy. 

We investigate the molecular gas content of z  ∼ 6 quasar host galaxies using the Institut de Radioastronomie Millimétrique Northern Extended Millimeter Array. We targeted the 3 mm dust continuum, and the line emission from CO(6–5), CO(7–6), and [C  I ] 2−1 in ten infrared–luminous quasars that have been previously studied in their 1 mm dust continuum and [C  II ] line emission. We detected CO(7–6) at various degrees of significance in all the targeted sources, thus doubling the number of such detections in z  ∼ 6 quasars. The 3 mm to 1 mm flux density ratios are consistent with a modified black body spectrum with a dust temperature T dust  ∼ 47 K and an optical depth τ ν  = 0.2 at the [C  II ] frequency. Our study provides us with four independent ways to estimate the molecular gas mass, M H2 , in the targeted quasars. This allows us to set constraints on various parameters used in the derivation of molecular gas mass estimates, such as the mass per luminosity ratios α CO and α [CII] , the gas-to-dust mass ratio δ g/d , and the carbon abundance [C]/H 2 . Leveraging either on the dust, CO, [C  I ], or [C  II ] emission yields mass estimates of the entire sample in the range M H2  ∼ 10 10 –10 11 M ⊙ . We compared the observed luminosities of dust, [C  II ], [C  I ], and CO(7–6) with predictions from photo-dissociation and X-ray dominated regions. We find that the former provide better model fits to our data, assuming that the bulk of the emission arises from dense ( n H  > 10 4 cm −3 ) clouds with a column density N H  ∼ 10 23 cm −2 , exposed to a radiation field with an intensity of G 0  ∼ 10 3 (in Habing units). Our analysis reiterates the presence of massive reservoirs of molecular gas fueling star formation and nuclear accretion in z  ∼ 6 quasar host galaxies. It also highlights the power of combined 3 mm and 1 mm observations for quantitative studies of the dense gas content in massive galaxies at cosmic dawn. 

ABSTRACT We study environmental quenching using the spatial distribution of current star formation and stellar population ages with the full SAMI Galaxy Survey. By using a star formation concentration index [C-index, defined as log10(r50, H α/r50, cont)], we separate our sample into regular galaxies (C-index ≥−0.2) and galaxies with centrally concentrated star formation (SF-concentrated; C-index <−0.2). Concentrated star formation is a potential indicator of galaxies currently undergoing ‘outside-in’ quenching. Our environments cover ungrouped galaxies, low-mass groups (M200 ≤ 1012.5M⊙), high-mass groups (M200 in the range 1012.5–14 M⊙) and clusters (M200 > 1014M⊙). We find the fraction of SF-concentrated galaxies increases as halo mass increases by 9 ± 2 per cent, 8 ± 3 per cent, 19 ± 4 per cent, and 29 ± 4 per cent for ungrouped galaxies, low-mass groups, high-mass groups, and clusters, respectively. We interpret these results as evidence for ‘outside-in’ quenching in groups and clusters. To investigate the quenching time-scale in SF-concentrated galaxies, we calculate light-weighted age (AgeL) and mass-weighted age (AgeM) using full spectral fitting, as well as the Dn4000 and HδA indices. We assume that the average galaxy age radial profile before entering a group or cluster is similar to ungrouped regular galaxies. At large radius (1–2 Re), SF-concentrated galaxies in high-mass groups have older ages than ungrouped regular galaxies with an age difference of 1.83 ± 0.38 Gyr for AgeL and 1.34 ± 0.56 Gyr for AgeM. This suggests that while ‘outside-in’ quenching can be effective in groups, the process will not quickly quench the entire galaxy. In contrast, the ages at 1–2 Re of cluster SF-concentrated galaxies and ungrouped regular galaxies are consistent (difference of 0.19 ± 0.21 Gyr for AgeL, 0.40 ± 0.61 Gyr for AgeM), suggesting the quenching process must be rapid. 

Abstract The CO-to-H 2 conversion factor ( α CO ) is critical to studying molecular gas and star formation in galaxies. The value of α CO has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on ~kiloparsec scales revealed low values of α CO in the centers of some barred spiral galaxies, including NGC 3351. We present new Atacama Large Millimeter/submillimeter Array Band 3, 6, and 7 observations of 12 CO, 13 CO, and C 18 O lines on 100 pc scales in the inner ∼2 kpc of NGC 3351. Using multiline radiative transfer modeling and a Bayesian likelihood analysis, we infer the H 2 density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of 2–3 × 10 3 cm −3 in the central ∼1 kpc and a high temperature of 30–60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H 2 abundance of 3 × 10 −4 , our analysis yields α CO ∼ 0.5–2.0 M ⊙ (K km s −1 pc 2 ) −1 with a decreasing trend with galactocentric radius in the central ∼1 kpc. The inflows show a substantially lower α CO ≲ 0.1 M ⊙ (K km s −1 pc 2 ) −1 , likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted α CO of ∼1.5 M ⊙ (K km s −1 pc 2 ) −1 , which is similar to previous dust-modeling-based results at kiloparsec scales. This suggests that low α CO on kiloparsec scales in the centers of some barred galaxies may be due to the contribution of low-optical-depth CO emission in bar-driven inflows. 

ABSTRACT The kinematic morphology–density relation of galaxies is normally attributed to a changing distribution of galaxy stellar masses with the local environment. However, earlier studies were largely focused on slow rotators; the dynamical properties of the overall population in relation to environment have received less attention. We use the SAMI Galaxy Survey to investigate the dynamical properties of ∼1800 early and late-type galaxies with log (M⋆/M⊙) > 9.5 as a function of mean environmental overdensity (Σ5) and their rank within a group or cluster. By classifying galaxies into fast and slow rotators, at fixed stellar mass above log (M⋆/M⊙) > 10.5, we detect a higher fraction (∼3.4σ) of slow rotators for group and cluster centrals and satellites as compared to isolated-central galaxies. We find similar results when using Σ5 as a tracer for environment. Focusing on the fast-rotator population, we also detect a significant correlation between galaxy kinematics and their stellar mass as well as the environment they are in. Specifically, by using inclination-corrected or intrinsic $$\lambda _{R_{\rm {e}}}$$ values, we find that, at fixed mass, satellite galaxies on average have the lowest $$\lambda _{\, R_{\rm {e}},\rm {intr}}$$, isolated-central galaxies have the highest $$\lambda _{\, R_{\rm {e}},\rm {intr}}$$, and group and cluster centrals lie in between. Similarly, galaxies in high-density environments have lower mean $$\lambda _{\, R_{\rm {e}},\rm {intr}}$$ values as compared to galaxies at low environmental density. However, at fixed Σ5, the mean $$\lambda _{\, R_{\rm {e}},\rm {intr}}$$ differences for low and high-mass galaxies are of similar magnitude as when varying Σ5 ($$\Delta \lambda _{\, R_{\rm {e}},\rm {intr}} \sim 0.05$$, with σrandom = 0.025, and σsyst < 0.03). Our results demonstrate that after stellar mass, environment plays a significant role in the creation of slow rotators, while for fast rotators we also detect an independent, albeit smaller, impact of mass and environment on their kinematic properties. 

ABSTRACT Understanding the spatial distribution of metals within galaxies allows us to study the processes of chemical enrichment and mixing in the interstellar medium. In this work, we map the 2D distribution of metals using a Gaussian Process Regression (GPR) for 19 star-forming galaxies observed with the Very Large Telescope/Multi Unit Spectroscopic Explorer (VLT–MUSE) as a part of the PHANGS–MUSE survey. We find that 12 of our 19 galaxies show significant 2D metallicity variation. Those without significant variations typically have fewer metallicity measurements, indicating this is due to the dearth of $${\rm H\, {\small II}}$$ regions in these galaxies, rather than a lack of higher-order variation. After subtracting a linear radial gradient, we see no enrichment in the spiral arms versus the disc. We measure the 50 per cent correlation scale from the two-point correlation function of these radially subtracted maps, finding it to typically be an order of magnitude smaller than the fitted GPR kernel scale length. We study the dependence of the two-point correlation scale length with a number of global galaxy properties. We find no relationship between the 50 per cent correlation scale and the overall gas turbulence, in tension with existing theoretical models. We also find more actively star-forming galaxies, and earlier type galaxies have a larger 50 per cent correlation scale. The size and stellar mass surface density do not appear to correlate with the 50 per cent correlation scale, indicating that perhaps the evolutionary state of the galaxy and its current star formation activity is the strongest indicator of the homogeneity of the metal distribution. 

Abstract We compare mid-infrared (mid-IR), extinction-corrected H α , and CO (2–1) emission at 70–160 pc resolution in the first four PHANGS–JWST targets. We report correlation strengths, intensity ratios, and power-law fits relating emission in JWST’s F770W, F1000W, F1130W, and F2100W bands to CO and H α . At these scales, CO and H α each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H α emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to a good match with the molecular-gas-tracing CO, and as a heating tracer, leading to a good match with the H α . By combining mid-IR, CO, and H α at scales where the overall correlation between cold gas and star formation begins to break down, we are able to separate these two effects. We model the mid-IR above I ν = 0.5 MJy sr −1 at F770W, a cut designed to select regions where the molecular gas dominates the interstellar medium (ISM) mass. This bright emission can be described to first order by a model that combines a CO-tracing component and an H α -tracing component. The best-fitting models imply that ∼50% of the mid-IR flux arises from molecular gas heated by the diffuse interstellar radiation field, with the remaining ∼50% associated with bright, dusty star-forming regions. We discuss differences between the F770W, F1000W, and F1130W bands and the continuum-dominated F2100W band and suggest next steps for using the mid-IR as an ISM tracer. 

ABSTRACT We investigate the mean locally measured velocity dispersions of ionized gas (σgas) and stars (σ*) for 1090 galaxies with stellar masses $$\log \, (M_{\!\ast }/M_{\odot }) \ge 9.5$$ from the SAMI Galaxy Survey. For star-forming galaxies, σ* tends to be larger than σgas, suggesting that stars are in general dynamically hotter than the ionized gas (asymmetric drift). The difference between σgas and σ* (Δσ) correlates with various galaxy properties. We establish that the strongest correlation of Δσ is with beam smearing, which inflates σgas more than σ*, introducing a dependence of Δσ on both the effective radius relative to the point spread function and velocity gradients. The second strongest correlation is with the contribution of active galactic nuclei (AGN) (or evolved stars) to the ionized gas emission, implying that the gas velocity dispersion is strongly affected by the power source. In contrast, using the velocity dispersion measured from integrated spectra (σap) results in less correlation between the aperture-based Δσ (Δσap) and the power source. This suggests that the AGN (or old stars) dynamically heat the gas without causing significant deviations from dynamical equilibrium. Although the variation of Δσap is much smaller than that of Δσ, a correlation between Δσap and gas velocity gradient is still detected, implying that there is a small bias in dynamical masses derived from stellar and ionized gas velocity dispersions.