Search

Over the past decade, several millimeter interferometer programs have mapped the nearby star-forming galaxy M51 at a spatial resolution of ${\le}170$ pc. This study combines observations from three major programs: the PdBI Arcsecond Whirlpool Survey (PAWS), the SMA M51 large program (SMA-PAWS), and the Surveying the Whirlpool at Arcseconds with NOEMA (SWAN). The dataset includes the (1-0) and (2-1) rotational transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O isotopologues. The observations cover the $r{<}\rm 3\,kpc$ region including center and part of the disk, thereby ensuring strong detections of the weaker $^{13}$CO and C$^{18}$O lines. All observations are convolved in this analysis to an angular resolution of 4$''$, corresponding to a physical scale of ${\sim}$170 pc. We investigate empirical line ratio relations and quantitatively evaluate molecular gas conditions such as temperature, density, and the CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$). We employ two approaches to study the molecular gas conditions: (i) assuming local thermal equilibrium (LTE) to analytically determine the CO column density and $\alpha_{\rm CO}$, and (ii) using non-LTE modeling with RADEX to fit physical conditions to observed CO isotopologue intensities. We find that the $\alpha_{\rm CO}$ values {in the center and along the inner spiral arm} are $\sim$0.5 dex (LTE) and ${\sim}$0.1 dex (non-LTE) below the Milky Way inner disk value. The average non-LTE $\alpha_{\rm CO}$ is $2.4{\pm}0.5$ M$_\odot$ pc$^{-2}$ (K km s$^{-1}$)$^{-1}$. While both methods show dispersion due to underlying assumptions, the scatter is larger for LTE-derived values. This study underscores the necessity for robust CO line modeling to accurately constrain the molecular ISM's physical and chemical conditions in nearby galaxies., Comment: accepted for publication in AJ; 31 pages, 16 figures, 7 tables

Classifying ionised nebulae in nearby galaxies is crucial to studying stellar feedback mechanisms and understanding the physical conditions of the interstellar medium. This classification task is generally performed by comparing observed line ratios with photoionisation simulations of different types of nebulae (HII regions, planetary nebulae, and supernova remnants). However, due to simplifying assumptions, such simulations are generally unable to fully reproduce the line ratios in observed nebulae. This discrepancy limits the performance of the classical machine-learning approach, where a model is trained on the simulated data and then used to classify real nebulae. In this study, we use a domain-adversarial neural network (DANN) to bridge the gap between photoionisation models (source domain) and observed ionised nebulae from the PHANGS-MUSE survey (target domain). The DANN is an example of a domain-adaptation algorithm, whose goal is to maximise the performance of a model trained on labelled data in the source domain on an unlabelled target domain by extracting domain-invariant features. Our results indicate a significant improvement in classification performance in the target domain when employing the DANN framework compared to a classical neural network (NN) classifier. Additionally, we investigate the impact of adding noise to the source dataset, finding that noise injection acts as a form of regularisation, further enhancing the performances of both the NN and DANN models on the observational data. The combined use of domain adaptation and noise injection improves the classification accuracy in the target domain by 24%. This study highlights the potential of domain adaptation methods in tackling the domain-shift challenge when using theoretical models to train machine-learning pipelines in astronomy., Comment: submitted for publication, comments welcome

Bars are remarkable stellar structures that can transport gas toward centers and drive the secular evolution of galaxies. In this context, it is important to locate dynamical resonances associated with bars. For this study, we used ${Spitzer}$ near-infrared images as a proxy for the stellar gravitational potential and the ALMA CO(J=2-1) gas distribution from the PHANGS survey to determine the position of the main dynamical resonances associated with the bars in the PHANGS sample of 74 nearby star-forming galaxies. We used the gravitational torque method to estimate the location of the bar corotation radius ($R_{\rm CR}$), where stars and gas rotate at the same angular velocity as the bar. Of the 46 barred galaxies in PHANGS, we have successfully determined the corotation (CR) for 38 of them. The mean ratio of the $R_{\rm CR}$ to the bar radius ($R_{\rm bar}$) is $\mathcal{R} = R_{\rm CR}/R_{\rm bar} = 1.12$, with a standard deviation of $0.39$. This is consistent with the average value expected from theory and suggests that bars are predominantly fast. We also compared our results with other bar CR measurements from the literature, which employ different methods, and find good agreement ($\rho = 0.64$). Finally, using rotation curves, we have estimated other relevant resonances such as the inner Lindblad resonance (ILR) and the outer Lindblad resonance (OLR), which are often associated with rings. This work provides a useful catalog of resonances for a large sample of nearby galaxies and emphasizes the clear connection between bar dynamics and morphology.

We present the first spatially-resolved infrared images of supernova remnants (SNRs) in M33 with the unprecedented sensitivity and resolution of JWST. We analyze 43 SNRs in four JWST fields: two covering central and southern M33 with separate NIRCam (F335M, F444W) and MIRI (F560W, F2100W) observations, one $\sim$5 kpc-long radial strip observed with MIRI F770W, and one covering the giant HII region NGC 604 with multiple NIRCam and MIRI broad/narrowband filters. Of the 21 SNRs in the MIRI field, we found three clear detections (i.e., identical infrared and \ha morphologies), and six partial-detections, implying a detection fraction of 43\% in these bands. In contrast, only one SNR (out of 16) is detectable in the NIRCam field. One of the SNRs, L10-080, is a potential candidate for having freshly-formed ejecta dust, based on its size and centrally-concentrated 21 \mum emission. Two SNRs near NGC 604 have strong evidence of molecular (H$_2$) emission at 4.7 \mum, making them the farthest known SNRs with visible molecular shocks. Five SNRs have F770W observations, with the smaller younger objects showing tentative signs of emission, while the older, larger ones have voids. Multi-wavelength data indicate that the clearly-detected SNRs are also among the smallest, brightest at other wavelengths (\ha, radio and X-ray), have the broadest line widths (H$\alpha$ FWHM$\sim$250-350 km/s), and the densest environments. No strong correlation with star-formation histories are seen, with the clearly-detected SNRs having both high-mass ($\sim$35 \Msun) and low-mass ($\lesssim$10 \Msun) progenitors., Comment: 33 pages, 17 figures, 2 tables, submitted to ApJ

Combining Atacama Large Millimeter/sub-millimeter Array CO(2-1) mapping and JWST near- and mid-infrared imaging, we characterize the relationship between CO(2-1) and polycyclic aromatic hydrocarbon (PAH) emission at ~100 pc resolution in 66 nearby star-forming galaxies, expanding the sample size from previous ~100 pc resolution studies by more than an order of magnitude. Focusing on regions of galaxies where most of the gas is likely to be molecular, we find strong correlations between CO(2-1) and 3.3 micron, 7.7 micron, and 11.3 micron PAH emission, estimated from JWST's F335M, F770W, and F1130W filters. We derive power law relations between CO(2-1) and PAH emission, which have indices in the range 0.8-1.2, implying relatively weak variations in the observed CO-to-PAH ratios across the regions that we study. We find that CO-to-PAH ratios and scaling relationships near HII regions are similar to those in diffuse sight lines. The main difference between the two types of regions is that sight lines near HII regions show higher intensities in all tracers. Galaxy centers, on the other hand, show higher overall intensities and enhanced CO-to-PAH ratios compared to galaxy disks. Individual galaxies show 0.19 dex scatter in the normalization of CO at fixed I_PAH and this normalization anti-correlates with specific star formation rate (SFR/M*) and correlates with stellar mass. We provide a prescription that accounts for these galaxy-to-galaxy variations and represents our best current empirical predictor to estimate CO(2-1) intensity from PAH emission, which allows one to take advantage of JWST's excellent sensitivity and resolution to trace cold gas., Comment: 21 pages, 4 figures, 3 tables. Submitted to ApJ

