Showing 1–31 of 31 results for author: Bešlić, I

A 260 pc resolution ALMA map of HCN(1-0) in the galaxy NGC 4321

Authors: Lukas Neumann, Frank Bigiel, Ashley T. Barnes, Molly J. Gallagher, Adam Leroy, Antonio Usero, Erik Rosolowsky, Ivana Bešlić, Médéric Boquien, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Cosima Eibensteiner, Kathryn Grasha, Jonathan D. Henshaw, María J. Jiménez-Donaire, Sharon Meidt, Shyam H. Menon, Eric J. Murphy, Hsi-An Pan, Miguel Querejeta, Toshiki Saito, Eva Schinnerer, Sophia K. Stuber , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 18 pages, 9 figures, accepted for pub in A&A, Jun 13, 2024

Do spiral arms enhance star formation efficiency?

Authors: Miguel Querejeta, Adam K. Leroy, Sharon E. Meidt, Eva Schinnerer, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Jiayi Sun, Mattia Sormani, Ivana Bešlic, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Santiago García-Burillo, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Eric. W. Koch, Lukas Neumann, Hsi-An Pan, Ismael Pessa, Jérôme Pety, Francesca Pinna, Lise Ramambason , et al. (10 additional authors not shown)

Abstract: 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 fr… ▽ More 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) △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 26 pages, 16 figures. Accepted for publication in A&A

Bias versus variance when fitting multi-species molecular lines with a non-LTE radiative transfer model

Authors: Antoine Roueff, Jérôme Pety, Maryvonne Gerin, Léontine Ségal, Javier Goicoechea, Harvey Liszt, Pierre Gratier, Ivana Bešlić, Lucas Einig, M. Gaudel, Jan Orkisz, Pierre Palud, Miriam Santa-Maria, Victor de Souza Magalhaes, Antoine Zakardjian, Sebastien Bardeau, Emeric E. Bron, Pierre Chainais, Simon Coudé, Karine Demyk, Viviana Guzman Veloso, Annie Hughes, David Languignon, François Levrier, Dariusz C Lis , et al. (6 additional authors not shown)

Abstract: Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity… ▽ More Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity of the problem introduce estimation biases and how to detect them.We focus on the CO and HCO+ isotopologues and analyze maps of a 50 square arcminutes field. We used the RADEX escape probability model to solve the statistical equilibrium equations and compute the emerging line profiles, assuming that all species coexist. Depending on the considered set of species, we also fixed the abundance ratio between some species and explored different values. We proposed a maximum likelihood estimator to infer the physical conditions and considered the effect of both the thermal noise and calibration uncertainty. We analyzed any potential biases induced by model misspecifications by comparing the results on the actual data for several sets of species and confirmed with Monte Carlo simulations. The variance of the estimations and the efficiency of the estimator were studied based on the Cram{é}r-Rao lower bound.Column densities can be estimated with 30% accuracy, while the best estimations of the volume density are found to be within a factor of two. Under the chosen model framework, the peak 12CO(1--0) is useful for constraining the kinetic temperature. The thermal pressure is better and more robustly estimated than the volume density and kinetic temperature separately. Analyzing CO and HCO+ isotopologues and fitting the full line profile are recommended practices with respect to detecting possible biases.Combining a non-local thermodynamic equilibrium model with a rigorous analysis of the accuracy allows us to obtain an efficient estimator and identify where the model is misspecified. We note that other combinations of molecular lines could be studied in the future. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: Astronomy and Astrophysics - A\&A, In press

The properties and kinematics of HCN emission across the closest starburst galaxy NGC 253 observed with ALMA

Authors: Ivana Beslic, Ashley T. Barnes, Frank Bigiel, Maria Jesus Jimenez-Donaire, Antonio Usero, Jonathan D. Henshaw, Christopher Faesi, Adam K. Leroy, Erik Rosolowsky, Jakob S. den Brok, Melanie Chevance, Cosima Eibensteiner, Kathryn Grasha, Ralf S. Klessen, J. M. Diedrerik Kruijssen, Daizhong Liu, Sharon Meidt, Justus Neumann, Lukas Neumann, Hsi-An Pan, Johannes Puschnig, Miguel Querejeta, Eva Schinnerer, Thomas G. Williams

Abstract: 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 - an… ▽ More 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. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: Accepted for publication to Astronomy and Astrophysics

The magnetic field in the Flame nebula

Authors: Ivana Bešlić, Simon Coudé, Dariusz C. Lis, Maryvonne Gerin, Paul F. Goldsmith, Jerome Pety, Antoine Roueff, Karine Demyk, Charles D. Dowell, Lucas Einig, Javier R. Goicoechea, Francois Levrier, Jan Orkisz, Nicolas Peretto, Miriam G. Santa-Maria, Nathalie Ysard, Antoine Zakardjian

Abstract: Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatia… ▽ More Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatial resolution observations are required to understand the transition from diffuse gas to gravitationally unstable filaments. NGC 2024, the Flame Nebula, in the nearby Orion B molecular cloud, contains a young, expanding HII region and a dense filament that harbors embedded protostellar objects. Therefore, NGC 2024 is an excellent opportunity to study the role of B-fields in the formation, evolution, and collapse of filaments, as well as the dynamics and effects of young HII regions on the surrounding molecular gas. We combine new 154 and 216 micron dust polarization measurements carried out using the HAWC+ instrument aboard SOFIA with molecular line observations of 12CN(1-0) and HCO+(1-0) from the IRAM 30-meter telescope to determine the B-field geometry and to estimate the plane of the sky magnetic field strength across the NGC 2024. The HAWC+ observations show an ordered B-field geometry in NGC 2024 that follows the morphology of the expanding HII region and the direction of the main filament. The derived plane of the sky B-field strength is moderate, ranging from 30 to 80 micro G. The strongest B-field is found at the northern-west edge of the HII region, characterized by the highest gas densities and molecular line widths. In contrast, the weakest field is found toward the filament in NGC 2024. The B-field has a non-negligible influence on the gas stability at the edges of the expanding HII shell (gas impacted by the stellar feedback) and the filament (site of the current star formation). △ Less

Submitted 7 February, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 36 pages, 26 figures Accepted for publication in Astronomy & Astrophysics

Surveying the Whirlpool at Arcseconds with NOEMA (SWAN)- I. Mapping the HCN and N$_2$H$^+$ 3mm lines

