Search

Connecting the gas in H ii regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how H ii regions evolve over time. With PHANGS–MUSE, we detect nearly 24 000 H ii regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, and density). We use catalogues of multiscale stellar associations from PHANGS–HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4177 H ii regions that are clearly linked to a single ionizing association. A weak anticorrelation is observed between the association ages and the $\mathrm{H}\, \alpha$ equivalent width $\mathrm{EW}(\mathrm{H}\, \alpha)$, the $\mathrm{H}\, \alpha/\mathrm{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, $\mathrm{EW}(\mathrm{H}\, \alpha)$ and log q show the most consistent trends and appear to be most reliable tracers for the age of an H ii region.

Nebular HeII$\lambda$4686\AA~line emission is useful to unveil active galactic nuclei (AGN) residing in actively star-forming (SF) galaxies, typically missed by the standard BPT classification. Here we adopt the HeII diagnostic to identify hidden AGN in the Local Universe using for the first time spatially-resolved data from the Data Release 15 of the Mapping Nearby Galaxies at APO survey (MaNGA DR15). By combining results from HeII and BPT diagnostics, we overall select 459 AGN host candidates ($\sim$10% in MaNGA DR15), out of which 27 are identified as AGN by the HeII diagram only. The HeII-only AGN population is hosted by massive (M$_*\gtrsim10^{10}$ M$_{\odot}$) SF Main Sequence galaxies, and on average less luminous than the BPT-selected AGN. Given the HeII line faintness, we revisit our census accounting for incompleteness effects due to the HeII sensitivity limit of MaNGA. We thus obtain an overall increased fraction (11%) of AGN in MaNGA compared to the BPT-only census (9%), which further increases to 14% for galaxies more massive than $10^{10}$ M$_{\odot}$; interestingly, on the SF Main Sequence the increase is by about a factor of 2. A substantial number of AGN in SF galaxies points to significant, coeval star formation and black hole accretion, consistently with results from hydrodynamical simulations and with important implications on quenching scenarios. In view of exploring unprecedented high redshifts with JWST and new ground-based facilities, revisiting the standard BPT classification through novel emission-line diagnostics is fundamental to discover AGN in highly SF environments., Comment: 13 pages, 12 figures. Accepted for publication in MNRAS

We present a comparison of the nitrogen-to-oxygen ratio (N/O) in 37 high-redshift galaxies at $z\sim$2 taken from the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey with a comparison sample of local galaxies, taken from the Sloan Digital Sky Survey (SDSS). The KLEVER sample shows only a mild enrichment in N/O of $+$0.1 dex when compared to local galaxies at a given gas-phase metallicity (O/H), but shows a depletion in N/O of $-$0.36 dex when compared at a fixed stellar mass (M$_*$). We find a strong anti-correlation in local galaxies between N/O and SFR in the M$_*$-N/O plane, similar to the anti-correlation between O/H and SFR found in the mass-metallicity relation (MZR). We use this anti-correlation to construct a fundamental nitrogen relation (FNR), analogous to the fundamental metallicity relation (FMR). We find that KLEVER galaxies are consistent with both the FMR and the FNR. This suggests that the depletion of N/O in high-$z$ galaxies when considered at a fixed M$_*$ is driven by the redshift-evolution of the mass-metallicity relation in combination with a near redshift-invariant N/O-O/H relation. Furthermore, the existence of an fundamental nitrogen relation suggests that the mechanisms governing the fundamental metallicity relation must be probed by not only O/H, but also N/O, suggesting pure-pristine gas inflows are not the primary driver of the FMR, and other properties such as variations in galaxy age and star formation efficiency must be important., Comment: 23 pages, 11 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

We investigate which physical properties are most predictive of the position of local star forming galaxies on the BPT diagrams, by means of different Machine Learning (ML) algorithms. Exploiting the large statistics from the Sloan Digital Sky Survey (SDSS), we define a framework in which the deviation of star-forming galaxies from their median sequence can be described in terms of the relative variations in a variety of observational parameters. We train artificial neural networks (ANN) and random forest (RF) trees to predict whether galaxies are offset above or below the sequence (via classification), and to estimate the exact magnitude of the offset itself (via regression). We find, with high significance, that parameters primarily associated to variations in the nitrogen-over-oxygen abundance ratio (N/O) are the most predictive for the [N II]-BPT diagram, whereas properties related to star formation (like variations in SFR or EW[H$\alpha$]) perform better in the [S II]-BPT diagram. We interpret the former as a reflection of the N/O-O/H relationship for local galaxies, while the latter as primarily tracing the variation in the effective size of the S$^{+}$ emitting region, which directly impacts the [S II]emission lines. This analysis paves the way to assess to what extent the physics shaping local BPT diagrams is also responsible for the offsets seen in high redshift galaxies or, instead, whether a different framework or even different mechanisms need to be invoked., Comment: Accepted for publication on MNRAS (31 pages, 20 figures)

Galaxy spectra encode in their continuum and absorption/emission features a wealth of information on galaxy physics, mass assembly and chemical enrichment history. The 4MOST-StePS survey will collect high-quality spectra (with a median signal-to-noise ratio of about 30 Å–1, and resolution R ~ 5000) for a sample of about 3300 galaxies brighter than IAB = 20.5 within the RA-Dec-z footprint of the WAVES-Deep survey. These spectra will provide a precise empirical description of the evolutionary path of massive galaxies in the intermediate redshift range (0.3 < z < 0.7) between the LEGA-C and SDSS surveys. The locations of the galaxies within the cosmic web, unveiled by WAVES-Deep, will disclose the connection between galaxy properties and environment, down to the scales of galaxy pairs., Published in The Messenger vol. 190, pp. 22-24, March 2023.

We use new HCN(1-0) data from the ALMOND (ACA Large-sample Mapping Of Nearby galaxies in Dense gas) survey to trace the kpc-scale molecular gas density structure and CO(2-1) data from PHANGS-ALMA to trace the bulk molecular gas across 25 nearby, star-forming galaxies. At 2.1 kpc scale, we measure the density-sensitive HCN/CO line ratio and the SFR/HCN ratio to trace the star formation efficiency in the denser molecular medium. At 150 pc scale, we measure structural and dynamical properties of the molecular gas via CO(2-1) line emission, which is linked to the lower resolution data using an intensity-weighted averaging method. We find positive correlations (negative) of HCN/CO (SFR/HCN) with the surface density, the velocity dispersion and the internal turbulent pressure of the molecular gas. These observed correlations agree with expected trends from turbulent models of star formation, which consider a single free-fall time gravitational collapse. Our results show that the kpc-scale HCN/CO line ratio is a powerful tool to trace the 150 pc scale average density distribution of the molecular clouds. Lastly, we find systematic variations of the SFR/HCN ratio with cloud-scale molecular gas properties, which are incompatible with a universal star formation efficiency. Overall, these findings show that mean molecular gas density, molecular cloud properties and star formation are closely linked in a coherent way, and observations of density-sensitive molecular gas tracers are a useful tool to analyse these variations, linking molecular gas physics to stellar output across galaxy discs., Comment: 48 pages, 40 figures

We present maps of the 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature in NGC 628, NGC 1365, and NGC 7496 as observed with the Near-Infrared Camera (NIRCam) imager on JWST from the PHANGS-JWST Cycle 1 Treasury project. We create maps that isolate the 3.3 micron PAH feature in the F335M filter (F335M$_{\rm PAH}$) using combinations of the F300M and F360M filters for removal of starlight continuum. This continuum removal is complicated by contamination of the F360M by PAH emission and variations in the stellar spectral energy distribution slopes between 3.0 and 3.6 micron. We modify the empirical prescription from Lai et al. (2020) to remove the starlight continuum in our highly resolved galaxies, which have a range of starlight- and PAH-dominated lines-of-sight. Analyzing radially binned profiles of the F335M$_{\rm PAH}$ emission, we find that between 5-65% of the F335M intensity comes from the 3.3 micron feature within the inner 0.5 $r_{25}$ of our targets. This percentage systematically varies from galaxy to galaxy, and shows radial trends within the galaxies related to each galaxy's distribution of stellar mass, interstellar medium, and star formation. The 3.3 micron emission is well correlated with the 11.3 micron PAH feature traced with the MIRI F1130W filter, as is expected, since both features arise from C-H vibrational modes. The average F335M$_{\rm PAH}$/F1130W ratio agrees with the predictions of recent models by Draine et al. (2021) for PAHs with size and charge distributions shifted towards larger grains with normal or higher ionization., Comment: submitted to AAS journals and revised according to referee comments, part of a PHANGS-JWST Focus Issue to appear in ApJ

We compare mid-infrared (mid-IR), extinction-corrected H$\alpha$, 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$\alpha$. At these scales, CO and H$\alpha$ each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H$\alpha$ emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to the good match with the molecular gas-tracing CO, and as a heating tracer, leading to the good match with the H$\alpha$. By combining mid-IR, CO, and H$\alpha$ 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_\nu = 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$\alpha$-tracing component. The best-fitting models imply that $\sim 50\%$ of the mid-IR flux arises from molecular gas heated by the diffuse interstellar radiation field, with the remaining $\sim 50\%$ associated with bright, dusty star forming regions. We discuss differences between the F770W, F1000W, F1130W bands and the continuum dominated F2100W band and suggest next steps for using the mid-IR as an ISM tracer., Comment: 49 pages, 17 figures, Section 8 provides a detailed summary, accepted for publication in ApJ Letters, part of a PHANGS-JWST Focus Issue to appear in ApJ

The PHANGS collaboration has been building a reference data set for the multiscale, multiphase study of star formation and the interstellar medium (ISM) in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds (∼5–50 pc). In Cycle 1, PHANGS is conducting an eight-band imaging survey from 2 to 21 μm of 19 nearby spiral galaxies. Optical integral field spectroscopy, CO(2–1) mapping, and UV-optical imaging for all 19 galaxies have been obtained through large programs with ALMA, VLT-MUSE, and Hubble. PHANGS–JWST enables a full inventory of star formation, accurate measurement of the mass and age of star clusters, identification of the youngest embedded stellar populations, and characterization of the physical state of small dust grains. When combined with Hubble catalogs of ∼10,000 star clusters, MUSE spectroscopic mapping of ∼20,000 H II regions, and ∼12,000 ALMA-identified molecular clouds, it becomes possible to measure the timescales and efficiencies of the earliest phases of star formation and feedback, build an empirical model of the dependence of small dust grain properties on local ISM conditions, and test our understanding of how dust-reprocessed starlight traces star formation activity, all across a diversity of galactic environments. Here we describe the PHANGS–JWST Treasury survey, present the remarkable imaging obtained in the first few months of science operations, and provide context for the initial results presented in the first series of PHANGS–JWST publications., La lista completa de autores que integran el documento puede consultarse en el archivo., Instituto de Astrofísica de La Plata

We study the relations between the mass of the central black hole (BH) $M_{\rm BH}$, the dark matter halo mass $M_{\rm h}$, and the stellar-to-halo mass fraction $f_\star\propto M_\star/M_{\rm h}$ in a sample of $55$ nearby galaxies with dynamically measured $M_{\rm BH}>10^6\,{\rm M}_\odot$ and $M_{\rm h}>5\times10^{11}\,{\rm M}_\odot$. The main improvement with respect to previous studies is that we consider both early- and late-type systems for which $M_{\rm h}$ is determined either from globular cluster dynamics or from spatially resolved rotation curves. Independently of their structural properties, galaxies in our sample build a well defined sequence in the $M_{\rm BH}$-$M_{\rm h}$-$f_\star$ space. We find that: (i) $M_{\rm h}$ and $M_{\rm BH}$ strongly correlate with each other and anti-correlate with $f_\star$; (ii) there is a break in the slope of the $M_{\rm BH}$-$M_{\rm h}$ relation at $M_{\rm h}$ of $10^{12}\,{\rm M}_\odot$, and in the $f_\star$-$M_{\rm BH}$ relation at $M_{\rm BH}$ of $\sim10^7\!-\!10^8\,{\rm M}_\odot$; (iii) at a fixed $M_{\rm BH}$, galaxies with a larger $f_\star$ tend to occupy lighter halos and to have later morphological types. We show that the observed trends can be reproduced by a simple equilibrium model in the $��$CDM framework where galaxies smoothly accrete dark and baryonic matter at a cosmological rate, having their stellar and black hole build-up regulated both by the cooling of the available gas reservoir and by the negative feedback from star formation and active galactic nuclei (AGN). Feature (ii) arises as the BH population transits from a rapidly accreting phase to a more gentle and self-regulated growth, while scatter in the AGN feedback efficiency can account for feature (iii)., 20 pages, 12 figures, accepted by MNRAS. Key figures are Fig.1 and Fig.6

The Diffuse Ionized Gas (DIG) contributes to the nebular emission of galaxies, resulting in emission line flux ratios that can be significantly different from those produced by HII regions. Comparing the emission of [SII]6717,31 between pointed observations of HII regions in nearby galaxies and integrated spectra of more distant galaxies, it has been recently claimed that the DIG can also deeply affect the emission of bright, star-forming galaxies, and that a large correction must be applied to observed line ratios to recover the genuine contribution from HII regions. Here we show instead that the effect of DIG on the integrated spectra of star-forming galaxies is lower than assumed in previous work. Indeed, aperture effects on the spectroscopy of nearby HII regions are largely responsible for the observed difference: when spectra of local HII regions are extracted using large enough apertures while still avoiding the DIG, the observed line ratios are the same as in more distant galaxies. This result is highly relevant for the use of strong-line methods to measure metallicity., Comment: MNRAS, in press

A fundamental assumption adopted in nearly every extragalactic emission-line study is that the attenuation of different emission lines can be described by a single attenuation curve. Here we show this assumption fails in many cases with important implications for derived results. We developed a new method to measure the differential nebular attenuation among three kinds of transitions: the Balmer lines of hydrogen, high-ionization transitions, and low-ionization transitions. This method bins the observed data in a multidimensional space spanned by attenuation-insensitive line ratios. Within each small bin, the variations in line ratios are mainly driven by the variations in the nebular attenuation. This allows us to measure the nebular attenuation using both forbidden lines and Balmer lines. We applied this method to a sample of 2.4 million star-forming spaxels from SDSS-IV MaNGA. We found that the attenuation of high ionization lines and Balmer lines can be well described by a single Fitzpatrick (1999) extinction curve with $R_V=3.1$. However, no single attenuation curve can simultaneously account for all three transitions. This strongly suggests that different lines have different effective attenuations, likely because spectroscopy at kiloparsec resolutions mixes multiple regions with different intrinsic line ratios and different levels of attenuation. As a result, the assumption that different lines follow the same attenuation curve breaks down. Using a single attenuation curve determined by Balmer lines to correct attenuation-sensitive forbidden line ratios could bias the nebular parameters derived by 0.06--0.25 dex at $A_V = 1$, depending on the details of the dust attenuation model. Observations of a statistically large sample of H II regions with high spatial resolutions and large spectral coverage are vital for improved modeling and deriving accurate corrections for this effect., 28 pages, 17 figures, published in A&A

We analyse the chemical properties of three z~8 galaxies behind the galaxy cluster SMACS J0723.3-7327, observed as part of the Early Release Observations programme of the James Webb Space Telescope (JWST). Exploiting [O III]4363 auroral line detections in NIRSpec spectra, we robustly apply the direct Te method for the very first time at such high redshift, measuring metallicities ranging from extremely metal poor (12+log(O/H)~7) to about one-third solar. We also discuss the excitation properties of these sources, and compare them with local strong-line metallicity calibrations. We find that none of the considered diagnostics match simultaneously the observed relations between metallicity and strong-line ratios for the three sources, implying that a proper re-assessment of the calibrations may be needed at these redshifts. On the mass-metallicity plane, the two galaxies at z~7.6 (log(M*/M_sun) = 8.1, 8.7) have metallicities that are consistent with the extrapolation of the mass-metallicity relation at z~2-3, while the least massive galaxy at z~8.5 (log(M*/M_sun) = 7.8) shows instead a significantly lower metallicity . The three galaxies show different level of offset relative to the Fundamental Metallicity Relation, with two of them (at z~7.6) being marginally consistent, while the z~8.5 source deviating significantly, being probably far from the smooth equilibrium between gas flows, star formation and metal enrichment in place at later epochs., 13 pages, 7 figures. Accepted for publication on MNRAS

We analyze spatially resolved and co-added SDSS-IV MaNGA spectra with signal-to-noise ~100 from 2200 passive central galaxies (z~0.05) to understand how central galaxy assembly depends on stellar mass (M*) and halo mass (Mh). We control for systematic errors in Mh by employing a new group catalog from Tinker (2020a,b) and the widely-used Yang et al. (2007) catalog. At fixed M*, the strength of several stellar absorption features varies systematically with Mh. Completely model-free, this is one of the first indications that the stellar populations of centrals with identical M* are affected by the properties of their host halos. To interpret these variations, we applied full spectral fitting with the code alf. At fixed M*, centrals in more massive halos are older, show lower [Fe/H], and have higher [Mg/Fe] with 3.5 sigma confidence. We conclude that halos not only dictate how much M* galaxies assemble, but also modulate their chemical enrichment histories. Turning to our analysis at fixed Mh, high-M* centrals are older, show lower [Fe/H], and have higher [Mg/Fe] for Mh>10^{12}Msun/h with confidence > 4 sigma. While massive passive galaxies are thought to form early and rapidly, our results are among the first to distinguish these trends at fixed Mh. They suggest that high-M* centrals experienced unique early formation histories, either through enhanced collapse and gas fueling, or because their halos were early-forming and highly concentrated, a possible signal of galaxy assembly bias., Accepted for publication in ApJ. 28 pages and 12 figures

We study quenching in seven green valley galaxies on kpc scales by resolving their molecular gas content using \textsuperscript{12}CO(1-0) observations obtained with NOEMA and ALMA, and their star-formation rate using spatially resolved optical spectroscopy from the MaNGA survey. We perform radial stacking of both datasets to increase the sensitivity to molecular gas and star formation, thereby avoiding biases against strongly quenched regions. We find that both spatially resolved gas fraction ($\rm {f_{gas}}$) and star formation efficiency ($\rm {SFE}$) are responsible for quenching green valley galaxies at all radii: both quantities are suppressed with respect to typical star-forming regions. $\rm {f_{gas}}$ and $\rm {SFE}$ have roughly equal influence in quenching the outer disc. We are, however, unable to identify the dominant mechanism in the strongly quenched central regions. We find that $\rm {f_{gas}}$ is reduced by $\rm \sim 1~dex$ in the central regions, but the star formation rate is too low to be measured, leading to upper limits for the $\rm {SFE}$. Moving from the outer disc to central regions, the reduction in $\rm {f_{gas}}$ is driven by an increasing $\rm \Sigma_{\star}$ profile rather than a decreasing $\rm \Sigma_{H_{2}}$ profile. The reduced $\rm {f_{gas}}$ may therefore be caused by a decrease in the gas supply rather than molecular gas ejection mechanisms, such as winds driven by active galactic nuclei. We warn more generally that studies investigating $\rm {f_{gas}}$ may be deceiving in inferring the cause of quenching, particularly in the central (bulge-dominated) regions of galaxies., Comment: 8 pages, 11 figures, Accepted for publication in MNRAS Letters

We utilize the ALMA-MaNGA QUEnch and STar formation (ALMaQUEST) survey to investigate the kpc-scale scaling relations, presented as the resolved star-forming MS (rSFMS: ΣSFR versus Σ*), the resolved Schmidt–Kennicutt relation (rSK: ΣSFR versus Σ H 2 ), and the resolved molecular gas MS (rMGMS: Σ H 2 versus Σ*), for 11,478 star-forming and 1414 retired spaxels (oversampled by a factor of ∼20) located in 22 GV and 12 MS galaxies. For a given galaxy type (MS or GV), the retired spaxels are found to be offset from the sequences formed by the star-forming spaxels on the rSFMS, rSK, and rMGMS planes, toward lower absolute values of sSFR, SFE, and f H 2 by ∼1.1, 0.6, and 0.5 dex. The scaling relations for GV galaxies are found to be distinct from that of the MS galaxies, even if the analyses are restricted to the star-forming spaxels only. It is found that, for star-forming spaxels, sSFR, SFE, and f H 2 in GV galaxies are reduced by ∼0.36, 0.14, and 0.21 dex, respectively, compared to those in MS galaxies. Therefore, the suppressed sSFR/SFE/f gas in GV galaxies is associated with not only an increased proportion of retired regions in GV galaxies but also a depletion of these quantities in star-forming regions. Finally, the reduction of SFE and f H 2 in GV galaxies relative to MS galaxies is seen in both bulge and disk regions (albeit with larger uncertainties), suggesting that, statistically, quenching in the GV population may persist from the inner to the outer regions.

We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2 10.8$, AGN-dominated galaxies, suggesting that AGNs may be the primary drivers of these gas flows. Surprisingly, at $\log(M_\star/M_{\odot})\leq 9.6$, the observed line profiles are fully consistent with a rotating disc model, indicating that ionised gas outflows in dwarf galaxies might play a negligible role even during the peak of cosmic star-formation activity. Finally, we find that the observed mass loading factor scales with stellar mass as expected from the TNG50 cosmological simulation, but the ionised gas mass accounts for only 2$\%$ of the predicted value. This suggests that either the bulk of the outflowing mass is in other gaseous phases or the current feedback models implemented in cosmological simulations need to be revised., Comment: Accepted for publication in MNRAS, 28 pages, 20 figures

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$\alpha$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$\alpha$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$\alpha$ (PHANGS-ALMA and PHANGS-H$\alpha$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$\alpha$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$\alpha$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$\alpha$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process., Comment: 60 pages, 22 figures, 9 tables, accepted for publication in The Astrophysical Journal

Aims: Temperature uncertainties plague our understanding of abundance variations within the ISM. Using the PHANGS-MUSE large program, we develop and apply a new technique to model the strong emission lines arising from HII regions in 19 nearby spiral galaxies at ~50 pc resolution and infer electron temperatures for the nebulae. Methods: Due to the charge-exchange coupling of the ionization fraction of the atomic oxygen to that of hydrogen, the emissivity of the observed [OI]6300/Ha line ratio can be modeled as a function of gas phase oxygen abundance (O/H), ionization fraction (f_ion) and electron temperature (T_e). We measure (O/H) using a strong line metallicity calibration, and identify a correlation between f_ion and [SIII]9069/[SII]6716,6730, tracing ionization parameter variations. Results: We solve for T_e, and test the method by reproducing direct measurements of T_e([NII]5755) based on auroral line detections to within ~600 K. We apply this charge-exchange method of calculating T_e to 4,129 HII regions across 19 PHANGS-MUSE galaxies. We uncover radial temperature gradients, increased homogeneity on small scales, and azimuthal temperature variations in the disks that correspond to established abundance patterns. This new technique for measuring electron temperatures leverages the growing availability of optical integral field unit spectroscopic maps across galaxy samples, increasing the statistics available compared to direct auroral line detections., Comment: 21 pages, 22 figures, accepted by A&A

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre., 48 pages, 28 figures, 4 tables, accepted for pub in A&A, Dec 8, 2021. Data products released publicly at this address https://archive.eso.org/scienceportal/home?data_collection=PHANGS

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.

We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies., Comment: 34 pages; 12 figures; accepted for publication in the Astrophysical Journal; data tables available http://www.astronomy.ohio-state.edu/~leroy.42/int_corat_tab.ecsv and http://www.astronomy.ohio-state.edu/~leroy.42/co_model.ecsv before publication

We investigate the nature of the scaling relations between the surface density of star formation rate ($\Sigma _{\rm SFR}$), stellar mass ($\Sigma _*$), and molecular gas mass ($\Sigma _{\rm H_2}$), aiming at distinguishing between the relations that are primary, i.e. more fundamental, and those which are instead an indirect by-product of the other relations. We use the ALMaQUEST survey and analyse the data by using both partial correlations and Random Forest regression techniques. We unambiguously find that the strongest intrinsic correlation is between $\Sigma _{\rm SFR}$ and $\Sigma_{\rm H_2}$ (i.e. the resolved Schmidt-Kennicutt relation), followed by the correlation between $\Sigma _{\rm H_2}$ and $\Sigma _*$ (resolved Molecular Gas Main Sequence, rMGMS). Once these two correlations are taken into account, we find that there is no evidence for any intrinsic correlation between $\Sigma _{\rm SFR}$ and $\Sigma _*$, implying that SFR is entirely driven by the amount of molecular gas, while its dependence on stellar mass (i.e. the resolved Star Forming Main Sequence, rSFMS) simply emerges as a consequence of the relationship between molecular gas and stellar mass., Comment: Accepted for publication in MNRAS, 7 pages, 4 figures

Mapping star-formation rates (SFR) within galaxies is key to unveiling their assembly and evolution. Calibrations exist for computing SFR from a combination of ultraviolet and infrared bands for galaxies as integrated systems, but their applicability to sub-galactic (kpc) scales remains largely untested. Here we use integral field spectroscopy of 19 nearby ($D, Comment: accepted in A&A

The SDSS-IV Mapping Nearby Galaxies at APO (MaNGA) program has been operating from 2014-2020, and has now observed a sample of 9,269 galaxies in the low redshift universe (z ~ 0.05) with integral-field spectroscopy. With rest-optical (\lambda\lambda 0.36 - 1.0 um) spectral resolution R ~ 2000 the instrumental spectral line-spread function (LSF) typically has 1sigma width of about 70 km/s, which poses a challenge for the study of the typically 20-30 km/s velocity dispersion of the ionized gas in present-day disk galaxies. In this contribution, we present a major revision of the MaNGA data pipeline architecture, focusing particularly on a variety of factors impacting the effective LSF (e.g., undersampling, spectral rectification, and data cube construction). Through comparison with external assessments of the MaNGA data provided by substantially higher-resolution R ~ 10,000 instruments we demonstrate that the revised MPL-10 pipeline measures the instrumental line spread function sufficiently accurately ( 50. Velocity dispersions derived from [O II], Hbeta, [O III], [N II], and [S II] are consistent with those derived from Halpha to within about 2% at sigma_Halpha > 30 km/s. Although the impact of these changes to the estimated LSF will be minimal at velocity dispersions greater than about 100 km/s, scientific results from previous data releases that are based on dispersions far below the instrumental resolution should be reevaulated., Comment: 26 pages, 23 figures. Accepted for publication in AJ

The feedback from young stars (i.e. pre-supernova) is thought to play a crucial role in molecular cloud destruction. In this paper, we assess the feedback mechanisms acting within a sample of 5810 HII regions identified from the PHANGS-MUSE survey of 19 nearby ($ 1$, and expanding, yet there is a small sample of compact HII regions with $P_\mathrm{tot,max}/P_\mathrm{de} < 1$ ($\sim$1% of the sample). These mostly reside in galaxy centres ($R_\mathrm{gal}, Comment: 28 pages total, 16+3 figures, and 3 tables. Accepted for publication in MNRAS

We present PHANGS-ALMA, the first survey to map CO J=2-1 line emission at ~1" ~ 100pc spatial resolution from a representative sample of 90 nearby (d, Comment: 76 pages, 33 figures. Accepted for publication in the Astrophysical Journal Supplement series. Full resolution version and the image atlas to appear as a figure set in the published version can be found https://sites.google.com/view/phangs/publications . Data release coming soon to the ALMA archive and CADC temporarily available at http://phangs.org/data

It is still poorly constrained how the densest phase of the interstellar medium varies across galactic environment. A large observing time is required to recover significant emission from dense molecular gas at high spatial resolution, and to cover a large dynamic range of extragalactic disc environments. We present new NOrthern Extended Millimeter Array (NOEMA) observations of a range of high critical density molecular tracers (HCN, HNC, HCO+) and CO isotopologues (13CO, C18O) towards the nearby (11.3 Mpc), strongly barred galaxy NGC 3627. These observations represent the current highest angular resolution (1.85"; 100 pc) map of dense gas tracers across a disc of a nearby spiral galaxy, which we use here to assess the properties of the dense molecular gas, and their variation as a function of galactocentric radius, molecular gas, and star formation. We find that the HCN(1-0)/CO(2-1) integrated intensity ratio does not correlate with the amount of recent star formation. Instead, the HCN(1-0)/CO(2-1) ratio depends on the galactic environment, with differences between the galaxy centre, bar, and bar end regions. The dense gas in the central 600 pc appears to produce stars less efficiently despite containing a higher fraction of dense molecular gas than the bar ends where the star formation is enhanced. In assessing the dynamics of the dense gas, we find the HCN(1-0) and HCO+(1-0) emission lines showing multiple components towards regions in the bar ends that correspond to previously identified features in CO emission. These features are co-spatial with peaks of Halpha emission, which highlights that the complex dynamics of this bar end region could be linked to local enhancements in the star formation., Comment: 25 pages, 12 figures; accepted for publication in MNRAS