We present resolved $3.6-250~\mu$m dust spectral energy distribution (SED) fitting for $\sim 800$ nearby galaxies. We measure the distribution of radiation field intensities heating the dust, the dust mass surface density ($\Sigma_{\rm d}$), and the fraction of dust in the form of polycyclic aromatic hydrocarbons (PAHs; $q_{\rm PAH}$). We find that the average interstellar radiation field ($\overline{U}$) is correlated both with stellar mass surface density ($\Sigma_{\star}$) and star formation rate surface density ($\Sigma_{\rm SFR}$), while more intense radiation fields are only correlated with $\Sigma_{\rm SFR}$. We show that $q_{\rm PAH}$ is a steeply decreasing function of $\Sigma_{\rm SFR}$, likely reflecting PAH destruction in H II regions. Galaxy integrated $q_{\rm PAH}$ is strongly, negatively correlated with specific star formation rate (sSFR) and offset from the star-forming ``main sequence'' ($\Delta$MS), suggesting that both metallicity and star formation intensity play a role in setting the global $q_{\rm PAH}$. We also find a nearly constant M$_{\rm d}$/M$_\star$ ratio for galaxies on the main sequence, with a lower ratio for more quiescent galaxies, likely due to their lower gas fractions. From these results, we construct prescriptions to estimate the radiation field distribution in both integrated and resolved galaxies. We test these prescriptions by comparing our predicted $\overline{U}$ to results of SED fitting for stacked "main sequence" galaxies at $0

We present a chemical analysis of selected HII regions from the PHANGS-MUSE nebular catalogue. Our intent is to empirically re-calibrate strong-line diagnostics of gas-phase metallicity, applicable across a wide range of metallicities within nearby star-forming galaxies. To ensure reliable measurements of auroral line fluxes, we carried out a new spectral fitting procedure whereby only restricted wavelength regions around the emission lines of interest are taken into account: this assures a better fit for the stellar continuum. No prior cuts to nebulae luminosity were applied to limit biases in auroral line detections. Ionic abundances of O+, O++, N+, S+, and S++ were estimated by applying the direct method. We integrated the selected PHANGS-MUSE sample with other existing auroral line catalogues, appropriately re-analysed to obtain a homogeneous dataset. This was used to derive strong-line diagnostic calibrations that span from 12+log(O/H) = 7.5 to 8.8. We investigate their dependence on the ionisation parameter and conclude that it is likely the primary cause of the significant scatter observed in these diagnostics. We apply our newly calibrated strong-line diagnostics to the total sample of HII regions from the PHANGS-MUSE nebular catalogue, and we exploit these indirect metallicity estimates to study the radial metallicity gradient within each of the 19 galaxies of the sample. We compare our results with the literature and find good agreement, validating our procedure and findings. With this paper, we release the full catalogue of auroral and nebular line fluxes for the selected HII regions from the PHANGS-MUSE nebular catalogue. This is the first catalogue of direct chemical abundance measurements carried out with PHANGS-MUSE data., Comment: 20 pages, 14 figures, 6 tables; accepted for publication in Astronomy & Astrophysics

Modelling complex line emission in the interstellar medium (ISM) is a degenerate, high-dimensional problem. Here, we present McFine, a tool for automated multi-component fitting of emission lines with complex hyperfine structure, in a fully automated way. We use Markov chain Monte Carlo (MCMC) to efficiently explore the complex parameter space, allowing for characterising model denegeracies. This tool allows for both local thermodynamic equilibrium (LTE) and radiative-transfer (RT) models. McFine can fit individual spectra and data cubes, and for cubes encourage spatial coherence between neighbouring pixels. It is also built to fit the minimum number of distinct components, to avoid overfitting. We have carried out tests on synthetic spectra, where in around 90~per~cent of cases it fits the correct number of components, otherwise slightly fewer components. Typically, $T_{\rm ex}$ is overestimated and $\tau$ underestimated, but accurate within the estimated uncertainties. The velocity and line widths are recovered with extremely high accuracy, however. We verify McFine by applying to a large Atacama Large Millimeter/submillimeter Array (ALMA) N$_2$H$^+$ mosaic of an high-mass star forming region, G316.75-00.00. We find a similar quality of fit to our synthetic tests, aside from in the active regions forming O-stars, where the assumptions of Gaussian line profiles or LTE may break down. To show the general applicability of this code, we fit CO(J = 2-1) observations of NGC 3627, a nearby star-forming galaxy, again obtaining excellent fit quality. McFine provides a fully automated way to analyse rich datasets from interferometric observations, is open source, and pip-installable., Comment: 15 pages, 2 Appendices, 22 Figures (7 in Appendices), 1 Table. Accepted to MNRAS

We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $\sim 10 \times 3~{\rm kpc}^2$ area of the lower circumgalactic medium (CGM) was mapped with MIRI F770W at 12 pc scales. Thanks to the sensitivity and resolution of JWST, we detect dust emission out to $\sim 4$ kpc from the disk, in the form of filaments, arcs, and super-bubbles. Some of these filaments can be traced back to regions with recent star formation activity, suggesting that feedback-driven galactic winds play an important role in regulating baryonic cycling. The presence of dust at these altitudes raises questions about the transport mechanisms at play and suggests that small dust grains are able to survive for several tens of million years after having been ejected by galactic winds in the disk-halo interface. We lay out several scenarios that could explain this emission: dust grains may be shielded in the outer layers of cool dense clouds expelled from the galaxy disk, and/or the emission comes from the mixing layers around these cool clumps where material from the hot gas is able to cool down and mix with these cool cloudlets. This first set of data and upcoming spectroscopy will be very helpful to understand the survival of dust grains in energetic environments, and their contribution to recycling baryonic material in the mid-plane of galaxies., Comment: Accepted for publication in Astronomy & Astrophysics; 16 pages, 8 figures

Measuring the properties of the cold neutral medium (CNM) in low-metallicity galaxies provides insight into heating and cooling mechanisms in early Universe-like environments. We report detections of two localized atomic neutral hydrogen (HI) absorption features in NGC 6822, a low-metallicity (0.2 Z$_{\odot}$) dwarf galaxy in the Local Group. These are the first unambiguous CNM detections in a low-metallicity dwarf galaxy outside the Magellanic Clouds. The Local Group L-Band Survey (LGLBS) enabled these detections due to its high spatial (15 pc for HI emission) and spectral (0.4 \kms) resolution. We introduce LGLBS and describe a custom pipeline to search for HI absorption at high angular resolution and extract associated HI emission. A detailed Gaussian decomposition and radiative transfer analysis of the NGC 6822 detections reveals five CNM components, with key properties: a mean spin temperature of 32$\pm$6 K, a mean CNM column density of 3.1$\times$10$^{20}$ cm$^{-2}$, and CNM mass fractions of 0.33 and 0.12 for the two sightlines. Stacking non-detections does not reveal low-level signals below our median optical depth sensitivity of 0.05. One detection intercepts a star-forming region, with the HI absorption profile encompassing the CO (2$-$1) emission, indicating coincident molecular gas and a depression in high-resolution HI emission. We also analyze a nearby sightline with deep, narrow HI self-absorption dips, where the background warm neutral medium is attenuated by intervening CNM. The association of CNM, CO, and H$\alpha$ emissions suggests a close link between the colder, denser HI phase and star formation in NGC 6822., Comment: 20 pages, 11 figures, accepted for publication in ApJ

Neutral atomic gas (HI) effectively traces galactic dynamics across mid to large galactocentric radii. However, its limitations in observing small-scale changes within the central few kiloparsecs, coupled with the often observed HI deficit in galactic centers, necessitates using molecular gas emission as a preferred tracer in these regions. Understanding the dynamics of both neutral atomic and molecular gas is crucial for a more complete understanding of how galaxies evolve, funnel gas from the outer disk into their central parts, and eventually form stars. In this work we aim to quantify the dynamics of both, the neutral atomic and molecular gas, in the nearby spiral galaxies NGC 1512, NGC 4535, and NGC 7496 using new MeerKAT-HI observations together with ALMA CO (2-1) observations from the PHANGS collaboration. We use the analysis tool 3D-Barolo to fit tilted ring models to the HI and CO observations. A combined approach of using the HI to constrain the true disk orientation parameters before applying these to the CO datasets is tested. This paper sets expectations for the results of the upcoming high-resolution HI coverage of many galaxies in the PHANGS-ALMA sample using MeerKAT or VLA, to establish a robust methodology for characterizing galaxy orientations and deriving dynamics from combining new HI with existing CO data., Comment: accepted for publication in A&A; 13 pages, 9 Figures (+2 appendix pages)