Authors: Sophia K. Stuber, Jerome Pety, Eva Schinnerer, Frank Bigiel, Antonio Usero, Ivana Beslić, Miguel Querejeta, María J. Jiménez-Donaire, Adam Leroy, Jakob den Brok, Lukas Neumann, Cosima Eibensteiner, Yu-Hsuan Teng, Ashley Barnes, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Simon C. O. Glover, Daizhong Liu, Hsi-An Pan

Abstract: We present the first results from "Surveying the Whirlpool at Arcseconds with NOEMA" (SWAN), an IRAM Northern Extended Millimetre Array (NOEMA)+30m large program that maps emission from several molecular lines at 90 and 110 GHz in the iconic nearby grand-design spiral galaxy M~51 at cloud-scale resolution ($\sim$3\arcsec=125\,pc). As part of this work, we have obtained the first sensitive cloud-sc… ▽ More We present the first results from "Surveying the Whirlpool at Arcseconds with NOEMA" (SWAN), an IRAM Northern Extended Millimetre Array (NOEMA)+30m large program that maps emission from several molecular lines at 90 and 110 GHz in the iconic nearby grand-design spiral galaxy M~51 at cloud-scale resolution ($\sim$3\arcsec=125\,pc). As part of this work, we have obtained the first sensitive cloud-scale map of N$_2$H$^+$(1-0) of the inner $\sim5\,\times 7\,$kpc of a normal star-forming galaxy, which we compare to HCN(1-0) and CO(1-0) emission to test their ability in tracing dense, star-forming gas. The average N$_2$H$^+$-to-HCN line ratio of our total FoV is $0.20\pm0.09$, with strong regional variations of a factor of $\gtrsim 2$ throughout the disk, including the south-western spiral arm and the center. The central $\sim1\,$kpc exhibits elevated HCN emission compared to N$_2$H$^+$, probably caused by AGN-driven excitation effects. We find that HCN and N$_2$H$^+$ are strongly super-linearily correlated in intensity ($ρ_\mathrm{Sp}\sim 0.8$), with an average scatter of $\sim0.14\,$dex over a span of $\gtrsim 1.5\,$dex in intensity. When excluding the central region, the data is best described by a power-law of exponent $1.2$, indicating that there is more N$_2$H$^+$ per unit HCN in brighter regions. Our observations demonstrate that the HCN-to-CO line ratio is a sensitive tracer of gas density in agreement with findings of recent Galactic studies which utilize N$_2$H$^+$. The peculiar line ratios present near the AGN and the scatter of the power-law fit in the disk suggest that in addition to a first-order correlation with gas density, second-order physics (such as optical depth, gas temperature) or chemistry (abundance variations) are encoded in the N$_2$H$^+$/CO, HCN/CO and N$_2$H$^+$/HCN ratios. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A. 6 pages, 3 figures (+ Appendix 4 pages, 2 figures)

A sensitive, high-resolution, wide-field IRAM NOEMA CO(1-0) survey of the very nearby spiral galaxy IC 342

Authors: M. Querejeta, J. Pety, A. Schruba, A. K. Leroy, C. N. Herrera, I-D. Chiang, S. E. Meidt, E. Rosolowsky, E. Schinnerer, K. Schuster, J. Sun, K. A. Herrmann, A. T. Barnes, I. Beslic, F. Bigiel, Y. Cao, M. Chevance, C. Eibensteiner, E. Emsellem, C. M. Faesi, A. Hughes, J. Kim, R. S. Klessen, K. Kreckel, J. M. D. Kruijssen , et al. (9 additional authors not shown)

Abstract: We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6" ~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into >600 individual giant molecular clouds and ass… ▽ More We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6" ~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into >600 individual giant molecular clouds and associations. We assessed their properties and found that overall the clouds show approximate virial balance, with typical virial parameters of alpha_vir=1-2. The typical surface density and line width of molecular gas increase from the inter-arm region to the arm and bar region, and they reach their highest values in the inner kiloparsec of the galaxy (median Sigma_mol~80, 140, 160, and 1100 M_sun/pc^2, sigma_CO~6.6, 7.6, 9.7, and 18.4 km/s for inter-arm, arm, bar, and center clouds, respectively). Clouds in the central part of the galaxy show an enhanced line width relative to their surface densities and evidence of additional sources of dynamical broadening. All of these results agree well with studies of clouds in more distant galaxies at a similar physical resolution. Leveraging our measurements to estimate the density and gravitational free-fall time at 90 pc resolution, averaged on 1.5 kpc hexagonal apertures, we estimate a typical star formation efficiency per free-fall time of 0.45% with a 16-84% variation of 0.33-0.71% among such 1.5 kpc regions. We speculate that bar-driven gas inflow could explain the large gas concentration in the central kiloparsec and the buildup of the massive nuclear star cluster. This wide-area CO map of the closest face-on massive spiral galaxy demonstrates the current mapping power of NOEMA and has many potential applications. The data and products are publicly available. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 16 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 680, A4 (2023)

HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30m maps of Orion B GMC: Revisiting its role as tracer of the dense gas reservoir for star formation

Authors: M. G. Santa-Maria, J. R. Goicoechea, J. Pety, M. Gerin, J. H. Orkisz, F. Le Petit, L. Einig, P. Palud, V. de Souza Magalhaes, I. Bešlić, L. Segal, S. Bardeau, E. Bron, P. Chainais, J. Chanussot, P. Gratier, V. V. Guzmán, A. Hughes, D. Languignon, F. Levrier, D. C. Lis, H. S. Liszt, J. Le Bourlot, Y. Oya, K. Öberg , et al. (6 additional authors not shown)

Abstract: We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V <… ▽ More We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V < 8 mag. The HCN/CO J=1-0 line intensity ratio shows a bimodal behavior with an inflection point at A_V < 3 mag typical of translucent gas and UV-illuminated cloud edges. We find that most of the HCN J=1-0 emission arises from extended gas with n(H2) ~< 10^4 cm^-3, even lower density gas if the ionization fraction is > 10^-5 and electron excitation dominates. This result explains the low-A_V branch of the HCN/CO J=1-0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~0.1) at A_V < 3 mag correspond to regions of high [CI] 492 GHz/CO J=1-0 intensity ratios (>1) characteristic of low-density PDRs. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps. The low surface brightness HCN and HCO+ J=1-0 emission scale with I_FIR (a proxy of the stellar FUV radiation field) in a similar way. Together with CO J=1-0, these lines respond to increasing I_FIR up to G0~20. On the other hand, the bright HCN J=1-0 emission from dense gas in star-forming clumps weakly responds to I_FIR once the FUV radiation field becomes too intense (G0>1500). The different power law scalings (produced by different chemistries, densities, and line excitation regimes) in a single but spatially resolved GMC resemble the variety of Kennicutt-Schmidt law indexes found in galaxy averages. As a corollary for extragalactic studies, we conclude that high HCN/CO J=1-0 line intensity ratios do not always imply the presence of dense gas. △ Less