The molecular-to-atomic gas ratio is crucial to the evolution of the interstellar medium in galaxies. We investigate the balance between the atomic ($\Sigma_{\rm HI}$) and molecular gas ($\Sigma_{\rm H2}$) surface densities in eight nearby star-forming galaxies using new high-quality observations from MeerKAT and ALMA (for HI and CO, respectively). We define the molecular gas ratio as $R_{\rm mol} = \Sigma_{\rm H2} / \Sigma_{\rm HI}$ and measure how it depends on local conditions in the galaxy disks using multi-wavelength observations. We find that, depending on the galaxy, HI is detected at $>3\sigma$ out to 20-120 kpc in galactocentric radius ($r_{\rm gal}$). The typical radius at which $\Sigma_{\rm HI}$ reaches 1~$\rm M_\odot~pc^{-2}$ is $r_{\rm HI}\approx22$~kpc, which corresponds to 1-3 times the optical radius ($r_{25}$). $R_{\rm mol}$ correlates best with the dynamical equilibrium pressure, P$_{\rm DE}$, among potential drivers studied, with a median correlation coefficient of $<\rho>=0.89$. Correlations between $R_{\rm mol}$ and star formation rate, total gas and stellar surface density, metallicity, and $\Sigma_{\rm SFR}$/P$_{\rm DE}$ are present but somewhat weaker. Our results also show a direct correlation between P$_{\rm DE}$ and $\Sigma_{\rm SFR}$, supporting self-regulation models. Quantitatively, we measure similar scalings as previous works and attribute the modest differences that we find to the effect of varying resolution and sensitivity. At $r_{\rm gal} {\gtrsim}0.4~r_{25}$, atomic gas dominates over molecular gas, and at the balance of these two gas phases, we find that the baryon mass is dominated by stars, with $\Sigma_{*} > 5~\Sigma_{\rm gas}$. Our study constitutes an important step in the statistical investigation of how local galaxy properties impact the conversion from atomic to molecular gas in nearby galaxies., Comment: accepted for publication in A&A; 20 pages, 12 Figures (+4 appendix pages)

The 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. In this work, we study dense gas and star formation in the nearby spiral galaxy NGC 4321 to test how the amount of dense gas and its ability to form stars varies with environmental properties at 260 pc scales. We present new ALMA observations of HCN(1-0) line emission. Combined with existing CO(2-1) observations from ALMA, and H-alpha from MUSE, as well as F2100W from JWST to trace the SFR, we measure 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. Towards the centre of the galaxy, HCN/CO systematically increases while SFR/HCN decreases, but these ratios stay roughly constant throughout the disc. 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. We conclude that the centres of galaxies show the strongest environmental influence on dense gas and star formation, suggesting either that clouds couple strongly to the surrounding pressure or that HCN is tracing more of the bulk molecular gas that is less efficiently converted into stars. On the contrary, across the disc of NGC 4321, where the ISM pressure is typically low, SFR/HCN does not show large variations (< 0.3 dex) in agreement with Galactic observations of molecular clouds. Despite the large variations across environments and physical conditions, HCN/CO is a good predictor of the mean molecular gas surface density at 260 pc scales., Comment: 18 pages, 9 figures, accepted for pub in A&A, Jun 13, 2024

Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observations of the \textsuperscript{12}CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionised-gas emission lines from SITELLE at the Canada-France-Hawaii Telescope. Although constrained by the resolution of the ALMA observations ($\approx37$~pc), we identify only $52$ GMCs with radii ranging from $30$ to $140$~pc, a low mean molecular gas mass surface density $\langle\Sigma_{\rm gas}\rangle\approx125$~M$_\odot$~pc$^{-2}$ and a high mean virial parameter $\langle\alpha_{\rm obs,vir}\rangle\approx5.3$. We measure spatially-resolved molecular gas depletion times ($\tau_{\rm dep}\equiv1/{\rm SFE}$) with a spatial resolution of $\approx100$~pc within a galactocentric distance of $1.5$~kpc. The global depletion time is $\approx2.0$~Gyr but $\tau_{\rm dep}$ increases toward the galaxy centre, with a maximum $\tau_{\rm dep,max}\approx5.2$~Gyr. However, no pure \ion{H}{II} region is identified in NGC~524 using ionised-gas emission-line ratio diagnostics, so the $\tau_{\rm dep}$ inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC~524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC~524., Comment: 17 pages, 10 figures. To be published in MNRAS, accepted on May 27

NGC~613 is a nearby barred spiral galaxy with a nuclear ring. Exploiting high spatial resolution ($\approx20$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(1-0) observations, we study the giant molecular clouds (GMCs) in the nuclear ring and its vicinity, identifying $158$ spatially- and spectrally-resolved GMCs. The GMC sizes ($R_{\mathrm{c}}$) are comparable to those of the clouds in the Milky Way (MW) disc, but their gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities are larger. The GMC size -- linewidth relation ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{0.77}$) is steeper than that of the clouds of the MW disc and centre, and the GMCs are on average only marginally gravitationally bound (with a mean virial parameter $\langle\alpha_{\mathrm{obs,vir}}\rangle\approx1.7$). We discuss the possible origins of the steep size -- linewidth relation and enhanced observed linewidths of the clouds and suggest that a combination of mechanisms such as stellar feedback, gas accretion and cloud-cloud collisions, as well as the gas inflows driven by the large-scale bar, may play a role., Comment: 15 pages, 8 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2304.10471

Polycyclic aromatic hydrocarbons (PAHs) are a ubiquitous component of the interstellar medium (ISM) in z~0 massive, star-forming galaxies and play key roles in ISM energy balance, chemistry, and shielding. Wide field of view, high resolution mid-infrared (MIR) images from JWST provides the ability to map the fraction of dust in the form of PAHs and the properties of these key dust grains at 10-50 pc resolution in galaxies outside the Local Group. We use MIR JWST photometric observations of a sample of 19 nearby galaxies from the "Physics at High Angular Resolution in Nearby GalaxieS" (PHANGS) survey to investigate the variations of the PAH fraction. By comparison to lower resolution far-IR mapping, we show that a combination of the MIRI filters (R$_{\rm{PAH}}$ = [F770W+F1130W]/F2100W) traces the fraction of dust by mass in the form of PAHs (i.e., the PAH fraction, or q$_{\rm{PAH}}$). Mapping R$_{\rm{PAH}}$ across the 19 PHANGS galaxies, we find that the PAH fraction steeply decreases in HII regions, revealing the destruction of these small grains in regions of ionized gas. Outside HII regions, we find R$_{\rm{PAH}}$ is constant across the PHANGS sample with an average value of 3.43$\pm$0.98, which, for an illuminating radiation field of intensity 2-5 times that of the radiation field in the solar neighborhood, corresponds to q$_{\rm{PAH}}$ values of 3-6%., Comment: Accepted at ApJ, 39 pages, 25 figures

We present the largest extragalactic survey of supernova remnant (SNR) candidates in nearby star-forming galaxies using exquisite spectroscopic maps from MUSE. Supernova remnants exhibit distinctive emission-line ratios and kinematic signatures, which are apparent in optical spectroscopy. Using optical integral field spectra from the PHANGS-MUSE project, we identify SNRs in 19 nearby galaxies at ~ 100~pc scales. We use five different optical diagnostics: (1) line ratio maps of [SII]/H$\alpha$; (2) line ratio maps of [OI]/H$\alpha$; (3) velocity dispersion map of the gas; (4) and (5) two line ratio diagnostic diagrams from BPT diagrams to identify and distinguish SNRs from other nebulae. Given that our SNRs are seen in projection against HII regions and diffuse ionized gas, in our line ratio maps we use a novel technique to search for objects with [SII]/H$\alpha$ or [OI]/H$\alpha$ in excess of what is expected at fixed H$\alpha$ surface brightness within photoionized gas. In total, we identify 2,233 objects using at least one of our diagnostics, and define a subsample of 1,166 high-confidence SNRs that have been detected with at least two diagnostics. The line ratios of these SNRs agree well with the MAPPINGS shock models, and we validate our technique using the well-studied nearby galaxy M83, where all SNRs we found are also identified in literature catalogs and we recover 51% of the known SNRs. The remaining 1,067 objects in our sample are detected with only one diagnostic and we classify them as SNR candidates. We find that ~ 35% of all our objects overlap with the boundaries of HII regions from literature catalogs, highlighting the importance of using indicators beyond line intensity morphology to select SNRs. [OI]/H$\alpha$ line ratio is responsible for selecting the most objects (1,368; 61%), (abridged)., Comment: 35 pages, 24 figures,6 tables, submitted to A&A

Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies from the PHANGS survey. We find that the arm/interarm contrast in stellar mass surface density (Sigma_*) is very modest, typically a few tens of percent. This is much smaller than the contrasts measured for molecular gas (Sigma_mol) or star formation rate (Sigma_SFR) surface density, which typically reach a factor of ~2-3. Yet, Sigma_mol and Sigma_SFR contrasts show a significant correlation with the enhancement in Sigma_*, suggesting that the small stellar contrast largely dictates the stronger accumulation of gas and star formation. All these contrasts increase for grand-design spirals compared to multi-armed and flocculent systems (and for galaxies with high stellar mass). The median star formation efficiency (SFE) of the molecular gas is 16% higher in spiral arms than in interarm regions, with a large scatter, and the contrast increases significantly (median SFE contrast 2.34) for regions of particularly enhanced stellar contrast (Sigma_* contrast >1.97). The molecular-to-atomic gas ratio (Sigma_mol/Sigma_atom) is higher in spiral arms, pointing to a transformation of atomic to molecular gas. In conclusion, the boost in the star formation efficiency of molecular gas in spiral arms is generally modest or absent, except for locations with exceptionally large stellar contrasts. (abridged), Comment: 26 pages, 16 figures. Accepted for publication in A&A