Submitted 18 September, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: accepted for publication in A&A. 24 pages, 18 figures, plus Appendix. Abridged Abstract

Journal ref: A&A 679, A4 (2023)

Neural network-based emulation of interstellar medium models

Authors: Pierre Palud, Lucas Einig, Franck Le Petit, Emeric Bron, Pierre Chainais, Jocelyn Chanussot, Jérôme Pety, Pierre-Antoine Thouvenin, David Languignon, Ivana Bešlić, Miriam G. Santa-Maria, Jan H. Orkisz, Léontine E. Ségal, Antoine Zakardjian, Sébastien Bardeau, Maryvonne Gerin, Javier R. Goicoechea, Pierre Gratier, Viviana V. Guzman, Annie Hughes, François Levrier, Harvey S. Liszt, Jacques Le Bourlot, Antoine Roueff, Albrecht Sievers

Abstract: The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolatio… ▽ More The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolation of a grid of precomputed models. We propose a new general method to derive faster, lighter, and more accurate approximations of the model from a grid of precomputed models. These emulators are defined with artificial neural networks (ANNs) designed and trained to address the specificities inherent in ISM models. Indeed, such models often predict many observables (e.g., line intensities) from just a few input physical parameters and can yield outliers due to numerical instabilities or physical bistabilities. We propose applying five strategies to address these characteristics: 1) an outlier removal procedure; 2) a clustering method that yields homogeneous subsets of lines that are simpler to predict with different ANNs; 3) a dimension reduction technique that enables to adequately size the network architecture; 4) the physical inputs are augmented with a polynomial transform to ease the learning of nonlinearities; and 5) a dense architecture to ease the learning of simple relations. We compare the proposed ANNs with standard classes of interpolation methods to emulate the Meudon PDR code, a representative ISM numerical model. Combinations of the proposed strategies outperform all interpolation methods by a factor of 2 on the average error, reaching 4.5% on the Meudon PDR code. These networks are also 1000 times faster than accurate interpolation methods and require ten to forty times less memory. This work will enable efficient inferences on wide-field multiline observations of the ISM. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Journal ref: A&A 678, A198 (2023)

A constant N$_2$H$^+$(1-0)-to-HCN(1-0) ratio on kiloparsec scales

Authors: M. J. Jiménez-Donaire, A. Usero, I. Bešlić, M. Tafalla, A. Chacón-Tanarro, Q. Salomé, C. Eibensteiner, A. García-Rodríguez, A. Hacar, A. T. Barnes, F. Bigiel, M. Chevance, D. Colombo, D. A. Dale, T. A. Davis, S. C. O. Glover, J. Kauffmann, R. S. Klessen, A. K. Leroy, L. Neumann, H. Pan, J. Pety, M. Querejeta, T. Saito, E. Schinnerer , et al. (2 additional authors not shown)

Abstract: Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral gal… ▽ More Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral galaxy, NGC6946. We find a strong correlation ($r>0.98,p<0.01$) between the HCN(1-0) and N$_2$H$^+$(1-0) intensities across the inner $\sim8\mathrm{kpc}$ of the galaxy, at kiloparsec scales. This correlation is equally strong between the ratios N$_2$H$^+$(1-0)/CO(1-0) and HCN(1-0)/CO(1-0), tracers of dense gas fractions ($f_\mathrm{dense}$). We measure an average intensity ratio of N$_2$H$^+$(1-0)/HCN(1-0)$=0.15\pm0.02$ over our set of five IRAM-30m pointings. These trends are further supported by existing measurements for Galactic and extragalactic sources. This narrow distribution in the average ratio suggests that the observed systematic trends found in kiloparsec-scale extragalactic studies of $f_\mathrm{dense}$ and the efficiency of dense gas (SFE$_\mathrm{dense}$) would not change if we employed N$_2$H$^+$(1-0) as a more direct tracer of dense gas. At kiloparsec scales our results indicate that the HCN(1-0) emission can be used to predict the expected N$_2$H$^+$(1-0) over those regions. Our results suggest that, even if HCN(1-0) and N$_2$H$^+$(1-0) trace different density regimes within molecular clouds, subcloud differences average out at kiloparsec scales, yielding the two tracers proportional to each other. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 676, L11 (2023)

Deep learning denoising by dimension reduction: Application to the ORION-B line cubes

Authors: Lucas Einig, Jérôme Pety, Antoine Roueff, Paul Vandame, Jocelyn Chanussot, Maryvonne Gerin, Jan H. Orkisz, Pierre Palud, Miriam Garcia Santa-Maria, Victor de Souza Magalhaes, Ivana Bešlić, Sébastien Bardeau, Emeric E. Bron, Pierre Chainais, Javier R Goicoechea, Pierre Gratier, Viviana Guzman Veloso, Annie Hughes, Jouni Kainulainen, David Languignon, Rosine Lallement, François Levrier, Dariuscz C. Lis, Harvey Liszt, Jacques Le Bourlot , et al. (7 additional authors not shown)

Abstract: Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR… ▽ More Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed analysis of the noise statistics with an innovative autoencoder architecture is a promising path to denoise radio-astronomy line data cubes. In the future, exploring whether a better use of the spatial correlations of the noise may further improve the denoising performances seems a promising avenue. In addition, △ Less

Submitted 24 July, 2023; originally announced July 2023.

Journal ref: A&A 677, A158 (2023)

Spectral Stacking of Radio-Interferometric Data

Authors: Lukas Neumann, Jakob S. den Brok, Frank Bigiel, Adam Leroy, Antonio Usero, Ashley T. Barnes, Ivana Bešlić, Cosima Eibensteiner, Malena Held, María J. Jiménez-Donaire, Jérôme Pety, Erik W. Rosolowsky, Eva Schinnerer, Thomas G. Williams

Abstract: 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 extern… ▽ More 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. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: 10 pages, 10 figures, accepted for pub in A&A, Apr 28, 2023

Journal ref: A&A 675, A104 (2023)

The impact of HII regions on Giant Molecular Cloud properties in nearby galaxies sampled by PHANGS ALMA and MUSE

Authors: Antoine Zakardjian, Jérôme Pety, Cinthya N. Herrera, Annie Hughes, Elias Oakes, Kathryn Kreckel, Chris Faesi, Simon C. O. Glover, Brent Groves, Ralf S. Klessen, Sharon Meidt, Ashley Barnes, Francesco Belfiore, Ivana Bešlić, Frank Bigiel, Guillermo A. Blanc, Mélanie Chevance, Daniel A. Dale, Jakob den Brok, Cosima Eibensteiner, Eric Emsellem, Axel García-Rodríguez, Kathryn Grasha, Eric W. Koch, Adam K. Leroy , et al. (14 additional authors not shown)

Abstract: We identify giant molecular clouds (GMCs) associated with HII regions for a sample of 19 nearby galaxies using catalogs of GMCs and H regions released by the PHANGS-ALMA and PHANGS-MUSE surveys, using the overlap of the CO and Hα emission as the key criterion for physical association. We compare the distributions of GMC and HII region properties for paired and non-paired objects. We investigate co… ▽ More We identify giant molecular clouds (GMCs) associated with HII regions for a sample of 19 nearby galaxies using catalogs of GMCs and H regions released by the PHANGS-ALMA and PHANGS-MUSE surveys, using the overlap of the CO and Hα emission as the key criterion for physical association. We compare the distributions of GMC and HII region properties for paired and non-paired objects. We investigate correlations between GMC and HII region properties among galaxies and across different galactic environments to determine whether GMCs that are associated with HII regions have significantly distinct physical properties to the parent GMC population. We identify trends between the Hα luminosity of an HII region and the CO peak brightness and the molecular mass of GMCs that we tentatively attribute to a direct physical connection between the matched objects, and which arise independently of underlying environmental variations of GMC and HII region properties within galaxies. The study of the full sample nevertheless hides a large variability galaxy by galaxy. Our results suggests that at the ~100 pc scales accessed by the PHANGS-ALMA and PHANGS-MUSE data, pre-supernova feedback mechanisms in HII regions have a subtle but measurable impact on the properties of the surrounding molecular gas, as inferred from CO observations. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Journal ref: A&A 678, A171 (2023)

Kinematic analysis of the super-extended HI disk of the nearby spiral galaxy M83

Authors: Cosima Eibensteiner, Frank Bigiel, Adam K. Leroy, Eric W. Koch, Erik Rosolowsky, Eva Schinnerer, Amy Sardone, Sharon Meidt, W. J. G de Blok, David Thilker, D. J. Pisano, Jürgen Ott, Ashley Barnes, Miguel Querejeta, Eric Emsellem, Johannes Puschnig, Dyas Utomo, Ivana Bešlic, Jakob den Brok, Shahram Faridani, Simon C. O. Glover, Kathryn Grasha, Hamid Hassani, Jonathan D. Henshaw, Maria J. Jiménez-Donaire , et al. (11 additional authors not shown)

Abstract: 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… ▽ More 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 ($σ_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 $σ_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. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: accepted for publication in A&A; 16 pages, 12 figures (+8 pages appendix)

Journal ref: A&A 675, A37 (2023)

Wide-field CO isotopologue emission and the CO-to-H$_2$ factor across the nearby spiral galaxy M101

Authors: Jakob S. den Brok, Frank Bigiel, Jérémy Chastenet, Karin Sandstrom, Adam Leroy, Antonio Usero, Eva Schinnerer, Erik W. Rosolowsky, Eric W. Koch, I-Da Chiang, Ashley T. Barnes, Johannes Puschnig, Toshiki Saito, Ivana Bešlić, Melanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Simon Glover, María J. Jiménez-Donaire, Yu-Hsuan Teng, Thomas G. Williams

Abstract: Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $α_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $α_{\rm CO}$ across all environments across galaxies as a function of me… ▽ More Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $α_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $α_{\rm CO}$ across all environments across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using resolved spectral line observations of CO and its isotopologues, we can constrain the molecular gas conditions and link them to a variation in the conversion factor. We present new IRAM 30-m 1mm and 3mm line observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O} across the nearby galaxy M101. Based on the CO isotopologue line ratios, we find that selective nucleosynthesis and opacity changes are the main drivers of the variation in the line emission across the galaxy. Furthermore, we estimated $α_{\rm CO(1-0)}$ using different approaches, including (i) the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) LTE-based measurements using the optically thin $^{13}$CO(1-0) intensity. We find an average value of $α_{\rm CO}=4.4{\pm}0.9\rm\,M_\odot\,pc^{-2}(K\,km\,s^{-1})^{-1}$ across the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based values are lower by a factor of 2-3 across the disk relative to the dust-based result. Accounting for $α_{\rm CO}$ variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy's center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if $α_{\rm CO}$ variations are not accounted for. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A, 25 pages, 15 figures

Journal ref: A&A 676, A93 (2023)

The ALMOND Survey: Molecular cloud properties and gas density tracers across 25 nearby spiral galaxies with ALMA

Authors: Lukas Neumann, Molly J. Gallagher, Frank Bigiel, Adam K. Leroy, Ashley T. Barnes, Antonio Usero, Jakob S. den Brok, Francesco Belfiore, Ivana Bešlić, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Simon C. O. Glover, Kathryn Grasha, Jonathan D. Henshaw, María J. Jiménez-Donaire, Ralf S. Klessen, J. M. Diederik Kruijssen, Daizhong Liu, Sharon Meidt, Jérôme Pety, Johannes Puschnig, Miguel Querejeta, Erik Rosolowsky , et al. (6 additional authors not shown)

Abstract: 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 i… ▽ More 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. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 48 pages, 40 figures

Sub-kiloparsec empirical relations and excitation conditions of HCN and HCO+ J=3-2 in nearby star-forming galaxies

Authors: Axel Garcia-Rodriguez, Antonio Usero, Adam K. Leroy, Frank Bigiel, Maria Jesus Jimenez-Donaire, Daizhong Liu, Miguel Querejeta, Toshiki Saito, Eva Schinnerer, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Yixian Cao, Melanie Chevance, Daniel A. Dale, Jakob S. den Brok, Cosima Eibensteiner, Santiago Garcia-Burillo, Simon C. O. Glover, Ralf S. Klessen, Jerome Pety, Johannes Puschnig, Erik Rosolowsky, Karin Sandstrom, Mattia C. Sormani , et al. (2 additional authors not shown)