The mass ($M_\mathrm{BH}$) of a supermassive black hole (SMBH) can be measured using spatially-resolved kinematics of the region where the SMBH dominates gravitationally. The most reliable measurements are those that resolve the smallest physical scales around the SMBHs. We consider here three metrics to compare the physical scales probed by kinematic tracers dominated by rotation: the radius of the innermost detected kinematic tracer $R_\mathrm{min}$ normalised by respectively the SMBH's Schwarzschild radius ($R_\mathrm{Schw}\equiv 2GM_\mathrm{BH}/c^2$, where $G$ is the gravitational constant and $c$ the speed of light), sphere-of-influence (SOI) radius ($R_\mathrm{SOI}\equiv GM_\mathrm{BH}/\sigma_\mathrm{e}^2$, where $\sigma_\mathrm{e}$ is the stellar velocity dispersion within the galaxy's effective radius) and equality radius [the radius $R_\mathrm{eq}$ at which the SMBH mass equals the enclosed stellar mass, $M_\mathrm{BH}=M_*(R_\mathrm{eq})$, where $M_*(R)$ is the stellar mass enclosed within the radius $R$]. All metrics lead to analogous simple relations between $R_\mathrm{min}$ and the highest circular velocity probed $V_\mathrm{c}$. Adopting these metrics to compare the SMBH mass measurements using molecular gas kinematics to those using megamaser kinematics, we demonstrate that the best molecular gas measurements resolve material that is physically closer to the SMBHs in terms of $R_\mathrm{Schw}$ but is slightly farther in terms of $R_\mathrm{SOI}$ and $R_\mathrm{eq}$. However, molecular gas observations of nearby galaxies using the most extended configurations of the Atacama Large Millimeter/sub-millimeter Array can resolve the SOI comparably well and thus enable SMBH mass measurements as precise as the best megamaser measurements., Comment: 13 pages, 9 figures. Accepted by MNRAS

Supermassive black hole (SMBH) masses can be measured by observing the impacts of the SMBHs on dynamical tracers around them. We present high angular resolution ($0.19$ arcsec or $\approx24$ pc) Atacama Large Millimeter/submillimeter Array observations of the $^{12}$CO(3-2) line emission of the early-type galaxy NGC 4751, which reveal a highly-inclined regularly-rotating molecular gas disc with clear central Keplerian motions. Using a Hubble Space Telescope image to constrain the stellar mass distribution, we forward model the molecular gas kinematics and data cube in a Bayesian framework using the Kinematic Molecular Simulation code. Assuming a constant mass-to-light ratio ($M/L$), we infer a SMBH mass $M_\text{BH}=3.43^{+0.45}_{-0.44}\times10^9$ $\text{M}_\odot$ and a F160W filter stellar $M/L$ $M/L_\text{F160W}=(2.68\pm0.11)$ $\text{M}_\odot/\text{L}_{\odot,\text{F160W}}$ (all quoted uncertainties are at $3\sigma$ confidence). Assuming a linearly spatially-varying $M/L$, we infer $M_\text{BH}=2.79_{-0.57}^{+0.75}\times10^9$ $\text{M}_\odot$ and $\left(M/L_\text{F160W}\right)/\left(\text{M}_\odot/\text{L}_{\odot,\text{F160W}}\right)=3.07^{+0.27}_{-0.35}-0.09^{+0.08}_{-0.06}\,\left(R/\text{arcsec}\right)$, where $R$ is the galactocentric radius. We also present alternative SMBH mass estimates using the Jeans Anisotropic Modelling (JAM) method and SINFONI stellar kinematics. Assuming a cylindrically-aligned velocity ellipsoid (JAM$_\text{cyl}$) we infer $M_\text{BH}=(2.52\pm 0.36)\times10^9$ $\text{M}_\odot$, while assuming a spherically-aligned velocity ellipsoid (JAM$_\text{sph}$) we infer $M_\text{BH}=(3.24\pm0.87)\times10^9$ $\text{M}_\odot$. Our derived masses are all consistent with one another, but they are larger than (and inconsistent with) one previous stellar dynamical measurement using Schwarzschil's method and the same SINFONI kinematics., Comment: 15 pages, 13 figures, submitted to MNRAS. arXiv admin note: text overlap with arXiv:2401.16376

We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII regions cover 32% of the area within 1 kpc of any recent SN. Contrasting this local covering fraction with the fraction of SNe coincident with HII regions, we find a statistical excess of 7.6% +/- 8.7% of all SNe to be associated with HII regions. This increases to an excess of 19.2% +/- 10.4% when considering only core-collapse SNe. These estimates appear to be in good agreement with qualitative results from new, higher resolution HST H-alpha imaging, which also suggest many CCSNe detonate near but not in HII regions. Our results appear consistent with the expectation that only a modest fraction of stars explode during the first 5 Myr of the life of a stellar population, when H-alpha emission is expected to be bright. Of the HII region associated SNe, 8% (11/13) also have associated detected CO(2-1) emission, indicating the presence of molecular gas. The HII region associated SNe have typical Av extinctions approximately equal to 1 mag, consistent with a significant amount of pre-clearing of gas from the region before the SNe explode., Comment: Accepted for publication in ApJ; 33 pages, 13 figures, 3 tables in two-column AASTEX63 format

Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the building and consumption of central gas reservoirs are still not well understood. We are building a suite of hydro-dynamical RAMSES simulations of isolated, low-redshift galaxies that mimic the properties of the PHANGS sample. The initial conditions of the models reproduce the observed stellar mass, disc scale length, or gas fraction, and this paper presents a first subset of these models. Most of our simulated galaxies develop a prominent bar structure, which itself triggers central gas fuelling and the building of an over-density with a typical scale of 100-1000 pc. We confirm that if the host galaxy features an ellipsoidal component, the formation of the bar and gas fuelling are delayed. We show that most of our simulations follow a common time evolution, when accounting for mass scaling and the bar formation time. In our simulations, the stellar mass of $10^{10}$~M$_{\odot}$ seems to mark a change in the phases describing the time evolution of the bar and its impact on the interstellar medium. In massive discs (M$_{\star} \geq 10^{10}$~M$_{\odot}$), we observe the formation of a central gas reservoir with star formation mostly occurring within a restricted starburst region, leading to a gas depletion phase. Lower-mass systems (M$_{\star} < 10^{10}$~M$_{\odot}$) do not exhibit such a depletion phase, and show a more homogeneous spread of star-forming regions along the bar structure, and do not appear to host inner bar-driven discs or rings., Comment: 22 pages, 17 figures

We present a quantitative and statistical analysis of the molecular gas morphology in 73 nearby galaxies using high spatial resolution CO (J=2-1) data obtained from the Atacama Large Millimeter/submillimeter Array (ALMA) by the PHANGS large program. We applied three model-independent parameters: Concentration (C), Asymmetry (A), and Clumpiness (S) which are commonly used to parameterize the optical and infrared morphology of galaxies. We find a significant correlation between A and S, with a Spearman's rank correlation coefficient of 0.62 with a p-value of $4 \times 10^{-9}$. This suggests a higher abundance of molecular clumps (i.e. giant molecular cloud associations) in galaxies that display stronger distortion or biased large-scale molecular gas distribution. In addition, the analysis of the C parameter suggests high central molecular concentration in most barred spiral galaxies investigated in this study. Furthermore, we found a positive correlation between the length of the bar structure ($R_\mathrm{bar}/R_{25}$) and the C parameter, with a Spearman's rank correlation coefficient of 0.63 with a p-value of $3.8 \times 10^{-5}$ suggesting that larger bar structure can facilitate overall molecular gas transport and yield higher central concentration than galaxies with shorter bars. Finally, we offer a possible classification scheme of nearby disk galaxies which is based on the CAS parameters of molecular gas., Comment: 23 pages, 18 figures, 4 Tables

We present new catalogs of likely globular clusters (GCs) in 17 nearby spiral galaxies studied as part of the PHANGS-HST Treasury Survey. The galaxies were imaged in five broad-band filters from the near-ultraviolet through the $I$ band. PHANGS-HST has produced catalogs of stellar clusters of all ages by selecting extended sources (from multiple concentration index measurements) followed by morphological classification (centrally concentrated and symmetric or asymmetric, multiple peaks, contaminant) by visually examining the V-band image and separately by a machine-learning algorithm which classified larger samples to reach fainter limits. From both cluster catalogs, we select an initial list of candidate GCs to have $B-V \geq 0.5$ and $V-I \geq 0.73$~mag, then remove likely contaminants (including reddened young clusters, background galaxies misclassified by the neural network, and chance superpositions/blends of stars) after a careful visual inspection. We find that $\approx86$ % of the color-selected candidates classified as spherically symmetric, and $\approx68$ of those classified as centrally concentrated but asymmetric are likely to be GCs. The luminosity functions of the GC candidates in 2 of our 17 galaxies, NGC 628 and NGC 3627, are atypical, and continue to rise at least 1~mag fainter than the expected turnover near $M_V \sim -7.4$. These faint candidate GCs have more extended spatial distributions than their bright counterparts, and may reside in the disk rather than the bulge/halo, similar to faint GCs previously discovered in M101. These faint clusters may be somewhat younger since the age-metallicity degeneracy makes it difficult to determine precise cluster ages from integrated colors once they reach $\approx1$~Gyr., Comment: 21 pages, 12 figures

Studying molecular gas in nearby galaxies using hydrogen cyanide (HCN) as a tracer for higher densities than CO emission still poses a significant challenge. Even though several galaxies have HCN maps on a few kpc scales, higher-resolution maps are still required. Our goal is to examine the contrast in intensity between two tracers that probe different density regimes - HCN(1-0)/CO(2-1) ratio - and their kinematics across NGC 253. By utilizing the advanced capabilities of the Atacama Large Millimeter/submillimeter Array (ALMA), we can map these features at high resolution across a large field of view and uncover the nature of such dense gas in extragalactic systems. We present new ALMA Atacama Compact Array and Total Power (ACA+TP) observations of the HCN emission across NGC 253, covering the inner 8.6' of the galaxy disk at 300 pc scales. We analyze the integrated intensity and mean velocity of HCN and CO along each line of sight and use SCOUSE software to perform spectral decomposition, which considers each velocity component separately. Molecular gas traced by HCN piles up in a ring-like structure at a radius of 2 kpc. The HCN emission is enhanced by 2 orders of magnitude in the central 2 kpc regions, beyond which its intensity decreases with increasing galactocentric distance. The number of components in the HCN spectra shows a robust environmental dependence, with multiple velocity features across the center and bar. We have identified an increase in the HCN/CO ratio in these regions, corresponding to a velocity component likely associated with a molecular outflow. We have also discovered that the ratio between the total infrared luminosity and dense gas mass, which indicates the star formation efficiency of dense gas, is anti-correlated with the molecular gas surface density up to approximately 200 Msul/pc^2. In contrast, beyond this point, the ratio starts to increase., Comment: Accepted for publication to Astronomy and Astrophysics

We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$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-ALMA. We explore new ways of understanding this large dataset in a multi-scale 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 (SED) fitting techniques., Comment: 48 pages, 26 figures, article in review at APJ

The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorithms to mine this information-rich dataset to identify novel patterns. We focus on three of the PHANGS-JWST galaxies, for which we extract properties pertaining to their stellar populations; warm ionized and cold molecular gas; and Polycyclic Aromatic Hydrocarbons (PAHs), as measured over 150 pc-scale regions. We show that we can divide the regions into groups with distinct multiphase gas and PAH properties. In the process, we identify previously-unknown galaxy-wide correlations between PAH band and optical line ratios and use our identified groups to interpret them. The correlations we measure can be naturally explained in a scenario where the PAHs and the ionized gas are exposed to different parts of the same radiation field that varies spatially across the galaxies. This scenario has several implications for nearby galaxies: (i) The uniform PAH ionized fraction on 150 pc scales suggests significant self-regulation in the ISM, (ii) the PAH 11.3/7.7 \mic~ band ratio may be used to constrain the shape of the non-ionizing far-ultraviolet to optical part of the radiation field, and (iii) the varying radiation field affects line ratios that are commonly used as PAH size diagnostics. Neglecting this effect leads to incorrect or biased PAH sizes., Comment: Main results in figures 6 and 12. Submitted to ApJ, and comments are welcome!

The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as part of the PHANGS-JWST Cycle 1 Treasury program at spatial scales of $\sim$5-50pc. Here, we describe pjpipe, an image processing pipeline developed for the PHANGS-JWST program that wraps around and extends the official JWST pipeline. We release this pipeline to the community as it contains a number of tools generally useful for JWST NIRCam and MIRI observations. Particularly for extended sources, pjpipe products provide significant improvements over mosaics from the MAST archive in terms of removing instrumental noise in NIRCam data, background flux matching, and calibration of relative and absolute astrometry. We show that slightly smoothing F2100W MIRI data to 0.9" (degrading the resolution by about 30 percent) reduces the noise by a factor of $\approx$3. We also present the first public release (DR1.1.0) of the pjpipe processed eight-band 2-21 $\mu$m imaging for all nineteen galaxies in the PHANGS-JWST Cycle 1 Treasury program. An additional 55 galaxies will soon follow from a new PHANGS-JWST Cycle 2 Treasury program., Comment: 49 pages (27 in Appendices), 54 Figures (39 in Appendices), 3 Tables. Accepted for publication in ApJS. Updated to match accepted version. Data available at https://archive.stsci.edu/hlsp/phangs/phangs-jwst

We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. 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 (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows., Comment: 27 pages, 12 figures; ApJ accepted

Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the James Webb Space Telescope was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present MIRI imaging mosaics at 5.6 and 21 microns that cover a significant portion of one of M33's spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep ALMA CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new DOLPHOT MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further shows YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass., Comment: 12 pages, 10 figures, 1 tables, accepted for publication at MNRAS

We present the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-AstroSat atlas, which contains ultraviolet imaging of 31 nearby star-forming galaxies captured by the Ultraviolet Imaging Telescope (UVIT) on the AstroSat satellite. The atlas provides a homogeneous data set of far- and near-ultraviolet maps of galaxies within a distance of 22 Mpc and a median angular resolution of 1.4 arcseconds (corresponding to a physical scale between 25 and 160 pc). After subtracting a uniform ultraviolet background and accounting for Milky Way extinction, we compare our estimated flux densities to GALEX observations, finding good agreement. We find candidate extended UV disks around the galaxies NGC 6744 and IC 5332. We present the first statistical measurements of the clumping of the UV emission and compare it to the clumping of molecular gas traced with ALMA. We find that bars and spiral arms exhibit the highest degree of clumping, and the molecular gas is even more clumped than the FUV emission in galaxies. We investigate the variation of the ratio of observed FUV to H$\alpha$ in different galactic environments and kpc-sized apertures. We report that $\sim 65 \%$ varation of the $\log_{10}$(FUV/H$\alpha$) can be described through a combination of dust attenuation with star formation history parameters. The PHANGS-AstroSat atlas enhances the multi-wavelength coverage of our sample, offering a detailed perspective on star formation. When integrated with PHANGS data sets from ALMA, VLT-MUSE, HST and JWST, it develops our comprehensive understanding of attenuation curves and dust attenuation in star-forming galaxies., Comment: 35 pages, 16 figures. The survey webpage is available at https://sites.google.com/view/phangs/home/data/astrosat and the data archive can be accessed at https://www.canfar.net/storage/vault/list/phangs/RELEASES/PHANGS-AstroSat/v1p0

We use high-resolution data from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) project to investigate the connection between circumnuclear gas reservoirs and nuclear activity in a sample of nearby galaxies. Our sample spans a wide range of nuclear activity types including radio galaxies, Seyfert galaxies, low-luminosity active galactic nuclei (AGN) and inactive galaxies. We use measurements of nuclear millimetre continuum emission along with other archival tracers of AGN accretion/activity to investigate previous claims that at, circumnuclear scales (<100 pc), these should correlate with the mass of the cold molecular gas. We find that the molecular gas mass does not correlate with any tracer of nuclear activity. This suggests the level of nuclear activity cannot solely be regulated by the amount of cold gas around the supermassive black hole (SMBH). This indicates that AGN fuelling, that drives gas from the large scale galaxy to the nuclear regions, is not a ubiquitous process and may vary between AGN type, with timescale variations likely to be very important. By studying the structure of the central molecular gas reservoirs, we find our galaxies have a range of nuclear molecular gas concentrations. This could indicate that some of our galaxies may have had their circumnuclear regions impacted by AGN feedback, even though they currently have low nuclear activity. On the other hand, the nuclear molecular gas concentrations in our galaxies could instead be set by secular processes., Comment: 15 pages plus 3 in the appendix, 8 figures plus 1 in the appendix, 3 tables plus 4 in the appendix