Abstract: We present new HCN and HCO$^+$ ($J$=3-2) images of the nearby star-forming galaxies (SFGs) NGC 3351, NGC 3627, and NGC 4321. The observations, obtained with the Morita ALMA Compact Array, have a spatial resolution of $\sim$290-440 pc and resolve the inner $R_\textrm{gal} \lesssim$ 0.6-1 kpc of the targets, as well as the southern bar end of NGC 3627. We complement this data set with publicly avail… ▽ More We present new HCN and HCO$^+$ ($J$=3-2) images of the nearby star-forming galaxies (SFGs) NGC 3351, NGC 3627, and NGC 4321. The observations, obtained with the Morita ALMA Compact Array, have a spatial resolution of $\sim$290-440 pc and resolve the inner $R_\textrm{gal} \lesssim$ 0.6-1 kpc of the targets, as well as the southern bar end of NGC 3627. We complement this data set with publicly available images of lower excitation lines of HCN, HCO$^+$, and CO and analyse the behaviour of a representative set of line ratios: HCN(3-2)/HCN(1-0), HCN(3-2)/HCO$^+$(3-2), HCN(1-0)/CO(2-1), and HCN(3-2)/CO(2-1). Most of these ratios peak at the galaxy centres and decrease outwards. We compare the HCN and HCO$^+$ observations with a grid of one-phase, non-local thermodynamic equilibrium (non-LTE) radiative transfer models and find them compatible with models that predict subthermally excited and optically thick lines. We study the systematic variations of the line ratios across the targets as a function of the stellar surface density ($Σ_\textrm{star}$), the intensity-weighted CO(2-1) ($\langle I_\text{CO}\rangle$), and the star formation rate surface density ($Σ_\text{SFR}$). We find no apparent correlation with $Σ_\text{SFR}$, but positive correlations with the other two parameters, which are stronger in the case of $\langle I_\text{CO}\rangle$. The HCN/CO-$\langle I_\text{CO}\rangle$ relations show $\lesssim$0.3 dex galaxy-to-galaxy offsets, with HCN(3-2)/CO(2-1)-$\langle I_\text{CO}\rangle$ being $\sim$2 times steeper than HCN(1-0)/CO(2-1). In contrast, the HCN(3-2)/HCN(1-0)-$\langle I_\text{CO}\rangle$ relation exhibits a tighter alignment between galaxies. We conclude that the overall behaviour of the line ratios cannot be ascribed to variations in a single excitation parameter (e.g. density or temperature). △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A. 14 pages, 8 figures

Journal ref: A&A 672, A96 (2023)

PHANGS-JWST First Results: Mid-infrared emission traces both gas column density and heating at 100 pc scales

Authors: Adam K. Leroy, Karin Sandstrom, Erik Rosolowsky, Francesco Belfiore, Alberto D. Bolatto, Yixian Cao, Eric W. Koch, Eva Schinnerer, Ashley. T. Barnes, Ivana Bešlić, F. Bigiel, Guillermo A. Blanc, Jérémy Chastenet, Ness Mayker Chen, Mélanie Chevance, Ryan Chown, Enrico Congiu, Daniel A. Dale, Oleg V. Egorov, Eric Emsellem, Cosima Eibensteiner, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Brent Groves , et al. (26 additional authors not shown)

Abstract: We compare mid-infrared (mid-IR), extinction-corrected H$α$, and CO (2-1) emission at 70--160 pc resolution in the first four PHANGS-JWST targets. We report correlation strengths, intensity ratios, and power law fits relating emission in JWST's F770W, F1000W, F1130W, and F2100W bands to CO and H$α$. At these scales, CO and H$α$ each correlate strongly with mid-IR emission, and these correlations a… ▽ More We compare mid-infrared (mid-IR), extinction-corrected H$α$, and CO (2-1) emission at 70--160 pc resolution in the first four PHANGS-JWST targets. We report correlation strengths, intensity ratios, and power law fits relating emission in JWST's F770W, F1000W, F1130W, and F2100W bands to CO and H$α$. At these scales, CO and H$α$ each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H$α$ emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to the good match with the molecular gas-tracing CO, and as a heating tracer, leading to the good match with the H$α$. By combining mid-IR, CO, and H$α$ at scales where the overall correlation between cold gas and star formation begins to break down, we are able to separate these two effects. We model the mid-IR above $I_ν= 0.5$~MJy sr$^{-1}$ at F770W, a cut designed to select regions where the molecular gas dominates the interstellar medium (ISM) mass. This bright emission can be described to first order by a model that combines a CO-tracing component and an H$α$-tracing component. The best-fitting models imply that $\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. △ Less

Submitted 6 January, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 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

CI and CO in Nearby Spiral Galaxies -- I. Line Ratio and Abundance Variations at ~ 200 pc Scales

Authors: Daizhong Liu, Eva Schinnerer, Toshiki Saito, Erik Rosolowsky, Adam Leroy, Antonio Usero, Karin Sandstrom, Ralf S. Klessen, Simon C. O. Glover, Yiping Ao, Ivana Bešlić, Frank Bigiel, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, Daniel A. Dale, Yu Gao, Annie Hughes, Kathryn Kreckel, J. M. Diederik Kruijssen, Hsi-An Pan, Jérôme Pety, Dragan Salak, Francesco Santoro, Andreas Schruba , et al. (3 additional authors not shown)

Abstract: We present new neutral atomic carbon [CI](3P1-3P0) mapping observations within the inner ~7 kpc and ~4 kpc of the disks of NGC3627 and NGC4321 at a spatial resolution of 190 pc and 270 pc, respectively, using the ALMA Atacama Compact Array (ACA). We combine these with the CO(2-1) data from PHANGS-ALMA, and literature [CI] and CO data for two other starburst and/or active galactic nucleus (AGN) gal… ▽ More We present new neutral atomic carbon [CI](3P1-3P0) mapping observations within the inner ~7 kpc and ~4 kpc of the disks of NGC3627 and NGC4321 at a spatial resolution of 190 pc and 270 pc, respectively, using the ALMA Atacama Compact Array (ACA). We combine these with the CO(2-1) data from PHANGS-ALMA, and literature [CI] and CO data for two other starburst and/or active galactic nucleus (AGN) galaxies (NGC1808, NGC7469), to study: a) the spatial distributions of CI and CO emission; b) the observed line ratio RCICO = I_[CI](1-0)/I_CO(2-1) as a function of various galactic properties; and c) the abundance ratio of [CI/CO]. We find excellent spatial correspondence between CI and CO emission and nearly uniform RCICO ~0.1 across the majority of the star-forming disks of NGC3627 and NGC4321. However, RCICO strongly varies from ~0.05 at the centre of NGC4321 to >0.2-0.5 in NGC1808's starburst centre and NGC7469's centre with an X-ray AGN. Meanwhile, RCICO does not obviously vary with $U$, similar to the prediction of PDR models. We also find a mildly decreasing RCICO with an increasing metallicity over 0.7-0.85 solar metallicity, consistent with the literature. Assuming various typical ISM conditions representing GMCs, active star-forming regions and strong starbursting environments, we calculate the LTE radiative transfer and estimate the [CI/CO] abundance ratio to be ~0.1 across the disks of NGC3627 and NGC4321, similar to previous large-scale findings in Galactic studies. However, this abundance ratio likely has a substantial increase to ~1 and >1-5 in NGC1808's starburst and NGC7469's strong AGN environments, respectively, in line with the expectations for cosmic-ray dominated region (CRDR) and X-ray dominated region (XDR) chemistry. Finally, we do not find a robust evidence for a generally CO-dark, CI-bright gas in the disk areas we probed. (abbreviated) △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 23 pages, 13 figures and one table in total (17 pages and 9 figures in main text). Accepted for publication in A&A. For associated data cubes and moment maps, see https://www.canfar.net/storage/vault/list/phangs/RELEASES/DZLIU_etal_2022