The co-evolution of galaxies and supermassive black holes (SMBHs) underpins our understanding of galaxy evolution, but different methods to measure SMBH masses have only infrequently been cross-checked. We attempt to identify targets to cross-check two of the most accurate methods, megamaser and cold molecular gas dynamics. Three promising galaxies are selected from all those with existing megamaser SMBH mass measurements. We present Atacama Large Millimeter/sub-millimeter Array (ALMA) 12CO(2-1) and 230-GHz continuum observations with angular resolutions of about 0.5". Every galaxy has an extended rotating molecular gas disc and 230-GHz continuum source(s), but all also have irregularities and/or non-axisymmetric features: NGC1194 is highly inclined and has disturbed and lopsided central 12CO(2-1) emission; NGC3393 has a nuclear disc with fairly regular but patchy 12CO(2-1) emission with little gas near the kinematic major axis, faint emission in the very centre and two brighter structures reminiscent of a nuclear ring and/or spiral; NGC5765B has a strong bar and very bright 12CO(2-1) emission concentrated along two bisymmetric offset dust lanes and two bisymmetric nuclear spiral arms. 12CO(2-1) and 12CO(3-2) observations with the James Clerk Maxwell Telescope are compared with the ALMA observations. Because of the disturbed gas kinematics and the impractically long integration times required for higher angular resolution observations, none of the three galaxies is suitable for a future SMBH mass measurement. Nonetheless, increasing the number of molecular gas observations of megamaser galaxies is valuable, and the ubiquitous disturbances suggest a link between large-scale gas properties and the existence of megamasers., Comment: 17 pages, 5 figures, accepted by MNRAS