Journal ref: A&A 672, A36 (2023)

PHANGS-JWST First Results: Stellar Feedback-Driven Excitation and Dissociation of Molecular Gas in the Starburst Ring of NGC 1365?

Authors: Daizhong Liu, Eva Schinnerer, Yixian Cao, Adam Leroy, Antonio Usero, Erik Rosolowsky, Eric Emsellem, J. M. Diederik Kruijssen, Mélanie Chevance, Simon C. O. Glover, Mattia C. Sormani, Alberto D. Bolatto, Jiayi Sun, Sophia K. Stuber, Yu-Hsuan Teng, Frank Bigiel, Ivana Bešlić, Kathryn Grasha, Jonathan D. Henshaw, Ashley. T. Barnes, Jakob S. den Brok, Toshiki Saito, Daniel A. Dale, Elizabeth J. Watkins, Hsi-An Pan , et al. (14 additional authors not shown)

Abstract: We compare embedded young massive star clusters (YMCs) to (sub-)millimeter line observations tracing the excitation and dissociation of molecular gas in the starburst ring of NGC 1365. This galaxy hosts one of the strongest nuclear starbursts and richest populations of YMCs within 20 Mpc. Here we combine near-/mid-IR PHANGS-JWST imaging with new ALMA multi-J CO (1-0, 2-1 and 4-3) and [CI](1-0) map… ▽ More We compare embedded young massive star clusters (YMCs) to (sub-)millimeter line observations tracing the excitation and dissociation of molecular gas in the starburst ring of NGC 1365. This galaxy hosts one of the strongest nuclear starbursts and richest populations of YMCs within 20 Mpc. Here we combine near-/mid-IR PHANGS-JWST imaging with new ALMA multi-J CO (1-0, 2-1 and 4-3) and [CI](1-0) mapping, which we use to trace CO excitation via R42 = I_CO(4-3)/I_CO(2-1) and R21 = I_CO(2-1)/I_CO(1-0) and dissociation via RCICO = I_[CI](1-0)/I_CO(2-1) at 330 pc resolution. We find that the gas flowing into the starburst ring from northeast to southwest appears strongly affected by stellar feedback, showing decreased excitation (lower R42) and increased signatures of dissociation (higher RCICO) in the downstream regions. There, radiative transfer modeling suggests that the molecular gas density decreases and temperature and [CI/CO] abundance ratio increase. We compare R42 and RCICO with local conditions across the regions and find that both correlate with near-IR 2 um emission tracing the YMCs and with both PAH (11.3 um) and dust continuum (21 um) emission. In general, RCICO exhibits ~ 0.1 dex tighter correlations than R42, suggesting CI to be a more sensitive tracer of changing physical conditions in the NGC 1365 starburst than CO (4-3). Our results are consistent with a scenario where gas flows into the two arm regions along the bar, becomes condensed/shocked, forms YMCs, and then these YMCs heat and dissociate the gas. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 19 pages, 7 figures and 2 tables in total (12 pages and 6 figures in main text). Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJL

The PHANGS-JWST Treasury Survey: Star Formation, Feedback, and Dust Physics at High Angular resolution in Nearby GalaxieS

Authors: Janice C. Lee, Karin M. Sandstrom, Adam K. Leroy, David A. Thilker, Eva Schinnerer, Erik Rosolowsky, Kirsten L. Larson, Oleg V. Egorov, Thomas G. Williams, Judy Schmidt, Eric Emsellem, Gagandeep S. Anand, Ashley T. Barnes, Francesco Belfiore, Ivana Beslic, Frank Bigiel, Guillermo A. Blanc, Alberto D. Bolatto, Mederic Boquien, Jakob den Brok, Yixian Cao, Rupali Chandar, Jeremy Chastenet, Melanie Chevance, I-Da Chiang , et al. (52 additional authors not shown)

Abstract: The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium 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 ($\sim$5-50 pc)… ▽ More The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium 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 ($\sim$5-50 pc). In Cycle 1, PHANGS is conducting an 8-band imaging survey from 2-21$μ$m of 19 nearby spiral galaxies. CO(2-1) mapping, optical integral field spectroscopy, 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 $\sim$10,000 star clusters, MUSE spectroscopic mapping of $\sim$20,000 HII regions, and $\sim$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. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: Re-submitted after addressing minor comments from referee. To be published as part of PHANGS-JWST ApJL Focus Issue

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

Authors: Jakob S. den Brok, Frank Bigiel, Kazimierz Sliwa, Toshiki Saito, Antonio Usero, Eva Schinnerer, Adam K. Leroy, María J. Jiménez-Donaire, Erik Rosolowsky, Ashley T. Barnes, Johannes Puschnig, Jérôme Pety, Andreas Schruba, Ivana Bešlić, Yixian Cao, Cosima Eibensteiner, Simon C. O. Glover, Ralf S. Klessen, J. M. Diederik Kruijssen, Sharon E. Meidt, Lukas Neumann, Neven Tomičić, Hsi-An Pan, Miguel Querejeta, Elizabeth Watkins , et al. (2 additional authors not shown)

Abstract: We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical c… ▽ More We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 29 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A 662, A89 (2022)

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

Authors: Cosima Eibensteiner, Ashley T. Barnes, Frank Bigiel, Eva Schinnerer, Daizhong Liu, David S. Meier, Antonio Usero, Adam K. Leroy, Erik Rosolowsky, Johannes Puschnig, Ilin Lazar, Jérôme Pety, Laura A. Lopez, Eric Emsellem, Ivana Bešlić, Miguel Querejeta, Eric J. Murphy, Jakob den Brok, Andreas Schruba, Mélanie Chevance, Simon C. O. Glover, Yu Gao, Kathryn Grasha, Hamid Hassani, Jonathan D. Henshaw , et al. (8 additional authors not shown)

Abstract: The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spec… ▽ More The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies. △ Less

Submitted 10 January, 2022; v1 submitted 6 January, 2022; originally announced January 2022.

Comments: 38 pages, 23 figures, accepted for publication in A&A

Journal ref: A&A 659, A173 (2022)

Comparing the pre-SNe feedback and environmental pressures for 6000 HII regions across 19 nearby spiral galaxies

Authors: A. T. Barnes, S. C. O. Glover, K. Kreckel, E. C. Ostriker, F. Bigiel, F. Belfiore, I. Bešlić, G. A. Blanc, M. Chevance, D. A. Dale, O. Egorov, C. Eibensteiner, E. Emsellem, K. Grasha, B. A. Groves, R. S. Klessen, J. M. D. Kruijssen, A. K. Leroy, S. N. Longmore, L. Lopez, R. McElroy, S. E. Meidt, E. J. Murphy, E. Rosolowsky, T. Saito , et al. (6 additional authors not shown)

Abstract: 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 ($<$ 20 Mpc) star-forming, main sequence spiral galaxies (log($M_\star$/M$_\odot$)= 9.4 $-$ 11). These optical spectroscopic maps are essen… ▽ More 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 ($<$ 20 Mpc) star-forming, main sequence spiral galaxies (log($M_\star$/M$_\odot$)= 9.4 $-$ 11). These optical spectroscopic maps are essential to constrain the physical properties of the HII regions, which we use to investigate their internal pressure terms. We estimate the photoionised gas ($P_\mathrm{therm}$), direct radiation ($P_\mathrm{rad}$), and mechanical wind pressure ($P_\mathrm{wind}$), which we compare to the confining pressure of their host environment ($P_\mathrm{de}$). The HII regions remain unresolved within our ${\sim}50{-}100$ pc resolution observations, so we place upper ($P_\mathrm{max}$) and lower ($P_\mathrm{min}$) limits on each of the pressures by using a minimum (i.e. clumpy structure) and maximum (i.e. smooth structure) size, respectively. We find that the $P_\mathrm{max}$ measurements are broadly similar, and for $P_\mathrm{min}$ the $P_\mathrm{therm}$ is mildly dominant. We find that the majority of HII regions are over-pressured, $P_\mathrm{tot}/P_\mathrm{de} = (P_\mathrm{therm}+P_\mathrm{wind}+P_\mathrm{rad})/P_\mathrm{de} > 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}<1$kpc), or, specifically, environments of high gas surface density; log($Σ_\mathrm{gas}/\mathrm{M_\odot} \mathrm{pc}^{-2}$)$\sim$2.5 (measured on kpc-scales). Lastly, we compare to a sample of literature measurements for $P_\mathrm{therm}$ and $P_\mathrm{rad}$ to investigate how dominant pressure term transitions over around 5dex in spatial dynamic range and 10 dex in pressure. △ Less

Submitted 25 October, 2021; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 28 pages total, 16+3 figures, and 3 tables. Accepted for publication in MNRAS

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

Authors: Eric Emsellem, Eva Schinnerer, Francesco Santoro, Francesco Belfiore, Ismael Pessa, Rebecca McElroy, Guillermo A. Blanc, Enrico Congiu, Brent Groves, I-Ting Ho, Kathryn Kreckel, Alessandro Razza, Patricia Sanchez-Blazquez, Oleg Egorov, Chris Faesi, Ralf S. Klessen, Adam K. Leroy, Sharon Meidt, Miguel Querejeta, Erik Rosolowsky, Fabian Scheuermann, Gagandeep S. Anand, Ashley T. Barnes, Ivana Bešlić, Frank Bigiel , et al. (23 additional authors not shown)

Abstract: 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 loca… ▽ More 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. △ Less

Submitted 5 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 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

Journal ref: A&A 659, A191 (2022)

Stellar structures, molecular gas, and star formation across the PHANGS sample of nearby galaxies

Authors: M. Querejeta, E. Schinnerer, S. Meidt, J. Sun, A. K. Leroy, E. Emsellem, R. S. Klessen, J. C. Munoz-Mateos, H. Salo, E. Laurikainen, I. Beslic, G. A. Blanc, M. Chevance, D. A. Dale, C. Eibensteiner, C. Faesi, A. Garcia-Rodriguez, S. C. O. Glover, K. Grasha, J. Henshaw, C. Herrera, A. Hughes, K. Kreckel, J. M. D. Kruijssen, D. Liu , et al. (10 additional authors not shown)

Abstract: We identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on Spitzer 3.6 micron images. At the simplest level, we distinguish centres, bars, spiral arms, interarm and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, publicly released but not explicitly used in this paper. W… ▽ More We identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on Spitzer 3.6 micron images. At the simplest level, we distinguish centres, bars, spiral arms, interarm and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, publicly released but not explicitly used in this paper. We examine trends using PHANGS-ALMA CO(2-1) intensity maps and tracers of star formation. The interarm regions and discs without strong spirals dominate in area, whereas molecular gas and star formation are quite evenly distributed among the five basic environments. We reproduce the molecular Kennicutt-Schmidt relation with a slope compatible with unity within the uncertainties, without significant slope differences among environments. In contrast to early studies, we find that bars are not always deserts devoid of gas and star formation, but instead they show large diversity. Similarly, spiral arms do not account for most of the gas and star formation in disc galaxies, and they do not have shorter depletion times than the interarm regions. Spiral arms accumulate gas and star formation, without systematically boosting the star formation efficiency. Centres harbour remarkably high surface densities and on average shorter depletion times than other environments. Centres of barred galaxies show higher surface densities and wider distributions compared to the outer disc; yet, depletion times are similar to unbarred galaxies, suggesting highly intermittent periods of star formation when bars episodically drive gas inflow, without enhancing the central star formation efficiency permanently. In conclusion, we provide quantitative evidence that stellar structures in galaxies strongly affect the organisation of molecular gas and star formation, but their impact on star formation efficiency is more subtle. △ Less

Submitted 7 October, 2021; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 28 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 656, A133 (2021)

Dense molecular gas properties on 100 pc scales across the disc of NGC 3627

Authors: I. Bešlić, A. T. Barnes, F. Bigiel, J. Puschnig, J. Pety, C. Herrera Contreras, A. K. Leroy, A. Usero, E. Schinnerer, S. E. Meidt, E. Emsellem, A. Hughes, C. Faesi, K. Kreckel, F. M. C. Belfiore, M. Chevance, J. S. den Brok, C. Eibensteiner, S. C. O. Glover, K. Grasha, M. J. Jimenez-Donaire, R. S. Klessen, J. M. D. Kruijssen, D. Liu, I. Pessa , et al. (7 additional authors not shown)

Abstract: 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 cri… ▽ More 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. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: 25 pages, 12 figures; accepted for publication in MNRAS

PHANGS-ALMA: Arcsecond CO(2-1) Imaging of Nearby Star-Forming Galaxies

Authors: Adam K. Leroy, Eva Schinnerer, Annie Hughes, Erik Rosolowsky, Jérôme Pety, Andreas Schruba, Antonio Usero, Guillermo A. Blanc, Mélanie Chevance, Eric Emsellem, Christopher M. Faesi, Cinthya N. Herrera, Daizhong Liu, Sharon E. Meidt, Miguel Querejeta, Toshiki Saito, Karin M. Sandstrom, Jiayi Sun, Thomas G. Williams, Gagandeep S. Anand, Ashley T. Barnes, Erica A. Behrens, Francesco Belfiore, Samantha M. Benincasa, Ivana Bešlić , et al. (47 additional authors not shown)

Abstract: 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<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PH… ▽ More 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<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PHANGS-ALMA, each beam reaches the size of a typical individual giant molecular cloud (GMC), so that these data can be used to measure the demographics, life-cycle, and physical state of molecular clouds across the population of galaxies where the majority of stars form at z=0. This paper describes the scientific motivation and background for the survey, sample selection, global properties of the targets, ALMA observations, and characteristics of the delivered ALMA data and derived data products. As the ALMA sample serves as the parent sample for parallel surveys with VLT/MUSE, HST, AstroSat, VLA, and other facilities, we include a detailed discussion of the sample selection. We detail the estimation of galaxy mass, size, star formation rate, CO luminosity, and other properties, compare estimates using different systems and provide best-estimate integrated measurements for each target. We also report the design and execution of the ALMA observations, which combine a Cycle~5 Large Program, a series of smaller programs, and archival observations. Finally, we present the first 1" resolution atlas of CO emission from nearby galaxies and describe the properties and contents of the first PHANGS-ALMA public data release. △ Less

Submitted 28 April, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 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

PHANGS-ALMA Data Processing and Pipeline

Authors: Adam K. Leroy, Annie Hughes, Daizhong Liu, Jerome Pety, Erik Rosolowsky, Toshiki Saito, Eva Schinnerer, Andreas Schruba, Antonio Usero, Christopher M. Faesi, Cinthya N. Herrera, Melanie Chevance, Alexander P. S. Hygate, Amanda A. Kepley, Eric W. Koch, Miguel Querejeta, Kazimierz Sliwa, David Will, Christine D. Wilson, Gagandeep S. Anand, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Frank Bigiel, Guillermo A. Blanc , et al. (43 additional authors not shown)

Abstract: We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and… ▽ More We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u-v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to VLA and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12-m+7-m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7-m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS--ALMA and test the performance of popular short-spacing correction algorithms. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Astrophysical Journal Supplement series. 65 pages, 33 figures. Software available at https://github.com/akleroy/phangs_imaging_scripts . For a full resolution version see https://sites.google.com/view/phangs/publications

Applying the Tremaine-Weinberg Method to Nearby Galaxies: Stellar Mass-Based Pattern Speeds, and Comparisons with ISM Kinematics

Authors: Thomas G. Williams, Eva Schinnerer, Eric Emsellem, Sharon Meidt, Miguel Querejeta, Francesco Belfiore, Ivana Bešlić, Frank Bigiel, Mélanie Chevance, Daniel A. Dale, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, J. M. Diederik Kruijssen, Adam K. Leroy, Hsi-An Pan, Jérôme Pety, Ismael Pessa, Erik Rosolowsky, Toshiki Saito, Francesco Santoro, Andreas Schruba, Mattia C. Sormani, Jiayi Sun, Elizabeth J. Watkins

Abstract: We apply the Tremaine-Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds ($Ω_{\rm P}$). After quality checks, we find that around half (10) of these stellar mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously publish… ▽ More We apply the Tremaine-Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds ($Ω_{\rm P}$). After quality checks, we find that around half (10) of these stellar mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously published pattern speeds, and use rotation curves to calculate major resonance locations (co-rotation radii and Lindblad resonances). We also compare these stellar-mass derived pattern speeds with H$α$ (from MUSE) and CO($J=2{-}1$) emission from the PHANGS-ALMA survey. We find that in the case of these clumpy ISM tracers, this method erroneously gives a signal that is simply the angular frequency at a representative radius set by the distribution of these clumps ($Ω_{\rm clump}$), and that this $Ω_{\rm clump}$ is significantly different to $Ω_{\rm P}$ ($\sim$20% in the case of H$α$, and $\sim$50% in the case of CO). Thus, we conclude that it is inadvisable to use "pattern speeds" derived from ISM kinematics. Finally, we compare our derived pattern speeds and co-rotation radii, along with bar properties, to the global parameters of these galaxies. Consistent with previous studies, we find that galaxies with a later Hubble type have a larger ratio of co-rotation radius to bar length, more molecular-gas rich galaxies have higher $Ω_{\rm P}$, and more bulge-dominated galaxies have lower $Ω_{\rm P}$. Unlike earlier works, however, there are no clear trends between the bar strength and $Ω_{\rm P}$, nor between the total stellar mass surface density and the pattern speed. △ Less

Submitted 15 March, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 22 pages, 13 figures, plus an additional 24 pages and 19 figures in appendices. Accepted for publication in AJ. Include DOI, as well as link to GitHub page containing .fits machine-readable version of Table 2

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

Authors: Erik Rosolowsky, Annie Hughes, Adam K. Leroy, Jiayi Sun, Miguel Querejeta, Andreas Schruba, Antonio Usero, Cinthya N. Herrera, Daizhong Liu, Jérôme Pety, Toshiki Saito, Ivana Bešlić, Frank Bigiel, Guillermo Blanc, Mélanie Chevance, Daniel A. Dale, Sinan Deger, Christopher M. Faesi, Simon C. O. Glover, Jonathan D. Henshaw, Ralf S. Klessen, J. M. Diederik Kruijssen, Kirsten Larson, Janice Lee, Sharon Meidt , et al. (4 additional authors not shown)

Abstract: We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data… ▽ More We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available python package, PYCPROPS. We apply these methods to ten of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4). △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: Accepted to MNRAS. Catalogues, code, and atlas available at http://phangs.org/data