We measure the CO-to-H$_2$ conversion factor ($\alpha_\mathrm{CO}$) in 37 galaxies at 2~kpc resolution, using dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metals ratio. In total, we have $\sim790$ and $\sim610$ independent measurements of $\alpha_\mathrm{CO}$ for CO (2-1) and (1-0), respectively. The mean values for $\alpha_\mathrm{CO~(2-1)}$ and $\alpha_\mathrm{CO~(1-0)}$ are $9.3^{+4.6}_{-5.4}$ and $4.2^{+1.9}_{-2.0}~M_\odot~pc^{-2}~(K~km~s^{-1})^{-1}$, respectively. The CO-intensity-weighted mean for $\alpha_\mathrm{CO~(2-1)}$ is 5.69, and 3.33 for $\alpha_\mathrm{CO~(1-0)}$. We examine how $\alpha_\mathrm{CO}$ scales with several physical quantities, e.g.\ star-formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ($\overline{U}$). Among them, $\overline{U}$, $\Sigma_{\rm SFR}$, and integrated CO intensity ($W_\mathrm{CO}$) have the strongest anti-correlation with spatially resolved $\alpha_\mathrm{CO}$. We provide linear regression results to \aco for all quantities tested. At galaxy integrated scales, we observe significant correlations between $\alpha_\mathrm{CO}$ and $W_\mathrm{CO}$, metallicity, $\overline{U}$, and $\Sigma_{\rm SFR}$. We also find that the normalized $\alpha_\mathrm{CO}$ decreases with stellar mass surface density ($\Sigma_\star$) in the high surface density regions ($\Sigma_\star\geq100~{\rm M_\odot~pc^{-2}}$), following the power-law relations $\alpha_\mathrm{CO~(2-1)}\propto\Sigma_\star^{-0.5}$ and $\alpha_\mathrm{CO~(1-0)}\propto\Sigma_\star^{-0.2}$. The power-law index is insensitive to the assumed dust-to-metals ratio. (abridged), Comment: 27 pages, 13 figures, 4 tables, submitted to ApJ

Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars $>$30 M$_{\odot}$, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction ($\sim$45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies., Comment: 34 pages, 14 figures. Submitted to ApJ. Comments welcome! The density distributions will be made publicly available after journal acceptance of manuscript. Please feel free to contact us in the meantime if you would like to use them

Determining how 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$_2$ conversion factor ($\alpha_\rm{CO}$) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence of $\alpha_\rm{CO}$ on local CO velocity dispersion at 150-pc scales using a new set of dust-based $\alpha_\rm{CO}$ measurements, and propose a new $\alpha_\rm{CO}$ prescription 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-ALMA survey. We find increasing SFE towards high surface density regions like galaxy centers, while using a constant or metallicity-based $\alpha_\rm{CO}$ results 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 3-4 times shorter than in non-barred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2-3 Gyr regardless of the choice of $\alpha_\rm{CO}$ prescription. 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 an enhanced SFE compared to non-barred centers or disk regions., Comment: 12 pages of main text, 4 figures, accepted for publication in ApJ

Observations indicate that the central gas discs are smoother in early-type galaxies than their late-type counterparts, while recent simulations predict that the dynamical suppression of star formation in spheroid-dominated galaxies is preceded by the suppression of fragmentation of their interstellar media. The mass surface density power spectrum is a powerful tool to constrain the degree of structure within a gas reservoir. Specifically here, we focus on the power spectrum slope and aim to constrain whether the shear induced by a dominant spheroidal potential can induce sufficient turbulence to suppress fragmentation, resulting in the smooth central gas discs observed. We compute surface density power spectra for the nuclear gas reservoirs of fourteen simulated isolated galaxies and twelve galaxies observed as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project. Both simulated and observed galaxies range from disc-dominated galaxies to spheroids, with central stellar mass surface densities, a measure of bulge dominance, varying by more than an order of magnitude. For the simulations, the power spectra steepen with increasing central stellar mass surface density, thereby clearly linking the suppression of fragmentation to the shear-driven turbulence induced by the spheroid. The WISDOM observations show a different (but potentially consistent) picture: while there is no correlation between the power spectrum slopes and the central stellar mass surface densities, the slopes scatter around a value of 2.6. This is similar to the behaviour of the slopes of the simulated galaxies with high central stellar mass surface densities, and could indicate that high shear eventually drives incompressible turbulence., Comment: 14 pages, 7 figures + appendices (7 pages, 9 figures); accepted by MNRAS (October 09th, 2023)

HII region electron temperatures are a critical ingredient in metallicity determinations and recent observations reveal systematic variations in the temperatures measured using different ions. We present electron temperatures ($T_e$) measured using the optical auroral lines ([NII]$\lambda5756$, [OII]$\lambda\lambda7320,7330$, [SII]$\lambda\lambda4069,4076$, [OIII]$\lambda4363$, and [SIII]$\lambda6312$) for a sample of HII regions in seven nearby galaxies. We use observations from the Physics at High Angular resolution in Nearby Galaxies survey (PHANGS) obtained with integral field spectrographs on Keck (Keck Cosmic Web Imager; KCWI) and the Very Large Telescope (Multi-Unit Spectroscopic Explorer; MUSE). We compare the different $T_e$ measurements with HII region and interstellar medium environmental properties such as electron density, ionization parameter, molecular gas velocity dispersion, and stellar association/cluster mass and age obtained from PHANGS. We find that the temperatures from [OII] and [SII] are likely over-estimated due to the presence of electron density inhomogeneities in HII regions. We observe that differences between [NII] and [SIII] temperatures are weakly correlated with stellar association mass and molecular gas velocity dispersion. We measure high [OIII] temperatures in a subset of regions with high molecular gas velocity dispersion and low ionization parameter, which may be explained by the presence of low-velocity shocks. In agreement with previous studies, the $T_{\rm{e}}$--$T_{\rm{e}}$ between [NII] and [SIII] temperatures have the lowest observed scatter and generally follow predictions from photoionization modeling, which suggests that these tracers reflect HII region temperatures across the various ionization zones better than [OII], [SII], and [OIII]., Comment: Accepted to ApJ, 66 pages, 36 figures

We present a study of the molecular gas of seven early-type galaxies with high angular resolution data obtained as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project with the Atacama Large Millimeter/submillimeter Array. Using a fixed spatial scale approach, we study the mass surface density ($\Sigma$) and velocity dispersion ($\sigma$) of the molecular gas on spatial scales ranging from $60$ to $120$pc. Given the spatial resolution of our data ($20$ - $70$pc), we characterise these properties across many thousands of individual sight lines ($\approx50,000$ at our highest physical resolution). The molecular gas along these sight lines has a large range ($\approx2$dex) of mass surface densities and velocity dispersions $\approx40\%$ higher than those of star-forming spiral galaxies. It has virial parameters $\alpha_\mathrm{vir}$ that depend weakly on the physical scale observed, likely due to beam smearing of the bulk galactic rotation, and is generally super-virial. Comparing the internal turbulent pressure ($P_\mathrm{turb}$) to the pressure required for dynamic equilibrium ($P_\mathrm{DE}$), the ratio $P_\mathrm{turb}$/$P_\mathrm{DE}$ is significantly less than unity in all galaxies, indicating that the gas is not in dynamic equilibrium and is strongly compressed, in apparent contradiction to the virial parameters. This may be due to our neglect of shear and tidal forces, and/or the combination of three-dimensional and vertical diagnostics. Both $\alpha_\mathrm{vir}$ and $P_\mathrm{turb}$ anti-correlate with the global star-formation rate of our galaxies. We therefore conclude that the molecular gas in early-type galaxies is likely unbound, and that large-scale dynamics likely plays a critical role in its regulation. This contrasts to the giant molecular clouds in the discs of late-type galaxies, that are much closer to dynamical equilibrium., Comment: 32 pages (16 of Appendices), 39 Figures (27 in Appendices). Accepted for publication in MNRAS

We report the discovery of the ``mm fundamental plane of black-hole accretion'', which is a tight correlation between the nuclear 1 mm luminosity ($L_{\rm \nu, mm}$), the intrinsic $2$ -- $10$~keV X-ray luminosity ($L_{\rm X,2-10}$) and the supermassive black hole (SMBH) mass ($M_{\rm BH}$) with an intrinsic scatter ($\sigma_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies, most of which are low-luminosity active galactic nuclei (LLAGN). Combining these sources with a sample of high-luminosity (quasar-like) nearby AGN, we find that the plane still holds. We also find that $M_{\rm BH}$ correlates with $L_{\rm \nu, mm}$ at a highly significant level, although such correlation is less tight than the mm fundamental plane ($\sigma_{\rm int}=0.51$ dex). Crucially, we show that spectral energy distribution (SED) models for both advection-dominated accretion flows (ADAFs) and compact jets can explain the existence of these relations, which are not reproduced by the standard torus-thin accretion disc models usually associated to quasar-like AGN. The ADAF models reproduces the observed relations somewhat better than those for compact jets, although neither provides a perfect prediction. Our findings thus suggest that radiatively-inefficient accretion processes such as those in ADAFs or compact (and thus possibly young) jets may play a key role in both low- and high-luminosity AGN. This mm fundamental plane also offers a new, rapid method to (indirectly) estimate SMBH masses., Comment: 9 pages, 5 figures

We investigate the ionised gas morphology, excitation properties, and kinematics in 19 nearby star-forming galaxies from the PHANGS-MUSE survey. We directly compare the kinetic energy of expanding superbubbles and the turbulent motions in the interstellar medium with the mechanical energy deposited by massive stars in the form of winds and supernovae, with the aim to answer whether the stellar feedback is responsible for the observed turbulent motions and to quantify the fraction of mechanical energy retained in the superbubbles. Based on the distribution of the flux and velocity dispersion in the H$\alpha$ line, we select 1484 regions of locally elevated velocity dispersion ($\sigma$(H$\alpha$)>45 km/s), including at least 171 expanding superbubbles. We analyse these regions and relate their properties to those of the young stellar associations and star clusters identified in PHANGS-HST data. We find a good correlation between the kinetic energy of the ionised gas and the total mechanical energy input from supernovae and stellar winds from the stellar associations, with a typical coupling efficiency of 10-20%. The contribution of mechanical energy by the supernovae alone is not sufficient to explain the measured kinetic energy of the ionised gas, which implies that pre-supernova feedback in the form of radiation/thermal pressure and winds is necessary. We find that the gas kinetic energy decreases with metallicity for our sample covering Z=0.5-1.0 Zsun, reflecting the lower impact of stellar feedback. For the sample of superbubbles, we find that about 40% of the young stellar associations are preferentially located in their rims. We also find a slightly higher (by ~15%) fraction of the youngest (<3 Myr) stellar associations in the rims of the superbubbles than in the centres, and the opposite for older associations, which implies possible propagation or triggering of star formation., Comment: Accepted by A&A. The abstract is abridged. 31 pages (including 8 pages in appendix), 20 figures, 2 tables

Pattern speeds are a fundamental parameter of the dynamical features (e.g. bars, spiral arms) of a galaxy, setting resonance locations. Pattern speeds are not directly observable, so the Tremaine-Weinberg (TW) method has become the most common method used to measure them in galaxies. However, it has not been tested properly whether this method can straightforwardly be applied to gas tracers, despite this being widely done in the literature. When applied to observations, the TW method may return invalid results, which are difficult to diagnose due to a lack of ground truth for comparison. Although some works applying the TW method to simulated galaxies exist, only stellar populations have been tested. Therefore, here we explore the applicability of the TW method for gas gracers, by applying it to hydrodynamical simulations of galaxies, where we know the true value of the bar pattern speed. We perform some simple tests to see if the TW method has a physically reasonable output. We add different kinds of uncertainties (e.g. in position angle or flux) to the data to mock observational errors based on the magnitude of uncertainty present in the observations. Second, we test the method on 3D simulations with chemical networks. We show that in general, applying TW to observations of gas will not recover the true pattern speed. These results have implications for many "pattern speeds" reported in the literature, and based on these tests we also give some best practices for measuring pattern speeds using gas tracers going forwards., Comment: 9 pages, 8 figures, submitted to MNRAS

Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dust attenuation is complicated by the combined effects of dust heating from old stellar populations and an uncertain dust geometry with respect to heating sources. We use images obtained with NIRCam and MIRI as part of the PHANGS--JWST survey to calibrate dust emission at 21$\rm \mu m$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm \mu m$ as tracers of obscured star formation. We analyse $\sim$ 20000 optically selected HII regions across 19 nearby star-forming galaxies, and benchmark their IR emission against dust attenuation measured from the Balmer decrement. We model the extinction-corrected H$\alpha$ flux as the sum of the observed H$\alpha$ emission and a term proportional to the IR emission, with $a_{IR}$ as the proportionality coefficient. A constant $a_{IR}$ leads to extinction-corrected H$\alpha$ estimates which agree with those obtained with the Balmer decrement with a scatter of $\sim$ 0.1 dex for all bands considered. Among these bands, 21$\rm \mu m$ emission is demonstrated to be the best tracer of dust attenuation. The PAH-tracing bands underestimate the correction for bright HII regions, since in these environments the ratio of PAH-tracing bands to 21$\rm \mu m$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background., Comment: accepted for publication in A&A

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 onto 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 timescales of $\sim$10 Myr. Most of the inflow appears to be consumed by star formation in the ring which is currently occurring at a rate of ${\rm SFR}\simeq~1.8$-$2 \rm M_\odot\, yr^{-1}$, suggesting that the inflow is causally controlling the star formation rate in the ring as a function of time., Comment: Accepted in MNRAS

Mapping molecular line emission beyond the bright low-J CO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)mm interferometers, such as ALMA and NOEMA. We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1-0), HCO+(1-0), and even fainter lines in external galaxies. The goal is to study the capabilities and limitations of the stacking technique when applied to imaged interferometric observations. The core idea of spectral stacking is to align spectra of the low S/N spectral lines to a known velocity field calculated from a higher S/N line expected to share the kinematics of the fainter line, e.g., CO(1-0) or 21-cm emission. Then these aligned spectra can be coherently averaged to produce potentially high S/N spectral stacks. Here, we use imaged simulated interferometric and total power observations at different signal-to-noise levels, based on real CO observations. For the combined interferometric and total power data, we find that the spectral stacking technique is capable of recovering the integrated intensities even at low S/N levels across most of the region where the high S/N prior is detected. However, when stacking interferometer-only data for low S/N emission, the stacks can miss up to 50% of the emission from the fainter line. A key result of this analysis is that the spectral stacking method is able to recover the true mean line intensities in low S/N cubes and to accurately measure the statistical significance of the recovered lines. To facilitate the application of this technique we provide a public Python package, called PyStacker., Comment: 10 pages, 10 figures, accepted for pub in A&A, Apr 28, 2023

We present high spatial resolution ($\approx24$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(2-1) observations of the central region of the nearby barred spiral galaxy NGC 5806. NGC 5806 has a highly structured molecular gas distribution with a clear nucleus, a nuclear ring and offset dust lanes. We identify $170$ spatially- and spectrally-resolved giant molecular clouds (GMCs). These clouds have comparable sizes ($R_{\mathrm{c}}$) and larger gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities than those of clouds in the Milky Way disc. The size -- linewidth relation of the clouds is one of the steepest reported so far ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{1.20}$), the clouds are on average only marginally bound (with a mean virial parameter $\langle\alpha_{\mathrm{vir}}\rangle\approx2$), and high velocity dispersions are observed in the nuclear ring. These behaviours are likely due to bar-driven gas shocks and inflows along the offset dust lanes, and we infer an inflow velocity of $\approx120$ kms$^{-1}$ and a total molecular gas mass inflow rate of $\approx5$ M$_\odot$ yr$^{-1}$ into the nuclear ring. The observed internal velocity gradients of the clouds are consistent with internal turbulence. The number of clouds in the nuclear ring decreases with azimuthal angle downstream from the dust lanes without clear variation of cloud properties. This is likely due to the estimated short lifetime of the clouds ($\approx6$ Myr), which appears to be mainly regulated by cloud-cloud collision and/or shear processes. Overall, it thus seems that the presence of the large-scale bar and gas inflows to the centre of NGC 5806 affect cloud properties., Comment: Accepted for publication in MNRAS, 20 pages, 16 figures

We present a CO dynamical estimate of the mass of the super-massive black hole (SMBH) in three nearby early-type galaxies: NGC0612, NGC1574 and NGC4261. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3-6 observations of the $^{12}$CO(2-1) emission line with spatial resolutions of $14-58$ pc ($0.01"-0.26"$). We detect disc-like CO distributions on scales from $\lesssim200$ pc (NGC1574 and NGC4261) to $\approx10$ kpc (NGC0612). In NGC0612 and NGC1574 the bulk of the gas is regularly rotating. The data also provide evidence for the presence of a massive dark object at the centre of NGC1574, allowing us to obtain the first measure of its mass, $M_{\rm BH}=(1.0\pm0.2)\times10^{8}$ M$_{\odot}$ (1$\sigma$ uncertainty). In NGC4261, the CO kinematics is clearly dominated by the SMBH gravitational influence, allowing us to determine an accurate black hole mass of $(1.62{\pm 0.04})\times10^{9}$ M$_{\odot}$ ($1\sigma$ uncertainty). This is fully consistent with a previous CO dynamical estimate obtained using a different modelling technique. Signs of non-circular gas motions (likely outflow) are also identified in the inner regions of NGC4261. In NGC0612, we are only able to obtain a (conservative) upper limit of $M_{\rm BH}\lesssim3.2\times10^{9}$ M$_{\odot}$. This has likely to be ascribed to the presence of a central CO hole (with a radius much larger than that of the SMBH sphere of influence), combined with the inability of obtaining a robust prediction for the CO velocity curve. The three SMBH mass estimates are overall in agreement with predictions from the $M_{\rm BH}-\sigma_{\star}$ relation., Comment: Main text: 20 pages, 14 Figures; Appendix: 7 pages, 6 Figures. Accepted for publication in MNRAS on 2023 March 28

The CO-to-H$_2$ conversion factor ($\alpha_\rm{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 $\alpha_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 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 multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $\alpha_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $\alpha_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $\alpha_\rm{CO}$ with low-J $^{12}$CO optical depths ($\tau_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $\tau_\rm{CO}$ correlation explains most of the $\alpha_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $\alpha_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $\tau_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $\alpha_\rm{CO}$ variation. We also test current simulation-based $\alpha_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations., Comment: Accepted for publication in ApJ; 30 pages of main text + 3 appendices

We present new HI observations of the nearby massive spiral galaxy M83, taken with the VLA at $21^{\prime\prime}$ angular resolution ($\approx500$ pc) of an extended ($\sim$1.5 deg$^2$) 10-point mosaic combined with GBT single dish data. We study the super-extended HI disk of M83 (${\sim}$50 kpc in radius), in particular disc kinematics, rotation and the turbulent nature of the atomic interstellar medium. We define distinct regions in the outer disk ($r_{\rm gal}>$central optical disk), including ring, southern area, and southern and northern arm. We examine HI gas surface density, velocity dispersion and non-circular motions in the outskirts, which we compare to the inner optical disk. We find an increase of velocity dispersion ($\sigma_v$) towards the pronounced HI ring, indicative of more turbulent HI gas. Additionally, we report over a large galactocentric radius range (until $r_{\rm gal}{\sim}$50 kpc) that $\sigma_v$ is slightly larger than thermal (i.e. $>8$km s$^{-1}$ ). We find that a higher star formation rate (as traced by FUV emission) is not always necessarily associated with a higher HI velocity dispersion, suggesting that radial transport could be a dominant driver for the enhanced velocity dispersion. We further find a possible branch that connects the extended HI disk to the dwarf irregular galaxy UGCA365, that deviates from the general direction of the northern arm. Lastly, we compare mass flow rate profiles (based on 2D and 3D tilted ring models) and find evidence for outflowing gas at r$_{\rm gal}$ $\sim$2 kpc, inflowing gas at r$_{\rm gal}$ $\sim$5.5 kpc and outflowing gas at r$_{\rm gal}$ $\sim$14 kpc. We caution that mass flow rates are highly sensitive to the assumed kinematic disk parameters, in particular, to the inclination., Comment: accepted for publication in A&A; 16 pages, 12 figures (+8 pages appendix)

Connecting the gas in HII regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how HII regions evolve over time. With PHANGS$\unicode{x2013}$MUSE we detect nearly 24,000 HII regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, density). We use catalogues of multi-scale stellar associations from PHANGS$\unicode{x2013}$HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4,177 HII regions that are clearly linked to a single ionizing association. A weak anti-correlation is observed between the association ages and the H$\alpha$ equivalent width EW(H$\alpha$), the H$\alpha$/FUV flux ratio and the ionization parameter, log q. As all three are expected to decrease as the stellar population ages, this could indicate that we observe an evolutionary sequence. This interpretation is further supported by correlations between all three properties. Interpreting these as evolutionary tracers, we find younger nebulae to be more attenuated by dust and closer to giant molecular clouds, in line with recent models of feedback-regulated star formation. We also observe strong correlations with the local metallicity variations and all three proposed age tracers, suggestive of star formation preferentially occurring in locations of locally enhanced metallicity. Overall, EW(H$\alpha$) and log q show the most consistent trends and appear to be most reliable tracers for the age of an HII region., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS